[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: DA Jacobian



Are you really using ADIC?

 Matt

On 12/9/06, Jianing Shi <js2615@xxxxxxxxxxxx> wrote:
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




--
"Failure has a thousand explanations. Success doesn't need one" -- Sir
Alec Guiness