Calling fortran from R via C does not link -

-

i have routine in fortran, wish use in r via c implementation. suppose fortran file bivnt.f. now, r cmd shlib c_binary.c , r cmd shlib bivnt.f not give problem when comes linking obtain

> dyn.load("bivnt.so") > dyn.load("c_binary.so") errore in dyn.load("c_binary.so") :    unable load shared object '/users/erlisruli/dropbox/wpapers/abc-cs/rcode/corbindata/c_binary.so':   dlopen(/users/erlisruli/dropbox/wpapers/abc-cs/rcode/corbindata/c_binary.so, 6): symbol not found: _smvbvt_   referenced from: /users/erlisruli/dropbox/wpapers/abc-cs/rcode/corbindata/c_binary.so   expected in: flat namespace  in /users/erlisruli/dropbox/wpapers/abc-cs/rcode/corbindata/c_binary.so

this seems pretty strange me, since nm obtain 

$ nm bivnt.so 0000000000003940 s _a.1868                  u _asin                  u _atan                  u _atan2                  u _exp 0000000000003540 t _mvbvn_ 00000000000027a0 t _mvbvt_ 0000000000002610 t _mvbvtc_ 0000000000002040 t _mvbvtl_ 0000000000002d10 t _mvbvu_ 0000000000001f20 t _mvphi_ 0000000000002b80 t _mvstdt_                  u _sin 00000000000038a0 t _smvbvt_ 0000000000003ba0 s _w.1368 0000000000003aa0 s _x.1369                  u dyld_stub_binder

any clue?

p.s. c_binary.c is

#include<r.h> /* calcolo ripartizione normale bivariata */ void f77_sub(smvbvt)( double* , int* , double*, double* , int* , double* ) ;  void biv_pmnorm( double *prob , int *nu, double *lower , double *upper , int *flag , double *corr ) {    f77_call(smvbvt)( prob , nu , lower , upper , flag , corr ) ; }  //pairwise log-likelihood-->output: units void bin_pairwise_llik( double *res , double *data , double *cov, double *theta , int *n , int *nc  ) /* n = #righe ; nc = #colonne*/ {     register int k, i, j, ite=0;     int nu=0, flag[2];     double res2=0.0 , res3=0.0, corr=theta[2], lower[2], upper[2];      for( k=0; k<*n; k++ )     {         for( i=0; i<(*nc-1); i++ )         {             for( j=(i+1); j<*nc; j++ )             {                  if(*(data + + k**nc)==1)                   {                      lower[0] = -*(theta+0)-*(theta+1)**(cov + + k**nc);                      upper[0] = 0.0;                      flag[0]  = 1;                  }                  else                  {                      lower[0] = 0.0;                      upper[0] = -*(theta+0)-*(theta+1)**(cov + + k**nc);                      flag[0]  = 0;                   }                   if(*(data + j + k**nc)==1)                  {                      lower[1] = -*(theta+0)-*(theta+1)**(cov + j + k**nc);                      upper[1] = 0.0;                      flag[1]  = 1;                  }                  else                  {                      lower[1] = 0.0;                      upper[1] = -*(theta+0)-*(theta+1)**(cov + j + k**nc);                      flag[1]  = 0;                  }                   biv_pmnorm( &res2, &nu, lower, upper, flag, &corr ) ;                  res3+=log(res2);                  ite += 1;             }          }          res[k]=res3;          res3=0.0;      } }   //pairwise log-likelihood-->output: 1 value void bin_pairwise_llik_sum( double *res , double *data , double *cov, double *theta , int *n , int *nc  ) /* n = #righe ; nc = #colonne*/ {     register int k, i, j, ite=0;     int nu=0, flag[2];     double res2=0.0 , corr=theta[2], lower[2], upper[2];     *res=0.0;      for( k=0; k<*n; k++ )     {         for( i=0; i<(*nc-1); i++ )         {             for( j=(i+1); j<*nc; j++ )             {                  if(*(data + + k**nc)==1)                   {                      lower[0] = -*(theta+0)-*(theta+1)**(cov + + k**nc);                      upper[0] = 0.0;                      flag[0]  = 1;                  }                  else                  {                      lower[0] = 0.0;                      upper[0] = -*(theta+0)-*(theta+1)**(cov + + k**nc);                      flag[0]  = 0;                   }                   if(*(data + j + k**nc)==1)                  {                      lower[1] = -*(theta+0)-*(theta+1)**(cov + j + k**nc);                      upper[1] = 0.0;                      flag[1]  = 1;                  }                  else                  {                      lower[1] = 0.0;                      upper[1] = -*(theta+0)-*(theta+1)**(cov + j + k**nc);                      flag[1]  = 0;                  }                   biv_pmnorm( &res2, &nu, lower, upper, flag, &corr ) ;                  *res+=log(res2);                    ite += 1;             }          }      } }


Comments

Post a Comment

Popular posts from this blog

c# - Unity IoC Lifetime per HttpRequest for UserStore -

-
i'm trying clean default implementation of accountcontroller.cs comes out of box in new mvc5/owin security implementation. have modified constructor this: private usermanager<applicationuser> usermanager; public accountcontroller(usermanager<applicationuser> usermanager) { this.usermanager = usermanager; } also, have created lifetime manager unity looks this: public class httpcontextlifetimemanager<t> : lifetimemanager, idisposable { private httpcontextbase _context = null; public httpcontextlifetimemanager() { _context = new httpcontextwrapper(httpcontext.current); } public httpcontextlifetimemanager(httpcontextbase context) { if (context == null) throw new argumentnullexception("context"); _context = context; } public void dispose() { this.removevalue(); } public override object get...
Read more

Change the color of an oval at click in Java AWT -

-
i have draw ovals in java, , @ click change color. beginning tried change color after 20 ms, doesn't work. my code is: public class mycomponentnew extends frame { public graphics2d g2d; public mycomponentnew(string title) { super(title); setsize(400, 550); } @override public void paint(graphics g) { this.g2d = (graphics2d) g; this.g2d.setcolor(color.red); this.g2d.filloval(10, 55, 50, 100); } public void changecolor () { this.g2d.setcolor(color.blue); this.g2d.filloval(10, 55, 50, 100); } } and in class main method have: mycomponentnew m; m = new mycomponentnew("fereastra cu baloane"); m.setvisible(true); m.addwindowlistener(new windowadapter() { @override public void windowclosing(windowevent we) { system.exit(0); } }); try { thread.sleep(20); } catch(interruptedexception e) {} m.changecolor(); the color of oval remains red. ...
Read more

I am trying to solve the error message 'incompatible ranks 0 and 1 in assignment' in a fortran 95 program. -

-
i writing program in fortran find velocity of parachuting person in relation time. keep getting error can't fix. new programming , appreciated. the error is v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) 1 error: incompatible ranks 0 , 1 in assignment @ (1) and program is program para integer :: real :: v(11) !velocity real :: q !initial velocity real :: h !time step real :: c !drag coefficient real :: m !mass ! gravity equal 32 ft/s^2 write (*,*)'enter time step' read(*,*)h write(*,*)'enter initial velocity' read(*,*)q write(*,*)'enter drag coefficient' read(*,*)c write(*,*)'enter mass' read(*,*)m i=1,10 ! 1 10, 1 being interval. end v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) q=v(1) end program you cannot use [] normal parenthesis in expressions. array constructor, [ items ] means array items elements. end do should after line.
Read more