Changeset 5772 for issm/trunk/src/c/solvers/solver_adjoint_linear.cpp

Timestamp:

09/13/10 13:18:07 (15 years ago)

Author:

Eric.Larour

Message:

New solution strategy, based on creation of Kff, Kfs directly,
instead of Kgg. Idem for pf instead of pg.
Kept backwards compatibility for now, and only applied to macayeal diagnostic
in 2d. Will switch to this new strategy in the near future.

File:

• issm/trunk/src/c/solvers/solver_adjoint_linear.cpp (modified) (1 diff)

issm/trunk/src/c/solvers/solver_adjoint_linear.cpp

@@ -17 +17 @@
         Vec pg  = NULL, pf=NULL;
 
-        SystemMatricesx(&Kgg, &pg,NULL,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
+        SystemMatricesx(&Kgg, NULL, NULL, &pg,NULL, NULL,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
 
         Reducematrixfromgtofx(&Kff,&Kfs,Kgg,femmodel->nodesets,femmodel->parameters); MatFree(&Kgg);