Return to assert1.rr CVS log

Up to [local] / OpenXM / src / ox_gsl

Annotation of OpenXM/src/ox_gsl/assert1.rr, Revision 1.2

1.1       takayama    1: #define myfail(s) printf("%a fails, Ans0=%a, Ans=%a\n",s,Ans0,Ans)
                      2: #define dcmp(a,b) (deval(abs(1-a/b))<1e-10?1:0)
                      3: Pid=ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_gsl/ox_gsl");
                      4: // Try 1.
                      5: ox_cmo_rpc(Pid,"gsl_sf_lngamma_complex_e",10^5,2); Ans=ox_pop_cmo(Pid);
                      6: Ans0=[1.05129e+06,-2.1069,0]$
                      7: if (dcmp(Ans[0],Ans0[0])) myfail("Try1 real"); else ;
                      8: if (dcmp(Ans[1],Ans0[1])) myfail("Try1 im"); else ;
                      9: // Try 2.
                     10: ox_cmo_rpc(Pid,"restart");
                     11: // Try 3.  error packet should be returned.
                     12: ox_cmo_rpc(Pid,"gsl_sf_lngamma_complex_e",10^10,2); Ans=ox_pop_cmo(Pid);
                     13: // Try 4.
                     14: ox_cmo_rpc(Pid,"gsl_sf_lngamma_complex_e",deval(10^10),2); Ans=ox_pop_cmo(Pid);
                     15: Ans0=[2.20259e+11,2.0694,0]$
                     16: if (dcmp(Ans[0],Ans0[0])) myfail("Try4 real"); else ;
                     17: if (dcmp(Ans[1],Ans0[1])) myfail("Try4 im"); else ;
1.2     ! takayama   18: // Try 5.
        !            19: ox_cmo_rpc(Pid,"gsl_integration_qags",quote(log(x)/x^(1/2)),0.01,1);
        !            20: Ans=ox_pop_cmo(Pid);
1.1       takayama   21: end$