|
From: James Lu <jam...@oe...> - 2006-09-06 16:17:15
|
hi Rainer: In looking at the adjoint code again, I saw that IntegratorInstance_reset calls IntegratorInstance_set, which consists of calls: CvodeData_initialize and IntegratorInstance_initializeSolver. I just found that the bug in running examples/adj_sensitivity.c is associated with the IntegratorInstance_set called after the forward phase. In fact, what is need is IntegratorInstance_initializeSolver: CvodeData_initialize shouldn't be called before the adjoint phase starts. Could you tell me circumstances in which IntegratorInstance_reset is used? And is it possible (perhaps after some reorganization) to have access to simply IntegratorInstance_initializeSolver...what is the best way to go about it? thanks, James |