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

tel: +64 (0)9 923 85353

 

 

 

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

 

On Mon, Dec 8, 2014 at 4:00 AM, John O'Sullivan <jp.osullivan@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
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