[petsc-users] VecGetArrayReadF90 with gfortran

[Satish Balay]

Attaching the modified version that works with complex [using complex build of petsc]

BTW: I'm using the following gfortran [but I don't think the version makes any difference] with petsc-3.8 (maint branch):

$ gfortran --version
GNU Fortran (GCC) 8.0.1 20180222 (Red Hat 8.0.1-0.16)

Satish

--------

balay at asterix /home/balay/download-pine
$ make ex2f
mpif90 -c -Wall -ffree-line-length-0 -Wno-unused-dummy-argument -g    -I/home/balay/petsc/include -I/home/balay/petsc/arch-cmplx/include    -o ex2f.o ex2f.F90
mpif90 -Wall -ffree-line-length-0 -Wno-unused-dummy-argument -g   -o ex2f ex2f.o  -Wl,-rpath,/home/balay/petsc/arch-cmplx/lib -L/home/balay/petsc/arch-cmplx/lib -Wl,-rpath,/home/balay/soft/mpich-3.3b1/lib -L/home/balay/soft/mpich-3.3b1/lib -Wl,-rpath,/usr/lib/gcc/x86_64-redhat-linux/8 -L/usr/lib/gcc/x86_64-redhat-linux/8 -lpetsc -llapack -lblas -lX11 -lpthread -lm -lmpifort -lgfortran -lm -lgfortran -lm -lquadmath -lmpicxx -lstdc++ -lm -Wl,-rpath,/home/balay/soft/mpich-3.3b1/lib -L/home/balay/soft/mpich-3.3b1/lib -Wl,-rpath,/usr/lib/gcc/x86_64-redhat-linux/8 -L/usr/lib/gcc/x86_64-redhat-linux/8 -ldl -Wl,-rpath,/home/balay/soft/mpich-3.3b1/lib -lmpi -lgcc_s -ldl 
/usr/bin/rm -f ex2f.o
balay at asterix /home/balay/download-pine
$ ./ex2f
 setting values
 assembling vector
 copying vector to array
 N =           10
 m =           10
 getting vecX size
 getting ownership range
 ownership range:            0          10
 copy values to global array
 vecX_pt:
 SIZE =           10
                  (1.0000000000000000,0.0000000000000000)
                  (2.0000000000000000,0.0000000000000000)
                  (3.0000000000000000,0.0000000000000000)
                  (4.0000000000000000,0.0000000000000000)
                  (5.0000000000000000,0.0000000000000000)
                  (6.0000000000000000,0.0000000000000000)
                  (7.0000000000000000,0.0000000000000000)
                  (8.0000000000000000,0.0000000000000000)
                  (9.0000000000000000,0.0000000000000000)
                  (10.000000000000000,0.0000000000000000)
 copy_buffer = 
                  (1.0000000000000000,0.0000000000000000)
                  (2.0000000000000000,0.0000000000000000)
                  (3.0000000000000000,0.0000000000000000)
                  (4.0000000000000000,0.0000000000000000)
                  (5.0000000000000000,0.0000000000000000)
                  (6.0000000000000000,0.0000000000000000)
                  (7.0000000000000000,0.0000000000000000)
                  (8.0000000000000000,0.0000000000000000)
                  (9.0000000000000000,0.0000000000000000)
                  (10.000000000000000,0.0000000000000000)
 calling allreduce
 buffer = 
                  (1.0000000000000000,0.0000000000000000)
                  (2.0000000000000000,0.0000000000000000)
                  (3.0000000000000000,0.0000000000000000)
                  (4.0000000000000000,0.0000000000000000)
                  (5.0000000000000000,0.0000000000000000)
                  (6.0000000000000000,0.0000000000000000)
                  (7.0000000000000000,0.0000000000000000)
                  (8.0000000000000000,0.0000000000000000)
                  (9.0000000000000000,0.0000000000000000)
                  (10.000000000000000,0.0000000000000000)
balay at asterix /home/balay/download-pine
$ 

On Fri, 23 Mar 2018, Satish Balay wrote:

> I'm attaching the modified code - that appears to work fine for me. My changes are:
> 
> 1. use .F90 suffix - so the compiler does the f90 freeform/preprocessing automatically
> 2. remove slepc references and use petsc only
> 3. compile with a petsc formatted makefile
> 4. replace complex with real*8 - thats compatible with PetscScalar
> 
> So perhaps you can modify your version of test code - to be closer to
> this one [or vice versa] to see what is triggering the difference in
> behavior.
> 
> Satish
> 
> On Fri, 23 Mar 2018, Samuel Lanthaler wrote:
> 
> > Hi all,
> > 
> > I am having trouble using the function VecGetArrayReadF90 under gfortran. I
> > have created a minimal example and put it in the attachment. Basically, it
> > appears that the input pointer is empty upon return from
> > VecGetArrayReadF90. The code runs fine, when compiled with ifort, so I
> > don't know what's going on.
> > 
> > I am attaching the code and makefile that I used, as well as the command
> > line output that I can see. Though the errors are just to do with the fact
> > that the output pointer is empty, so the command line output is not all
> > that helpful...
> > 
> > Thanks in advance for your help!
> > Sam
> > 
> 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: ex2f.F90
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180323/a21dbc36/attachment-0001.ksh>
-------------- next part --------------
CFLAGS	         =
FFLAGS	         =
CPPFLAGS         =
FPPFLAGS         =

include ${PETSC_DIR}/lib/petsc/conf/variables
include ${PETSC_DIR}/lib/petsc/conf/rules

ex1f: ex1f.o  chkopts
	-${FLINKER} -o ex1f ex1f.o  ${PETSC_KSP_LIB}
	${RM} ex1f.o

ex2f: ex2f.o  chkopts
	-${FLINKER} -o ex2f ex2f.o  ${PETSC_KSP_LIB}
	${RM} ex2f.o