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