OpenXM/src/asir-port/cgi/r-fd.rr - diff

Return to r-fd.rr CVS log

Up to [local] / OpenXM / src / asir-port / cgi

version 1.1, 2014/07/27 03:18:44

version 1.3, 2014/07/31 07:00:41

Line 1

/* $OpenXM: OpenXM/src/asir-port/cgi/r-fd.rr,v 1.2 2014/07/31 06:24:11 takayama Exp $ */

load("tk_fd.rr")$

import("tk_r.rr")$

import("oh_number.rr")$

/* r_d2rat(0.3) --> precision loss in truncation if not ctrl("bigfloat",1) */

ctrl("bigfloat",1)$

def r_d2rat(Y) {

if ((type(Y) ==4)||(type(Y)==5)||(type(Y)==6)) return map(r_d2rat,Y);

if ((type(Y) == 1) && (ntype(Y) >= 1)) {

S = rtostr(Y); Y = "eval(("+S+")*exp(0));";

/* print(Y); */

Y = eval_str(Y);

return oh_number.rats(Y);

}else return Y;

}

def r_ahvec(A,B,C,Y) {

Y = r_d2rat(Y);

R=tk_fd.ygahvec(A,B,C,Y|all=1);

Gamma=R[1];

Der=R[0];

Line 10 def r_ahvec(A,B,C,Y) {

Line 26 def r_ahvec(A,B,C,Y) {

Ans=cons(Z,Der2);

/* Fans=map(rtostr,map(tk_fd.tk_number_rattofloat,Ans)); */

Fans=map(deval,Ans);

Fans = tk_r.asir2r_c(Fans);

return Fans;

}

end$