Displaying 1 result from an estimated 1 matches for "ydot".
Did you mean:
dot
2004 Oct 06
1
Foreign code problem
..., hmin=0, hmax=Inf)
I get the error message:
Error in lsoda(y, times, "mond", parms, rtol, atol, tcrit = NULL,
jacfunc = NULL, :
Unable to find mond in mond
The C code is:
#include <math.h>
void mond(long int *neq, double *t, double *y, double *ydot)
{
double a = sqrt(y[2]*y[2] + y[4]*y[4]);
a = -1.0/(a*a*a);
ydot[0] = a*y[1];
ydot[1] = y[0];
ydot[2] = a*y[3];
ydot[3] = y[2];
}
and mond.so is build by:
gcc -O3 -shared mond.c -lm -o mond.so
I have set LD_LIBRARY_PATH to the correct directory and also tried to use the full path of mond...