[f2py] (no subject)
Razibul Islam <tamal105-/[email protected]> Thu, 24 Jan 2013 14:19:39 -0800 (PST)
| Newsgroups | gmane.comp.python.f2py.user |
|---|---|
| Message-ID | <[email protected]> |
--===============5109665736910653735==
Content-Type: multipart/alternative; boundary="1234169369-2058904177-1359065979=:90075"
--1234169369-2058904177-1359065979=:90075
Content-Type: text/plain; charset=iso-8859-1
Content-Transfer-Encoding: quoted-printable
Hi, I need some information about wrapping MPI/Fortran code into Python mod=
ule. So, suppose I have code below, which is MPI/Fortran code.=0AI wonder w=
hat procedure I need to use to wrap it and execute. I really would like to =
know how to execute this code in Python, and I would like to =0Aknow what t=
he set up of the Python should be : regular python, mpi4py, or some other =
Python distribution. I prefer to not change Fortran code and =0AI wonder if=
it is possible. if so, I would very much appreciate some advice on how to =
do this. I have search the Internet and found info regarding =0Adifferent p=
arallel version to compile with f2py. So far I understood that the problem =
in my path of "mpif.h" file in f2py flag. It would be very helpful if someo=
ne=0A=A0could give me some information on this. I have asked for help in he=
re and this is really urgent as it is the last step of my master thesis. I =
would highly =0Aappreciate if anyone can give me some advice in this regard=
. I have tried the following code to compile using f2py with the option men=
tioned below and got =0Athe konsole output also mentioned below. =0A=0A#---=
------------------------------------------Fortran code --------------------=
-------------------------------------------=0A=0Aprogram hello_parallel ! I=
nclude the MPI library definitons: include 'mpif.h' integer numtasks, rank,=
ierr, rc, len, i character*(MPI_MAX_PROCESSOR_NAME) name ! Initialize the =
MPI library: call MPI_INIT(ierr) if (ierr .ne. MPI_SUCCESS) then print *,'E=
rror starting MPI program. Terminating.' call MPI_ABORT(MPI_COMM_WORLD, rc,=
ierr) end if ! Get the number of processors this job is using: call MPI_CO=
MM_SIZE(MPI_COMM_WORLD, numtasks, ierr) ! Get the rank of the processor thi=
s thread is running on. (Each ! processor has a unique rank.) call MPI_COM=
M_RANK(MPI_COMM_WORLD, rank, ierr) ! Get the name of this processor (usuall=
y the hostname) call MPI_GET_PROCESSOR_NAME(name, len, ierr) if (ierr .ne. =
MPI_SUCCESS) then print *,'Error getting processor name. Terminating.' call=
MPI_ABORT(MPI_COMM_WORLD, rc, ierr) end if print "('hello_parallel.f: Numb=
er of tasks=3D',I3,' My rank=3D',I3,' My name=3D',A,'')",& numtasks, rank, =
trim(name) ! Tell the MPI library to release all
resources it is using: call MPI_FINALIZE(ierr) end program hello_parallel=
=0A=0A=0A#--------------------------------Used f2py command=A0 a python ext=
ension module---------------------------=0A=0Af2py -m partest01 -c hello_pa=
r_sub.f -lmpi -L/opt/mpi/lib -L/usr/lib/mpi/gcc -lgm=0A=0A=0A=0A#----------=
-----------------------------Konsole output begins-------------------------=
-------------=0A=0A=0Arunning build=0Arunning config_cc=0Aunifing config_cc=
, config,=0A build_clib, build_ext, build commands --compiler options=0Arun=
ning config_fc=0Aunifing config_fc, config, build_clib, build_ext, build co=
mmands=0A --fcompiler options=0Arunning build_src=0Abuild_src=0Abuilding ex=
tension "partest01" sources=0Af2py options: []=0Af2py:> /tmp/tmpSO6Ywv/src.=
linux-i686-2.7/partest01module.c=0Acreating /tmp/tmpSO6Ywv=0Acreating /tmp/=
tmpSO6Ywv/src.linux-i686-2.7=0AReading fortran codes...=0A=A0=A0=A0=A0=A0=
=A0=A0 Reading file 'hello_par_sub.f' (format:fix,strict)=0ALine #5 in hell=
o_par_sub.f:"=A0=A0=A0=A0=A0 include 'mpif.h'"=0A=A0=A0=A0=A0=A0=A0=A0 read=
fortrancode: could not find include file 'mpif.h' in . Ignoring.=0Armbadnam=
e1: Replacing "rank" with "rank_bn".=0Armbadname1: Replacing "len" with "le=
n_bn".=0Armbadname1: Replacing "rank" with "rank_bn".=0Armbadname1: Replaci=
ng "len" with "len_bn".=0APost-processing...=0A=A0=A0=A0=A0=A0=A0=A0 Block:=
=0A partest01=0A=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0 Block: hello_parallel=0APost-processing (stage 2)...=0ABuildin=
g modules...=0A=A0=A0=A0=A0=A0=A0=A0 Building module "partest01"...=0A=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Constructing wrapper function "h=
ello_parallel"...=0A=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 hel=
lo_parallel()=0A=A0=A0=A0=A0=A0=A0=A0 Wrote C/API module "partest01" to fil=
e "/tmp/tmpSO6Ywv/src.linux-i686-2.7/partest01module.c"=0A=A0 adding '/tmp/=
tmpSO6Ywv/src.linux-i686-2.7/fortranobject.c' to sources.=0A=A0 adding '/tm=
p/tmpSO6Ywv/src.linux-i686-2.7' to include_dirs.=0Acopying /usr/lib/python2=
.7/site-packages/numpy/f2py/src/fortranobject.c ->=0A /tmp/tmpSO6Ywv/src.li=
nux-i686-2.7=0Acopying /usr/lib/python2.7/site-packages/numpy/f2py/src/fort=
ranobject.h -> /tmp/tmpSO6Ywv/src.linux-i686-2.7=0Abuild_src: building npy-=
pkg config files=0Arunning build_ext=0Acustomize UnixCCompiler=0Acustomize =
UnixCCompiler using build_ext=0Acustomize GnuFCompiler=0ACould not locate e=
xecutable g77=0ACould not locate executable=0A f77=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =0Acustomize=0A IntelFCompile=
r=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0 =0ACould not locate=0A executable ifort=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =0ACould not locate execut=
able=0A ifc=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0 =0Acustomize=0A LaheyFCompiler=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =0ACould not locate=
=0A executable lf95=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0 =0Acustomize=0A PGroupFCompiler=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =0ACould not loca=
te=0A executable pgfortran=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0 =0Acustomize=0A AbsoftFCompiler=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =0ACould not locate=0A exe=
cutable f90=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0 =0Acustomize=0A NAGFCompiler=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =0ACould=0A not l=
ocate executable f95=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0 =0Acustomize VastFCompiler=0Acustomize=0A CompaqFCompiler=
=0ACould not locate executable fort=0Acustomize IntelItaniumFCompiler=0ACou=
ld not locate executable efort=0ACould not locate executable efc=0Acustomiz=
e IntelEM64TFCompiler=0Acustomize Gnu95FCompiler=0AFound executable /usr/bi=
n/gfortran=0Acustomize Gnu95FCompiler=0Acustomize Gnu95FCompiler using buil=
d_ext=0Abuilding 'partest01' extension=0Acompiling C sources=0AC compiler: =
gcc -pthread -fno-strict-aliasing -fomit-frame-pointer -fmessage-length=3D0=
-O2 -Wall -D_FORTIFY_SOURCE=3D2 -fstack-protector -funwind-tables -fasynch=
ronous-unwind-tables -g -DNDEBUG -fomit-frame-pointer -fmessage-length=3D0 =
-O2 -Wall -D_FORTIFY_SOURCE=3D2 -fstack-protector -funwind-tables -fasynchr=
onous-unwind-tables -g -fPIC=0A=0Acreating /tmp/tmpSO6Ywv/tmp=0Acreating /t=
mp/tmpSO6Ywv/tmp/tmpSO6Ywv=0Acreating /tmp/tmpSO6Ywv/tmp/tmpSO6Ywv/src.linu=
x-i686-2.7=0Acompile options: '-I/tmp/tmpSO6Ywv/src.linux-i686-2.7=0A -I/us=
r/lib/python2.7/site-packages/numpy/core/include -I/usr/include/python2.7 -=
c'=0Agcc: /tmp/tmpSO6Ywv/src.linux-i686-2.7/partest01module.c=0Agcc: /tmp/t=
mpSO6Ywv/src.linux-i686-2.7/fortranobject.c=0Acompiling Fortran sources=0AF=
ortran f77 compiler: /usr/bin/gfortran -Wall -ffixed-form -fno-second-under=
score -fPIC -O3 -funroll-loops=0AFortran f90 compiler: /usr/bin/gfortran -W=
all -fno-second-underscore -fPIC -O3 -funroll-loops=0AFortran fix compiler:=
/usr/bin/gfortran -Wall -ffixed-form -fno-second-underscore -Wall -fno-sec=
ond-underscore -fPIC -O3 -funroll-loops=0Acompile options: '-I/tmp/tmpSO6Yw=
v/src.linux-i686-2.7 -I/usr/lib/python2.7/site-packages/numpy/core/include =
-I/usr/include/python2.7 -c'=0Agfortran:f77: hello_par_sub.f=0Ahello_par_su=
b.f:4: Error: Can't open included file 'mpif.h'=0Ahello_par_sub.f:4: Error:=
Can't open included file 'mpif.h'=0Aerror: Command "/usr/bin/gfortran -Wal=
l -ffixed-form -fno-second-underscore -fPIC -O3=0A -funroll-loops -I/tmp/tm=
pSO6Ywv/src.linux-i686-2.7 -I/usr/lib/python2.7/site-packages/numpy/core/in=
clude -I/usr/include/python2.7 -c -c hello_par_sub.f -o /tmp/tmpSO6Ywv/hell=
o_par_sub.o" failed with exit status 1=0A=0A=0A=0A=0A=0A#------------------=
------------------------ End of Konsole Output=0A -------------------------=
--
--1234169369-2058904177-1359065979=:90075
Content-Type: text/html; charset=iso-8859-1
Content-Transfer-Encoding: quoted-printable
<html><body><div style=3D"color:#000; background-color:#fff; font-family:ti=
mes new roman, new york, times, serif;font-size:10pt"><div id=3D"yiv4408169=
82"><div><div style=3D"color:#000;background-color:#fff;font-family:times n=
ew roman, new york, times, serif;font-size:10pt;"><div class=3D"yiv44081698=
2yui_3_7_2_17_1359034494621_59" id=3D"yiv440816982yui_3_7_2_17_135903449462=
1_60" style=3D"font-family:times new roman, new york, times, serif;font-siz=
e:10pt;"><div id=3D"yiv440816982yui_3_7_2_17_1359034494621_93" class=3D"yiv=
440816982yui_3_7_2_17_1359034494621_63" style=3D"font-family:times new roma=
n, new york, times, serif;font-size:12pt;"><div id=3D"yiv440816982"><div id=
=3D"yiv440816982yui_3_7_2_17_1359034494621_114"><div id=3D"yiv440816982yui_=
3_7_2_17_1359034494621_113" class=3D"yiv440816982yui_3_7_2_17_1359034494621=
_69" style=3D"color:#000;background-color:#fff;font-family:times new roman,=
new york, times, serif;font-size:12pt;"><pre id=3D"yiv97336293yui_3_7_2_16=
_1358775566147_68">Hi,=0A=0AI need some information about wrapping MPI/Fort=
ran code into Python module. So, suppose I have code below, which is MPI/Fo=
rtran code.<br>I wonder what procedure I need to use to wrap it and execute=
. I really would like to know how to execute this code in Python, and I wou=
ld like to <br>know what the set up of the Python should be : regular pytho=
n, mpi4py, or some other Python distribution. I prefer to not change Fortr=
an code and <br>I wonder if it is possible. if so, I would very much apprec=
iate some advice on how to do this. I have search the Internet and found in=
fo regarding <br>different parallel version to compile with f2py. So far I =
understood that the problem in my path of "mpif.h" file in f2py flag. It wo=
uld be very helpful if someone<br> could give me some information on t=
his. I have asked for help in here and this is really urgent as it is the l=
ast step of my master thesis. I would highly <br>appreciate if anyone can g=
ive me some advice in
this regard. I have tried the following code to compile using f2py with th=
e option mentioned below and got <br>the konsole output also mentioned belo=
w. </pre><div class=3D"yiv440816982yui_3_7_2_17_1359034494621_74" style=3D"=
color:rgb(0, 0, 0);font-size:16px;font-family:times new roman, new york, ti=
mes, serif;background-color:transparent;font-style:normal;"><br></div><div =
id=3D"yiv440816982yui_3_7_2_17_1359034494621_134" class=3D"yiv440816982yui_=
3_7_2_17_1359034494621_75" style=3D"color:rgb(0, 0,=0A 0);font-size:16px;fo=
nt-family:times new roman, new york, times, serif;background-color:transpar=
ent;font-style:normal;">#---------------------------------------------Fortr=
an code=0A ---------------------------------------------------------------<=
br></div><pre><span style=3D"font-size:13px;">program hello_parallel=0A=0A =
! Include the MPI library definitons:=0A include 'mpif.h'=0A=0A integer =
numtasks, rank, ierr, rc, len, i=0A character*(MPI_MAX_PROCESSOR_NAME) nam=
e=0A=0A ! Initialize the MPI library:=0A call MPI_INIT(ierr)=0A if (ierr=
.ne. MPI_SUCCESS) then=0A print *,'Error starting MPI program. Termina=
ting.'=0A call MPI_ABORT(MPI_COMM_WORLD, rc, ierr)=0A end if=0A=0A ! =
Get the number of processors this job is using:=0A call MPI_COMM_SIZE(MPI_=
COMM_WORLD, numtasks, ierr)=0A=0A ! Get the rank of the processor this thr=
ead is running on. (Each=0A ! processor has a unique rank.)=0A call MPI_=
COMM_RANK(MPI_COMM_WORLD, rank, ierr)=0A=0A ! Get the name of this process=
or (usually the hostname)=0A call MPI_GET_PROCESSOR_NAME(name, len, ierr)=
=0A if (ierr .ne. MPI_SUCCESS) then=0A print *,'Error getting processo=
r name. Terminating.'=0A call MPI_ABORT(MPI_COMM_WORLD, rc, ierr)=0A e=
nd if=0A=0A print "('hello_parallel.f: Number of tasks=3D',I3,' My rank=3D=
',I3,' My name=3D',A,'')",&=0A numtasks, rank, trim(name)=0A=0A =
! Tell the MPI library to release all resources it is using:=0A call MPI_F=
INALIZE(ierr)=0A=0Aend program hello_parallel</span></pre><div class=3D"yiv=
440816982yui_3_7_2_17_1359034494621_77" style=3D"color:rgb(0, 0, 0);font-si=
ze:16px;font-family:times new roman, new york, times, serif;background-colo=
r:transparent;font-style:normal;"><br></div><div class=3D"yiv440816982yui_3=
_7_2_17_1359034494621_78" style=3D"color:rgb(0, 0, 0);font-size:16px;font-f=
amily:times new roman, new york, times, serif;background-color:transparent;=
font-style:normal;"><br></div><div class=3D"yiv440816982yui_3_7_2_17_135903=
4494621_79" style=3D"color:rgb(0, 0, 0);font-size:16px;font-family:times ne=
w roman, new york, times, serif;background-color:transparent;font-style:nor=
mal;">#--------------------------------Used f2py command a python ext=
ension module---------------------------<br></div><div class=3D"yiv44081698=
2yui_3_7_2_17_1359034494621_80" style=3D"color:rgb(0, 0, 0);font-size:16px;=
font-family:times new roman, new york, times,
serif;background-color:transparent;font-style:normal;">f2py -m partest01 -=
c hello_par_sub.f -lmpi -L/opt/mpi/lib -L/usr/lib/mpi/gcc -lgm<br><br></div=
><div class=3D"yiv440816982yui_3_7_2_17_1359034494621_81" style=3D"color:rg=
b(0, 0, 0);font-size:16px;font-family:times new roman, new york, times, ser=
if;background-color:transparent;font-style:normal;"><br></div><div class=3D=
"yiv440816982yui_3_7_2_17_1359034494621_82" style=3D"color:rgb(0, 0, 0);fon=
t-size:16px;font-family:times new roman, new york, times, serif;background-=
color:transparent;font-style:normal;">#------------------------------------=
---Konsole output begins--------------------------------------</div><div cl=
ass=3D"yiv440816982yui_3_7_2_17_1359034494621_83" style=3D"color:rgb(0, 0, =
0);font-size:16px;font-family:times new roman, new york, times, serif;backg=
round-color:transparent;font-style:normal;"><br></div><div class=3D"yiv4408=
16982yui_3_7_2_17_1359034494621_84" style=3D"color:rgb(0, 0,
0);font-size:16px;font-family:times new roman, new york, times, serif;back=
ground-color:transparent;font-style:normal;"><br></div><div class=3D"yiv440=
816982yui_3_7_2_17_1359034494621_85" style=3D"color:rgb(0, 0, 0);font-size:=
16px;font-family:times new roman, new york, times, serif;background-color:t=
ransparent;font-style:normal;">running build<br>running config_cc<br>unifin=
g config_cc, config,=0A build_clib, build_ext, build commands --compiler op=
tions<br>running config_fc<br>unifing config_fc, config, build_clib, build_=
ext, build commands=0A --fcompiler options<br>running build_src<br>build_sr=
c<br>building extension "partest01" sources<br>f2py options: []<br>f2py:>=
; /tmp/tmpSO6Ywv/src.linux-i686-2.7/partest01module.c<br>creating /tmp/tmpS=
O6Ywv<br>creating /tmp/tmpSO6Ywv/src.linux-i686-2.7<br>Reading fortran code=
s...<br> Reading file 'hello_par_=
sub.f' (format:fix,strict)<br>Line #5 in hello_par_sub.f:"  =
; include 'mpif.h'"<br> &nbs=
p; readfortrancode: could not find include file 'mpif.h' in . Ignoring.<br>=
rmbadname1: Replacing "rank" with "rank_bn".<br>rmbadname1: Replacing "len"=
with "len_bn".<br>rmbadname1: Replacing "rank" with "rank_bn".<br>rmbadnam=
e1: Replacing "len" with "len_bn".<br>Post-processing...<br> &nb=
sp; Block:=0A partest01<br> =
&nb=
sp; Block: hello_parallel<br>Post-proce=
ssing (stage 2)...<br>Building modules...<br> =
Building module "partest01"...<br> &nbs=
p; Constructing=
wrapper function "hello_parallel"...<br> &nbs=
p; hello_=
parallel()<br> Wrote C/API module=
"partest01" to file "/tmp/tmpSO6Ywv/src.linux-i686-2.7/partest01module.c"<=
br> adding '/tmp/tmpSO6Ywv/src.linux-i686-2.7/fortranobject.c' to sou=
rces.<br> adding '/tmp/tmpSO6Ywv/src.linux-i686-2.7' to include_dirs.=
<br>copying /usr/lib/python2.7/site-packages/numpy/f2py/src/fortranobject.c=
->=0A /tmp/tmpSO6Ywv/src.linux-i686-2.7<br>copying /usr/lib/python2.7/s=
ite-packages/numpy/f2py/src/fortranobject.h -> /tmp/tmpSO6Ywv/src.linux-=
i686-2.7<br>build_src: building npy-pkg config files<br>running build_ext<b=
r>customize UnixCCompiler<br>customize UnixCCompiler using build_ext<br>cus=
tomize GnuFCompiler<br>Could not locate executable g77<br>Could not locate =
executable=0A f77 &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; <br>customize=0A IntelFCompiler &n=
bsp;  =
; &n=
bsp;  =
; &n=
bsp;  =
; &n=
bsp;  =
; &n=
bsp;  =
; &n=
bsp;  =
; <br>Cou=
ld not locate=0A executable ifort =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; <br>Could not locate executable=0A ifc &nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
<br>customize=0A LaheyFCompiler &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; &nbs=
p; &=
nbsp; <br>Could not locate=0A executable lf95 =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; <br>customize=0A PGroupFCompi=
ler =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; <br>Could not locate=0A executable pgfortran&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
<br>customize=0A AbsoftFCompiler =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
<br>Could not locate=0A executable f90 =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; <br>customize=0A NAGFCompiler =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
<br>Could=0A not locate executable f95 =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; =
&nb=
sp; <br>customize Vas=
tFCompiler<br>customize=0A CompaqFCompiler<br>Could not locate executable f=
ort<br>customize IntelItaniumFCompiler<br>Could not locate executable efort=
<br>Could not locate executable efc<br>customize IntelEM64TFCompiler<br>cus=
tomize Gnu95FCompiler<br>Found executable /usr/bin/gfortran<br>customize Gn=
u95FCompiler<br>customize Gnu95FCompiler using build_ext<br>building 'parte=
st01' extension<br>compiling C sources<br>C compiler: gcc -pthread -fno-str=
ict-aliasing -fomit-frame-pointer -fmessage-length=3D0 -O2 -Wall -D_FORTIFY=
_SOURCE=3D2 -fstack-protector -funwind-tables -fasynchronous-unwind-tables =
-g -DNDEBUG -fomit-frame-pointer -fmessage-length=3D0 -O2 -Wall -D_FORTIFY_=
SOURCE=3D2 -fstack-protector -funwind-tables -fasynchronous-unwind-tables -=
g -fPIC<br><br>creating /tmp/tmpSO6Ywv/tmp<br>creating /tmp/tmpSO6Ywv/tmp/t=
mpSO6Ywv<br>creating /tmp/tmpSO6Ywv/tmp/tmpSO6Ywv/src.linux-i686-2.7<br>com=
pile options: '-I/tmp/tmpSO6Ywv/src.linux-i686-2.7=0A -I/usr/lib/python2.7/=
site-packages/numpy/core/include -I/usr/include/python2.7 -c'<br>gcc: /tmp/=
tmpSO6Ywv/src.linux-i686-2.7/partest01module.c<br>gcc: /tmp/tmpSO6Ywv/src.l=
inux-i686-2.7/fortranobject.c<br>compiling Fortran sources<br>Fortran f77 c=
ompiler: /usr/bin/gfortran -Wall -ffixed-form -fno-second-underscore -fPIC =
-O3 -funroll-loops<br>Fortran f90 compiler: /usr/bin/gfortran -Wall -fno-se=
cond-underscore -fPIC -O3 -funroll-loops<br>Fortran fix compiler: /usr/bin/=
gfortran -Wall -ffixed-form -fno-second-underscore -Wall -fno-second-unders=
core -fPIC -O3 -funroll-loops<br>compile options: '-I/tmp/tmpSO6Ywv/src.lin=
ux-i686-2.7 -I/usr/lib/python2.7/site-packages/numpy/core/include -I/usr/in=
clude/python2.7 -c'<br>gfortran:f77: hello_par_sub.f<br>hello_par_sub.f:4: =
Error: Can't open included file 'mpif.h'<br>hello_par_sub.f:4: Error: Can't=
open included file 'mpif.h'<br>error: Command "/usr/bin/gfortran -Wall -ff=
ixed-form -fno-second-underscore -fPIC -O3=0A -funroll-loops -I/tmp/tmpSO6Y=
wv/src.linux-i686-2.7 -I/usr/lib/python2.7/site-packages/numpy/core/include=
-I/usr/include/python2.7 -c -c hello_par_sub.f -o /tmp/tmpSO6Ywv/hello_par=
_sub.o" failed with exit status 1<br></div><div class=3D"yiv440816982yui_3_=
7_2_17_1359034494621_86" style=3D"color:rgb(0, 0, 0);font-size:16px;font-fa=
mily:times new roman, new york, times, serif;background-color:transparent;f=
ont-style:normal;"><br></div><div class=3D"yiv440816982yui_3_7_2_17_1359034=
494621_87" style=3D"color:rgb(0, 0, 0);font-size:16px;font-family:times new=
roman, new york, times, serif;background-color:transparent;font-style:norm=
al;"><br></div><div class=3D"yiv440816982yui_3_7_2_17_1359034494621_88" sty=
le=3D"color:rgb(0, 0, 0);font-size:16px;font-family:times new roman, new yo=
rk, times, serif;background-color:transparent;font-style:normal;"><br></div=
><div class=3D"yiv440816982yui_3_7_2_17_1359034494621_89" style=3D"color:rg=
b(0, 0, 0);font-size:16px;font-family:times new roman,
new york, times, serif;background-color:transparent;font-style:normal;"> <=
br>#------------------------------------------ End of Konsole Output=0A ---=
------------------------<br></div><div></div></div></div></div><br><br> </d=
iv> </div> </div></div></div></div></body></html>
--1234169369-2058904177-1359065979=:90075--
--===============5109665736910653735==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline
_______________________________________________
f2py-users mailing list
f2py-users-Y4l6ocDipWCuvFJfX82//[email protected]
http://cens.ioc.ee/mailman/listinfo/f2py-users
--===============5109665736910653735==--