Changeset 322 for issm/trunk/src/c/objects/Beam.cpp

Timestamp:

05/08/09 14:41:55 (16 years ago)

Author:

seroussi

Message:

fixed problems with Beam element

File:

• issm/trunk/src/c/objects/Beam.cpp (modified) (2 diffs)

issm/trunk/src/c/objects/Beam.cpp

@@ -362 +362 @@
                 GetJacobianDeterminant(&Jdet, &z_list[0],gauss_coord);
                 
-                //Get nodal functions
-                GetNodalFunctions(l1l2, gauss_coord);
-
                 for(j=0;j<NDOF2;j++){
                         pe_g_gaussian[NDOF2+j]=constant_part*pow((s[0]-z_g)/B,n)*slope[j]*Jdet*gauss_weight;
                 }
-
+                
                 //add pe_gaussian vector to pe:
                 for(j=0;j<numdofs;j++){
@@ -387 +384 @@
                 pe_g[1]+=vb;
         }
+
+        /*Add pe_g to global vector pg: */
+        VecSetValues(pg,numdofs,doflist,(const double*)pe_g,ADD_VALUES);
+
+        cleanup_and_return: 
+        xfree((void**)&segment_gauss_coord);
+        xfree((void**)&gauss_weights);
 }