| version 1.3, 2000/08/22 05:04:04 |
version 1.4, 2001/02/21 07:10:18 |
|
|
| * 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/engine/PUM.c,v 1.2 2000/08/21 08:31:26 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/PUM.c,v 1.3 2000/08/22 05:04:04 noro Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| |
|
|
|
| |
|
| for ( i = n1; i >= n2; i-- ) { |
for ( i = n1; i >= n2; i-- ) { |
| if ( pw[i] ) { |
if ( pw[i] ) { |
| m = pw[i]; |
chsgnmp(mod,pw[i],&m); |
| for ( j = i; j >= 0; j-- ) { |
for ( j = i; j >= 0; j-- ) { |
| mulmp(vl,mod,pw[j],LC(p2),&m1); pw[j] = m1; |
mulmp(vl,mod,pw[j],LC(p2),&m1); pw[j] = m1; |
| } |
} |
| |
|
| for ( dc = DC(p2), d = i - n2; dc; dc = NEXT(dc) ) { |
for ( dc = DC(p2), d = i - n2; dc; dc = NEXT(dc) ) { |
| mulmp(vl,mod,COEF(dc),m,&m1); |
mulmp(vl,mod,COEF(dc),m,&m1); |
| submp(vl,mod,pw[QTOS(DEG(dc))+d],m1,&m2); |
addmp(vl,mod,pw[QTOS(DEG(dc))+d],m1,&m2); |
| pw[QTOS(DEG(dc))+d] = m2; |
pw[QTOS(DEG(dc))+d] = m2; |
| } |
} |
| } else |
} else |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to PUM.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / engine