/[ascend]/trunk/base/generic/integrator/integrator.c
ViewVC logotype

Diff of /trunk/base/generic/integrator/integrator.c

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

revision 1386 by jpye, Sat Mar 10 13:57:01 2007 UTC revision 1387 by jpye, Sat Apr 7 14:43:31 2007 UTC
# Line 40  Line 40 
40  #include "ida.h"  #include "ida.h"
41  #include "samplelist.h"  #include "samplelist.h"
42    
43  /* #define ANALYSE_DEBUG */  #define ANALYSE_DEBUG
44  /* #define SOLVE_DEBUG */  /* #define SOLVE_DEBUG */
45  /* #define CLASSIFY_DEBUG */  /* #define CLASSIFY_DEBUG */
46  /* #define DESTROY_DEBUG */  /* #define DESTROY_DEBUG */
47    /* #define ATOL_DEBUG */
48    
49  /*------------------------------------------------------------------------------  /*------------------------------------------------------------------------------
50     The following names are of solver_var children or attributes     The following names are of solver_var children or attributes
# Line 1171  void integrator_set_ydot(IntegratorSyste Line 1172  void integrator_set_ydot(IntegratorSyste
1172  */  */
1173  double *integrator_get_atol(IntegratorSystem *sys, double *atol){  double *integrator_get_atol(IntegratorSystem *sys, double *atol){
1174      long i;      long i;
1175    #ifdef ATOL_DEBUG
1176      char *varname;      char *varname;
1177    #endif
1178    
1179      if (atol==NULL) {      if (atol==NULL) {
1180          atol = ASC_NEW_ARRAY_CLEAR(double, sys->n_y);          atol = ASC_NEW_ARRAY_CLEAR(double, sys->n_y);
# Line 1181  double *integrator_get_atol(IntegratorSy Line 1184  double *integrator_get_atol(IntegratorSy
1184          asc_assert(sys->y[i]!=NULL);          asc_assert(sys->y[i]!=NULL);
1185          atol[i] = var_odeatol(sys->y[i]);          atol[i] = var_odeatol(sys->y[i]);
1186          asc_assert(atol[i]!=-1);          asc_assert(atol[i]!=-1);
1187    #ifdef ATOL_DEBUG
1188          varname = var_make_name(sys->system,sys->y[i]);          varname = var_make_name(sys->system,sys->y[i]);
1189          CONSOLE_DEBUG("%s.ode_atol = %8.2e",varname,atol[i]);          CONSOLE_DEBUG("%s.ode_atol = %8.2e",varname,atol[i]);
1190          ASC_FREE(varname);          ASC_FREE(varname);
1191    #endif
1192      }      }
1193      return atol;      return atol;
1194  }  }

Legend:
Removed from v.1386  
changed lines
  Added in v.1387

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