Use Fortran subroutine in R? Error: Return type mismatch -


i'm trying learn how use fortran code inside r. able follow this tutorial. now, i'm trying use base plus this fortran program calculate pi. create file fpi.f90 code:

subroutine pi(avepi, darts, rounds) double precision, intent(out)   ::  avepi integer, intent(in)             ::  darts, rounds integer                         ::  master, rank, i, n integer, allocatable            ::  seed(:) double precision                ::  pi_est, homepi, pirecv, pisum  ! set 0 in sequential run rank = 0 ! initialize random number generator ! make sure seed different each task call random_seed() call random_seed(size = n) allocate(seed(n)) seed = 12 + rank*11 call random_seed(put=seed(1:n)) deallocate(seed)  avepi = 0 = 0, rounds-1   pi_est = dboard(darts)   ! calculate average value of pi on iterations   avepi = ((avepi*i) + pi_est)/(i + 1) end end subroutine pi   double precision function dboard(darts) integer, intent(in)           :: darts double precision              :: x_coord, y_coord integer                       :: score, n  score = 0 n = 1, darts   call random_number(x_coord)   call random_number(y_coord)    if ((x_coord**2 + y_coord**2) <= 1.0d0)   score = score + 1   end if end dboard = 4.0d0*score/darts  end function

when

$ r cmd shlib ./fortran/fpi.f90 gfortran  -fpic -g -o2 -fstack-protector-strong  -c  fortran/fpi.f90 -o fortran/fpi.o fortran/fpi.f90:22.15:  pi_est = dboard(darts)                1 error: return type mismatch of function 'dboard' @ (1) (real(4)/real(8)) /usr/lib/r/etc/makeconf:161: recipe target 'fortran/fpi.o' failed make: *** [fortran/fpi.o] error 1

what doing wrong?

after adding double precision :: dboard pi different error. 

r cmd shlib ./fortran/fpi.f90 gfortran -shared -l/usr/lib/r/lib -wl,-z,relro -o fortran/fpi.so ./fortran/fpi.o -l/usr/lib/r/lib -lr /usr/bin/ld: ./fortran/fpi.o: relocation r_x86_64_32 against `.rodata' can not used when making shared object; recompile -fpic ./fortran/fpi.o: error adding symbols: bad value collect2: error: ld returned 1 exit status /usr/share/r/share/make/shlib.mk:6: recipe target 'fortran/fpi.so' failed make: *** [fortran/fpi.so] error 1

you not use implicit none. bad! due implicit typing dboard thought default real inside pi.

declare double precision, or if possible r, use modules. interface block can used declare dboard inside pi.


Comments

Post a Comment

Popular posts from this blog

javascript - Using jquery append to add option values into a select element not working -

i'm looking dynamically populate select element names of items stored array. for reason it's not working , i'm not sure why. html: <div class="col-md-4 col-md-offset-4"> <select class="select1"> </select> , <select class="select2"> </select> </div> note i'm using bootstrap layout. here's jquery i'm using try , populate ".select1" jquery: select: function() { $.each(state.greens, function(index, value) { $(".select1").append("<option value =' " + index + " '>" + state.greens[index].name + "</option>"); }); } note function 'select' stored within object called 'display'. state.greens name of array i'm trying access. so @ end of html page call display.select(); call function. no error messages showing in console. also, saw question: appending
Read more

Android soft keyboard reverts to default keyboard on orientation change -

i'm building android soft keyboard , can't seem fix bug - have arabic , qwerty keyboard , when rotate device on qwerty keyboard (or arabic shift), it's if program has "restarted" , becomes arabic keyboard without shift. the onsaveinstancestate(bundle savedinstancestate) not work because application not extend activity inputmethodservice . i put following in android manifest android:configchanges="keyboard|keyboardhidden|orientation" android:windowsoftinputmode="stateunchanged|adjustresize"> i tried using @override public void onconfigurationchanged(configuration newconfig) { super.onconfigurationchanged(newconfig); log.i(mydebug, "config changed " + currentkeyboard.equals(qwerty)); } however, currentkeyboard.equals(qwerty)) results false , made sure true before orientation change. any appreciated. i think applications fault. if application restarting on orientation change
Read more

jquery - javascript onscroll fade same class but with different div -

i new @ javascript , want learn new. there wrong code : why div fades @ same time. i want ask how fade different div s same class name.can me , explain me. thanks in advance. javascript: $(window).scroll(function(){ if ($( "div.fade" ).offset().top - $(window).scrolltop() <= 100){ $('.fade').addclass( "fade-in"); } else { $('.fade').removeclass("fade-in"); } }); here jsfiddle onscroll = function () { var scrolltop = $(this).scrolltop() ; $("div.fade").each(function(){ if(($(this).offset().top - scrolltop) <= 100){ $(this).toggleclass("fade-in"); } }) }
Read more