DA Jacobian

Jianing Shi js2615 at columbia.edu
Sat Dec 9 01:03:33 CST 2006


I have a question about creating the jacobian here.
Here is the code I use to create it.

-------------------------------------------------------
  ierr = DAGetMatrix(user.da,MATAIJ,&J);CHKERRQ(ierr);
  A = J;

  ierr =
DASetLocalAdicMFFunction(user.da,admf_FormFunctionLocal);CHKERRQ(ierr);
  ierr = MatRegisterDAAD();CHKERRQ(ierr);
  ierr = MatCreateDAAD(user.da,&A);CHKERRQ(ierr);
  ierr = MatDAADSetSNES(A,snes);CHKERRQ(ierr);
  ierr = MatDAADSetCtx(A,&user);CHKERRQ(ierr);
------------------------------------------------------

I included the head files
petscsnes.h    petscda.h

When I compile it, I got the following linking error:

ld: Undefined symbols:
_MatCreateDAAD
_MatDAADSetCtx
_MatDAADSetSNES
_MatRegisterDAAD
make: [paraImplicitV1] Error 1 (ignored)

Not sure what went wrong.

Jianing




More information about the petsc-users mailing list