| version 1.55, 2009/10/09 04:02:11 |
version 1.56, 2009/10/15 07:08:40 |
|
|
| * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
| * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
| * |
* |
| * $OpenXM: OpenXM_contrib2/asir2000/builtin/dp-supp.c,v 1.54 2009/06/01 07:31:54 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp-supp.c,v 1.55 2009/10/09 04:02:11 noro Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include "base.h" |
#include "base.h" |
| Line 2365 void qltozl(Q *w,int n,Q *dvr) |
|
| Line 2365 void qltozl(Q *w,int n,Q *dvr) |
|
| v.id = O_VECT; v.len = n; v.body = (pointer *)w; |
v.id = O_VECT; v.len = n; v.body = (pointer *)w; |
| igcdv(&v,dvr); return; |
igcdv(&v,dvr); return; |
| } |
} |
| c = w[0]; nm = NM(c); dn = INT(c) ? ONEN : DN(c); |
for ( i = 0; !w[i]; i++ ); |
| for ( i = 1; i < n; i++ ) { |
c = w[i]; nm = NM(c); dn = INT(c) ? ONEN : DN(c); |
| c = w[i]; l1 = INT(c) ? ONEN : DN(c); |
for ( i++; i < n; i++ ) { |
| |
c = w[i]; |
| |
if ( !c ) continue; |
| |
l1 = INT(c) ? ONEN : DN(c); |
| gcdn(nm,NM(c),&g); nm = g; |
gcdn(nm,NM(c),&g); nm = g; |
| gcdn(dn,l1,&l2); muln(dn,l1,&l3); divsn(l3,l2,&dn); |
gcdn(dn,l1,&l2); muln(dn,l1,&l3); divsn(l3,l2,&dn); |
| } |
} |