version 1.1, 1999/12/03 07:39:07 |
version 1.5, 2000/11/15 00:07:20 |
|
|
/* $OpenXM: OpenXM/src/asir99/builtin/parif.c,v 1.2 1999/11/18 05:42:01 noro Exp $ */ |
/* |
|
* Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED |
|
* All rights reserved. |
|
* |
|
* FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited, |
|
* non-exclusive and royalty-free license to use, copy, modify and |
|
* redistribute, solely for non-commercial and non-profit purposes, the |
|
* computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and |
|
* conditions of this Agreement. For the avoidance of doubt, you acquire |
|
* only a limited right to use the SOFTWARE hereunder, and FLL or any |
|
* third party developer retains all rights, including but not limited to |
|
* copyrights, in and to the SOFTWARE. |
|
* |
|
* (1) FLL does not grant you a license in any way for commercial |
|
* purposes. You may use the SOFTWARE only for non-commercial and |
|
* non-profit purposes only, such as academic, research and internal |
|
* business use. |
|
* (2) The SOFTWARE is protected by the Copyright Law of Japan and |
|
* international copyright treaties. If you make copies of the SOFTWARE, |
|
* with or without modification, as permitted hereunder, you shall affix |
|
* to all such copies of the SOFTWARE the above copyright notice. |
|
* (3) An explicit reference to this SOFTWARE and its copyright owner |
|
* shall be made on your publication or presentation in any form of the |
|
* results obtained by use of the SOFTWARE. |
|
* (4) In the event that you modify the SOFTWARE, you shall notify FLL by |
|
* e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification |
|
* for such modification or the source code of the modified part of the |
|
* SOFTWARE. |
|
* |
|
* THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL |
|
* MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND |
|
* EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS |
|
* FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES' |
|
* RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY |
|
* MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY. |
|
* UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT, |
|
* OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY |
|
* DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL |
|
* DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES |
|
* ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES |
|
* FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY |
|
* DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF |
|
* SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART |
|
* OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY |
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
|
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
|
* |
|
* $OpenXM: OpenXM_contrib2/asir2000/builtin/parif.c,v 1.4 2000/08/22 05:03:59 noro Exp $ |
|
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
|
|
|
|
#endif |
#endif |
|
|
void Peval(),Psetprec(),p_pi(),p_e(),p_mul(),p_gcd(); |
void Peval(),Psetprec(),p_pi(),p_e(),p_mul(),p_gcd(); |
|
void asir_cgiv(GEN); |
|
|
struct ftab pari_tab[] = { |
struct ftab pari_tab[] = { |
{"eval",Peval,-2}, {"setprec",Psetprec,-1}, {0,0,0}, |
{"eval",Peval,-2}, {"setprec",Psetprec,-1}, {0,0,0}, |
Line 28 struct ftab pari_tab[] = { |
|
Line 77 struct ftab pari_tab[] = { |
|
#define MKPREC(a,i,b) (argc(a)==(i)?mkprec(QTOS((Q)(b))):prec) |
#define MKPREC(a,i,b) (argc(a)==(i)?mkprec(QTOS((Q)(b))):prec) |
|
|
#define CALLPARI1(f,a,p,r)\ |
#define CALLPARI1(f,a,p,r)\ |
ritopa((Obj)a,&_pt1_); _pt2_ = f(_pt1_,p); patori(_pt2_,r); cgiv(_pt2_); cgiv(_pt1_) |
ritopa((Obj)a,&_pt1_); _pt2_ = f(_pt1_,p); patori(_pt2_,r); asir_cgiv(_pt2_); asir_cgiv(_pt1_) |
#define CALLPARI2(f,a,b,p,r)\ |
#define CALLPARI2(f,a,b,p,r)\ |
ritopa((Obj)a,&_pt1_); ritopa((Obj)b,&_pt2_); _pt3_ = f(_pt1_,_pt2_,p); patori(_pt3_,r); cgiv(_pt3_); cgiv(_pt2_); cgiv(_pt1_) |
ritopa((Obj)a,&_pt1_); ritopa((Obj)b,&_pt2_); _pt3_ = f(_pt1_,_pt2_,p); patori(_pt3_,r); asir_cgiv(_pt3_); asir_cgiv(_pt2_); asir_cgiv(_pt1_) |
|
|
#define PARIF1P(f,pf)\ |
#define PARIF1P(f,pf)\ |
void f(NODE,Obj *);\ |
void f(NODE,Obj *);\ |
Line 48 void f(ar,rp) NODE ar; Obj *rp;\ |
|
Line 97 void f(ar,rp) NODE ar; Obj *rp;\ |
|
#define PREC_CONV 0.051905126 |
#define PREC_CONV 0.051905126 |
#endif |
#endif |
|
|
|
/* XXX : we should be more careful when we free PARI pointers. */ |
|
|
|
void asir_cgiv(ptr) |
|
GEN ptr; |
|
{ |
|
if ( ptr != gzero && ptr != gun |
|
&& ptr != gdeux && ptr != ghalf |
|
&& ptr != polvar && ptr != gi ) |
|
cgiv(ptr); |
|
} |
|
|
mkprec(p) |
mkprec(p) |
int p; |
int p; |
{ |
{ |
|
|
GEN x; |
GEN x; |
|
|
x = mppi(MKPREC(arg,1,ARG0(arg))); |
x = mppi(MKPREC(arg,1,ARG0(arg))); |
patori(x,rp); cgiv(x); |
patori(x,rp); asir_cgiv(x); |
} |
} |
|
|
void p_e(arg,rp) |
void p_e(arg,rp) |
|
|
{ |
{ |
GEN x; |
GEN x; |
|
|
x = gexp(gun,MKPREC(arg,1,ARG0(arg))); patori(x,rp); cgiv(x); |
x = gexp(gun,MKPREC(arg,1,ARG0(arg))); patori(x,rp); asir_cgiv(x); |
} |
} |
|
|
void p_mul(a,b,r) |
void p_mul(a,b,r) |
|
|
|
|
ritopa((Obj)a,&p1); ritopa((Obj)b,&p2); |
ritopa((Obj)a,&p1); ritopa((Obj)b,&p2); |
p3 = mulii(p1,p2); |
p3 = mulii(p1,p2); |
patori(p3,r); cgiv(p3); cgiv(p2); cgiv(p1); |
patori(p3,r); asir_cgiv(p3); asir_cgiv(p2); asir_cgiv(p1); |
} |
} |
|
|
void p_gcd(a,b,r) |
void p_gcd(a,b,r) |
|
|
|
|
ritopa_i(a,1,&p1); ritopa_i(b,1,&p2); |
ritopa_i(a,1,&p1); ritopa_i(b,1,&p2); |
p3 = mppgcd(p1,p2); |
p3 = mppgcd(p1,p2); |
patori_i(p3,r); cgiv(p3); cgiv(p2); cgiv(p1); |
patori_i(p3,r); asir_cgiv(p3); asir_cgiv(p2); asir_cgiv(p1); |
} |
} |
|
|
PARIF1P(p_sin,gsin) PARIF1P(p_cos,gcos) PARIF1P(p_tan,gtan) |
PARIF1P(p_sin,gsin) PARIF1P(p_cos,gcos) PARIF1P(p_tan,gtan) |
Line 314 struct ftab pari_tab[] = { |
|
Line 374 struct ftab pari_tab[] = { |
|
|
|
void parif_init() {} |
void parif_init() {} |
|
|
pointer evalparif() { |
pointer evalparif(f,arg) |
|
FUNC f; |
|
NODE arg; |
|
{ |
error("evalparif : PARI is not combined."); |
error("evalparif : PARI is not combined."); |
} |
} |
#endif /*PARI */ |
#endif /*PARI */ |