diff --git a/test/newfmt/Makefile b/test/newfmt/Makefile deleted file mode 100644 index 6e641d18..00000000 --- a/test/newfmt/Makefile +++ /dev/null @@ -1,33 +0,0 @@ -include ../../Make.inc -# -# Libraries used -# -LIBDIR=../../lib/ -OPTDIR=../../opt -PSBLAS_LIB= -L$(LIBDIR) -lpsb_util -lpsb_krylov -lpsb_prec -lpsb_base -L$(OPTDIR) -lopt\ - -lspgpu -L/usr/local/spgpu/lib -lcudart -lcublas -L/usr/local/cuda/lib64 -LDLIBS=$(PSBLDLIBS) -# -# Compilers and such -# -CCOPT= -g -FINCLUDES=$(FMFLAG)$(LIBDIR) $(FMFLAG). $(FMFLAG)$(OPTDIR) - - -EXEDIR=./runs - -all: ppde - -ppde: ppde.o - $(F90LINK) ppde.o -o ppde $(PSBLAS_LIB) $(LDLIBS) -lgomp - /bin/mv ppde $(EXEDIR) - -clean: - /bin/rm -f ppde.o spde.o $(EXEDIR)/ppde -verycleanlib: - (cd ../..; make veryclean) -lib: - (cd ../../; make library) - - - diff --git a/test/newfmt/ppde.F90 b/test/newfmt/ppde.F90 deleted file mode 100644 index c2443250..00000000 --- a/test/newfmt/ppde.F90 +++ /dev/null @@ -1,403 +0,0 @@ -!!$ -!!$ Parallel Sparse BLAS version 3.1 -!!$ (C) Copyright 2006, 2007, 2008, 2009, 2010, 2012, 2013 -!!$ Salvatore Filippone University of Rome Tor Vergata -!!$ Alfredo Buttari CNRS-IRIT, Toulouse -!!$ -!!$ Redistribution and use in source and binary forms, with or without -!!$ modification, are permitted provided that the following conditions -!!$ are met: -!!$ 1. Redistributions of source code must retain the above copyright -!!$ notice, this list of conditions and the following disclaimer. -!!$ 2. Redistributions in binary form must reproduce the above copyright -!!$ notice, this list of conditions, and the following disclaimer in the -!!$ documentation and/or other materials provided with the distribution. -!!$ 3. The name of the PSBLAS group or the names of its contributors may -!!$ not be used to endorse or promote products derived from this -!!$ software without specific written permission. -!!$ -!!$ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS -!!$ ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED -!!$ TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR -!!$ PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS -!!$ BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR -!!$ CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF -!!$ SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS -!!$ INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN -!!$ CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) -!!$ ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE -!!$ POSSIBILITY OF SUCH DAMAGE. -!!$ -!!$ -! File: ppde.f90 -! -! Program: ppde -! This sample program solves a linear system obtained by discretizing a -! PDE with Dirichlet BCs. -! -! -! The PDE is a general second order equation in 3d -! -! a1 dd(u) a2 dd(u) a3 dd(u) b1 d(u) b2 d(u) b3 d(u) -! - ------ - ------ - ------ + ----- + ------ + ------ + c u = f -! dxdx dydy dzdz dx dy dz -! -! with Dirichlet boundary conditions -! u = g -! -! on the unit cube 0<=x,y,z<=1. -! -! -! Note that if b1=b2=b3=c=0., the PDE is the Laplace equation. -! -! In this sample program the index space of the discretized -! computational domain is first numbered sequentially in a standard way, -! then the corresponding vector is distributed according to a BLOCK -! data distribution. -! -program ppde3d - use psb_base_mod - use psb_prec_mod - use psb_krylov_mod - use psb_util_mod - use psb_d_ell_mat_mod - use psb_d_rsb_mat_mod - implicit none - - ! input parameters - character(len=20) :: kmethd, ptype - character(len=5) :: afmt - integer(psb_ipk_) :: idim - - ! miscellaneous - real(psb_dpk_), parameter :: one = 1.d0 - real(psb_dpk_) :: t1, t2, tprec - - ! sparse matrix and preconditioner - type(psb_dspmat_type) :: a,bm - type(psb_dprec_type) :: prec - type(psb_d_ell_sparse_mat) :: aell - type(psb_d_rsb_sparse_mat) :: arsb - ! descriptor - type(psb_desc_type) :: desc_a - ! dense vectors - type(psb_d_vect_type) :: xxv,bv - ! parallel environment - integer(psb_ipk_) :: ictxt, iam, np - - ! solver parameters - integer(psb_ipk_) :: iter, itmax,itrace, istopc, irst - integer(psb_long_int_k_) :: amatsize, precsize, descsize - real(psb_dpk_) :: err, eps - - ! other variables - integer(psb_ipk_) :: info, i - character(len=20) :: name,ch_err, fname - - info=psb_success_ - - - call psb_init(ictxt) - call psb_info(ictxt,iam,np) - info=psb_rsb_matmod_init() - if(info.ne.psb_success_)then - info=psb_err_from_subroutine_ - ch_err='create_matrix: rsb module initialization failed' - call psb_errpush(info,name,a_err=ch_err) - endif - - if (iam < 0) then - ! This should not happen, but just in case - call psb_exit(ictxt) - stop - endif - if(psb_get_errstatus() /= 0) goto 9999 - name='pde3d90' - call psb_set_errverbosity(itwo) - ! - ! Hello world - ! - if (iam == psb_root_) then - write(*,*) 'Welcome to PSBLAS version: ',psb_version_string_ - write(*,*) 'This is the ',trim(name),' sample program' - end if - ! - ! get parameters - ! - call get_parms(ictxt,kmethd,ptype,afmt,idim,istopc,itmax,itrace,irst) - - ! - ! allocate and fill in the coefficient matrix, rhs and initial guess - ! - call psb_barrier(ictxt) - t1 = psb_wtime() - call psb_gen_pde3d(ictxt,idim,a,bv,xxv,desc_a,afmt,& - & a1,a2,a3,b1,b2,b3,c,g,info,amold=arsb) - call psb_barrier(ictxt) - t2 = psb_wtime() - t1 - if(info /= psb_success_) then - info=psb_err_from_subroutine_ - ch_err='psb_gen_pde3d' - call psb_errpush(info,name,a_err=ch_err) - goto 9999 - end if - if (iam == psb_root_) write(psb_out_unit,'("Overall matrix creation time : ",es12.5)')t2 - if (iam == psb_root_) write(psb_out_unit,'(" ")') - ! - ! prepare the preconditioner. - ! - if(iam == psb_root_) write(psb_out_unit,'("Setting preconditioner to : ",a)')ptype - call psb_precinit(prec,ptype,info) - - call psb_barrier(ictxt) - t1 = psb_wtime() - call psb_precbld(a,desc_a,prec,info,amold=arsb) - if(info /= psb_success_) then - info=psb_err_from_subroutine_ - ch_err='psb_precbld' - call psb_errpush(info,name,a_err=ch_err) - goto 9999 - end if - - tprec = psb_wtime()-t1 - - call psb_amx(ictxt,tprec) - - if (iam == psb_root_) write(psb_out_unit,'("Preconditioner time : ",es12.5)')tprec - if (iam == psb_root_) write(psb_out_unit,'(" ")') - ! - ! iterative method parameters - ! - if(iam == psb_root_) write(psb_out_unit,'("Calling iterative method ",a)')kmethd - call psb_barrier(ictxt) - t1 = psb_wtime() - eps = 1.d-9 - call psb_krylov(kmethd,a,prec,bv,xxv,eps,desc_a,info,& - & itmax=itmax,iter=iter,err=err,itrace=itrace,istop=istopc,irst=irst) - - if(info /= psb_success_) then - info=psb_err_from_subroutine_ - ch_err='solver routine' - call psb_errpush(info,name,a_err=ch_err) - goto 9999 - end if - - call psb_barrier(ictxt) - t2 = psb_wtime() - t1 - call psb_amx(ictxt,t2) - amatsize = a%sizeof() - descsize = desc_a%sizeof() - precsize = prec%sizeof() - call psb_sum(ictxt,amatsize) - call psb_sum(ictxt,descsize) - call psb_sum(ictxt,precsize) - - if (iam == psb_root_) then - write(psb_out_unit,'(" ")') - write(psb_out_unit,'("Time to solve matrix : ",es12.5)')t2 - write(psb_out_unit,'("Time per iteration : ",es12.5)')t2/iter - write(psb_out_unit,'("Number of iterations : ",i0)')iter - write(psb_out_unit,'("Convergence indicator on exit : ",es12.5)')err - write(psb_out_unit,'("Info on exit : ",i0)')info - write(psb_out_unit,'("Total memory occupation for A: ",i12)')amatsize - write(psb_out_unit,'("Total memory occupation for PREC: ",i12)')precsize - write(psb_out_unit,'("Total memory occupation for DESC_A: ",i12)')descsize - write(psb_out_unit,'("Storage type for DESC_A: ",a)') desc_a%indxmap%get_fmt() - end if - - ! - ! cleanup storage and exit - ! - call psb_gefree(bv,desc_a,info) - call psb_gefree(xxv,desc_a,info) - call psb_spfree(a,desc_a,info) - call psb_precfree(prec,info) - call psb_cdfree(desc_a,info) - if(info /= psb_success_) then - info=psb_err_from_subroutine_ - ch_err='free routine' - call psb_errpush(info,name,a_err=ch_err) - goto 9999 - end if - -9999 continue - if(info /= psb_success_) then - call psb_error(ictxt) - end if - - info=psb_rsb_matmod_exit() - if(info.ne.psb_success_)then - info=psb_err_from_subroutine_ - ch_err='create_matrix: rsb module finalization failed' - call psb_errpush(info,name,a_err=ch_err) - endif - - call psb_exit(ictxt) - stop - -contains - ! - ! get iteration parameters from standard input - ! - subroutine get_parms(ictxt,kmethd,ptype,afmt,idim,istopc,itmax,itrace,irst) - integer(psb_ipk_) :: ictxt - character(len=*) :: kmethd, ptype, afmt - integer(psb_ipk_) :: idim, istopc,itmax,itrace,irst - integer(psb_ipk_) :: np, iam - integer(psb_ipk_) :: intbuf(10), ip - - call psb_info(ictxt, iam, np) - - if (iam == 0) then - read(psb_inp_unit,*) ip - if (ip >= 3) then - read(psb_inp_unit,*) kmethd - read(psb_inp_unit,*) ptype - read(psb_inp_unit,*) afmt - - ! broadcast parameters to all processors - call psb_bcast(ictxt,kmethd) - call psb_bcast(ictxt,afmt) - call psb_bcast(ictxt,ptype) - - - read(psb_inp_unit,*) idim - if (ip >= 4) then - read(psb_inp_unit,*) istopc - else - istopc=1 - endif - if (ip >= 5) then - read(psb_inp_unit,*) itmax - else - itmax=500 - endif - if (ip >= 6) then - read(psb_inp_unit,*) itrace - else - itrace=-1 - endif - if (ip >= 7) then - read(psb_inp_unit,*) irst - else - irst=1 - endif - ! broadcast parameters to all processors - - intbuf(1) = idim - intbuf(2) = istopc - intbuf(3) = itmax - intbuf(4) = itrace - intbuf(5) = irst - call psb_bcast(ictxt,intbuf(1:5)) - - write(psb_out_unit,'("Solving matrix : ell1")') - write(psb_out_unit,& - & '("Grid dimensions : ",i4," x ",i4," x ",i4)') & - & idim,idim,idim - write(psb_out_unit,'("Number of processors : ",i0)')np - write(psb_out_unit,'("Data distribution : BLOCK")') - write(psb_out_unit,'("Preconditioner : ",a)') ptype - write(psb_out_unit,'("Iterative method : ",a)') kmethd - write(psb_out_unit,'(" ")') - else - ! wrong number of parameter, print an error message and exit - call pr_usage(izero) - call psb_abort(ictxt) - stop 1 - endif - else - call psb_bcast(ictxt,kmethd) - call psb_bcast(ictxt,afmt) - call psb_bcast(ictxt,ptype) - call psb_bcast(ictxt,intbuf(1:5)) - idim = intbuf(1) - istopc = intbuf(2) - itmax = intbuf(3) - itrace = intbuf(4) - irst = intbuf(5) - end if - return - - end subroutine get_parms - ! - ! print an error message - ! - subroutine pr_usage(iout) - integer(psb_ipk_) :: iout - write(iout,*)'incorrect parameter(s) found' - write(iout,*)' usage: pde3d90 methd prec dim & - &[istop itmax itrace]' - write(iout,*)' where:' - write(iout,*)' methd: cgstab cgs rgmres bicgstabl' - write(iout,*)' prec : bjac diag none' - write(iout,*)' dim number of points along each axis' - write(iout,*)' the size of the resulting linear ' - write(iout,*)' system is dim**3' - write(iout,*)' istop stopping criterion 1, 2 ' - write(iout,*)' itmax maximum number of iterations [500] ' - write(iout,*)' itrace <=0 (no tracing, default) or ' - write(iout,*)' >= 1 do tracing every itrace' - write(iout,*)' iterations ' - end subroutine pr_usage - - ! - ! functions parametrizing the differential equation - ! - function b1(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: b1 - real(psb_dpk_), intent(in) :: x,y,z - b1=1.d0/sqrt(3.d0) - end function b1 - function b2(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: b2 - real(psb_dpk_), intent(in) :: x,y,z - b2=1.d0/sqrt(3.d0) - end function b2 - function b3(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: b3 - real(psb_dpk_), intent(in) :: x,y,z - b3=1.d0/sqrt(3.d0) - end function b3 - function c(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: c - real(psb_dpk_), intent(in) :: x,y,z - c=0.d0 - end function c - function a1(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: a1 - real(psb_dpk_), intent(in) :: x,y,z - a1=1.d0/80 - end function a1 - function a2(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: a2 - real(psb_dpk_), intent(in) :: x,y,z - a2=1.d0/80 - end function a2 - function a3(x,y,z) - use psb_base_mod, only : psb_dpk_ - real(psb_dpk_) :: a3 - real(psb_dpk_), intent(in) :: x,y,z - a3=1.d0/80 - end function a3 - function g(x,y,z) - use psb_base_mod, only : psb_dpk_, done - real(psb_dpk_) :: g - real(psb_dpk_), intent(in) :: x,y,z - g = dzero - if (x == done) then - g = done - else if (x == dzero) then - g = exp(y**2-z**2) - end if - end function g - -end program ppde3d - - diff --git a/test/newfmt/runs/Makefile b/test/newfmt/runs/Makefile deleted file mode 100644 index d8fa43ac..00000000 --- a/test/newfmt/runs/Makefile +++ /dev/null @@ -1,12 +0,0 @@ -lib: - (cd ..; $(MAKE) lib) -clean: - (cd ..; $(MAKE) clean) -verycleanlib: - (cd ..; $(MAKE) verycleanlib) -ppde90: - (cd ..; $(MAKE) ppde90) -pp2d: - (cd ..; $(MAKE) pp2d) - -.PHONY: pp2d ppde90 diff --git a/test/newfmt/runs/ppde.inp b/test/newfmt/runs/ppde.inp deleted file mode 100644 index 762c4095..00000000 --- a/test/newfmt/runs/ppde.inp +++ /dev/null @@ -1,11 +0,0 @@ -7 Number of entries below this -BICGSTAB Iterative method BICGSTAB CGS BICG BICGSTABL RGMRES -BJAC Preconditioner NONE DIAG BJAC -CSR Storage format for matrix A: CSR COO JAD -040 Domain size (acutal system is this**3) -2 Stopping criterion -0100 MAXIT -01 ITRACE -20 IRST restart for RGMRES and BiCGSTABL - -