Changeset 13596


Ignore:
Timestamp:
10/10/12 22:17:11 (12 years ago)
Author:
Eric.Larour
Message:

CHG: fixed segfault in NewMat call

Location:
issm/trunk-jpl/src/c
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/Makefile.am

    r13589 r13596  
    771771                                        ./toolkits/petsc/patches/KSPFree.cpp\
    772772                                        ./toolkits/petsc/patches/MatFree.cpp\
    773                                         ./toolkits/petsc/patches/GetOwnershipBoundariesFromRange.cpp\
    774773                                        ./toolkits/petsc/patches/VecPartition.cpp\
    775774                                        ./toolkits/petsc/patches/MatInvert.cpp\

  • issm/trunk-jpl/src/c/toolkits/petsc/objects/PetscMat.cpp

    r13595 r13596  
    3737PetscMat::PetscMat(int M,int N, IssmDouble sparsity){
    3838
    39         this->matrix=NewMat(M,N,sparsity);
     39        this->matrix=NewMat(M,N,sparsity,IssmComm::GetComm());
    4040}
    4141/*}}}*/
    … …  
    5454
    5555
    56         this->matrix=NewMat(M,N,sparsity);
     56        this->matrix=NewMat(M,N,sparsity,IssmComm::GetComm());
    5757        MatSetValues(this->matrix,M,idxm,N,idxn,serial_mat,INSERT_VALUES);
    5858        MatAssemblyBegin(this->matrix,MAT_FINAL_ASSEMBLY);
    … …  
    6767PetscMat::PetscMat(int M,int N, int connectivity,int numberofdofspernode){
    6868       
    69         this->matrix=NewMat(M,N,connectivity,numberofdofspernode);
     69        this->matrix=NewMat(M,N,connectivity,numberofdofspernode,IssmComm::GetComm());
    7070
    7171}

  • issm/trunk-jpl/src/c/toolkits/petsc/patches/NewMat.cpp

    r13595 r13596  
    4848}
    4949/*}}}*/
    50 /*NewMat(int M,int N,double sparsity){{{*/
     50/*NewMat(int M,int N,double sparsity,COMM comm){{{*/
    5151Mat NewMat(int M,int N,double sparsity,COMM comm){
    5252

  • issm/trunk-jpl/src/c/toolkits/petsc/patches/petscpatches.h

    r13595 r13596  
    2424int VecTranspose(Vec* tvector,Vec vector);
    2525int VecToMPISerial(double** pgathered_vector, Vec vector,COMM comm);
    26 void GetOwnershipBoundariesFromRange(int* plower_row,int* pupper_row,int range,COMM comm);
    2726void MatFree(Mat* pmat);
    2827void ISFree(IS* pis);
Note: See TracChangeset for help on using the changeset viewer.