PETSc version 3.15.5

DMCreateMatrix

Gets empty Jacobian for a DM

Synopsis

#include "petscdm.h"          
#include "petscdmlabel.h"     
#include "petscds.h"     
PetscErrorCode  DMCreateMatrix(DM dm,Mat *mat)

Collective on dm

Input Parameter

- the DM object

Output Parameter

mat

- the empty Jacobian

Options Database Keys

-dm_preallocate_only

- Only preallocate the matrix for DMCreateMatrix(), but do not fill it with zeros

Notes

This properly preallocates the number of nonzeros in the sparse matrix so you do not need to do it yourself.

By default it also sets the nonzero structure and puts in the zero entries. To prevent setting the nonzero pattern call DMSetMatrixPreallocateOnly()

For structured grid problems, when you call MatView() on this matrix it is displayed using the global natural ordering, NOT in the ordering used internally by PETSc.

For structured grid problems, in general it is easiest to use MatSetValuesStencil() or MatSetValuesLocal() to put values into the matrix because MatSetValues() requires the indices for the global numbering for DMDAs which is complicated.

Level

beginner

DMCreateMatrix

Synopsis

Input Parameter

Output Parameter

Options Database Keys

Notes

See Also

Level

Location

Examples

Implementations