diff --git a/Makefile b/Makefile index 7d76992b..05d90238 100644 --- a/Makefile +++ b/Makefile @@ -32,7 +32,7 @@ lib: mpobjs $(OBJS) $(AR) $(HERE)/$(LIBNAME) $(OBJS) $(RANLIB) $(HERE)/$(LIBNAME) /bin/cp -p $(HERE)/$(LIBNAME) $(LIBDIR) - /bin/cp -p $(LIBMOD) $(LIBDIR) + /bin/cp -p $(LIBMOD) $(LOCAL_MODS) $(LIBDIR) $(F90OBJS) $(MPFOBJS): $(MODOBJS) psb_prec_mod.o: psb_prec_type.o diff --git a/README b/README index de4ad17a..79205a2c 100644 --- a/README +++ b/README @@ -1,129 +1,18 @@ -This directory contains the PSBLAS library, version 2.1.0 +This directory contains the MLD2P4 set of preconditioners. - -Version 1.0 of the library was described in: -S. Filippone, M. Colajanni -PSBLAS: A library for parallel linear algebra computation on sparse matrices -ACM Trans. on Math. Software, 26(4), Dec. 2000, pp. 527-550. - -PLATFORMS: - -The compilation process relies on the choice of an appropriate -Make.inc file; we have tested with AIX XLF, Intel ifc/Linux, Lahey -F95/Linux, Nag f95/Linux, GNU Fortran/Linux. If you succeed in compiling with -other compiler/operating systems please let us know. - - -LINUX: - -There finally exist a GNU Fortran 95 implementation: we are using the -development snapshots from GCC 3.5.0, later 4.1 and 4.2 since July -2004, and it appears to work. The 4.2 version of GNU compilers is now -our reference platform; the official 4.2.0 release is due pretty -soon. It now includes support for ALLOCATABLES. - - -The Lahey version we got access to (6.0 and 6.1) seems to suffer from -spurious extra copies problem; this is most apparent in the matrix -build process. - -For the Intel compilers, we recommend moving to version 9; previous -versions of the library have been compiled with version 7 and 8 of -ifort. - -IBM SP. -The library has been tested on an IBM SP2, SP4 and SP5, with XLC and XLF -compilers, and a version of the BLACS based on MPI. -The setting - F90=xlf90 -qsuffix=f=f90 -in Make.inc.rs6k takes care of the f90 extension. -WARNING: xlf 8.1 introduced a performance bug, whereas a Fortan 90 -code calling a Fortan 77 code would incur spurious array copies; -please make sure your system has the PTF xlf 8102 installed. - - - -UTILITIES -The TEST/Fileread directory contains some utilities to convert to/from -Harwell-Boeing and MatrixMarket file formats. - - -DOCUMENTATION - -See userguide.pdf -Please consult the sample programs, especially TEST/pargen/ppde90.f90. - - -OTHER SOFTWARE CREDITS - -We include our modified implementation of some of the Sparker (serial -sparse BLAS) material, e.g. Jagged diagonal, plus a number of -extensions of our own design. The original file spblas.f can be -downloaded from matisa.cc.rl.ac.uk; of course any bugs in our -implementation are our own to fix. The main reference for the serial -sparse BLAS is: -Duff, I., Marrone, M., Radicati, G., and Vittoli, C. -Level 3 basic linear algebra subprograms for sparse matrices: a user -level interface -ACM Trans. Math. Softw., 23(3), 379-401, 1997. - -In the multilevel preconditioners we use SMMP by Randolph E. Bank and -Craig C. Douglas na.bank@na-net.ornl.gov and -na.cdouglas@na-net.ornl.gov; we wrapped it in a Fortran 95 interface -with dynamic memory allocation. - -To compile and run our software you will need - -1. A working version of MPI - -2. The MPI version of the BLACS from - http://www.netlib.org/blacs/index.html - -3. A version of the BLAS; if you don't have a specific version for - your platform you may try ATLAS available from - http://math-atlas.sourceforge.net/ - -4. We have had good results with the METIS library, from - http://www-users.cs.umn.edu/~karypis/metis/metis/main.html - This is not necessary to compile our library, but the test program - in test/Fileread assumes you have it installed. - -5. For our preconditioners we include interfaces to the following - software packages: - -- SuperLU 3.0 http://crd.lbl.gov/~xiaoye/SuperLU/ - -- UMFPACK 4.4 http://www.cise.ufl.edu/research/sparse/umfpack/ - These are optional, you only need to install them if you actually - want to use them. - - -TODO: -- As soon as TR 15581 and MOVE_ALLOC become available in GNU Fortran - we shall move from POINTERS to ALLOCATABLE for both vectors and - derived data types components. -- The GLOB_TO_LOC array should be changed for large test cases. - - -The PSBLAS team. - - -Contact: Salvatore Filippone salvatore.filippone@uniroma2.it - -Credits for version 2.0: -Salvatore Filippone -Alfredo Buttari - -The MD2P4 multilevel parallel preconditioners contained in directory -src/prec were developed with the contribution of: +developed by +Salvatore Filippone Pasqua D'Ambra Daniela di Serafino -They are still in an early experimental stage, use at your own risk! -Credits for version 1.0: -Salvatore Filippone -Michele Colajanni -Fabio Cerioni -Stefano Maiolatesi -Dario Pascucci + +To compile: +1. Edit Makefile and point PSBLASDIR to the main directory containing + the PSBLAS source code +2. Edit the Makefile in the main PSBLAS dir modifying the PREC + variable to make it point to this directory, e.g. + PREC=../mld2p4 +3. Run make from the main PSBLAS directory. diff --git a/psb_dasmatbld.f90 b/psb_dasmatbld.f90 index a15740f2..65b8fda3 100644 --- a/psb_dasmatbld.f90 +++ b/psb_dasmatbld.f90 @@ -138,7 +138,7 @@ Subroutine psb_dasmatbld(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt) goto 9999 endif - if (novr == 0) then + if ((novr == 0).or.(np==1)) then ! ! This is really just Block Jacobi..... ! @@ -176,16 +176,16 @@ Subroutine psb_dasmatbld(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt) ! ! Build the auiliary descriptor',desc_p%matrix_data(psb_n_row_) ! - call psb_cdbldovr(a,desc_data,novr,desc_p,info) + call psb_cdbldext(a,desc_data,novr,desc_p,info,extype=psb_ovt_asov_) if(info /= 0) then info=4010 - ch_err='psb_cdbldovr' + ch_err='psb_cdbldext' call psb_errpush(info,name,a_err=ch_err) goto 9999 end if Endif - if(debug) write(0,*) me,' From cdbldovr _:',desc_p%matrix_data(psb_n_row_),& + if(debug) write(0,*) me,' From cdbldext _:',desc_p%matrix_data(psb_n_row_),& & desc_p%matrix_data(psb_n_col_) @@ -202,11 +202,11 @@ Subroutine psb_dasmatbld(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt) !!$ blk%k = n_row if (present(outfmt)) then - if(debug) write(0,*) me,': Calling SPHALO with ',size(blk%ia2) - Call psb_sphalo(a,desc_p,blk,info,outfmt=outfmt) + if(debug) write(0,*) me,': Calling outfmt SPHALO with ',size(blk%ia2) + Call psb_sphalo(a,desc_p,blk,info,outfmt=outfmt,data=psb_comm_ext_) else if(debug) write(0,*) me,': Calling SPHALO with ',size(blk%ia2) - Call psb_sphalo(a,desc_p,blk,info) + Call psb_sphalo(a,desc_p,blk,info,data=psb_comm_ext_) end if @@ -235,7 +235,7 @@ Subroutine psb_dasmatbld(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dbaseprc_aply.f90 b/psb_dbaseprc_aply.f90 index 9bf9012b..fcc220a5 100644 --- a/psb_dbaseprc_aply.f90 +++ b/psb_dbaseprc_aply.f90 @@ -191,7 +191,7 @@ subroutine psb_dbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) tx(desc_data%matrix_data(psb_n_row_)+1:isz) = dzero if (prec%iprcparm(restr_)==psb_halo_) then - call psb_halo(tx,prec%desc_data,info,work=aux) + call psb_halo(tx,prec%desc_data,info,work=aux,data=psb_comm_ext_) if(info /=0) then info=4010 ch_err='psb_halo' @@ -271,7 +271,7 @@ subroutine psb_dbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) 9999 continue call psb_errpush(info,name,i_err=int_err,a_err=ch_err) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dbaseprc_bld.f90 b/psb_dbaseprc_bld.f90 index 6b5c326a..55e24695 100644 --- a/psb_dbaseprc_bld.f90 +++ b/psb_dbaseprc_bld.f90 @@ -257,7 +257,7 @@ subroutine psb_dbaseprc_bld(a,desc_a,p,info,upd) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dbjac_aply.f90 b/psb_dbjac_aply.f90 index f77dcafd..ab351b93 100644 --- a/psb_dbjac_aply.f90 +++ b/psb_dbjac_aply.f90 @@ -260,7 +260,7 @@ subroutine psb_dbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) 9999 continue call psb_errpush(info,name,i_err=int_err,a_err=ch_err) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dbldaggrmat.f90 b/psb_dbldaggrmat.f90 index 9f599385..0801bafd 100644 --- a/psb_dbldaggrmat.f90 +++ b/psb_dbldaggrmat.f90 @@ -89,7 +89,7 @@ subroutine psb_dbldaggrmat(a,desc_a,ac,desc_ac,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -265,8 +265,9 @@ contains !if(.not.associated(p%av(ap_nd_)%aspk)) p%iprcparm(jac_sweeps_) = 1 !------------------------------------------------------------------ ! Split AC=M+N N off-diagonal part + nzl = psb_sp_get_nnzeros(ac) call psb_sp_all(ac%m,ac%k,p%av(ap_nd_),nzl,info) - if(info /= 0) then + if (info /= 0) then call psb_errpush(4010,name,a_err='psb_sp_all') goto 9999 end if @@ -324,7 +325,7 @@ contains 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -1031,7 +1032,7 @@ contains 9999 continue call psb_errpush(info,name) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_ddiagsc_bld.f90 b/psb_ddiagsc_bld.f90 index 7c3c6eb0..0e01356d 100644 --- a/psb_ddiagsc_bld.f90 +++ b/psb_ddiagsc_bld.f90 @@ -158,7 +158,7 @@ subroutine psb_ddiagsc_bld(a,desc_a,p,upd,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dgenaggrmap.f90 b/psb_dgenaggrmap.f90 index 245bb453..adc6dc7b 100644 --- a/psb_dgenaggrmap.f90 +++ b/psb_dgenaggrmap.f90 @@ -283,7 +283,7 @@ subroutine psb_dgenaggrmap(aggr_type,a,desc_a,nlaggr,ilaggr,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dilu_bld.f90 b/psb_dilu_bld.f90 index 6744155f..7e22984d 100644 --- a/psb_dilu_bld.f90 +++ b/psb_dilu_bld.f90 @@ -339,11 +339,17 @@ subroutine psb_dilu_bld(a,desc_a,p,upd,info) if (psb_sp_getifld(psb_upd_,p%av(u_pr_),info) /= psb_upd_perm_) then call psb_sp_trimsize(p%av(u_pr_),i1,i2,ia,info) if (info == 0) call psb_sp_reall(p%av(u_pr_),i1,i2,ia,info) + if (info /=0) then + write(0,*) 'Error from trimsize 1',info + endif endif if (psb_sp_getifld(psb_upd_,p%av(l_pr_),info) /= psb_upd_perm_) then call psb_sp_trimsize(p%av(l_pr_),i1,i2,ia,info) if (info == 0) call psb_sp_reall(p%av(l_pr_),i1,i2,ia,info) + if (info /=0) then + write(0,*) 'Error from trimsize 2',info + endif endif @@ -354,7 +360,7 @@ subroutine psb_dilu_bld(a,desc_a,p,upd,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dilu_fct.f90 b/psb_dilu_fct.f90 index 4b965702..950a9b1c 100644 --- a/psb_dilu_fct.f90 +++ b/psb_dilu_fct.f90 @@ -122,7 +122,7 @@ subroutine psb_dilu_fct(a,l,u,d,info,blck) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -466,7 +466,7 @@ contains 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dmlprc_aply.f90 b/psb_dmlprc_aply.f90 index b228130d..257eba21 100644 --- a/psb_dmlprc_aply.f90 +++ b/psb_dmlprc_aply.f90 @@ -753,7 +753,7 @@ subroutine psb_dmlprc_aply(alpha,baseprecv,x,beta,y,desc_data,trans,work,info) 9999 continue call psb_errpush(info,name) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dmlprc_bld.f90 b/psb_dmlprc_bld.f90 index 56ca6bd4..d26582a8 100644 --- a/psb_dmlprc_bld.f90 +++ b/psb_dmlprc_bld.f90 @@ -189,7 +189,7 @@ subroutine psb_dmlprc_bld(a,desc_a,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dprc_aply.f90 b/psb_dprc_aply.f90 index 0ec31ae6..c64e4fe7 100644 --- a/psb_dprc_aply.f90 +++ b/psb_dprc_aply.f90 @@ -134,7 +134,7 @@ subroutine psb_dprc_aply(prec,x,y,desc_data,info,trans, work) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -242,7 +242,7 @@ subroutine psb_dprc_aply1(prec,x,desc_data,info,trans) 9999 continue call psb_errpush(info,name) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dprecbld.f90 b/psb_dprecbld.f90 index 53debd43..7032b6c9 100644 --- a/psb_dprecbld.f90 +++ b/psb_dprecbld.f90 @@ -144,7 +144,7 @@ subroutine psb_dprecbld(a,desc_a,p,info,upd) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dprecfree.f90 b/psb_dprecfree.f90 index e1dd3264..c1497caf 100644 --- a/psb_dprecfree.f90 +++ b/psb_dprecfree.f90 @@ -63,7 +63,7 @@ subroutine psb_dprecfree(p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dslu_bld.f90 b/psb_dslu_bld.f90 index b1d824fd..daea1dfa 100644 --- a/psb_dslu_bld.f90 +++ b/psb_dslu_bld.f90 @@ -196,7 +196,7 @@ subroutine psb_dslu_bld(a,desc_a,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dsp_renum.f90 b/psb_dsp_renum.f90 index 347f2582..532d0a5c 100644 --- a/psb_dsp_renum.f90 +++ b/psb_dsp_renum.f90 @@ -353,7 +353,7 @@ subroutine psb_dsp_renum(a,desc_a,blck,p,atmp,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -448,7 +448,7 @@ contains 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_dumf_bld.f90 b/psb_dumf_bld.f90 index cb3c9008..1569f979 100644 --- a/psb_dumf_bld.f90 +++ b/psb_dumf_bld.f90 @@ -200,7 +200,7 @@ subroutine psb_dumf_bld(a,desc_a,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zasmatbld.f90 b/psb_zasmatbld.f90 index dbcf1e98..ff176981 100644 --- a/psb_zasmatbld.f90 +++ b/psb_zasmatbld.f90 @@ -176,16 +176,16 @@ Subroutine psb_zasmatbld(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt) ! ! Build the auiliary descriptor',desc_p%matrix_data(psb_n_row_) ! - call psb_cdbldovr(a,desc_data,novr,desc_p,info) + call psb_cdbldext(a,desc_data,novr,desc_p,info,extype=psb_ovt_asov_) if(info /= 0) then info=4010 - ch_err='psb_cdbldovr' + ch_err='psb_cdbldext' call psb_errpush(info,name,a_err=ch_err) goto 9999 end if Endif - if(debug) write(0,*) me,' From cdovr _:',desc_p%matrix_data(psb_n_row_),desc_p%matrix_data(psb_n_col_) + if(debug) write(0,*) me,' From bldext:',desc_p%matrix_data(psb_n_row_),desc_p%matrix_data(psb_n_col_) n_row = desc_p%matrix_data(psb_n_row_) @@ -226,7 +226,7 @@ Subroutine psb_zasmatbld(ptype,novr,a,blk,desc_data,upd,desc_p,info,outfmt) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zbaseprc_aply.f90 b/psb_zbaseprc_aply.f90 index 951aa47e..8490b6af 100644 --- a/psb_zbaseprc_aply.f90 +++ b/psb_zbaseprc_aply.f90 @@ -270,7 +270,7 @@ subroutine psb_zbaseprc_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) 9999 continue call psb_errpush(info,name,i_err=int_err,a_err=ch_err) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zbaseprc_bld.f90 b/psb_zbaseprc_bld.f90 index 2c493f61..c7ad01f5 100644 --- a/psb_zbaseprc_bld.f90 +++ b/psb_zbaseprc_bld.f90 @@ -252,7 +252,7 @@ subroutine psb_zbaseprc_bld(a,desc_a,p,info,upd) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zbjac_aply.f90 b/psb_zbjac_aply.f90 index 9ef190e7..aa228e0d 100644 --- a/psb_zbjac_aply.f90 +++ b/psb_zbjac_aply.f90 @@ -260,7 +260,7 @@ subroutine psb_zbjac_aply(alpha,prec,x,beta,y,desc_data,trans,work,info) 9999 continue call psb_errpush(info,name,i_err=int_err,a_err=ch_err) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zbldaggrmat.f90 b/psb_zbldaggrmat.f90 index 65c461fd..3d9d2c5d 100644 --- a/psb_zbldaggrmat.f90 +++ b/psb_zbldaggrmat.f90 @@ -88,7 +88,7 @@ subroutine psb_zbldaggrmat(a,desc_a,ac,desc_ac,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -323,7 +323,7 @@ contains 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -1027,7 +1027,7 @@ contains 9999 continue call psb_errpush(info,name) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zdiagsc_bld.f90 b/psb_zdiagsc_bld.f90 index b925954f..077d3179 100644 --- a/psb_zdiagsc_bld.f90 +++ b/psb_zdiagsc_bld.f90 @@ -154,7 +154,7 @@ subroutine psb_zdiagsc_bld(a,desc_a,p,upd,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zgenaggrmap.f90 b/psb_zgenaggrmap.f90 index f89b00ca..c203616a 100644 --- a/psb_zgenaggrmap.f90 +++ b/psb_zgenaggrmap.f90 @@ -283,7 +283,7 @@ subroutine psb_zgenaggrmap(aggr_type,a,desc_a,nlaggr,ilaggr,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zilu_bld.f90 b/psb_zilu_bld.f90 index 47085381..f7a2673f 100644 --- a/psb_zilu_bld.f90 +++ b/psb_zilu_bld.f90 @@ -352,7 +352,7 @@ subroutine psb_zilu_bld(a,desc_a,p,upd,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zilu_fct.f90 b/psb_zilu_fct.f90 index 22c96ab8..0c00375f 100644 --- a/psb_zilu_fct.f90 +++ b/psb_zilu_fct.f90 @@ -119,7 +119,7 @@ subroutine psb_zilu_fct(a,l,u,d,info,blck) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -463,7 +463,7 @@ contains 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zmlprc_aply.f90 b/psb_zmlprc_aply.f90 index d55bf707..01d0bcdc 100644 --- a/psb_zmlprc_aply.f90 +++ b/psb_zmlprc_aply.f90 @@ -750,7 +750,7 @@ subroutine psb_zmlprc_aply(alpha,baseprecv,x,beta,y,desc_data,trans,work,info) 9999 continue call psb_errpush(info,name) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zmlprc_bld.f90 b/psb_zmlprc_bld.f90 index 6fb737c7..b53bc401 100644 --- a/psb_zmlprc_bld.f90 +++ b/psb_zmlprc_bld.f90 @@ -189,7 +189,7 @@ subroutine psb_zmlprc_bld(a,desc_a,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zprc_aply.f90 b/psb_zprc_aply.f90 index ca48ce96..46b825ff 100644 --- a/psb_zprc_aply.f90 +++ b/psb_zprc_aply.f90 @@ -134,7 +134,7 @@ subroutine psb_zprc_aply(prec,x,y,desc_data,info,trans, work) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -239,7 +239,7 @@ subroutine psb_zprc_aply1(prec,x,desc_data,info,trans) 9999 continue call psb_errpush(info,name) call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zprecbld.f90 b/psb_zprecbld.f90 index 2a51df83..7417838e 100644 --- a/psb_zprecbld.f90 +++ b/psb_zprecbld.f90 @@ -144,7 +144,7 @@ subroutine psb_zprecbld(a,desc_a,p,info,upd) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zprecfree.f90 b/psb_zprecfree.f90 index 7db0b54d..7dbe944a 100644 --- a/psb_zprecfree.f90 +++ b/psb_zprecfree.f90 @@ -65,7 +65,7 @@ subroutine psb_zprecfree(p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zslu_bld.f90 b/psb_zslu_bld.f90 index 1b51d6d1..dd597402 100644 --- a/psb_zslu_bld.f90 +++ b/psb_zslu_bld.f90 @@ -194,7 +194,7 @@ subroutine psb_zslu_bld(a,desc_a,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zsp_renum.f90 b/psb_zsp_renum.f90 index 7cdf5f83..0f1c1638 100644 --- a/psb_zsp_renum.f90 +++ b/psb_zsp_renum.f90 @@ -352,7 +352,7 @@ subroutine psb_zsp_renum(a,desc_a,blck,p,atmp,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if @@ -447,7 +447,7 @@ contains 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if diff --git a/psb_zumf_bld.f90 b/psb_zumf_bld.f90 index c82ce39f..9379dbb9 100644 --- a/psb_zumf_bld.f90 +++ b/psb_zumf_bld.f90 @@ -199,7 +199,7 @@ subroutine psb_zumf_bld(a,desc_a,p,info) 9999 continue call psb_erractionrestore(err_act) - if (err_act.eq.act_abort) then + if (err_act.eq.psb_act_abort_) then call psb_error() return end if