odd bug catched; in italian:
Usavo dichiarare e inizializzare "flags" come:
integer :: flags=c_def_flags
e da una chiamata all'altra, rimaneva il flag della chiamata precedente.
Poiche' nel modulo fortran, "aggiungo" is flag usando + anziche' un
operatore OR, il flag (un bit) rimanendo dalla chiamata precedente,
(la costruzione di L: qui siamo nella costruzione di U) veniva "spostato"
di una posizione in su, risultando nell'assenza effettiva del bit
"triangolare", a fronte della sua aggiunta.
Ho rimediato con:
integer :: flags
flags=c_def_flags
test/newfmt/runs/ppde.inp
test/serial/Makefile
test/serial/psb_d_rsb_mat_mod.F03
test/serial/rsb_mod.f03
Take out rsb stuff from test/serial, since it's now in opt.
opt/Makefile
opt/psb_d_rsb_mat_mod.F03
opt/rsb_mod.f03
test/newfmt/Makefile
test/newfmt/ppde.f90
test/newfmt/runs/ppde.inp
test/serial/d_matgen.f03
test/serial/psb_d_rsb_mat_mod.F03
Minor fix for interface of rsb_init.
Copied RSB interface to OPT/
Linked into TEST/NEWFMT: start of debug.
Status: TO/FROM COO/FMT (CSR), SPMV, GET_DIAG seem to be working
(tested from ppde with DIAG preconditioner).
GETBLK segfaults.
base/internals/Makefile
base/internals/psi_desc_impl.f90
base/internals/psi_ovrl_restr.f90
base/internals/psi_ovrl_save.f90
base/internals/psi_ovrl_upd.f90
base/serial/Makefile
base/serial/psi_desc_impl.f90
base/serial/psi_ovrl_restr.f90
base/serial/psi_ovrl_save.f90
base/serial/psi_ovrl_upd.f90
Moved routines from serial to internals where they really belong.
prec/psb_c_base_prec_mod.f03
prec/psb_c_bjacprec.f03
prec/psb_c_diagprec.f03
prec/psb_c_nullprec.f03
prec/psb_cprc_aply.f90
prec/psb_cprecbld.f90
prec/psb_d_base_prec_mod.f03
prec/psb_d_bjacprec.f03
prec/psb_d_diagprec.f03
prec/psb_d_nullprec.f03
prec/psb_dprc_aply.f90
prec/psb_dprecbld.f90
prec/psb_prec_mod.f90
prec/psb_s_base_prec_mod.f03
prec/psb_s_bjacprec.f03
prec/psb_s_diagprec.f03
prec/psb_s_nullprec.f03
prec/psb_sprc_aply.f90
prec/psb_sprecbld.f90
prec/psb_z_base_prec_mod.f03
prec/psb_z_bjacprec.f03
prec/psb_z_diagprec.f03
prec/psb_z_nullprec.f03
prec/psb_zprc_aply.f90
prec/psb_zprecbld.f90
test/newfmt
test/newfmt/Makefile
test/newfmt/ppde.f90
test/newfmt/runs
test/newfmt/runs/Makefile
test/newfmt/runs/ppde.inp
test/newfmt/spde.f90
1. New test/newfmt directory in which to test for new storage
versions.
2. New MOLD/AFMT arguments to PRECBLD for storing the matrices.
In preconditioners such as DIAG and NULL they are ignored.
psb_c_csc_impl.f03
psb_c_csr_impl.f03
psb_d_csc_impl.f03
psb_d_csr_impl.f03
psb_s_csc_impl.f03
psb_s_csr_impl.f03
psb_z_csc_impl.f03
psb_z_csr_impl.f03
Fixed use module for missing initialization of debug_unit in some methods.
--with-rsb Specify Recursive Sparse BLAS library linkage info
(that is, the output of librsb-config --static
--ldflags, or a directory where the usual
bin/include/lib subdirs with a regular RSB
installation resides, or nothing to make the
configure script invoke librsb-config)
interfaced the following operations: psb_d_rsb_rowsum, psb_d_rsb_colsum, d_rsb_get_nz_row, psb_d_rsb_reinit, psb_d_rsb_mold.
the colsum and rowsum operations were implemented brand new today in rsb.
Salvatore Filippone
Michele Martone