DA Jacobian

Matthew Knepley knepley at gmail.com
Sat Dec 9 08:19:57 CST 2006


Are you really using ADIC?

  Matt

On 12/9/06, Jianing Shi <js2615 at columbia.edu> 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




More information about the petsc-users mailing list