version 1.3, 2004/10/06 11:58:51 |
version 1.5, 2018/03/29 01:32:50 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/builtin/user.c,v 1.2 2004/09/29 08:50:23 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/builtin/user.c,v 1.4 2015/09/03 23:05:36 noro Exp $ */ |
|
|
/* a sample file for adding builtin functions */ |
/* a sample file for adding builtin functions */ |
|
|
|
|
void Ppartial_derivative(); |
void Ppartial_derivative(); |
void partial_derivative(VL vl,P p,V v,P *r); |
void partial_derivative(VL vl,P p,V v,P *r); |
void Pzadd(),Pzsub(),Pzmul(); |
void Pzadd(),Pzsub(),Pzmul(); |
|
void Pcomp_f(); |
|
|
struct ftab user_tab[] = { |
struct ftab user_tab[] = { |
/* |
/* |
Line 16 struct ftab user_tab[] = { |
|
Line 17 struct ftab user_tab[] = { |
|
{"zadd",Pzadd,2}, |
{"zadd",Pzadd,2}, |
{"zsub",Pzsub,2}, |
{"zsub",Pzsub,2}, |
{"zmul",Pzmul,2}, |
{"zmul",Pzmul,2}, |
|
{"comp_f",Pcomp_f,2}, |
{0,0,0}, |
{0,0,0}, |
}; |
}; |
|
|
|
/* compare two [[F,M],...] */ |
|
|
|
void Pcomp_f(NODE arg,Q *rp) |
|
{ |
|
NODE l1,l2,e1,e2; |
|
int m1,m2,r; |
|
|
|
l1 = BDY((LIST)ARG0(arg)); |
|
l2 = BDY((LIST)ARG1(arg)); |
|
for ( ; l1 && l2; l1 = NEXT(l1), l2 = NEXT(l2) ) { |
|
e1 = BDY((LIST)BDY(l1)); |
|
e2 = BDY((LIST)BDY(l2)); |
|
r = compp(CO,(P)ARG0(e1),(P)ARG0(e2)); |
|
if ( r ) { STOQ(r,*rp); return; } |
|
m1 = QTOS((Q)ARG1(e1)); |
|
m2 = QTOS((Q)ARG1(e2)); |
|
r = m1>m2?1:(m1<m2?-1:0); |
|
if ( r ) { STOQ(r,*rp); return; } |
|
} |
|
r = l1?1:(l2?-1:0); |
|
STOQ(r,*rp); |
|
} |
|
|
void Pzadd(NODE arg,Q *rp) |
void Pzadd(NODE arg,Q *rp) |
{ |
{ |
Z z0,z1,z2; |
Z z0,z1,z2; |
|
|
z0 = qtoz((Q)ARG0(arg)); |
z0 = qtoz((Q)ARG0(arg)); |
z1 = qtoz((Q)ARG1(arg)); |
z1 = qtoz((Q)ARG1(arg)); |
z2 = addz(z0,z1); |
z2 = addz(z0,z1); |
printz(z2); printf(" "); |
printz(z2); printf(" "); |
*rp = ztoq(z2); |
*rp = ztoq(z2); |
} |
} |
|
|
void Pzsub(NODE arg,Q *rp) |
void Pzsub(NODE arg,Q *rp) |
{ |
{ |
Z z0,z1,z2; |
Z z0,z1,z2; |
|
|
z0 = qtoz((Q)ARG0(arg)); |
z0 = qtoz((Q)ARG0(arg)); |
z1 = qtoz((Q)ARG1(arg)); |
z1 = qtoz((Q)ARG1(arg)); |
z2 = subz(z0,z1); |
z2 = subz(z0,z1); |
printz(z2); printf(" "); |
printz(z2); printf(" "); |
*rp = ztoq(z2); |
*rp = ztoq(z2); |
} |
} |
|
|
void Pzmul(NODE arg,Q *rp) |
void Pzmul(NODE arg,Q *rp) |
{ |
{ |
Z z0,z1,z2; |
Z z0,z1,z2; |
|
|
z0 = qtoz((Q)ARG0(arg)); |
z0 = qtoz((Q)ARG0(arg)); |
z1 = qtoz((Q)ARG1(arg)); |
z1 = qtoz((Q)ARG1(arg)); |
z2 = mulz(z0,z1); |
z2 = mulz(z0,z1); |
printz(z2); printf(" "); |
printz(z2); printf(" "); |
*rp = ztoq(z2); |
*rp = ztoq(z2); |
} |
} |
|
|
/* |
/* |