Changeset 13545
- Timestamp:
- 10/04/12 18:25:22 (12 years ago)
- Location:
- issm/trunk-jpl/src/c/solutions
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/solutions/ad_core.cpp
r13543 r13545 26 26 27 27 /*diverse: */ 28 int i; 29 int dummy; 30 31 bool isautodiff = false; 32 int num_dependents; 33 int num_independents; 34 char* driver=NULL; 35 size_t tape_stats[11]; 28 int i; 29 int dummy; 30 int num_dependents; 31 int num_independents; 32 bool isautodiff = false; 33 char *driver = NULL; 34 size_t tape_stats[11]; 36 35 37 36 /*state variables: */ 38 IssmDouble *axp 39 double *xp 37 IssmDouble *axp = NULL; 38 double *xp = NULL; 40 39 41 40 /*AD mode on?: */ … … 72 71 73 72 /* these are always needed regardless of the interpreter */ 74 75 73 anEDF_for_solverx_p->dp_x=xNew<double>(anEDF_for_solverx_p->max_n); 74 anEDF_for_solverx_p->dp_y=xNew<double>(anEDF_for_solverx_p->max_m); 76 75 77 76 if (strcmp(driver,"fos_forward")==0){ 78 77 79 78 int anIndepIndex; 80 double *tangentDir =NULL;81 double *jacTimesTangentDir =NULL;82 double *theOutput =NULL;79 double *tangentDir = NULL; 80 double *jacTimesTangentDir = NULL; 81 double *theOutput = NULL; 83 82 84 83 /*retrieve direction index: */ … … 113 112 else if ((strcmp(driver,"fov_forward")==0) || (strcmp(driver,"fov_forward_all")==0)){ 114 113 115 int * indepIndices=NULL;116 int tangentDirNum;117 int dummy;118 double **jacTimesSeed =NULL;119 double **seed =NULL;120 double *theOutput=NULL;114 int tangentDirNum; 115 int dummy; 116 int *indepIndices = NULL; 117 double **jacTimesSeed = NULL; 118 double **seed = NULL; 119 double *theOutput = NULL; 121 120 std::set<unsigned int> anIndexSet; 122 123 121 124 122 /*retrieve directions:*/ … … 302 300 } 303 301 if(VerboseAutodiff())_pprintLine_(" end AD core"); 302 /*}}}*/ 304 303 305 304 /*Free resources: */ -
issm/trunk-jpl/src/c/solutions/issm.cpp
r13543 r13545 75 75 /*call cores: */ 76 76 _pprintLine_("call computational core:"); 77 profiler->Tag(StartCore); solutioncore(femmodel); profiler->Tag(FinishCore); 78 profiler->Tag(StartAdCore); ad_core(femmodel); profiler->Tag(FinishAdCore); 79 77 profiler->Tag(StartCore); solutioncore(femmodel); profiler->Tag(FinishCore); 78 profiler->Tag(StartAdCore); ad_core(femmodel); profiler->Tag(FinishAdCore); 80 79 ProfilerEnd(profiler,femmodel->results,femmodel->parameters); 81 80 -
issm/trunk-jpl/src/c/solutions/solutions.h
r13540 r13545 57 57 void AdjointCorePointerFromSolutionEnum(void (**padjointcore)(FemModel*),int solutiontype); 58 58 59 60 59 #endif
Note:
See TracChangeset
for help on using the changeset viewer.