Changeset 11001 for issm/trunk-jpl/src/c/objects/ElementResults/TriaVertexElementResult.cpp

Timestamp:

12/01/11 10:14:29 (13 years ago)

Author:

Mathieu Morlighem

Message:

Results on vertices are now averaged using node connectivity as in PatchToVec

File:

• issm/trunk-jpl/src/c/objects/ElementResults/TriaVertexElementResult.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/objects/ElementResults/TriaVertexElementResult.cpp

@@ -184 +184 @@
 /*}}}*/
 /*FUNCTION TriaVertexElementResult::GetVectorFromResults{{{1*/
-void TriaVertexElementResult::GetVectorFromResults(Vec vector,int* doflist,int numdofs){
+void TriaVertexElementResult::GetVectorFromResults(Vec vector,int* doflist,int* connectivitylist,int numdofs){
+
+        double data[3];
 
         if(numdofs!=3)_error_("Result %s is a TriaVertexElementResult and cannot write vector of %i dofs",numdofs);
-        VecSetValues(vector,numdofs,doflist,(const double*)this->values,INSERT_VALUES);
+        for(int i=0;i<3;i++) data[i]=this->values[i]/(double)connectivitylist[i];
+        VecSetValues(vector,numdofs,doflist,(const double*)&data,ADD_VALUES);
 
 } /*}}}*/