MatShellSetOperation

Synopsis

#include "petscmat.h" 
PetscErrorCode MatShellSetOperation(Mat mat,MatOperation op,void (*f)(void))

Input Parameters

Usage

     extern PetscErrorCode usermult(Mat,Vec,Vec);

     ierr = MatCreateShell(comm,m,n,M,N,ctx,&A);

     ierr = MatShellSetOperation(A,MATOP_MULT,(void(*)(void))usermult);

Notes

All user-provided functions (except for MATOP_DESTROY) should have the same calling sequence as the usual matrix interface routines, since they are intended to be accessed via the usual matrix interface routines, e.g.,

      MatMult(Mat,Vec,Vec) -> usermult(Mat,Vec,Vec)

In particular each function MUST return an error code of 0 on success and nonzero on failure.

Within each user-defined routine, the user should call MatShellGetContext() to obtain the user-defined context that was set by MatCreateShell().

Fortran Notes

Use MatSetOperation() to set an operation for any matrix type

Keywords

See Also

Level

Location

Examples