Changeset 28259

Timestamp:

05/07/24 10:52:50 (12 months ago)

Author:

Mathieu Morlighem

Message:

CHG: implementing sum for IssmMpiVec

File:

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

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

@@ -588 +588 @@
                 /*}}}*/
                 void Pow(doubletype scale_factor){/*{{{*/
-
-                        int i;
-                        for(i=0;i<this->M;i++)this->vector[i]=pow(this->vector[i],scale_factor);
-
+                        for(int i=0;i<this->M;i++)this->vector[i]=pow(this->vector[i],scale_factor);
                 }
                 /*}}}*/
                 void Sum(doubletype* pvalue){/*{{{*/
-                        _error_("not support yet!");
+
+                        doubletype local_sum=0;
+                        doubletype sum;
+
+                        for(int i=0;i<this->m;i++) local_value+=this->vector[i];
+                        ISSM_MPI_Reduce(&local_sum, &sum, 1, ISSM_MPI_DOUBLE, ISSM_MPI_SUM, 0, IssmComm::GetComm());
+                        ISSM_MPI_Bcast(&sum,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm());
+                        *pvalue=sum;
+
                 }
                 /*}}}*/