Asir-contrib 2018 ( ) asir contrib update() import("names.rr"); asir_contrib_update( update=1); Win: AppData/Roaming/OpenXM/lib/asir-contrib. :.OpenXM

Size: px

Start display at page:

Download "Asir-contrib 2018 ( ) asir contrib update() import("names.rr"); asir_contrib_update( update=1); Win: AppData/Roaming/OpenXM/lib/asir-contrib. :.OpenXM"

自磐闻人
4 years ago
Views:

1 Asir-contrib 2018 ( ) asir contrib update() import("names.rr"); asir_contrib_update( update=1); Win: AppData/Roaming/OpenXM/lib/asir-contrib. :.OpenXM/OpenXM/lib/asir-contrib. update=2 : helph() document. update=4 : win ox pari, ox sm1. [check] 6:30AM. 1 / 15

2 module module sample; static A; localf init; localf foo; def init() { A=2; } def foo(n) { return(a^n); } endmodule; sample.init()$ sample.foo(8); module localf : openxm oxlistlocalf : http: // html-ja/man/man_21.html#g_t_00e3_0083_00a2_00e3_0082_ 00b8_00e3_0083_00a5_00e3_0083_00bc_00e3_0083_00ab 2 / 15

3 4, 5 testing/noro/grcrt.rr ( ). chinese remainder Th GB? [1828] load("grcrt.rr"); [1885] noro_grcrt.f4_chr([x^2+y^2-4,x*y-1],[x,y],0); [0][1][2] [CRT,2.5e-05,IR, ,#P,1] [x+y^3-4*y,x^2+y^2-4,y*x-1] noro_grcrt.f4_chr(cyclic(5),vars(cyclic(5)),0); 3 / 15

4 6, 7 ox gsl (tk, oh, ). GNU scientific library ox. : (π 3 (1 cos(x 0 ) cos(x 1 ) cos(x 2 ))) 1 [0, π] 3. Pid=ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_gsl/ox_gsl"); ox_cmo_rpc(pid,"gsl_monte_vegas_integrate", quote(1/(3.14^3*(1-cos(x0)*cos(x1)*cos(x2)))), [0.0,0.0,0.0], // [3.14,3.14,3.14], // 10000); // Ans=ox_pop_cmo(Pid); ( ). 4 / 15

5 6, 7 ox pari. ox pari ( ) gtt ds.rr (gtt=goto,tachibana,takayama, ( )), tk ds ahg.rr [1822] load("gtt_ds.rr"); [2720] gtt_ds.direct_sampler([[4,14,3],[10,6,5]], [[1,9/10,11/10],[1,13/10,99/100],[1,1,1]]); [ ] [ ] [ ] [2721] gtt_ds.direct_sampler([[4,14,3],[10,6,5]], [[1,9/10,11/10],[1,13/10,99/100],[1,1,1]]); [ ] [ ] [ ] openxm listusage gtt ds.rr 5 / 15

6 8, 9 os muldif.rr (os=oshima, commit). other-docs/muldif/changelog.pdf, xypoch ( ) mkdir -p ~/asir/tex /asir/tex/risaout.tex package.. \usepackage[all]{xy} \usepackage{tikz} import("os_muldif.rr"); os_md.dviout("\\begin{xy} "+os_md.xypoch(70,3,7,10 ar=0) +"\\end{xy}"); // ~/asir/tex/out.tex ( ) os_md.dviout(""); 6 / 15

7 /asir/tex/out.tex 7 / 15

8 8, 9 os muldif.rr ( commit). taylorode ( ) os_md.show(os_md.taylorode(2 series=h,dif=1,runge=2,c1=1))$ 1 2 h(((2f y b 2 a 21 f y )h + 2b 1 + 2b 2 2)f + (2f x b 1 c 1 + 2f x b 2 c 2 f x )h) (1) /home/web/openxm/current/doc/other-docs/cron-other-docs.sh (10:30AM). ox python ( ) Pid=ox_launch(0,getenv("OpenXM_HOME")+"/src/ox_python/ox_python")$ ox_cmo_rpc(pid,"pyrun_string","str(eval( 1+2 ))")$ printf("result of PyRun_String str(eval( 1+2 )) = %a\n", ox_pop_cmo(pid))$ ox_cmo_rpc(pid,"eval","1+2")$ printf("result of eval 1+2= %a\n",ox_pop_cmo(pid))$ 8 / 15

9 10, 11 os muldif ( ), OpenXM/src/sage/asir.py for SAGE ( ), 9 / 15

10 10, 11, tk approx-r.rr ( ) cont frac, tk rk-mfac.rr ( ), ( Runge-Kutta, for ec project), dp monomial hilbert poincare ( , , asir2018) [1830] poly_hilbert_polynomial([x^2*y^2,y^2*z^2] v=[x,y,z]); [not_implemeted_yet,2*n+5,[ ],h^6-2*h^4+1,3] [1831] poly_hilbert_polynomial([x^2*y^2,y^2*z^2] v=[x,y,z],sm1=1); [h^2+6*h-7,2*h+5] 10 / 15

11 12, 1 noro pd.rr ( ), I J nd f4 trace(). asir-contrib/ja/noro_pd-html/noro_pd-ja.html A=[j*h*g*f*e*d*b,j*i*g*d*c*b,j*i*h*g*d*b,j*i*h*e*b,i*e*c*b, B=[a*d-j*c,b*c,d*e-f*g*h]$ V=[a,b,c,d,e,f,g,h,i,j,z]$ noro_pd.gbf4(1)$ noro_pd.ideal_intersection(a,b,v,0); cgiasir.sm1 ( ), 11 / 15

12 2 3 os muldif, solveeq, baseode ( ), c2csp --> L=os_md.baseODE([y-x^2-t/2,2*x*y+b]); [x2-2*x^3-t*x-b, ] ODE, ODE. x2 = x x = y x 2 t/2, y = 2 x y + b x 2x 3 tx b (P 2 ). gtt ekn3.rr ( ). 2. contiguity, igcdctrl asir2018, binary splitting / 15

13 2 3 gtt ekn3.rr, 3 3 binary splitting. [4617] gtt_ekn3.assert2(15 bs=1); Marginal=[[1950,2550,5295],[1350,1785,6660]] P=[[17/100,1,10],[7/50,1,33/10],[1,1,1]] gtt_ekn3: Timing (int) =4.644 (CPU) + 0 (GC) = (total), real time= C-c C-cinterrupt?(q/t/c/d/u/w/?) t Abort this computation? (y or n) y return to toplevel [4618] gtt_ekn3.assert2(15); Marginal=[[1950,2550,5295],[1350,1785,6660]] P=[[17/100,1,10],[7/50,1,33/10],[1,1,1]] gtt_ekn3: Timing (int) =6.736 (CPU) + 0 (GC) = (total), real time= C-c C-cinterrupt?(q/t/c/d/u/w/?) t 13 / 15

14 2 3 gtt ekn3.rr Time(s) ,, : [[4*N,4*N,4*N,4*N,4*N], [2*N,3*N,5*N,5*N,5*N]]. N = 20 asir2000, igcdcntl(0) 101s asir2000, igcdcntl(3) 91.6s asir2018, bs, newpath 119s asir2018, newpath 50.1s N [1828] load("gtt_ekn3.rr"); [3757] gtt_ekn3.test5x5(20 bs=1); [[1,1/2,1/3,1/5,1/7],[1,1/11,1/13,1/17,1/19], [1,1/23,1/29,1/31,1/37],[1,1/37,1/41,1/43,1/47],[1,1,1,1,1]] Binary splitting: M(t + n 1)M(t + n 2) M(t). M(t + i). 14 / 15

15 2 3 gtt ekn3.rr ( ). Time(s) ,, : [[4*N,4*N,4*N,4*N,4*N], [2*N,3*N,5*N,5*N,5*N]] N [1828] load("gtt_ekn3.rr"); [3757] gtt_ekn3.test5x5(20); // bs [[1,1/2,1/3,1/5,1/7],[1,1/11,1/13,1/17,1/19], [1,1/23,1/29,1/31,1/37],[1,1/37,1/41,1/43,1/47],[1,1,1,1,1]] 15 / 15

哈尔滨应用职业技术学院2012年教育教学质量报告

哈尔滨应用职业技术学院2012年教育教学质量报告哈尔滨应用职业技术学院高等职业教育质量年度报告 (2015) 哈尔滨应用职业技术学院 2015 年 1 月目录一学院简介... 1 二办学条件... 1 1. 办学资源... 1 2. 办学经费... 2 3. 师资队伍... 3 三学院教学... 3 1. 专业建设... 3 2.