[petsc-dev] DMViewFromOptions() in fortran

[Matthew Knepley]

On Tue, Mar 25, 2014 at 8:52 PM, Adrian Croucher
<a.croucher at auckland.ac.nz>wrote:

>
>
> PETSC_STATIC_INLINE PetscErrorCode DMViewFromOptions(DM A,const char
> prefix[],const char name[]) {return
> PetscObjectViewFromOptions((PetscObject)A,prefix,name);}
>
>  so I think it might be better to just call the PetscObject version from
> Fortran.
>
>
> I tried that but got a similar linker error from
> PetscObjectViewFromOptions().
>
>
> Can you run an example, like SNES ex12, and give -dm_view in parallel?
>
>
> Yes, that works.  The C version of my test code also works- that's why I
> figured it might be just a missing fortran interface. My fortran test code
> is below.
>

Below, you should not need an if, and can just use DMView(). If that does
not work in parallel, something is wrong.
Likely on the interpretation of the default viewer. Let me try it.

   Matt


> Cheers, Adrian
>
> --
> program main
>
> ! distributing a DMPlex
>
>   implicit none
>
> #include <finclude/petsc.h90>
>
>   DM :: dm, interp_dm, dist_dm
>   PetscInt, parameter :: dim = 3
>   PetscInt :: depth = 1
>   PetscErrorCode :: ierr
>   PetscInt :: rank, size
>   PetscInt, allocatable :: numPoints(:)
>   PetscInt, allocatable :: coneSize(:)
>   PetscInt, allocatable :: cones(:)
>   PetscInt, allocatable :: coneOrientations(:)
>   PetscScalar, allocatable :: vertexCoords(:)
>   character(5) :: partitioner = "chaco"
>
>   call PetscInitialize(PETSC_NULL_CHARACTER,ierr); CHKERRQ(ierr)
>   call MPI_Comm_rank(PETSC_COMM_WORLD, rank, ierr); CHKERRQ(ierr);
>   call MPI_Comm_size(PETSC_COMM_WORLD, size, ierr); CHKERRQ(ierr);
>
>   call DMPlexCreate(PETSC_COMM_WORLD, dm, ierr); CHKERRQ(ierr)
>   call PetscObjectSetName(dm, "testplex", ierr); CHKERRQ(ierr)
>   call DMPlexSetDimension(dm, dim, ierr); CHKERRQ(ierr)
>
>   if (rank == 0) then
>      numPoints = [12, 2]
>      coneSize  = [8, 8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
>      cones = [2,5,4,3,6,7,8,9,  3,4,11,10,7,12,13,8]
>      coneOrientations = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
>      vertexCoords = [ &
>           -0.5,0.0,0.0, 0.0,0.0,0.0, 0.0,1.0,0.0, -0.5,1.0,0.0, &
>           -0.5,0.0,1.0, 0.0,0.0,1.0, 0.0,1.0,1.0, -0.5,1.0,1.0, &
>           0.5,0.0,0.0, 0.5,1.0,0.0, 0.5,0.0,1.0,  0.5,1.0,1.0 ]
>      call DMPlexCreateFromDAG(dm, depth, numPoints, coneSize, cones, &
>           coneOrientations, vertexCoords, ierr)
>   else
>      numPoints = [0, 0, 0]
>      coneSize = [0]
>      cones = [0]
>      coneOrientations = [0]
>      vertexCoords = [0]
>      call DMPlexCreateFromDAG(dm, 1, numPoints, coneSize, &
>           cones, coneOrientations, vertexCoords, ierr);
>   end if
>   CHKERRQ(ierr)
>
>   call DMPlexInterpolate(dm, interp_dm, ierr); CHKERRQ(ierr)
>   call DMPlexCopyCoordinates(dm, interp_dm, ierr); CHKERRQ(ierr)
>   call DMDestroy(dm, ierr); CHKERRQ(ierr)
>   dm = interp_dm
>
>   call DMPlexDistribute(dm, partitioner, 0, PETSC_NULL_OBJECT, dist_dm,
> ierr)
>   CHKERRQ(ierr)
>   if (size > 1) then
>     call DMViewFromOptions(dist_dm, PETSC_NULL_OBJECT, "-dm_view", ierr)
>     CHKERRQ(ierr)
>   else ! one processor
>     call DMView(dm, PETSC_VIEWER_STDOUT_WORLD, ierr); CHKERRQ(ierr)
>   end if
>   call DMDestroy(dm, ierr); CHKERRQ(ierr)
>   dm = dist_dm
>
>   call DMDestroy(dm, ierr); CHKERRQ(ierr)
>   call PetscFinalize(ierr); CHKERRQ(ierr)
>   deallocate(numPoints, coneSize, cones, coneOrientations, vertexCoords)
>
> end program main
>
>
> --
> Dr Adrian Croucher
> Department of Engineering Science
> University of Auckland
> New Zealand
> tel 64-9-373-7599 ext 84611
>
>


-- 
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/20140325/5d7ccc49/attachment.html>