[petsc-dev] More Fortran Interfaces

[John O'Sullivan]

Hi Matt,

Merry Christmas!

I was wondering if you had a chance to put in the Fortran interfaces? No worries if you haven’t. We’re breaking up from work on Tuesday and won’t be back until the 7th of Jan so I won’t be working on it for a couple of weeks. It would be great to be able to get back to working on out Petsc code once we’re back though. Is there anyway I can tell that I need to pull a new dev version or should I just wait for you to email me?

Anyway I hope you have a good break and its not too freezing where you are. We’ll be heading to the beach down here.

Cheers
John

--

Dr John O'Sullivan
Lecturer
Department of Engineering Science
University of Auckland, New Zealand
email: jp.osullivan at auckland.ac.nz<https://lists.mcs.anl.gov/mailman/listinfo/petsc-dev>
tel: +64 (0)9 923 85353



From: Matthew Knepley [mailto:knepley at gmail.com]
Sent: Sunday, 14 December 2014 4:52 p.m.
To: John O'Sullivan
Cc: petsc-dev at mcs.anl.gov
Subject: Re: [petsc-dev] More Fortran Interfaces

On Mon, Dec 8, 2014 at 4:00 AM, John O'Sullivan <jp.osullivan at auckland.ac.nz<mailto:jp.osullivan at auckland.ac.nz>> wrote:
Hi,

I'm working with Adrian on our PetSc project at the University of Auckland. I think I've found a couple more fortran interfaces that are missing:

DMPlexGetFaceFields
DMPlexGefFaceGeometry

I get the error:

Undefined symbols for architecture x86_64:
  "_dmplexgetfacefields_", referenced from:
      _rhs_darcy_local.2542 in darcy.o
  "_dmplexgetfacegeometry_", referenced from:
      _rhs_darcy_local.2542 in darcy.o

I haven't got to the restore stage yet but I think probably these ones are missing too:

DMPlexRestoreFaceFields
DMPlexRestoreFaceGeometry

Also I'm having trouble declaring a pointer to a PetscFVFaceGeom type to use with the DMPlexGetFaceGeometry. Is there a trick to this or maybe its definition in fortran is missing too?

Sorry, this is a really bad time of year for me with family and work travel. I can put in the Fortran bindings this week I think.
I do not know how to handle structs in Fortran (will ask Fortran gurus). Right now, I am just going to make a PetscScalar
array with one more dimension since the struct is a bunch of Scalars, and you will have to index into it treating the struct
as another dim.

  Thanks,

    Matt

Thanks very much!

Cheers
John

--

Dr John O'Sullivan

Lecturer

Department of Engineering Science

University of Auckland, New Zealand

email: jp.osullivan at auckland.ac.nz<https://lists.mcs.anl.gov/mailman/listinfo/petsc-dev>

tel: +64 (0)9 923 85353



--
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/20141221/6d2abb4b/attachment.html>