On Thu, Jun 1, 2017 at 2:59 PM, Justin Pogacnik <j.pogacnik@auckland.ac.nz> wrote:

​All good. Thanks Matt. Will keep an eye out for that update. :)

I have checked, and I was wrong before. I did write the code for PetscQuadratureGet/RestoreData()
in Fortran. Since it uses array arguments, I used F90. Thus you have to use F90 pointers, in the
same style as DMPlexGet/RestoreCone() in this example:

    https://bitbucket.org/petsc/petsc/src/a19fbe4d52f99f875359274419a2d40a87edfba3/src/dm/impls/plex/examples/tutorials/ex1f90.F90?at=master&fileviewer=file-view-default

Let me know if that is not understandable

  Thanks,

    Matt

-Justin

From: Matthew Knepley <knepley@gmail.com>
Sent: Friday, June 2, 2017 12:52 AM

To: Justin Pogacnik
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PetscFECreateDefault in Fortran
 
On Wed, May 31, 2017 at 10:00 PM, Justin Pogacnik <j.pogacnik@auckland.ac.nz> wrote:

Thanks Matt! That works perfectly now. I have another question regarding accessing the quadrature information.


When I use PetscFEGetQuadrature(), then PetscQuadratureView(), I see what I expect regarding point locations, weights.


However, when I try to use PetscQuadratureGetData() the pointers seem to point to random memory locations. 


The exact line from my test problem is: call PetscQuadratureGetData(quad,q_nc,q_dim,q_num,pq_points,pq_weights,ierr);

where the pq_* are the pointers giving strange output. The q_nc, q_dim, and q_num are all giving what I would expect to see. 

You are clearly the first Fortran user interested in this stuff ;) Handling of arrays in Fortran demands some more work
from us. I need to write a wrapper for that function. I will do it as soon as I can.

  Thanks,

    Matt 

Happy to send along the file if that helps.


Thanks again,


Justin

From: Matthew Knepley <knepley@gmail.com>
Sent: Thursday, June 1, 2017 1:34 AM
To: Justin Pogacnik
Cc: petsc-users@mcs.anl.gov
Subject: Re: [petsc-users] PetscFECreateDefault in Fortran
 
On Wed, May 31, 2017 at 7:53 AM, Matthew Knepley <knepley@gmail.com> wrote:
On Tue, May 30, 2017 at 10:19 PM, Justin Pogacnik <j.pogacnik@auckland.ac.nz> wrote:

Hello,

I'm developing a finite element code in fortran 90. I recently updated my PETSc and am now getting the following error during compile/linking on an existing application:

Undefined symbols for architecture x86_64:

  "_petscfecreatedefault_", referenced from:

      _MAIN__ in fe_test.o

ld: symbol(s) not found for architecture x86_64

collect2: error: ld returned 1 exit status

make: *** [dist/fe_test] Error 1


I'm running Mac OS X Yosemite (10.10.5). I've created a "minimum working example" (attached) that re-creates the problem. It's basically just dm/impls/plex/examples/tutorials/ex3f90, but tries to create a PetscFE object. Everything goes fine and the DM looks like what is expected if PetscFECreateDefault is commented out. Any idea what am I missing?
Yes, I had not made a Fortran binding for this function. I will do it now.

I have merged it to the 'next' branch, and it will be in 'master' soon.

  Thanks,

     Matt
 
  Thanks,

     Matt
 
Many thanks!

Justin







--
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




--
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




--
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




--
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

http://www.caam.rice.edu/~mk51/