Function opensn::CreateSquareMatrix(int64_t, int64_t)

Defined in File petsc_utils.h

Function Documentation

Mat opensn::CreateSquareMatrix(int64_t local_size, int64_t global_size)

Creates a general square matrix.

This is a function for:

Mat A;
MatCreate(opensn::mpi_comm,&A);
MatSetType(A,MATMPIAIJ);
MatSetSizes(A,local_size, local_size,
              global_size, global_size);

MatMPIAIJSetPreallocation(A,1, nullptr,
                          0, nullptr);
MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);
MatSetOption(A, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE);

return A;