[petsc-dev] Jed!!!

[Jed Brown]

Barry Smith <bsmith at mcs.anl.gov> writes:
> cmake/Findcutil.cmake:  add_subdirectory("${_moddir}/nested_targets/cutil")
> FEM/rf_pcs_TH.cpp:              ierr = MatCreateNest(PETSC_COMM_WORLD, dof, NULL, dof, NULL, &eqs_new->vec_subA[0], &eqs_new->A);CHKERRCONTINUE(ierr);
> FEM/rf_pcs_TH.cpp:              ierr = MatNestGetISs(eqs_new->A, &eqs_new->vec_isg[0], NULL);CHKERRCONTINUE(ierr);
> FEM/rf_pcs_TH.cpp:              MatNestGetSubMats(mat,NULL,NULL,&mats);

This is a "developer" function.  They have already shot themselves in
the foot because they can't build a monolithic matrix (so can't use LU,
SVD, etc).  Dmitry wrote MatConvert_Nest_AIJ, but I don't recall if it
was ever finished.  We can finish that, but I really mean it when I tell
people NOT to use MatNestGetSubMats.  The whole point of MatNest (versus
Daves' earlier generation) was to avoid writing assembly in that way.

> FEM/rf_pcs_TH.cpp~:             ierr = MatCreateNest(PETSC_COMM_WORLD, dof, NULL, dof, NULL, &eqs_new->vec_subA[0], &eqs_new->A);CHKERRCONTINUE(ierr);
> FEM/rf_pcs_TH.cpp~:             ierr = MatNestGetISs(eqs_new->A, &eqs_new->vec_isg[0], NULL);CHKERRCONTINUE(ierr);
> FEM/rf_pcs_TH.cpp~:             MatNestGetSubMats(*mat,NULL,NULL,&mats);

-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140427/1880a142/attachment.sig>