[petsc-dev] EM solver in XGC

[Mark Adams]

I inserted this and get this error.  I will push this to master as it is
trivial:

    ierr =
PetscObjectComposeFunction((PetscObject)A,"MatNestSetSubMats_C",
MatNestSetSubMats_Nest);CHKERRQ(ierr);
+  ierr =
PetscObjectComposeFunction((PetscObject)A,"MatConvert_nest_aij_C",MatConvert_Nest_AIJ);CHKERRQ(ierr);

I will work on this now.  Any suggestions welcome.

Mark

[3]PETSC ERROR: [11]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[11]PETSC ERROR: Argument out of range
[11]PETSC ERROR: Only local rows
[11]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
for trouble shooting.
[11]PETSC ERROR: Petsc Development GIT revision: v3.5.3-3422-g47edbfd  GIT
Date: 2015-04-14 16:28:49 -0700
[11]PETSC ERROR: ../../XGC1_3/xgc-em on a arch-macosx-gnu-g named
MacBook-Pro.local by markadams Wed Apr 15 09:37:54 2015
[11]PETSC ERROR: Configure options --download-hypre=1 --download-p4est
--download-superlu=1 --download-superlu_dist=1 --download-parmetis=1
--download-metis=1 --download-triangle=1
--with-hdf5-dir=/Users/markadams/Codes/hdf5 --with-x=0 --with-debugging=1
PETSC_ARCH=arch-macosx-gnu-g
[11]PETSC ERROR: #1 MatGetRow_MPIAIJ() line 1811 in
/Users/markadams/Codes/petsc/src/mat/impls/aij/mpi/mpiaij.c
[11]PETSC ERROR: #2 MatGetRow() line 425 in
/Users/markadams/Codes/petsc/src/mat/interface/matrix.c
[11]PETSC ERROR: #3 MatConvert_Nest_AIJ() line 1461 in
/Users/markadams/Codes/petsc/src/mat/impls/nest/matnest.c


On Tue, Apr 14, 2015 at 11:05 PM, Matthew Knepley <knepley at gmail.com> wrote:

> On Tue, Apr 14, 2015 at 9:59 PM, Mark Adams <mfadams at lbl.gov> wrote:
>
>> I am trying to see why MatConvert is failing for MatNest to MPIAIJ and
>> added print statements in to see what it is checking for.
>>
>> Where is the function list created for MatNest?  I am grepping on
>> MatConvert_[n/N] and not finding anything other than:
>>
>
> The person who added that code declined to hook it into the function
> table, set at the bottom of matnest.c. This does not bode well.
>
> MatNest should never ever ever be used as a dependency. Its an
> optimization. At least this is my view.
>
>    Matt
>
>
>> 21:54 next ~/Codes/petsc$ git grep MatConvert_N
>> src/mat/impls/nest/matnest.c:#define __FUNCT__ "MatConvert_Nest_AIJ"
>> src/mat/impls/nest/matnest.c:PETSC_EXTERN PetscErrorCode
>> MatConvert_Nest_AIJ(Mat A,MatType newtype,MatReuse reuse,Mat *newmat)
>> 21:57 next ~/Codes/petsc$ git grep MatConvert_n
>> 21:57 1 next ~/Codes/petsc$
>>
>> Mark
>>
>>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150415/6c06c755/attachment.html>