===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/io/ox.c,v
retrieving revision 1.21
retrieving revision 1.25
diff -u -p -r1.21 -r1.25
--- OpenXM_contrib2/asir2000/io/ox.c	2003/12/09 03:07:45	1.21
+++ OpenXM_contrib2/asir2000/io/ox.c	2004/08/28 12:50:30	1.25
@@ -44,7 +44,7 @@
  * 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/io/ox.c,v 1.20 2003/03/07 03:12:28 noro Exp $
+ * $OpenXM: OpenXM_contrib2/asir2000/io/ox.c,v 1.24 2003/12/11 05:48:04 noro Exp $
 */
 #include "ca.h"
 #include "parse.h"
@@ -56,6 +56,7 @@
 
 extern Obj VOIDobj;
 
+extern int nserver_102, myrank_102;
 extern int ox_need_conv;
 int ox_usr1_sent, ox_int_received, critical_when_signal;
 unsigned int ox_serial;
@@ -119,6 +120,7 @@ static int available_cmo[] = {
 	CMO_ZZ, CMO_QQ, CMO_ZERO,
 	CMO_DMS_GENERIC, CMO_DMS_OF_N_VARIABLES,
 	CMO_RING_BY_NAME, CMO_DISTRIBUTED_POLYNOMIAL,
+    CMO_RATIONAL,
 	CMO_RECURSIVE_POLYNOMIAL, CMO_UNIVARIATE_POLYNOMIAL,
 	CMO_INDETERMINATE,
 	CMO_TREE,
@@ -145,6 +147,7 @@ static int ox_asir_available_sm[] = {
 	SM_mathcap, SM_setMathcap, SM_nop,
 	SM_beginBlock, SM_endBlock,
 	SM_set_rank_102, SM_tcp_accept_102, SM_tcp_connect_102, SM_reset_102,
+	SM_bcast_102, SM_reduce_102,
 	0
 };
 
@@ -529,6 +532,60 @@ void ox_send_data_102(int rank,pointer p)
 	ox_write_cmo_102(rank,p);
 	ox_flush_stream_102(rank);
 	end_critical();
+}
+
+void ox_bcast_102(int root)
+{
+	Obj data;
+	int r,mask,id,src,dst;
+
+	r = myrank_102-root;
+	if ( r == 0 )
+		data = (Obj)asir_pop_one();
+		
+	if ( r < 0 ) r += nserver_102;
+	for ( mask = 1; mask < nserver_102; mask <<= 1 )
+		if ( r&mask ) {
+			src = myrank_102-mask;
+			if ( src < 0 ) src += nserver_102;
+			ox_recv_102(src,&id,&data);
+			break;
+		}
+	for ( mask >>= 1; mask > 0; mask >>= 1 )
+		if ( (r+mask) < nserver_102 ) {
+			dst = myrank_102+mask;
+			if ( dst >= nserver_102 ) dst -= nserver_102;
+			ox_send_data_102(dst,data);
+		}
+	asir_push_one(data);
+}
+
+/* func : an arithmetic funcion func(vl,a,b,*c) */
+
+void ox_reduce_102(int root,void (*func)())
+{
+	Obj data,data0,t;
+	int r,mask,id,src,dst;
+
+	r = myrank_102-root;
+	if ( r < 0 ) r += nserver_102;
+	data = (Obj)asir_pop_one();
+	for ( mask = 1; mask < nserver_102; mask <<= 1 )
+		if ( r&mask ) {
+			dst = (r-mask)+root;
+			if ( dst >= nserver_102 ) dst -= nserver_102;
+			ox_send_data_102(dst,data);
+			break;
+		} else {
+			src = r+mask;
+			if ( src < nserver_102 ) {
+				src += root;
+				if ( src >= nserver_102 ) src -= nserver_102;
+				ox_recv_102(src,&id,&data0);
+				(*func)(CO,data,data0,&t); data = t;
+			}
+		}
+	asir_push_one(r?0:data);
 }
 
 void ox_send_cmd(int s,int id)