Changeset 11737 for issm/trunk-jpl/src/c/modules

Timestamp:

03/19/12 20:53:40 (13 years ago)

Author:

Eric.Larour

Message:

Vector object does not have an M size.

File:

• issm/trunk-jpl/src/c/modules/Solverx/Solverx.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/modules/Solverx/Solverx.cpp

@@ -20 +20 @@
 
         /*output: */
-        Vector *uf=new Vector();
+        Vector *uf=NULL;
 
         /*In debugging mode, check that stiffness matrix and load vectors are not NULL (they can be empty)*/
@@ -29 +29 @@
                 Vec uf0_vector = NULL;
                 Vec df_vector  = NULL;
+                Vec uf_vector  = NULL;
                 if(uf0) uf0_vector = uf0->vector;
                 if(df)  df_vector  = df->vector;
@@ -38 +39 @@
                 #endif
 
-                SolverxPetsc(&uf->vector,Kff->matrix,pf->vector,uf0_vector,df_vector,parameters);
-                if(uf->vector == NULL){
-                        uf->M = 0;
-                }
-                else{
-                        VecGetSize(uf->vector,&uf->M);
-                }
+                SolverxPetsc(&uf_vector,Kff->matrix,pf->vector,uf0_vector,df_vector,parameters);
+
+                /*Create vector out of petsc vector: */
+                uf=new Vector(uf_vector);
         #else
                 #ifdef _HAVE_GSL_