C------------------------------------------------------------------------------ C Example parallel program to show the use of the simple interface c for solving symmetric systems in the PWSMP library in 0-master mode. C------------------------------------------------------------------------------ C.. This program can be obtained from: C C http://www.cs.umn.edu/~agupta/wsmp C C (C) IBM Corporation, 2002. C c Acceptance and use of this program constitutes the user's understanding c that he/she will have no recourse to IBM for any actual or consequential c damages, including, but not limited to, lost profits or savings, arising c out of the use or inability to use this program. C------------------------------------------------------------------------------ program pzssmp_ex2 implicit none include 'mpif.h' integer mperr, myid call mpi_init(mperr) if (mperr .ne. 0) then print *,'MPI initialization error' stop end if call mpi_comm_rank(MPI_COMM_WORLD,myid,mperr) if (myid .eq. 0) then call master (myid) else call slave (myid) end if call mpi_finalize(mperr) stop end C------------------------------------------------------------------------------ subroutine master (myid) implicit none integer myid C.. This subroutine is called only on node 0. C.. It is very similar to program wssmp_ex2. integer n, ldb, nrhs, naux, nnzl, wspace, niter, info integer options(5) integer ia(10) integer ja(29) integer perm(9), invp(9) complex*16 avals(29) complex*16 b(9) double precision rcond complex*16 aux ! just placeholders in this program integer i double precision waltime, wsmprtc C.. Fill all arrays containing matrix data. data n /9/, ldb /9/, nrhs /1/, naux /0/ data ia /1,5,9,13,17,21,25,27,29,30/ data ja 1 /1, 3, 7, 8, 2 2, 3, 8, 9, 3 3, 7, 8, 9, 4 4, 6, 7, 8, 5 5, 6, 8, 9, 6 6, 7, 8, 9, 7 7, 8, 8 8, 9, 9 9/ data avals 1/(14.d0,0.d0), (-1.d0,0.d0), (-1.d0,0.d0), (-3.d0,0.d0), 2 (14.d0,0.d0), (-1.d0,0.d0), (-3.d0,0.d0), (-1.d0,0.d0), 3 (16.d0,0.d0), (-2.d0,0.d0), (-4.d0,0.d0), (-2.d0,0.d0), 4 (14.d0,0.d0), (-1.d0,0.d0), (-1.d0,0.d0), (-3.d0,0.d0), 5 (14.d0,0.d0), (-1.d0,0.d0), (-3.d0,0.d0), (-1.d0,0.d0), 6 (16.d0,0.d0), (-2.d0,0.d0), (-4.d0,0.d0), (-2.d0,0.d0), 7 (16.d0,0.d0), (-4.d0,0.d0), 8 (71.d0,0.d0), (-4.d0,0.d0), 9 (16.d0,0.d0)/ data options /5 * 0/ C.. Ordering and Symbolic. waltime = wsmprtc() call pzscalz (n, ia, ja, options, perm, invp, nnzl, wspace, + aux, naux, info) print *,'Analyze step complete in time - ',wsmprtc()-waltime if (info .ne. 0) then print *,'The following ERROR was detected: ',info return end if print *,'Number of nonzeros in factor L = ',nnzl print *,'Double words needed to factor on',myid,'=',wspace C.. Perform LDL' factorization. waltime = wsmprtc() call pzscldl (n, ia, ja, avals, perm, invp, aux, naux, info) print *,'LDL^T factorization done in time - ',wsmprtc()-waltime if (info .ne. 0) then print *,'The following ERROR was detected: ',info return end if C.. Back substitution. do i = 1, n b(i) = (1.d0,0.d0) end do waltime = wsmprtc() niter = 1 call pzsslv (n, perm, invp, b, ldb, nrhs, niter, aux, naux) print *,'Back substitution done in time - ',wsmprtc()-waltime print *,'The solution of the system is as follows:' do i = 1, n print *,i,' : ',b(i) end do C.. Solve the same system using the expert driver. waltime = wsmprtc() do i = 1, n b(i) = (1.d0,0.d0) end do rcond = 0.d0 call pzscsvx (n, ia, ja, avals, perm, invp, b, ldb, nrhs, + aux, naux, rcond, info) print *,'Run time for expert driver = ',wsmprtc() - waltime if (info .ne. 0) then print *,'The following ERROR was detected: ',info return end if print *,'The solution of the system is as follows:' do i = 1, n print *,i,' : ',b(i) end do return end C------------------------------------------------------------------------------ subroutine slave (myid) implicit none integer myid integer n, ldb, nrhs, naux, niter, wspace, nnzl, info integer ia, ja integer options(5) integer perm(9), invp(9) double precision rcond complex*16 avals, b ! just placeholders in this routine complex*16 aux ! just placeholders in this program data n /0/, ldb /1/, nrhs /1/, naux /0/ data options /5 * 0/ C.. Ordering and Symbolic. call pzscalz (n, ia, ja, options, perm, invp, nnzl, wspace, + aux, naux, info) if (info .ne. 0) return print *,'Double words needed to factor on',myid,'=',wspace C.. Perform LDL' factorization call pzscldl (n, ia, ja, avals, perm, invp, aux, naux, + info) if (info .ne. 0) return C.. Back substitution. niter = 1 call pzsslv (n, perm, invp, b, ldb, nrhs, niter, aux, naux) C.. Solve the same system using the expert driver. call pzscsvx (n, ia, ja, avals, perm, invp, b, ldb, nrhs, + aux, naux, rcond, info) if (info .ne. 0) return return end