OpenXM/src/ox_gsl/assert4.rr - view

Return to assert4.rr CVS log

Up to [local] / OpenXM / src / ox_gsl

File: [local] / OpenXM / src / ox_gsl / assert4.rr (download)

Revision 1.1, Wed Oct 23 07:00:43 2019 UTC (4 years, 6 months ago) by takayama
Branch: MAIN
CVS Tags: HEAD

gsl.rr is the lapper of ox_gsl.
Example.
 Ans=gsl.eigen_nonsymmv(M=matrix_list_to_matrix([[1,0],[3,4]]));

calls gsl_eigen_nonsymmv() of GSL.
It returns a list of (complex) eigenvalues and eigenvectors.

It is still in a testing mode and ox_gsl is loaded from OpenXM/src/ox_gsl/ox_gsl

/* $OpenXM: OpenXM/src/ox_gsl/assert4.rr,v 1.1 2019/10/23 07:00:43 takayama Exp $ */
#define myfail(s) printf("%a fails, Ans0=%a, Ans=%a\n",s,Ans0,Ans)
#define dcmp2(a,b,e) (deval(abs(1-a/b))>e?1:0)
Pid=ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_gsl/ox_gsl");
// Try 1.
ox_cmo_rpc(Pid,"gsl_eigen_nonsymmv",[1.0, 2.0, 3.0, 4.0]);
Ans2=ox_pop_cmo(Pid);
Ans1=ox_pop_cmo(Pid);
end$