1. petsc
  2. PETSc
  3. petsc
  4. Issues
Issue #66 new

Fortran support for multicomponent DMDAVecGetArrayF90()

BarryFSmith
created an issue

Add support for using u(i,j)%p,v,x etc in Fortran. this means handling each dof size with custom DMDAVecGetArrayF90() code for each dimension. Also add

type DMDACoor2dType real*8 :: x,y end type DMDACoor2dType

define DMDACoor2d type(DMDACoor2dType)

DMDACoor2d ,pointer :: coor(:,:)
DMDAVecGetArrayF90(da,vec,coor,ierr)

Note that because of the way F90 interfaces work the fieldnames need to the same in the fortran interface file as in the user code! Ok for coordinates but a pain in general!

Comments (0)

  1. Log in to comment