Powered by Apache

CVS log for OpenXM/src/ox_gsl/ox_gsl.c

[BACK] Up to [local] / OpenXM / src / ox_gsl

Request diff between arbitrary revisions


Default branch: MAIN


Revision 1.18 / (download) - annotate - [select for diffs], Wed Mar 10 06:51:57 2021 UTC (3 years, 1 month ago) by takayama
Branch: MAIN
CVS Tags: HEAD
Changes since 1.17: +2 -2 lines
Diff to previous 1.17 (colored)

gsl-config is used to set flags.

Revision 1.17 / (download) - annotate - [select for diffs], Wed Oct 23 07:00:43 2019 UTC (4 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.16: +3 -1 lines
Diff to previous 1.16 (colored)

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

Revision 1.16 / (download) - annotate - [select for diffs], Mon Oct 21 05:37:20 2019 UTC (4 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.15: +4 -1 lines
Diff to previous 1.15 (colored)

Added the incomplete gamma function \Gamma(a,x)
[1856] ox_cmo_rpc(Pid,"gsl_sf_gamma_inc",deval(1/2),deval(2));
0
[1857] ox_pop_cmo(Pid);
0.0806471

Revision 1.15 / (download) - annotate - [select for diffs], Fri Jun 8 00:03:43 2018 UTC (5 years, 10 months ago) by takayama
Branch: MAIN
Changes since 1.14: +10 -1 lines
Diff to previous 1.14 (colored)

ode solver: Example of y0'=y1,y1'=-4*x^2*y0+y1/x with y0(3.14/2)=1, y1(3.14/2)=0, h=1e-6,  x in [3.14/2,10]
Pid=ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_gsl/ox_gsl");
ox_cmo_rpc(Pid,"gsl_odeiv_step_rk4",[quote(y1),quote(-4*x^2*y0+y1/x)],
   [1.0,0], 3.14/2, 10.0, 1e-6,"rk4");
Ans=ox_pop_cmo(Pid);

Revision 1.14 / (download) - annotate - [select for diffs], Thu Jun 7 11:49:51 2018 UTC (5 years, 10 months ago) by takayama
Branch: MAIN
Changes since 1.13: +63 -1 lines
Diff to previous 1.13 (colored)

cmo2double() and cmo2double_list() are added.

Revision 1.13 / (download) - annotate - [select for diffs], Thu Jun 7 11:13:05 2018 UTC (5 years, 10 months ago) by takayama
Branch: MAIN
Changes since 1.12: +32 -1 lines
Diff to previous 1.12 (colored)

element_of_at(cmo *list,int k) returns the k-th element of the list.

Revision 1.12 / (download) - annotate - [select for diffs], Thu Jun 7 01:53:33 2018 UTC (5 years, 10 months ago) by takayama
Branch: MAIN
Changes since 1.11: +6 -2 lines
Diff to previous 1.11 (colored)

added monte_miser_integrate and monte_vegas_integrate.

Revision 1.11 / (download) - annotate - [select for diffs], Wed Jun 6 07:40:32 2018 UTC (5 years, 10 months ago) by takayama
Branch: MAIN
Changes since 1.10: +3 -1 lines
Diff to previous 1.10 (colored)

gsl_monte_plain_integrate is added, but it is still buggy.

Revision 1.10 / (download) - annotate - [select for diffs], Mon Jun 4 06:39:26 2018 UTC (5 years, 10 months ago) by ohara
Branch: MAIN
Changes since 1.9: +7 -1 lines
Diff to previous 1.9 (colored)

Since signal(3) of cygwin has a bug, our own implementation is used in cygwin environment.

Revision 1.9 / (download) - annotate - [select for diffs], Wed Apr 18 02:20:51 2018 UTC (6 years ago) by takayama
Branch: MAIN
Changes since 1.8: +2 -1 lines
Diff to previous 1.8 (colored)

Added some codes for debugging.

Revision 1.8 / (download) - annotate - [select for diffs], Tue Apr 17 00:56:38 2018 UTC (6 years ago) by takayama
Branch: MAIN
Changes since 1.7: +3 -1 lines
Diff to previous 1.7 (colored)

gsl_integration_qags (adaptive numerical integration) is added.
ox_cmo_rpc(Pid,"gsl_integration_qags", f(x),a,b) returns an approximate value of
int_a^b f(x) dx.

Example:
ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_gsl/ox_gsl");
0
[1817] ox_cmo_rpc(0,"gsl_integration_qags",quote(x^2*log(x)),0.0,1.0);
0
[1818] ox_pop_cmo(0);
[-0.111111,9.27977e-09,0]

Revision 1.7 / (download) - annotate - [select for diffs], Fri Apr 6 01:56:49 2018 UTC (6 years ago) by takayama
Branch: MAIN
Changes since 1.6: +5 -25 lines
Diff to previous 1.6 (colored)

Clean unnecessary codes.

Revision 1.6 / (download) - annotate - [select for diffs], Thu Apr 5 10:50:17 2018 UTC (6 years ago) by ohara
Branch: MAIN
Changes since 1.5: +5 -1 lines
Diff to previous 1.5 (colored)

Invoked init_dic().

Revision 1.5 / (download) - annotate - [select for diffs], Wed Apr 4 01:03:59 2018 UTC (6 years ago) by takayama
Branch: MAIN
Changes since 1.4: +35 -1 lines
Diff to previous 1.4 (colored)

ox_eval.c is linked to ox_gsl.
Pid=ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_gsl/ox_gsl");
ox_cmo_rpc(Pid,"test_ox_eval",quote(x+1));
ox_pop_cmo(Pid);
does not return the correct value. We need debugs.

Revision 1.4 / (download) - annotate - [select for diffs], Fri Mar 30 08:48:23 2018 UTC (6 years, 1 month ago) by takayama
Branch: MAIN
Changes since 1.3: +16 -14 lines
Diff to previous 1.3 (colored)

Fixed a bug on creating an error cmo.

Revision 1.3 / (download) - annotate - [select for diffs], Fri Mar 30 04:43:16 2018 UTC (6 years, 1 month ago) by takayama
Branch: MAIN
Changes since 1.2: +108 -21 lines
Diff to previous 1.2 (colored)

An error handling function myhander() is added.

Revision 1.2 / (download) - annotate - [select for diffs], Thu Mar 29 11:52:18 2018 UTC (6 years, 1 month ago) by takayama
Branch: MAIN
Changes since 1.1: +173 -19 lines
Diff to previous 1.1 (colored)

The complex gamma function implemented in GSL is installed.
Ex.
ox_launch(0,"/Users/nobuki/OX4/OpenXM/src/ox_gsl/ox_gsl");
0
[1916] ox_cmo_rpc(0,"gsl_sf_lngamma_complex_e",3,0); ox_pop_cmo(0);
0
[1917] [0.693147,0,0]   [log( |Gamma(z)| ), arg( Gamma(z)), status]
[1918] ox_cmo_rpc(0,"gsl_sf_lngamma_complex_e",3,2); ox_pop_cmo(0);
0
[1919] [4.44406e+15,2.32392,0]
CVS ----------------------------------------------------------------------

Revision 1.1 / (download) - annotate - [select for diffs], Thu Mar 29 05:47:11 2018 UTC (6 years, 1 month ago) by takayama
Branch: MAIN

The initial skelton for ox_gsl
ex.
Pid=ox_launch(0,Path_of_ox_gsl+"/ox_gsl");
ox_cmo_rpc(Pid,"add_int32",2,3);
ox_pop_cmo(Pid);

This form allows you to request diff's between any two revisions of a file. You may select a symbolic revision name using the selection box or you may type in a numeric name using the type-in text box.




FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>