Changeset 865

Timestamp:

06/09/09 11:46:06 (16 years ago)

Author:

seroussi

Message:

fixed problem with velocity in inputs but not elegant

File:

• issm/trunk/src/c/parallel/transient_core_2d.cpp (modified) (3 diffs)

issm/trunk/src/c/parallel/transient_core_2d.cpp

@@ -29 +29 @@
 
         /*solutions: */
-        int step;
+        int step,i;
         double time;
         
         Vec u_g=NULL;
         double* u_g_serial=NULL;
+        double* u_g_serial3d=NULL;
         Vec p_g=NULL;
         double* p_g_serial=NULL;
@@ -82 +83 @@
         time=0;
         
-        fem_p->parameters->FindParam((void*)&u_g_serial,"u_g");
+        fem_p->parameters->FindParam((void*)&u_g_serial3d,"u_g");
+        u_g_serial=(double*)xmalloc(2*numberofnodes*sizeof(double));
+        for(i=0;i<numberofnodes;i++){
+                u_g_serial[2*i+0]=u_g_serial3d[3*i+0];
+                u_g_serial[2*i+1]=u_g_serial3d[3*i+1];
+        }
         u_g=SerialToVec(u_g_serial,3*numberofnodes); xfree((void**)&u_g_serial);
 
@@ -153 +159 @@
         /*Free ressources:*/
         xfree((void**)&u_g_serial);
+        xfree((void**)&u_g_serial3d);
         xfree((void**)&p_g_serial);
         xfree((void**)&h_g_serial);