===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/N.c,v
retrieving revision 1.6
retrieving revision 1.9
diff -u -p -r1.6 -r1.9
--- OpenXM_contrib2/asir2000/engine/N.c	2003/09/02 07:00:51	1.6
+++ OpenXM_contrib2/asir2000/engine/N.c	2007/09/15 10:17:08	1.9
@@ -45,7 +45,7 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2000/engine/N.c,v 1.5 2001/10/09 01:36:10 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/engine/N.c,v 1.8 2007/02/18 05:36:27 ohara Exp $ 
 */
 #include "ca.h"
 #include "base.h"
@@ -98,6 +98,7 @@ void addn(N n1,N n2,N *nr)
 		}
 #else
 		asm volatile("\
+		pushl	%%ebx;\
 		movl	%1,%%esi;\
 		movl	%2,%%edi;\
 		movl	%3,%%ebx;\
@@ -115,10 +116,11 @@ void addn(N n1,N n2,N *nr)
 		jnz Lstart_addn;\
 		movl	$0,%%eax;\
 		adcl	%%eax,%%eax;\
-		movl	%%eax,%0"\
+		movl	%%eax,%0;\
+		popl	%%ebx"\
 		:"=m"(c)\
 		:"m"(m1),"m"(m2),"m"(mr),"m"(d2)\
-		:"eax","ebx","ecx","edx","esi","edi");
+		:"eax","ecx","edx","esi","edi");
 #endif
 		for ( i = d2, m1 += d2, mr += d2; (i < d1) && c ; i++ ) {
 			tmp = *m1++ + c;
@@ -202,6 +204,7 @@ int subn(N n1,N n2,N *nr)
 		}
 #else
 		asm volatile("\
+		pushl	%%ebx;\
 		movl	%1,%%esi;\
 		movl	%2,%%edi;\
 		movl	%3,%%ebx;\
@@ -219,10 +222,11 @@ int subn(N n1,N n2,N *nr)
 		jnz Lstart_subn;\
 		movl	$0,%%eax;\
 		adcl	%%eax,%%eax;\
-		movl	%%eax,%0"\
+		movl	%%eax,%0;\
+		popl	%%ebx"\
 		:"=m"(br)\
 		:"m"(m1),"m"(m2),"m"(mr),"m"(d2)\
-		:"eax","ebx","ecx","edx","esi","edi");
+		:"eax","ecx","edx","esi","edi");
 #endif
 		for ( i = d2, m1 += d2, mr += d2; (i < d1) && br; i++ ) {
 			t = *m1++;
@@ -284,6 +288,7 @@ void _addn(N n1,N n2,N nr)
 		}
 #else
 		asm volatile("\
+		pushl	%%ebx;\
 		movl	%1,%%esi;\
 		movl	%2,%%edi;\
 		movl	%3,%%ebx;\
@@ -301,10 +306,11 @@ void _addn(N n1,N n2,N nr)
 		jnz Lstart__addn;\
 		movl	$0,%%eax;\
 		adcl	%%eax,%%eax;\
-		movl	%%eax,%0"\
+		movl	%%eax,%0;\
+		popl	%%ebx"\
 		:"=m"(c)\
 		:"m"(m1),"m"(m2),"m"(mr),"m"(d2)\
-		:"eax","ebx","ecx","edx","esi","edi");
+		:"eax","ecx","edx","esi","edi");
 #endif
 		for ( i = d2, m1 += d2, mr += d2; (i < d1) && c ; i++ ) {
 			tmp = *m1++ + c;
@@ -387,6 +393,7 @@ int _subn(N n1,N n2,N nr)
 		}
 #else
 		asm volatile("\
+		pushl	%%ebx;\
 		movl	%1,%%esi;\
 		movl	%2,%%edi;\
 		movl	%3,%%ebx;\
@@ -404,10 +411,11 @@ int _subn(N n1,N n2,N nr)
 		jnz Lstart__subn;\
 		movl	$0,%%eax;\
 		adcl	%%eax,%%eax;\
-		movl	%%eax,%0"\
+		movl	%%eax,%0;\
+		popl	%%ebx"\
 		:"=m"(br)\
 		:"m"(m1),"m"(m2),"m"(mr),"m"(d2)\
-		:"eax","ebx","ecx","edx","esi","edi");
+		:"eax","ecx","edx","esi","edi");
 #endif
 		for ( i = d2, m1 += d2, mr += d2; (i < d1) && br; i++ ) {
 			t = *m1++;
@@ -669,6 +677,15 @@ void pwrn(N n,int e,N *nr)
 extern int igcd_algorithm;
 
 void gcdEuclidn(), gcdn_HMEXT();
+
+void lcmn(N n1,N n2,N *nr)
+{
+	N g,t;
+
+	gcdn(n1,n2,&g);
+	divsn(n1,g,&t);
+	muln(t,n2,nr);
+}
 
 void gcdn(N n1,N n2,N *nr)
 {