[petsc-dev] Question: Regarding MatGetRowIJF90()

[Barry Smith]

  We'll need more information about the crash or sample code the reproduces the problem. petsc-maint at mcs.anl.gov

   Note this stuff will not work with older version of the gfortran compiler so make sure you are using the latest version of gfortran if you are using it.

    Barry


On Jul 22, 2011, at 10:38 AM, Bisht, Gautam wrote:

> Hi,
> 
> I have a sparse MPI Matrix and am trying to access the values for a FORTRAN subroutine.
> 
> 1 Create the matrix: MatCreateMPIAIJ( ...., A)
> 2 Get diagonal and off-diagonal data: Get the MatMPIAIJGetSeqAIJ (A, Ad, Ao, ...)
> 3 MatGetRowIJ( Ad, ... )
> 4 MatGetArray( Ad, ... )
> 
> This approach works.
> 
> But when I try to use MatGetRowIJF90( Ad, ... ), the code crashes (error related to F90Array1dCreate ). Any ideas why this could be happening?
> 
> I have seen that MatGetRowIJF90() works successfully on a matrix created by MatPartitioningSetAdjacency().
> 
> I would appreciate any input on this.
> 
> Thanks,
> -Gautam.
> 
>