/[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 777 by johnpye, Tue Jul 18 04:40:08 2006 UTC revision 815 by johnpye, Fri Aug 11 04:51:27 2006 UTC
# Line 1072  double integrator_get_t(IntegratorSystem Line 1072  double integrator_get_t(IntegratorSystem
1072  */  */
1073  void integrator_set_t(IntegratorSystem *blsys, double value){  void integrator_set_t(IntegratorSystem *blsys, double value){
1074    var_set_value(blsys->x, value);    var_set_value(blsys->x, value);
1075    CONSOLE_DEBUG("set_t = %g", value);    /* CONSOLE_DEBUG("set_t = %g", value); */
1076  }  }
1077    
1078  /*------------------------------------------------------------------------------  /*------------------------------------------------------------------------------
# Line 1096  double *integrator_get_y(IntegratorSyste Line 1096  double *integrator_get_y(IntegratorSyste
1096    for (i=0; i< blsys->n_y; i++) {    for (i=0; i< blsys->n_y; i++) {
1097      assert(blsys->y[i]!=NULL);      assert(blsys->y[i]!=NULL);
1098      y[i] = var_value(blsys->y[i]);      y[i] = var_value(blsys->y[i]);
1099      CONSOLE_DEBUG("ASCEND --> y[%ld] = %g", i+1, y[i]);      /* CONSOLE_DEBUG("ASCEND --> y[%ld] = %g", i+1, y[i]); */
1100    }    }
1101    return y;    return y;
1102  }  }
# Line 1117  void integrator_set_y(IntegratorSystem * Line 1117  void integrator_set_y(IntegratorSystem *
1117      var_set_value(blsys->y[i],y[i]);      var_set_value(blsys->y[i],y[i]);
1118  #ifndef NDEBUG  #ifndef NDEBUG
1119      varname = var_make_name(blsys->system, blsys->y[i]);      varname = var_make_name(blsys->system, blsys->y[i]);
1120      CONSOLE_DEBUG("y[%ld] = \"%s\" = %g --> ASCEND", i+1, varname, y[i]);      /* CONSOLE_DEBUG("y[%ld] = \"%s\" = %g --> ASCEND", i+1, varname, y[i]); */
1121      ASC_FREE(varname);      ASC_FREE(varname);
1122  #endif  #endif
1123    }    }
# Line 1142  double *integrator_get_ydot(IntegratorSy Line 1142  double *integrator_get_ydot(IntegratorSy
1142      if(blsys->ydot[i]!=NULL){      if(blsys->ydot[i]!=NULL){
1143          dydx[i] = var_value(blsys->ydot[i]);          dydx[i] = var_value(blsys->ydot[i]);
1144      }      }
1145      CONSOLE_DEBUG("ASCEND --> ydot[%ld] = %g", i+1, dydx[i]);      /* CONSOLE_DEBUG("ASCEND --> ydot[%ld] = %g", i+1, dydx[i]); */
1146    }    }
1147    return dydx;    return dydx;
1148  }  }

Legend:
Removed from v.777  
changed lines
  Added in v.815

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