Parallel Programming
on the Macintosh
Timothy H. Kaiser
tkaiser@sdsc.edu
University of New Mexico
High Performance Computing
Education and Research Center
Albuquerque Resource Center
#include#include #include #include "p4.h" double f(a) double a; { return (4.0 / (1.0 + a*a)); } main(int argc, char ** argv) { double h,sum,x,mypi; int ijk,type,from,size,myid,n,numprocs,jid; ijk=p4_initenv(&argc,argv) ijk=p4_create_procgroup() myid=p4_get_my_id(); numprocs=p4_num_total_ids(); if(myid == 0){ n=1000; printf("enter the number of iterations: "); scanf("%d",&n); for (ijk=1;ijk<=numprocs-1;ijk++) p4_send(1,ijk,(char*)&n,sizeof(int)); } else { type=-1; from=-1; size=sizeof(int); p4_recv(&type,&from,(char**)&n,&size); } h = 1.0 / (double) n; sum = 0.0; for (ijk = myid + 1; ijk <= n; ijk += numprocs) { x = h * ((double)ijk - 0.5); sum += f(x); } mypi = h * sum; type=2; size=sizeof(double); if(myid == 0){ for (ijk=1;ijk<=numprocs-1;ijk++){ jid=ijk; p4_recv(&type,&jid,(char**)&x,&size); mypi=mypi+x; } printf("pi = %g\n",mypi); } else { p4_send(type,0,(char*)&mypi,size); } p4_wait_for_end(); }
#include "mpi.h" #include#include double f(a) double a; { return (4.0 / (1.0 + a*a)); } int main(argc,argv) int argc; char *argv[]; { int n, myid, numprocs, i, rc; double mypi, pi, h, sum, x; MPI_Init(&argc,&argv); MPI_Comm_size(MPI_COMM_WORLD,&numprocs); MPI_Comm_rank(MPI_COMM_WORLD,&myid); if (myid == 0) n=10000; MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD); h = 1.0 / (double) n; sum = 0.0; for (i = myid + 1; i <= n; i += numprocs){ x = h * ((double)i - 0.5); sum += f(x); } mypi = h * sum; MPI_Reduce(&mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD); if (myid == 0) printf("pi is approximately %.16f\n",pi); MPI_Finalize(); }
mpich/configure remove top blank lines mpich/src/env/Makefile.in $(CC) $(CFLAGS) -c -DCONFIGURE_ARGS_CLEAN \ mpich/mpid/ch_p4/p4-1.4/lib/p4_secure.c /* extern char *sys_errlist[]; */ mpich/mpid/ch_p4/p4-1.4/lib/p4_debug.c #define VPRINTF mpich/mpid/ch_p4/p4init.c /* #include*/ #if defined(solaris) || defined(HAVE_SYSINFO) sysinfo( SI_SRPC_DOMAIN,name+l,nlen-l); #elif defined(HAVE_GETDOMAINNAME) getdomainname( name+l, nlen - l ); #endif becomes getdomainname( name+l, nlen - l );
P4_HOME_DIR =/home/tkaiser/p4-1.4 INCLUDEDIR = $(P4_HOME_DIR)/include LIBDIR = $(P4_HOME_DIR)/lib LIBS = $(LIBDIR)/libp4.a $(MDEP_LIBS) CLINKER = cc CFLAGS = -g -I$(INCLUDEDIR) $(MDEP_CFLAGS) $(USER_CFLAGS) MDEP_LIBS = -lrpc /home/tkaiser/mpich/lib/NeXT/ch_p4/xdr_float.o .c.o:$(P) $(CC) $(CFLAGS) -c $*.c bcast:$(P) bcast.o $(LIBDIR)/libp4.a $(CLINKER) $(CFLAGS) -o bcast bcast.o $(LIBS)
SHELL = /bin/sh ARCH = NeXT COMM = ch_p4 MPIR_HOME = /home/tkaiser/mpich CC = cc CLINKER = cc LIB_PATH = -L$(MPIR_HOME)/lib/$(ARCH)/$(COMM) LIB_LIST = -lmpi -lrpc $(MPIR_HOME)/lib/$(ARCH)/$(COMM)/*.o INCLUDE_DIR = -I$(MPIR_HOME)/include CFLAGS = -DMPID_NO_FORTRAN -DHAS_XDR=1 -DSTDC_HEADERS=1 -DHAVE_STDLIB_H=1 -DMALLOC_RET_VOID=1 -DHAVE_SYSTEM=1 -DHAVE_NICE=1 -DHAVE_LONG_DOUBLE=1 -DHAVE_LONG_LONG_INT=1 $(OPTFLAGS) $(INCLUDE_DIR) -DMPI_$(ARCH) MPILIB = $(MPIR_HOME)/lib/$(ARCH)/$(COMM)/libmpi.a LIBS = $(LIB_PATH) $(LIB_LIST) cpi: cpi.o $(MPIR_HOME)/include/mpir.h $(MPILIB) $(CLINKER) $(OPTFLAGS) -o cpi cpi.o \ $(LIB_PATH) $(LIB_LIST) -lm .c.o: $(CC) $(CFLAGS) -c $*.c
local 0 osage.arc.unm.edu 1 /home/tkaiser/p4-1.4/contrib/bcast cochiti.arc.unm.edu 1 /home/tkaiser/p4/SGI/bcast acoma.arc.unm.edu 1 /home/tkaiser/p4/RS6000/bcast
alice classic.arc.unm.edu hd:pi.68k ozark ozark.arc.unm.edu pb520hd:pi.68k osage osage.arc.unm.edu hd:pi.PPC
-- script to send a preference file and launch a remote job
set filespec to choose file with prompt "Select Process Group File"
set p5_file to (open for access filespec without write permission)
set mac_list to ""
set ftp_list to ""
set app_list to {}
set mlist to ""
set the_count to 0
set theline to "continue"
repeat while theline ""
try
set mac_name to " "
repeat while mac_name = " "
set mac_name to read p5_file until " "
end repeat
set ftp_name to " "
repeat while ftp_name = " "
set ftp_name to read p5_file until " "
end repeat
set app_name to read p5_file before return as string
--display dialog app_name
repeat while character 1 of app_name is " "
set z to number of characters of app_name
set app_name to subString(2, z, app_name)
end repeat
set mlist to mlist & mac_name & " " & ftp_name & " " & app_name & return
if theline "" then
set the_count to the_count + 1
set mac_list to mac_list & " " & mac_name
set ftp_list to ftp_list & " " & ftp_name
set app_list to app_list & app_name
end if
on error
set theline to ""
end try
end repeat
set launched to 0
repeat with i from 1 to the_count
set machine_name to word i of mac_list
set app_name to "" & item i of app_list
set launched to launched + open_remote(machine_name, app_name, mlist)
end repeat
if launched the_count then
repeat with i from 1 to the_count - launched
beep
end repeat
end if
on open_remote(machine_name, app_path, machine_list)
tell application "Finder" of machine machine_name
activate
try
set the_path to path to preferences folder as string
set the_path to the_path & "p4.pg"
set p4_file to (open for access file the_path with write permission)
set eof p4_file to 0
write machine_list to p4_file
close access p4_file
open file app_path
set worked to 1
on error
set worked to 0
end try
end tell
return worked
end open_remote
on subString(x, y, theString)
return (characters x thru y of theString) as string
end subString
p4_initenv(&argc,argv); Init p4 p4_create_procgroup(); Init p4 p4_get_my_id(); Get my id p4_num_total_ids(); Number of processors p4_send(1,ijk,(char*)&n,sizeof(int)); Nonblocking send p4_recv(&type,&from,(char**)&n,&size); Blocking Receive p4_wait_for_end(); Termination