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);CHKE
Shit, I may need to do interlaced but let me try. I added
ierr = PetscObjectComposeFunction((PetscObject)A,"MatConvert_Nest_AIJ_C",
MatNestSetSubMats_Nest);CHKERRQ(ierr);
But I still don't find it. I added print statements for what it is looking
for. I will keep trying.
Thanks,
[3] * ma
On Tue, Apr 14, 2015 at 9:59 PM, Mark Adams 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
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:
21:54 next ~/Codes/petsc$ git grep MatConvert_N
OK, I use:
call
MatConvert(a_ts%FJacobian,MATAIJ,MAT_INITIAL_MATRIX,a_ts%FJacobian2,ierr);CHKERRQ(ierr)
I added code to src/snes/examples/tutorials/ex73f90t.F90 that reproduces
this error, in next.
I'll look at it.
Mark
On Tue, Apr 14, 2015 at 6:33 PM, Matthew Knepley wrote:
> On Tue,
On Tue, Apr 14, 2015 at 5:19 PM, Mark Adams wrote:
> Matt, in looking at MatConvert_Nest_AIJ, the code looks like it is 2 years
> old.
>
> I need a MatConvert_Nest_MPIAIJ, should I write this or does it exist
> someplace?
>
This is MPIAIJ. Look at the code.
Matt
>
> Mark
>
>
>
> On Tue, Apr
Matt, in looking at MatConvert_Nest_AIJ, the code looks like it is 2 years
old.
I need a MatConvert_Nest_MPIAIJ, should I write this or does it exist
someplace?
Mark
On Tue, Apr 14, 2015 at 2:54 PM, Mark Adams wrote:
> Also, this is with MPI so I need MPIAIJ. Is that supposed to work?
>
> O
Also, this is with MPI so I need MPIAIJ. Is that supposed to work?
On Tue, Apr 14, 2015 at 2:53 PM, Mark Adams wrote:
> It looks like MatConvert is looking for "MatConvert_nest_aij_C"
>
> but I see:
>
> > git grep MatConvert_Nest
> impls/nest/matnest.c:#define __FUNCT__ "MatConvert_Nest_AIJ"
>
It looks like MatConvert is looking for "MatConvert_nest_aij_C"
but I see:
> git grep MatConvert_Nest
impls/nest/matnest.c:#define __FUNCT__ "MatConvert_Nest_AIJ"
impls/nest/matnest.c:PETSC_EXTERN PetscErrorCode MatConvert_Nest_AIJ(Mat
A,MatType newtype,MatReuse reuse,Mat *newmat)
Is there a pro
I have tried
call
MatConvert(a_ts%FJacobian,MATMPIAIJ,MAT_INITIAL_MATRIX,a_ts%FJacobian2,ierr);CHKERRQ(ierr)
and
call
MatConvert(a_ts%FJacobian,MATAIJ,MAT_INITIAL_MATRIX,a_ts%FJacobian2,ierr);CHKERRQ(ierr)
in 'next' and it is not finding the specialized routine.
I will look into it.
Mark
10 matches
Mail list logo
