/[ascend]/trunk/base/generic/solver/lsode.c
ViewVC logotype

Diff of /trunk/base/generic/solver/lsode.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1133 by johnpye, Sun Jan 14 11:51:48 2007 UTC revision 1134 by johnpye, Sun Jan 14 11:56:30 2007 UTC
# Line 946  int integrator_lsode_solve(IntegratorSys Line 946  int integrator_lsode_solve(IntegratorSys
946      /* read METH and MITER parameters, create MF value */      /* read METH and MITER parameters, create MF value */
947      method = SLV_PARAM_CHAR(&(blsys->params),LSODE_PARAM_METH);      method = SLV_PARAM_CHAR(&(blsys->params),LSODE_PARAM_METH);
948      miter = SLV_PARAM_INT(&(blsys->params),LSODE_PARAM_MITER);      miter = SLV_PARAM_INT(&(blsys->params),LSODE_PARAM_MITER);
949      maxord = SLV_PARAM_INT(&(blsys->params),LSODE_PARAM_MITER);      maxord = SLV_PARAM_INT(&(blsys->params),LSODE_PARAM_MAXORD);
950      if(miter < 0 || miter > 3){      if(miter < 0 || miter > 3){
951          ERROR_REPORTER_HERE(ASC_USER_ERROR,"Unacceptable value '%d' of parameter 'miter'",miter);          ERROR_REPORTER_HERE(ASC_USER_ERROR,"Unacceptable value '%d' of parameter 'miter'",miter);
952          return 5;          return 5;
# Line 1029  int integrator_lsode_solve(IntegratorSys Line 1029  int integrator_lsode_solve(IntegratorSys
1029    rwork[6] = integrator_get_minstep(blsys);    rwork[6] = integrator_get_minstep(blsys);
1030    iwork[5] = integrator_get_maxsubsteps(blsys);    iwork[5] = integrator_get_maxsubsteps(blsys);
1031      iwork[4] = maxord;      iwork[4] = maxord;
1032        CONSOLE_DEBUG("MAXORD = %d",maxord);
1033    
1034    if(x[0] > integrator_getsample(blsys, 2)){    if(x[0] > integrator_getsample(blsys, 2)){
1035      ERROR_REPORTER_HERE(ASC_USER_ERROR,"Invalid initialisation time: exceeds second timestep value");      ERROR_REPORTER_HERE(ASC_USER_ERROR,"Invalid initialisation time: exceeds second timestep value");

Legend:
Removed from v.1133  
changed lines
  Added in v.1134

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22