From: Harshad S. <hsa...@pu...> - 2015-06-08 21:51:20
|
Hi Paul, Following up on our discussion on FD Jacobian through PetscNonlinearSolver. It is not the coloring you need to generate, just the nonzero structure. >> (From that PETSc computes the coloring) so do as I suggest in my other >> email. >> > FYI, this information is contained in the libMesh::CouplingMatrix which > you can get from the libMesh::DofMap. > I tried to use CouplingMatrix in DofMap for building a "wrong" Jacobian just for the non-zero structure. Apparently the DofMap class contains a NULL pointer to the Coupling matrix. Do I need to call any function to initialize that? Thanks, Harshad On Fri, Jun 5, 2015 at 2:49 PM, Paul T. Bauman <ptb...@gm...> wrote: > > > On Fri, Jun 5, 2015 at 2:42 PM, Harshad Sahasrabudhe <hsa...@pu...> > wrote: > >> It is not the coloring you need to generate, just the nonzero >>> structure. (From that PETSc computes the coloring) so do as I suggest in my >>> other email. >> >> >> Thanks, I'll try that out. >> > > FYI, this information is contained in the libMesh::CouplingMatrix which > you can get from the libMesh::DofMap. > |