Changeset 23606

Timestamp:

01/06/19 09:38:48 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added GetLocalVector to IssmSeqVec

File:

• issm/trunk-jpl/src/c/toolkits/issm/IssmSeqVec.h (modified) (1 diff)

issm/trunk-jpl/src/c/toolkits/issm/IssmSeqVec.h

@@ -144 +144 @@
                 }/*}}}*/
                 void GetLocalVector(doubletype** pvector,int** pindices){/*{{{*/
-                        _error_("not implemented");
+
+                        /*First, check that vector size is not 0*/
+                        int vector_size;
+                        this->GetSize(&vector_size);
+                        if(vector_size==0){
+                                *pvector=NULL;
+                                *pindices=NULL;
+                                return;
+                        }
+
+                        /*Build indices*/
+                        int* indices=xNew<int>(vector_size); 
+                        for(int i=0;i<vector_size;i++) indices[i]=i;
+
+                        /*Get vector*/
+                        doubletype* values =xNew<doubletype>(vector_size);
+                        xMemCpy<doubletype>(values,this->vector,vector_size);
+
+                        *pvector  = values;
+                        *pindices = indices;
                 } /*}}}*/
                 IssmSeqVec<doubletype>* Duplicate(void){/*{{{*/