| version 1.29, 2015/08/19 05:29:23 |
version 1.36, 2017/03/31 06:10:13 |
|
|
| /* $OpenXM: OpenXM_contrib2/asir2000/builtin/parif.c,v 1.28 2015/08/18 05:35:17 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/builtin/parif.c,v 1.35 2017/03/31 04:15:34 ohara Exp $ */ |
| #include "ca.h" |
#include "ca.h" |
| #include "parse.h" |
#include "parse.h" |
| #include "ox.h" |
#include "ox.h" |
|
|
| Q ox_pari_stream; |
Q ox_pari_stream; |
| int ox_pari_stream_initialized = 0; |
int ox_pari_stream_initialized = 0; |
| int ox_get_pari_result = 0; |
int ox_get_pari_result = 0; |
| |
P ox_pari_starting_function = 0; |
| |
|
| typedef void (*mpfr_func)(NODE,Obj *); |
typedef void (*mpfr_func)(NODE,Obj *); |
| |
|
| Line 63 Obj list_to_vect(Obj a) |
|
| Line 64 Obj list_to_vect(Obj a) |
|
| return (Obj)v; |
return (Obj)v; |
| } |
} |
| |
|
| |
Obj vect_to_mat(VECT v) |
| |
{ |
| |
MAT m; |
| |
int len,col,i,j; |
| |
|
| |
len = v->len; |
| |
if ( v->body[0] && OID((Obj)v->body[0]) == O_VECT ) { |
| |
col = ((VECT)v->body[0])->len; |
| |
for ( i = 1; i < len; i++ ) |
| |
if ( !v->body[i] || OID((Obj)v->body[i]) != O_VECT |
| |
|| ((VECT)v->body[i])->len != col ) |
| |
break; |
| |
if ( i == len ) { |
| |
/* convert to a matrix */ |
| |
MKMAT(m,len,col); |
| |
for ( i = 0; i < len; i++ ) |
| |
for ( j = 0; j < col; j++ ) |
| |
m->body[i][j] = ((VECT)v->body[i])->body[j]; |
| |
return (Obj)m; |
| |
} |
| |
} |
| |
return (Obj)v; |
| |
} |
| |
|
| void reset_ox_pari() |
void reset_ox_pari() |
| { |
{ |
| NODE nd; |
NODE nd; |
| Line 85 pointer evalparif(FUNC f,NODE arg) |
|
| Line 110 pointer evalparif(FUNC f,NODE arg) |
|
| USINT ui; |
USINT ui; |
| Obj ret,dmy; |
Obj ret,dmy; |
| mpfr_func mpfr_function; |
mpfr_func mpfr_function; |
| |
V v; |
| |
|
| if ( mpfr_function = mpfr_search(f->name) ) { |
if ( arg && ARG0(arg) && NID((Num)ARG0(arg)) != N_C |
| (*mpfr_function)(arg,&ret); |
&& (mpfr_function = mpfr_search(f->name)) ) { |
| return (pointer) ret; |
(*mpfr_function)(arg,&ret); |
| |
return (pointer) ret; |
| } |
} |
| |
|
| if ( !ox_pari_stream_initialized ) { |
if ( !ox_pari_stream_initialized ) { |
| |
if ( ox_pari_starting_function && OID(ox_pari_starting_function) == O_P ) { |
| |
v = VR(ox_pari_starting_function); |
| |
if ( (int)v->attr != V_SR ) { |
| |
error("pari : no handler."); |
| |
} |
| |
MKNODE(nd,0,0); |
| |
r = (Q)bevalf((FUNC)v->priv,0); |
| |
}else { |
| |
#if !defined(VISUAL) |
| MKSTR(name,"ox_pari"); |
MKSTR(name,"ox_pari"); |
| nd = mknode(2,NULL,name); |
nd = mknode(2,NULL,name); |
| Pox_launch_nox(nd,&r); |
Pox_launch_nox(nd,&r); |
| |
#else |
| |
error("Please load names.rr from latest asir-contrib library before using pari functions."); |
| |
#endif |
| |
} |
| ox_pari_stream = r; |
ox_pari_stream = r; |
| ox_pari_stream_initialized = 1; |
ox_pari_stream_initialized = 1; |
| } |
} |
| Line 117 pointer evalparif(FUNC f,NODE arg) |
|
| Line 157 pointer evalparif(FUNC f,NODE arg) |
|
| ox_get_pari_result = 1; |
ox_get_pari_result = 1; |
| Pox_pop_cmo(oxarg,&ret); |
Pox_pop_cmo(oxarg,&ret); |
| ox_get_pari_result = 0; |
ox_get_pari_result = 0; |
| if ( ret && OID(ret) == O_LIST ) |
if ( ret && OID(ret) == O_ERR ) { |
| |
char buf[BUFSIZ]; |
| |
soutput_init(buf); |
| |
sprintexpr(CO,((ERR)ret)->body); |
| |
error(buf); |
| |
} |
| |
if ( ret && OID(ret) == O_LIST ) { |
| ret = list_to_vect(ret); |
ret = list_to_vect(ret); |
| |
ret = vect_to_mat((VECT)ret); |
| |
} |
| return ret; |
return ret; |
| } |
} |
| |
|
| Line 130 struct pariftab { |
|
| Line 178 struct pariftab { |
|
| |
|
| /* |
/* |
| * type = 1 => argc = 1, second arg = precision |
* type = 1 => argc = 1, second arg = precision |
| * type = 2 => argc = 1, second arg = optional (long int) |
* type = 2 => argc = 1, second arg = (long int)0 |
| * |
* |
| */ |
*/ |
| /* |
/* |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to parif.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / builtin