[petsc-users] TAO setup with modules in Fortran 90

[Barry Smith]

> On Oct 6, 2017, at 11:18 PM, Randy Michael Churchill <rchurchi at pppl.gov> wrote:
> 
> So if I'm limited to petsc 3.7.6 for reasons of eventually using within an existing, larger codebase that depends on 3.7.6, is it possible to use TAO with a user-defined module in Fortran90 using 3.7.6?

  You should really pus this "existing, larger codebase" to transition to PETSc 3.8 sooner, rather than later. Especially for developments in Fortran it will make life better for everyone. We are always willing to help users, once they have read the changes information, with information to make transitioning to the latest PETSc version easy. For any parts of the code in C, transitioning from 3.7 to 3.8 should be really simple.

   Barry

> 
> I had tried the various forms of includes listed in the documentation, e.g. see below. I think I now realize this is an issue with the petsc installation on Edison, it does not seem to have the petsctao module in the library file (confirmed using nm -D on the library file). If I do the same include and use statement but with, for example, petscmat, it compiles fine. 
> 
> I built v3.8 from source, and the petsctao module is in the library file, and now the make works. 
> 
> commondat.F90
> module commondat
> #include <petsc/finclude/petscdef.h>
>       use petsc
>       PetscReal    ::    alpha
>       PetscInt     ::    n
> end module commondat
> 
> program rosenbrock1f
> #include <petsc/finclude/petsctaodef.h>
>       use petsctao
>       use commondat
> 
> 
> 
> 
> On Fri, Oct 6, 2017 at 7:54 AM, Matthew Knepley <knepley at gmail.com> wrote:
> On Fri, Oct 6, 2017 at 7:36 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> 
>    Randy,
> 
>     First you absolutely must use version 3.8 or the master development copy. We improved and simplified dramatically how Fortran (90) is utilized from PETSc.
> 
>     Note that there is only one simple set of include files and modules for Fortran; see the newest documentation.
> 
> 
> http://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/Sys/UsingFortran.html
> 
>    Matt
>  
> 
>    Barry
> 
> 
> > On Oct 5, 2017, at 11:48 PM, Randy Michael Churchill <rchurchi at pppl.gov> wrote:
> >
> > A simple setup question with TAO: if I were to convert the rosenbrock1f.F90 example to use a module instead of common structures, how would I setup the include statements? I've tried various combinations (using petscXXXdef.h, petscXXX.h, petscXXX.h90, along with use petscXXX), but seem to get errors with all.
> >
> > file:rosenbrock1f.h:
> > module commondat
> >       PetscReal    ::    alpha
> >       PetscInt     ::    n
> > end module commondat
> >
> > file:rosenbrock1f.90:
> > program rosenbrock1f
> > !!include statements??? which and where???!!!
> > use commondat
> > ...
> >
> > subroutine FormFunctionGradient(tao, X, f, G, dummy, ierr)
> >      use commondat
> >      implicit none
> >      ...
> >
> > (https://www.mcs.anl.gov/petsc/petsc-dev/src/tao/unconstrained/examples/tutorials/rosenbrock1f.F90.html)
> 
> 
> 
> 
> -- 
> 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
> 
> https://www.cse.buffalo.edu/~knepley/
> 
> 
> 
> -- 
> R. Michael Churchill