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

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

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

revision 1067 by johnpye, Sun Jan 7 07:54:06 2007 UTC revision 1068 by johnpye, Mon Jan 8 04:45:47 2007 UTC
# Line 594  int integrator_analyse_dae(IntegratorSys Line 594  int integrator_analyse_dae(IntegratorSys
594          if(info->derivative){          if(info->derivative){
595              sys->y[yindex] = info->i;              sys->y[yindex] = info->i;
596              var_set_deriv(info->i,FALSE);              var_set_deriv(info->i,FALSE);
             var_set_diff(info->i,TRUE);  
597    
598              sys->ydot[yindex] = info->derivative->i;              sys->ydot[yindex] = info->derivative->i;
599              var_set_deriv(info->derivative->i,TRUE);              var_set_deriv(info->derivative->i,TRUE);
             var_set_diff(info->derivative->i,FALSE);  
600    
601              if(info->varindx >= 0){              if(info->varindx >= 0){
602                  ASC_ASSERT_RANGE(yindex, -1e7L, 1e7L);                  ASC_ASSERT_RANGE(yindex, -1e7L, 1e7L);
# Line 617  int integrator_analyse_dae(IntegratorSys Line 615  int integrator_analyse_dae(IntegratorSys
615          }else{          }else{
616              sys->y[yindex] = info ->i;              sys->y[yindex] = info ->i;
617              var_set_deriv(info->i,FALSE);              var_set_deriv(info->i,FALSE);
             var_set_diff(info->i,FALSE);  
618    
619              sys->ydot[yindex] = NULL;              sys->ydot[yindex] = NULL;
620              /* in this case, ydot is not present in the system so we can't set flags on it! */              /* in this case, ydot is not present in the system so we can't set flags on it! */

Legend:
Removed from v.1067  
changed lines
  Added in v.1068

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