Changeset 13599
- Timestamp:
- 10/10/12 22:39:12 (12 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/OutputRiftsx/OutputRiftsx.cpp
r4214 r13599 17 17 18 18 /*Allocate grad_g: */ 19 riftproperties=NewVec(numrifts );19 riftproperties=NewVec(numrifts,IssmComm::GetComm()); 20 20 21 21 /*Compute rift properties : */ -
issm/trunk-jpl/src/c/modules/SmearFunctionx/SmearFunctionx.cpp
r11259 r13599 28 28 /*allocate: */ 29 29 gsize=nodes->NumberOfDofs(configuration_type,GsetEnum); 30 smearedvector=NewVec(gsize );30 smearedvector=NewVec(gsize,IssmComm::GetComm()); 31 31 32 x=NewVec(gsize );33 y=NewVec(gsize );34 z=NewVec(gsize );32 x=NewVec(gsize,IssmComm::GetComm()); 33 y=NewVec(gsize,IssmComm::GetComm()); 34 z=NewVec(gsize,IssmComm::GetComm()); 35 35 36 36 for(i=0;i<nodes->Size();i++){ -
issm/trunk-jpl/src/c/modules/Solverx/SolverxPetsc.cpp
r13590 r13599 70 70 MatGetSize(Kff,&global_m,&global_n); _assert_(global_m==global_m); 71 71 if(!global_n){ 72 *puf=NewVec(0 ); return;72 *puf=NewVec(0,IssmComm::GetComm()); return; 73 73 } 74 74 … … 86 86 } 87 87 else{ 88 MatGetLocalSize(Kff,&local_m,&local_n);uf=NewVec(local_n, fromlocalsize);88 MatGetLocalSize(Kff,&local_m,&local_n);uf=NewVec(local_n,IssmComm::GetComm(),fromlocalsize); 89 89 } 90 90 -
issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscVec.cpp
r13595 r13599 30 30 PetscVec::PetscVec(int M,bool fromlocalsize){ 31 31 32 this->vector=NewVec(M, fromlocalsize);32 this->vector=NewVec(M,IssmComm::GetComm(),fromlocalsize); 33 33 34 34 } … … 38 38 39 39 if(petsc_vec==NULL){ 40 this->vector=NewVec(0 );40 this->vector=NewVec(0,IssmComm::GetComm()); 41 41 } 42 42 else{ … … 55 55 for(int i=0;i<M;i++) idxm[i]=i; 56 56 57 this->vector=NewVec(M );57 this->vector=NewVec(M,IssmComm::GetComm()); 58 58 VecSetValues(this->vector,M,idxm,serial_vec,INSERT_VALUES); 59 59 VecAssemblyBegin(this->vector); -
issm/trunk-jpl/src/c/toolkits/petsc/patches/NewVec.cpp
r11695 r13599 17 17 #include "../../mpi/patches/mpipatches.h" 18 18 19 Vec NewVec(int size, bool fromlocalsize){19 Vec NewVec(int size,COMM comm,bool fromlocalsize){ 20 20 21 21 int ierr; … … 33 33 } 34 34 35 VecCreate( PETSC_COMM_WORLD,&vector);35 VecCreate(comm,&vector); 36 36 37 37 VecSetSizes(vector,local_size,PETSC_DECIDE); -
issm/trunk-jpl/src/c/toolkits/petsc/patches/SerialToVec.cpp
r12431 r13599 30 30 31 31 /*Create parallel vector: */ 32 outvector=NewVec(vector_size );32 outvector=NewVec(vector_size,MPI_COMM_WORLD); 33 33 34 34 /*plug values from local vector into new parallel vector: */ -
issm/trunk-jpl/src/c/toolkits/petsc/patches/VecPartition.cpp
r12431 r13599 100 100 if(kffpartition){ 101 101 VecFree(&outvector); 102 outvector=NewVec(row_partition_vector_size );102 outvector=NewVec(row_partition_vector_size,MPI_COMM_WORLD); 103 103 } 104 104 -
issm/trunk-jpl/src/c/toolkits/petsc/patches/VecTranspose.cpp
r12431 r13599 25 25 26 26 /*Create new vector of same size: */ 27 tvector=NewVec(size );27 tvector=NewVec(size,MPI_COMM_WORLD); 28 28 29 29 /*Extract values locally from input vector: */ -
issm/trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h
r13597 r13599 17 17 class Parameters; 18 18 19 Vec NewVec(int size, bool fromlocalsize=false);19 Vec NewVec(int size,COMM comm,bool fromlocalsize=false); 20 20 Mat NewMat(int M,int N,COMM comm); 21 21 Mat NewMat(int M,int N,double sparsity,COMM comm);
Note:
See TracChangeset
for help on using the changeset viewer.