combining python with fortran, trouble with tutorial -

-

i'm following tutorial http://www.sam.math.ethz.ch/~raoulb/teaching/pythontutorial/combining.html

i use same code

      program hwtest       real*8 r1, r2 ,s       r1 = 1.0       r2 = 0.0       s = hw1(r1, r2)       write(*,*) 'hw1, result:',s       write(*,*) 'hw2, result:'       call hw2(r1, r2)       call hw3(r1, r2, s)       write(*,*) 'hw3, result:', s       end        real*8 function hw1(r1, r2)       real*8 r1, r2       hw1 = sin(r1 + r2)       return       end        subroutine hw2(r1, r2)       real*8 r1, r2, s       s = sin(r1 + r2)       write(*,1000) 'hello, world! sin(',r1+r2,')=',s  1000 format(a,f6.3,a,f8.6)       return       end  c     special version of hw1 result c     returned argument:        subroutine hw3_v1(r1, r2, s)       real*8 r1, r2, s       s = sin(r1 + r2)       return       end  c     f2py treats s input arg. in hw3_v1; fix this:        subroutine hw3(r1, r2, s)       real*8 r1, r2, s cf2py intent(out) s       s = sin(r1 + r2)       return       end  c     test case sensitivity:        subroutine hw4(r1, r2, s)       real*8 r1, r2, s cf2py intent(out) s       s = sin(r1 + r2)       return       end  c end of f77 file

but when use command

f2py -m hw -c ../hw.f

i following error:

gfortran:f77: ../hw.f ../hw.f:5.13:        s = hw1(r1, r2)                                                                 1 error: return type mismatch of function 'hw1' @ (1) (real(4)/real(8)) ../hw.f:5.13:        s = hw1(r1, r2)                                                                 1 error: return type mismatch of function 'hw1' @ (1) (real(4)/real(8)) error: command "/usr/bin/gfortran -wall -ffixed-form -fno-second-underscore -fpic -o3 -funroll-loops -i/tmp/tmprghzqj/src.linux-x86_64-2.7 -i/usr/lib/python2.7/dist-packages/numpy/core/include -i/usr/include/python2.7 -c -c ../hw.f -o /tmp/tmprghzqj/fortran/hw.o" failed exit status 1

i have no experience fortran, need implement existing fortran code use in python.

the interface hw1 missing in code. either specify explicitly:

program hwtest   real*8 r1, r2 ,s   real*8 hw1 c [...]

or, better, put functions module, , use module.

btw: change *8 better suited, like: 

program hwtest   integer,parameter :: dp = kind(1.d0)   real(dp) r1, r2 ,s   real(dp) hw1

or, use iso_fortran_env, , pre-defined parameters real64 , real32.


Comments

Post a Comment

Popular posts from this blog

java - Plugin org.apache.maven.plugins:maven-install-plugin:2.4 or one of its dependencies could not be resolved -

-
i getting error message , further 1 here: failed read artifact descriptor org.apache.maven.plugins:maven-install-plugin:jar:2.4: even though when go .m2 repository path c:\users\a591024.m2\repository\org\apache\maven\plugins\maven-resources-plugin\2.4 the maven-resources-plugin-2.4.jar file there , maven-resources-plugin-2.4.pom there. have tried delete .m2 repository no luck. far maven not build project desperately need to. any ideas? it seems proxy issue please aad following in setting.xml grabbing lan setting , proxy information. <proxies> <proxy> <id>proxy</id> <active>true</active> <protocol>http</protocol> <host>hostname</host> <port>port</port> </proxy> </proxies> this frequent issue, please see below link well. plugin org.apache.maven.plugins:maven-clean-plugin:2.5 or 1 of dependencies not resolved maven: failed ret...
Read more

Round ImageView Android -

-
Image
i working on android application in doing rounded image view. working fine images, images 160x120 resolution shows oval shaped. code custom imageview given below, please me out here: public class roundimage extends drawable { private final bitmap mbitmap; private final paint mpaint; private final rectf mrectf; private final int mbitmapwidth; private final int mbitmapheight; public roundimage(bitmap bitmap) { mbitmap = bitmap; mrectf = new rectf(); mpaint = new paint(); mpaint.setantialias(true); mpaint.setdither(true); final bitmapshader shader = new bitmapshader(bitmap, shader.tilemode.clamp, shader.tilemode.clamp); mpaint.setshader(shader); mbitmapwidth = mbitmap.getwidth(); mbitmapheight = mbitmap.getheight(); } @override public void draw(canvas canvas) { canvas.drawoval(mrectf, mpaint); } @override protected void onboundschange(rect bound...
Read more

How can I utilize Yahoo Weather API in android -

-
Image
i want develop app using yahoo weather api android. want current , forecast weather details. i have tried different sources (e.g.: https://github.com/survivingwithandroid/weatherlib , getting error weatherlib api android , unable execute activity in android weatherlib ) , getting different errors m unable resplve. there can proper tutorial or guide! error: -17 08:50:48.961 3168-3168/com.survivingwithandroid.weather e/androidruntime﹕ fatal exception: main process: com.survivingwithandroid.weather, pid: 3168 java.lang.runtimeexception: unable start activity componentinfo{com.survivingwithandroid.weather/com.survivingwithandroid.weather.mainactivity}: java.lang.nullpointerexception: attempt invoke interface method 'boolean java.util.list.isempty()' on null object reference @ android.app.activitythread.performlaunchactivity(activitythread.java:2298) @ android.app.activitythread.handlelaunchactivity(activitythread.java:2360) @ android.app.activi...
Read more