diff --git a/config/pac.m4 b/config/pac.m4 index 5dc4e0c5..c9f55d94 100644 --- a/config/pac.m4 +++ b/config/pac.m4 @@ -1792,6 +1792,9 @@ AC_DEFUN(PAC_CHECK_METIS, Default: "-lmetis"]), [psblas_cv_metis=$withval], [psblas_cv_metis='-lmetis']) +AC_ARG_WITH(metisincfile, AC_HELP_STRING([--with-metisincfile=DIR], [Specify the name for METIS include file.]), + [psblas_cv_metisincfile=$withval], + [psblas_cv_metisincfile='metis.h']) AC_ARG_WITH(metisdir, AC_HELP_STRING([--with-metisdir=DIR], [Specify the directory for METIS library and includes.]), [psblas_cv_metisdir=$withval], [psblas_cv_metisdir='']) @@ -1821,7 +1824,7 @@ if test "x$psblas_cv_metislibdir" != "x"; then fi AC_MSG_NOTICE([metis dir $psblas_cv_metisdir]) -AC_CHECK_HEADERS([limits.h metis.h], +AC_CHECK_HEADERS([limits.h $psblas_cv_metisincfile], [pac_metis_header_ok=yes], [pac_metis_header_ok=no; METIS_INCLUDES=""]) if test "x$pac_metis_header_ok" == "xno" ; then @@ -1831,7 +1834,7 @@ dnl Maybe Include or include subdirs? CPPFLAGS="$METIS_INCLUDES $SAVE_CPPFLAGS" AC_MSG_CHECKING([for metis_h in $METIS_INCLUDES]) - AC_CHECK_HEADERS([limits.h metis.h], + AC_CHECK_HEADERS([limits.h $psblas_cv_metisincfile], [pac_metis_header_ok=yes], [pac_metis_header_ok=no; METIS_INCLUDES=""]) fi @@ -1840,11 +1843,27 @@ dnl Maybe new structure with METIS UFconfig METIS? unset ac_cv_header_metis_h METIS_INCLUDES="-I$psblas_cv_metisdir/UFconfig -I$psblas_cv_metisdir/METIS/Include -I$psblas_cv_metisdir/METIS/Include" CPPFLAGS="$METIS_INCLUDES $SAVE_CPPFLAGS" - AC_CHECK_HEADERS([limits.h metis.h], + AC_CHECK_HEADERS([limits.h $psblas_cv_metisincfile], [pac_metis_header_ok=yes], [pac_metis_header_ok=no; METIS_INCLUDES=""]) fi +if test "x$pac_metis_header_ok" == "xyes" ; then + AC_LANG_PUSH([C]) + AC_MSG_CHECKING([for METIS idx size]) + AC_LINK_IFELSE([AC_LANG_SOURCE( + #include + #include "$psblas_cv_metisincfile" + void main(){ + printf("%d\n",IDXTYPEWIDTH); + } + )], + [pac_cv_metis_idx=`./conftest${ac_exeext} | sed 's/^ *//'`], + [pac_cv_metis_idx="unknown"]) + AC_MSG_RESULT($pac_cv_metis_idx) + + AC_LANG_POP() +fi if test "x$pac_metis_header_ok" == "xyes" ; then psblas_cv_metis_includes="$METIS_INCLUDES" diff --git a/configure b/configure index 727f56c3..ccefc2c0 100755 --- a/configure +++ b/configure @@ -1,6 +1,6 @@ #! /bin/sh # Guess values for system-dependent variables and create Makefiles. -# Generated by GNU Autoconf 2.69 for PSBLAS 3.5. +# Generated by GNU Autoconf 2.69 for PSBLAS 3.7.0. # # Report bugs to . # @@ -580,8 +580,8 @@ MAKEFLAGS= # Identity of this package. PACKAGE_NAME='PSBLAS' PACKAGE_TARNAME='psblas' -PACKAGE_VERSION='3.5' -PACKAGE_STRING='PSBLAS 3.5' +PACKAGE_VERSION='3.7.0' +PACKAGE_STRING='PSBLAS 3.7.0' PACKAGE_BUGREPORT='https://github.com/sfilippone/psblas3/issues' PACKAGE_URL='' @@ -626,6 +626,7 @@ ac_subst_vars='am__EXEEXT_FALSE am__EXEEXT_TRUE LTLIBOBJS LIBOBJS +METISINCFILE UTILLIBNAME METHDLIBNAME PRECLIBNAME @@ -676,7 +677,6 @@ am__nodep AMDEPBACKSLASH AMDEP_FALSE AMDEP_TRUE -am__quote am__include DEPDIR am__untar @@ -753,7 +753,8 @@ PACKAGE_VERSION PACKAGE_TARNAME PACKAGE_NAME PATH_SEPARATOR -SHELL' +SHELL +am__quote' ac_subst_files='' ac_user_opts=' enable_option_checking @@ -775,6 +776,7 @@ with_blasdir with_lapack with_rsb with_metis +with_metisincfile with_metisdir with_metisincdir with_metislibdir @@ -1336,7 +1338,7 @@ if test "$ac_init_help" = "long"; then # Omit some internal or obsolete options to make the list less imposing. # This message is too long to be a string in the A/UX 3.1 sh. cat <<_ACEOF -\`configure' configures PSBLAS 3.5 to adapt to many kinds of systems. +\`configure' configures PSBLAS 3.7.0 to adapt to many kinds of systems. Usage: $0 [OPTION]... [VAR=VALUE]... @@ -1402,7 +1404,7 @@ fi if test -n "$ac_init_help"; then case $ac_init_help in - short | recursive ) echo "Configuration of PSBLAS 3.5:";; + short | recursive ) echo "Configuration of PSBLAS 3.7.0:";; esac cat <<\_ACEOF @@ -1454,6 +1456,7 @@ Optional Packages: configure script invoke librsb-config) --with-metis=LIBNAME Specify the library name for METIS library. Default: "-lmetis" + --with-metisincfile=DIR Specify the name for METIS include file. --with-metisdir=DIR Specify the directory for METIS library and includes. --with-metisincdir=DIR Specify the directory for METIS includes. @@ -1544,7 +1547,7 @@ fi test -n "$ac_init_help" && exit $ac_status if $ac_init_version; then cat <<\_ACEOF -PSBLAS configure 3.5 +PSBLAS configure 3.7.0 generated by GNU Autoconf 2.69 Copyright (C) 2012 Free Software Foundation, Inc. @@ -2180,7 +2183,7 @@ cat >config.log <<_ACEOF This file contains any messages produced by compilers while running configure, to aid debugging if configure makes a mistake. -It was created by PSBLAS $as_me 3.5, which was +It was created by PSBLAS $as_me 3.7.0, which was generated by GNU Autoconf 2.69. Invocation command line was $ $0 $@ @@ -2531,7 +2534,7 @@ ac_compiler_gnu=$ac_cv_c_compiler_gnu # VERSION is the file containing the PSBLAS version code # FIXME -psblas_cv_version="3.5" +psblas_cv_version="3.7.0" # A sample source file @@ -3515,6 +3518,68 @@ ac_compile='$CC -c $CFLAGS $CPPFLAGS conftest.$ac_ext >&5' ac_link='$CC -o conftest$ac_exeext $CFLAGS $CPPFLAGS $LDFLAGS conftest.$ac_ext $LIBS >&5' ac_compiler_gnu=$ac_cv_c_compiler_gnu +# Expand $ac_aux_dir to an absolute path. +am_aux_dir=`cd "$ac_aux_dir" && pwd` + +ac_ext=c +ac_cpp='$CPP $CPPFLAGS' +ac_compile='$CC -c $CFLAGS $CPPFLAGS conftest.$ac_ext >&5' +ac_link='$CC -o conftest$ac_exeext $CFLAGS $CPPFLAGS $LDFLAGS conftest.$ac_ext $LIBS >&5' +ac_compiler_gnu=$ac_cv_c_compiler_gnu +{ $as_echo "$as_me:${as_lineno-$LINENO}: checking whether $CC understands -c and -o together" >&5 +$as_echo_n "checking whether $CC understands -c and -o together... " >&6; } +if ${am_cv_prog_cc_c_o+:} false; then : + $as_echo_n "(cached) " >&6 +else + cat confdefs.h - <<_ACEOF >conftest.$ac_ext +/* end confdefs.h. */ + +int +main () +{ + + ; + return 0; +} +_ACEOF + # Make sure it works both with $CC and with simple cc. + # Following AC_PROG_CC_C_O, we do the test twice because some + # compilers refuse to overwrite an existing .o file with -o, + # though they will create one. + am_cv_prog_cc_c_o=yes + for am_i in 1 2; do + if { echo "$as_me:$LINENO: $CC -c conftest.$ac_ext -o conftest2.$ac_objext" >&5 + ($CC -c conftest.$ac_ext -o conftest2.$ac_objext) >&5 2>&5 + ac_status=$? + echo "$as_me:$LINENO: \$? = $ac_status" >&5 + (exit $ac_status); } \ + && test -f conftest2.$ac_objext; then + : OK + else + am_cv_prog_cc_c_o=no + break + fi + done + rm -f core conftest* + unset am_i +fi +{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $am_cv_prog_cc_c_o" >&5 +$as_echo "$am_cv_prog_cc_c_o" >&6; } +if test "$am_cv_prog_cc_c_o" != yes; then + # Losing compiler, so override with the script. + # FIXME: It is wrong to rewrite CC. + # But if we don't then we get into trouble of one sort or another. + # A longer-term fix would be to have automake use am__CC in this case, + # and then we could set am__CC="\$(top_srcdir)/compile \$(CC)" + CC="$am_aux_dir/compile $CC" +fi +ac_ext=c +ac_cpp='$CPP $CPPFLAGS' +ac_compile='$CC -c $CFLAGS $CPPFLAGS conftest.$ac_ext >&5' +ac_link='$CC -o conftest$ac_exeext $CFLAGS $CPPFLAGS $LDFLAGS conftest.$ac_ext $LIBS >&5' +ac_compiler_gnu=$ac_cv_c_compiler_gnu + + CFLAGS="$save_CFLAGS"; @@ -4942,7 +5007,7 @@ else fi -am__api_version='1.13' +am__api_version='1.16' { $as_echo "$as_me:${as_lineno-$LINENO}: checking whether build environment is sane" >&5 $as_echo_n "checking whether build environment is sane... " >&6; } @@ -5021,9 +5086,6 @@ test "$program_suffix" != NONE && ac_script='s/[\\$]/&&/g;s/;s,x,x,$//' program_transform_name=`$as_echo "$program_transform_name" | sed "$ac_script"` -# expand $ac_aux_dir to an absolute path -am_aux_dir=`cd $ac_aux_dir && pwd` - if test x"${MISSING+set}" != xset; then case $am_aux_dir in *\ * | *\ *) @@ -5041,7 +5103,7 @@ else $as_echo "$as_me: WARNING: 'missing' script is too old or missing" >&2;} fi -if test x"${install_sh}" != xset; then +if test x"${install_sh+set}" != xset; then case $am_aux_dir in *\ * | *\ *) install_sh="\${SHELL} '$am_aux_dir/install-sh'" ;; @@ -5278,45 +5340,45 @@ DEPDIR="${am__leading_dot}deps" ac_config_commands="$ac_config_commands depfiles" - -am_make=${MAKE-make} -cat > confinc << 'END' +{ $as_echo "$as_me:${as_lineno-$LINENO}: checking whether ${MAKE-make} supports the include directive" >&5 +$as_echo_n "checking whether ${MAKE-make} supports the include directive... " >&6; } +cat > confinc.mk << 'END' am__doit: - @echo this is the am__doit target + @echo this is the am__doit target >confinc.out .PHONY: am__doit END -# If we don't find an include directive, just comment out the code. -{ $as_echo "$as_me:${as_lineno-$LINENO}: checking for style of include used by $am_make" >&5 -$as_echo_n "checking for style of include used by $am_make... " >&6; } am__include="#" am__quote= -_am_result=none -# First try GNU make style include. -echo "include confinc" > confmf -# Ignore all kinds of additional output from 'make'. -case `$am_make -s -f confmf 2> /dev/null` in #( -*the\ am__doit\ target*) - am__include=include - am__quote= - _am_result=GNU - ;; -esac -# Now try BSD make style include. -if test "$am__include" = "#"; then - echo '.include "confinc"' > confmf - case `$am_make -s -f confmf 2> /dev/null` in #( - *the\ am__doit\ target*) - am__include=.include - am__quote="\"" - _am_result=BSD +# BSD make does it like this. +echo '.include "confinc.mk" # ignored' > confmf.BSD +# Other make implementations (GNU, Solaris 10, AIX) do it like this. +echo 'include confinc.mk # ignored' > confmf.GNU +_am_result=no +for s in GNU BSD; do + { echo "$as_me:$LINENO: ${MAKE-make} -f confmf.$s && cat confinc.out" >&5 + (${MAKE-make} -f confmf.$s && cat confinc.out) >&5 2>&5 + ac_status=$? + echo "$as_me:$LINENO: \$? = $ac_status" >&5 + (exit $ac_status); } + case $?:`cat confinc.out 2>/dev/null` in #( + '0:this is the am__doit target') : + case $s in #( + BSD) : + am__include='.include' am__quote='"' ;; #( + *) : + am__include='include' am__quote='' ;; +esac ;; #( + *) : ;; - esac -fi - - -{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $_am_result" >&5 -$as_echo "$_am_result" >&6; } -rm -f confinc confmf +esac + if test "$am__include" != "#"; then + _am_result="yes ($s style)" + break + fi +done +rm -f confinc.* confmf.* +{ $as_echo "$as_me:${as_lineno-$LINENO}: result: ${_am_result}" >&5 +$as_echo "${_am_result}" >&6; } # Check whether --enable-dependency-tracking was given. if test "${enable_dependency_tracking+set}" = set; then : @@ -5398,7 +5460,7 @@ fi # Define the identity of the package. PACKAGE='psblas' - VERSION='3.5' + VERSION='3.7.0' cat >>confdefs.h <<_ACEOF @@ -5428,12 +5490,12 @@ MAKEINFO=${MAKEINFO-"${am_missing_run}makeinfo"} # For better backward compatibility. To be removed once Automake 1.9.x # dies out for good. For more background, see: -# -# +# +# mkdir_p='$(MKDIR_P)' -# We need awk for the "check" target. The system "awk" is bad on -# some platforms. +# We need awk for the "check" target (and possibly the TAP driver). The +# system "awk" is bad on some platforms. # Always define AMTAR for backward compatibility. Yes, it's still used # in the wild :-( We should find a proper way to deprecate it ... AMTAR='$${TAR-tar}' @@ -5577,6 +5639,48 @@ fi +# POSIX will say in a future version that running "rm -f" with no argument +# is OK; and we want to be able to make that assumption in our Makefile +# recipes. So use an aggressive probe to check that the usage we want is +# actually supported "in the wild" to an acceptable degree. +# See automake bug#10828. +# To make any issue more visible, cause the running configure to be aborted +# by default if the 'rm' program in use doesn't match our expectations; the +# user can still override this though. +if rm -f && rm -fr && rm -rf; then : OK; else + cat >&2 <<'END' +Oops! + +Your 'rm' program seems unable to run without file operands specified +on the command line, even when the '-f' option is present. This is contrary +to the behaviour of most rm programs out there, and not conforming with +the upcoming POSIX standard: + +Please tell bug-automake@gnu.org about your system, including the value +of your $PATH and any error possibly output before this message. This +can help us improve future automake versions. + +END + if test x"$ACCEPT_INFERIOR_RM_PROGRAM" = x"yes"; then + echo 'Configuration will proceed anyway, since you have set the' >&2 + echo 'ACCEPT_INFERIOR_RM_PROGRAM variable to "yes"' >&2 + echo >&2 + else + cat >&2 <<'END' +Aborting the configuration process, to ensure you take notice of the issue. + +You can download and install GNU coreutils to get an 'rm' implementation +that behaves properly: . + +If you want to complete the configuration process using your problematic +'rm' anyway, export the environment variable ACCEPT_INFERIOR_RM_PROGRAM +to "yes", and re-run configure. + +END + as_fn_error $? "Your 'rm' program is bad, sorry." "$LINENO" 5 + fi +fi + @@ -7260,10 +7364,6 @@ ac_compiler_gnu=$ac_cv_c_compiler_gnu -# -# Optional features -# - { $as_echo "$as_me:${as_lineno-$LINENO}: checking support for Fortran VOLATILE" >&5 $as_echo_n "checking support for Fortran VOLATILE... " >&6; } ac_ext=${ac_fc_srcext-f} @@ -7283,12 +7383,13 @@ _ACEOF if ac_fn_fc_try_compile "$LINENO"; then : { $as_echo "$as_me:${as_lineno-$LINENO}: result: yes" >&5 $as_echo "yes" >&6; } - FDEFINES="$psblas_cv_define_prepend-DHAVE_VOLATILE $FDEFINES" + : else { $as_echo "$as_me:${as_lineno-$LINENO}: result: no" >&5 $as_echo "no" >&6; } echo "configure: failed program was:" >&5 cat conftest.$ac_ext >&5 + as_fn_error $? "Sorry, cannot build PSBLAS without support for VOLATILE" "$LINENO" 5 fi rm -f core conftest.err conftest.$ac_objext conftest.$ac_ext @@ -7300,31 +7401,21 @@ ac_compiler_gnu=$ac_cv_c_compiler_gnu -{ $as_echo "$as_me:${as_lineno-$LINENO}: checking test GENERIC interfaces" >&5 -$as_echo_n "checking test GENERIC interfaces... " >&6; } -ac_ext=${ac_fc_srcext-f} +{ $as_echo "$as_me:${as_lineno-$LINENO}: checking support for ISO_FORTRAN_ENV" >&5 +$as_echo_n "checking support for ISO_FORTRAN_ENV... " >&6; } + ac_ext=${ac_fc_srcext-f} ac_compile='$FC -c $FCFLAGS $ac_fcflags_srcext conftest.$ac_ext >&5' ac_link='$FC -o conftest$ac_exeext $FCFLAGS $LDFLAGS $ac_fcflags_srcext conftest.$ac_ext $LIBS >&5' ac_compiler_gnu=$ac_cv_fc_compiler_gnu ac_exeext='' - ac_ext='F90' + ac_ext='f90' ac_fc=${MPIFC-$FC}; cat > conftest.$ac_ext <<_ACEOF -module conftest - - interface foo - subroutine i_sub_foo(v) - integer, intent(inout) :: v(:) - end subroutine i_sub_foo - end interface foo - - interface bar - procedure i_sub_foo - end interface bar - -end module conftest + program test + use iso_fortran_env + end program test _ACEOF if ac_fn_fc_try_compile "$LINENO"; then : { $as_echo "$as_me:${as_lineno-$LINENO}: result: yes" >&5 @@ -7335,7 +7426,7 @@ else $as_echo "no" >&6; } echo "configure: failed program was:" >&5 cat conftest.$ac_ext >&5 - FDEFINES="$psblas_cv_define_prepend-DHAVE_BUGGY_GENERICS $FDEFINES" + as_fn_error $? "Sorry, cannot build PSBLAS without support for ISO_FORTRAN_ENV" "$LINENO" 5 fi rm -f core conftest.err conftest.$ac_objext conftest.$ac_ext @@ -7347,8 +7438,8 @@ ac_compiler_gnu=$ac_cv_c_compiler_gnu -{ $as_echo "$as_me:${as_lineno-$LINENO}: checking support for Fortran FLUSH statement" >&5 -$as_echo_n "checking support for Fortran FLUSH statement... " >&6; } +{ $as_echo "$as_me:${as_lineno-$LINENO}: checking support for Fortran FINAL" >&5 +$as_echo_n "checking support for Fortran FINAL... " >&6; } ac_ext=${ac_fc_srcext-f} ac_compile='$FC -c $FCFLAGS $ac_fcflags_srcext conftest.$ac_ext >&5' ac_link='$FC -o conftest$ac_exeext $FCFLAGS $LDFLAGS $ac_fcflags_srcext conftest.$ac_ext $LIBS >&5' @@ -7359,23 +7450,35 @@ ac_compiler_gnu=$ac_cv_fc_compiler_gnu ac_fc=${MPIFC-$FC}; cat > conftest.$ac_ext <<_ACEOF +module conftest_mod + type foo + integer :: i + contains + final :: destroy_foo + end type foo + + private destroy_foo +contains + subroutine destroy_foo(a) + type(foo) :: a + ! Just a test + end subroutine destroy_foo +end module conftest_mod program conftest - integer :: iunit=10 - open(10) - write(10,*) 'Test ' - flush(10) - close(10) + use conftest_mod + type(foo) :: foovar end program conftest _ACEOF if ac_fn_fc_try_compile "$LINENO"; then : { $as_echo "$as_me:${as_lineno-$LINENO}: result: yes" >&5 $as_echo "yes" >&6; } - FDEFINES="$psblas_cv_define_prepend-DHAVE_FLUSH_STMT $FDEFINES" + : else { $as_echo "$as_me:${as_lineno-$LINENO}: result: no" >&5 $as_echo "no" >&6; } echo "configure: failed program was:" >&5 cat conftest.$ac_ext >&5 + as_fn_error $? "Sorry, cannot build PSBLAS without support for FINAL" "$LINENO" 5 fi rm -f core conftest.err conftest.$ac_objext conftest.$ac_ext @@ -7387,31 +7490,48 @@ ac_compiler_gnu=$ac_cv_c_compiler_gnu -{ $as_echo "$as_me:${as_lineno-$LINENO}: checking support for ISO_FORTRAN_ENV" >&5 -$as_echo_n "checking support for ISO_FORTRAN_ENV... " >&6; } - ac_ext=${ac_fc_srcext-f} + +# +# Optional features +# + + +{ $as_echo "$as_me:${as_lineno-$LINENO}: checking test GENERIC interfaces" >&5 +$as_echo_n "checking test GENERIC interfaces... " >&6; } +ac_ext=${ac_fc_srcext-f} ac_compile='$FC -c $FCFLAGS $ac_fcflags_srcext conftest.$ac_ext >&5' ac_link='$FC -o conftest$ac_exeext $FCFLAGS $LDFLAGS $ac_fcflags_srcext conftest.$ac_ext $LIBS >&5' ac_compiler_gnu=$ac_cv_fc_compiler_gnu ac_exeext='' - ac_ext='f90' + ac_ext='F90' ac_fc=${MPIFC-$FC}; cat > conftest.$ac_ext <<_ACEOF - program test - use iso_fortran_env - end program test +module conftest + + interface foo + subroutine i_sub_foo(v) + integer, intent(inout) :: v(:) + end subroutine i_sub_foo + end interface foo + + interface bar + procedure i_sub_foo + end interface bar + +end module conftest _ACEOF if ac_fn_fc_try_compile "$LINENO"; then : { $as_echo "$as_me:${as_lineno-$LINENO}: result: yes" >&5 $as_echo "yes" >&6; } - FDEFINES="$psblas_cv_define_prepend-DHAVE_ISO_FORTRAN_ENV $FDEFINES" + : else { $as_echo "$as_me:${as_lineno-$LINENO}: result: no" >&5 $as_echo "no" >&6; } echo "configure: failed program was:" >&5 cat conftest.$ac_ext >&5 + FDEFINES="$psblas_cv_define_prepend-DHAVE_BUGGY_GENERICS $FDEFINES" fi rm -f core conftest.err conftest.$ac_objext conftest.$ac_ext @@ -7423,8 +7543,8 @@ ac_compiler_gnu=$ac_cv_c_compiler_gnu -{ $as_echo "$as_me:${as_lineno-$LINENO}: checking support for Fortran FINAL" >&5 -$as_echo_n "checking support for Fortran FINAL... " >&6; } +{ $as_echo "$as_me:${as_lineno-$LINENO}: checking support for Fortran FLUSH statement" >&5 +$as_echo_n "checking support for Fortran FLUSH statement... " >&6; } ac_ext=${ac_fc_srcext-f} ac_compile='$FC -c $FCFLAGS $ac_fcflags_srcext conftest.$ac_ext >&5' ac_link='$FC -o conftest$ac_exeext $FCFLAGS $LDFLAGS $ac_fcflags_srcext conftest.$ac_ext $LIBS >&5' @@ -7435,29 +7555,18 @@ ac_compiler_gnu=$ac_cv_fc_compiler_gnu ac_fc=${MPIFC-$FC}; cat > conftest.$ac_ext <<_ACEOF -module conftest_mod - type foo - integer :: i - contains - final :: destroy_foo - end type foo - - private destroy_foo -contains - subroutine destroy_foo(a) - type(foo) :: a - ! Just a test - end subroutine destroy_foo -end module conftest_mod program conftest - use conftest_mod - type(foo) :: foovar + integer :: iunit=10 + open(10) + write(10,*) 'Test ' + flush(10) + close(10) end program conftest _ACEOF if ac_fn_fc_try_compile "$LINENO"; then : { $as_echo "$as_me:${as_lineno-$LINENO}: result: yes" >&5 $as_echo "yes" >&6; } - FDEFINES="$psblas_cv_define_prepend-DHAVE_FINAL $FDEFINES" + FDEFINES="$psblas_cv_define_prepend-DHAVE_FLUSH_STMT $FDEFINES" else { $as_echo "$as_me:${as_lineno-$LINENO}: result: no" >&5 $as_echo "no" >&6; } @@ -8756,6 +8865,14 @@ else fi +# Check whether --with-metisincfile was given. +if test "${with_metisincfile+set}" = set; then : + withval=$with_metisincfile; psblas_cv_metisincfile=$withval +else + psblas_cv_metisincfile='metis.h' +fi + + # Check whether --with-metisdir was given. if test "${with_metisdir+set}" = set; then : withval=$with_metisdir; psblas_cv_metisdir=$withval @@ -8805,7 +8922,7 @@ fi { $as_echo "$as_me:${as_lineno-$LINENO}: metis dir $psblas_cv_metisdir" >&5 $as_echo "$as_me: metis dir $psblas_cv_metisdir" >&6;} -for ac_header in limits.h metis.h +for ac_header in limits.h $psblas_cv_metisincfile do : as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh` ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default" @@ -8827,7 +8944,7 @@ if test "x$pac_metis_header_ok" == "xno" ; then { $as_echo "$as_me:${as_lineno-$LINENO}: checking for metis_h in $METIS_INCLUDES" >&5 $as_echo_n "checking for metis_h in $METIS_INCLUDES... " >&6; } - for ac_header in limits.h metis.h + for ac_header in limits.h $psblas_cv_metisincfile do : as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh` ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default" @@ -8847,7 +8964,7 @@ if test "x$pac_metis_header_ok" == "xno" ; then unset ac_cv_header_metis_h METIS_INCLUDES="-I$psblas_cv_metisdir/UFconfig -I$psblas_cv_metisdir/METIS/Include -I$psblas_cv_metisdir/METIS/Include" CPPFLAGS="$METIS_INCLUDES $SAVE_CPPFLAGS" - for ac_header in limits.h metis.h + for ac_header in limits.h $psblas_cv_metisincfile do : as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh` ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default" @@ -8864,6 +8981,41 @@ done fi +if test "x$pac_metis_header_ok" == "xyes" ; then + ac_ext=c +ac_cpp='$CPP $CPPFLAGS' +ac_compile='$CC -c $CFLAGS $CPPFLAGS conftest.$ac_ext >&5' +ac_link='$CC -o conftest$ac_exeext $CFLAGS $CPPFLAGS $LDFLAGS conftest.$ac_ext $LIBS >&5' +ac_compiler_gnu=$ac_cv_c_compiler_gnu + + { $as_echo "$as_me:${as_lineno-$LINENO}: checking for METIS idx size" >&5 +$as_echo_n "checking for METIS idx size... " >&6; } + cat confdefs.h - <<_ACEOF >conftest.$ac_ext +/* end confdefs.h. */ +#include + #include "$psblas_cv_metisincfile" + void main(){ + printf("%d\n",IDXTYPEWIDTH); + } + +_ACEOF +if ac_fn_c_try_link "$LINENO"; then : + pac_cv_metis_idx=`./conftest${ac_exeext} | sed 's/^ *//'` +else + pac_cv_metis_idx="unknown" +fi +rm -f core conftest.err conftest.$ac_objext \ + conftest$ac_exeext conftest.$ac_ext + { $as_echo "$as_me:${as_lineno-$LINENO}: result: $pac_cv_metis_idx" >&5 +$as_echo "$pac_cv_metis_idx" >&6; } + + ac_ext=c +ac_cpp='$CPP $CPPFLAGS' +ac_compile='$CC -c $CFLAGS $CPPFLAGS conftest.$ac_ext >&5' +ac_link='$CC -o conftest$ac_exeext $CFLAGS $CPPFLAGS $LDFLAGS conftest.$ac_ext $LIBS >&5' +ac_compiler_gnu=$ac_cv_c_compiler_gnu + +fi if test "x$pac_metis_header_ok" == "xyes" ; then psblas_cv_metis_includes="$METIS_INCLUDES" @@ -9003,9 +9155,12 @@ fi LIBS="$SAVE_LIBS"; CPPFLAGS="$SAVE_CPPFLAGS"; -if test "x$psblas_cv_have_metis" == "xyes" ; then - FDEFINES="$psblas_cv_define_prepend-DHAVE_METIS $FDEFINES" - CDEFINES="-DHAVE_METIS_ $psblas_cv_metis_includes $CDEFINES" +if test "x$pac_cv_metis_idx" == "xunknown" ; then + true ; # do nothing +elif test "x$psblas_cv_have_metis" == "xyes" ; then + FDEFINES="$psblas_cv_define_prepend-DHAVE_METIS $psblas_cv_define_prepend-DMETIS_$pac_cv_metis_idx $FDEFINES" + CDEFINES="-DHAVE_METIS_ $psblas_cv_metis_includes $CDEFINES -DMETIS_$pac_cv_metis_idx" + METISINCFILE=$psblas_cv_metisincfile fi @@ -9218,11 +9373,6 @@ fi -# Note : We would like to detect PSBLAS, but this is complicated by the -# module symbols mangling rules, which are compiler specific ! -# -# Moreover, the PSBLAS doesn't have an installer, currently. - ############################################################################### # Library target directory and archive files. ############################################################################### @@ -9275,7 +9425,6 @@ UTILLIBNAME=libpsb_util.a - PSBLASRULES=' @@ -9301,11 +9450,14 @@ FDEFINES=$(PSBFDEFINES) + ############################################################################### # the following files will be created by Automake ac_config_files="$ac_config_files Make.inc" +ac_config_files="$ac_config_files util/psb_metis_int.h" + cat >confcache <<\_ACEOF # This file is a shell script that caches the results of configure # tests run on this system so they can be shared between configure @@ -9872,7 +10024,7 @@ cat >>$CONFIG_STATUS <<\_ACEOF || ac_write_fail=1 # report actual input values of CONFIG_FILES etc. instead of their # values after options handling. ac_log=" -This file was extended by PSBLAS $as_me 3.5, which was +This file was extended by PSBLAS $as_me 3.7.0, which was generated by GNU Autoconf 2.69. Invocation command line was CONFIG_FILES = $CONFIG_FILES @@ -9929,7 +10081,7 @@ _ACEOF cat >>$CONFIG_STATUS <<_ACEOF || ac_write_fail=1 ac_cs_config="`$as_echo "$ac_configure_args" | sed 's/^ //; s/[\\""\`\$]/\\\\&/g'`" ac_cs_version="\\ -PSBLAS config.status 3.5 +PSBLAS config.status 3.7.0 configured by $0, generated by GNU Autoconf 2.69, with options \\"\$ac_cs_config\\" @@ -10037,7 +10189,7 @@ cat >>$CONFIG_STATUS <<_ACEOF || ac_write_fail=1 # # INIT-COMMANDS # -AMDEP_TRUE="$AMDEP_TRUE" ac_aux_dir="$ac_aux_dir" +AMDEP_TRUE="$AMDEP_TRUE" MAKE="${MAKE-make}" _ACEOF @@ -10049,6 +10201,7 @@ do case $ac_config_target in "depfiles") CONFIG_COMMANDS="$CONFIG_COMMANDS depfiles" ;; "Make.inc") CONFIG_FILES="$CONFIG_FILES Make.inc" ;; + "util/psb_metis_int.h") CONFIG_FILES="$CONFIG_FILES util/psb_metis_int.h" ;; *) as_fn_error $? "invalid argument: \`$ac_config_target'" "$LINENO" 5;; esac @@ -10482,29 +10635,35 @@ $as_echo "$as_me: executing $ac_file commands" >&6;} # Older Autoconf quotes --file arguments for eval, but not when files # are listed without --file. Let's play safe and only enable the eval # if we detect the quoting. - case $CONFIG_FILES in - *\'*) eval set x "$CONFIG_FILES" ;; - *) set x $CONFIG_FILES ;; - esac + # TODO: see whether this extra hack can be removed once we start + # requiring Autoconf 2.70 or later. + case $CONFIG_FILES in #( + *\'*) : + eval set x "$CONFIG_FILES" ;; #( + *) : + set x $CONFIG_FILES ;; #( + *) : + ;; +esac shift - for mf + # Used to flag and report bootstrapping failures. + am_rc=0 + for am_mf do # Strip MF so we end up with the name of the file. - mf=`echo "$mf" | sed -e 's/:.*$//'` - # Check whether this is an Automake generated Makefile or not. - # We used to match only the files named 'Makefile.in', but - # some people rename them; so instead we look at the file content. - # Grep'ing the first line is not enough: some people post-process - # each Makefile.in and add a new line on top of each file to say so. - # Grep'ing the whole file is not good either: AIX grep has a line + am_mf=`$as_echo "$am_mf" | sed -e 's/:.*$//'` + # Check whether this is an Automake generated Makefile which includes + # dependency-tracking related rules and includes. + # Grep'ing the whole file directly is not great: AIX grep has a line # limit of 2048, but all sed's we know have understand at least 4000. - if sed -n 's,^#.*generated by automake.*,X,p' "$mf" | grep X >/dev/null 2>&1; then - dirpart=`$as_dirname -- "$mf" || -$as_expr X"$mf" : 'X\(.*[^/]\)//*[^/][^/]*/*$' \| \ - X"$mf" : 'X\(//\)[^/]' \| \ - X"$mf" : 'X\(//\)$' \| \ - X"$mf" : 'X\(/\)' \| . 2>/dev/null || -$as_echo X"$mf" | + sed -n 's,^am--depfiles:.*,X,p' "$am_mf" | grep X >/dev/null 2>&1 \ + || continue + am_dirpart=`$as_dirname -- "$am_mf" || +$as_expr X"$am_mf" : 'X\(.*[^/]\)//*[^/][^/]*/*$' \| \ + X"$am_mf" : 'X\(//\)[^/]' \| \ + X"$am_mf" : 'X\(//\)$' \| \ + X"$am_mf" : 'X\(/\)' \| . 2>/dev/null || +$as_echo X"$am_mf" | sed '/^X\(.*[^/]\)\/\/*[^/][^/]*\/*$/{ s//\1/ q @@ -10522,53 +10681,48 @@ $as_echo X"$mf" | q } s/.*/./; q'` - else - continue - fi - # Extract the definition of DEPDIR, am__include, and am__quote - # from the Makefile without running 'make'. - DEPDIR=`sed -n 's/^DEPDIR = //p' < "$mf"` - test -z "$DEPDIR" && continue - am__include=`sed -n 's/^am__include = //p' < "$mf"` - test -z "$am__include" && continue - am__quote=`sed -n 's/^am__quote = //p' < "$mf"` - # Find all dependency output files, they are included files with - # $(DEPDIR) in their names. We invoke sed twice because it is the - # simplest approach to changing $(DEPDIR) to its actual value in the - # expansion. - for file in `sed -n " - s/^$am__include $am__quote\(.*(DEPDIR).*\)$am__quote"'$/\1/p' <"$mf" | \ - sed -e 's/\$(DEPDIR)/'"$DEPDIR"'/g'`; do - # Make sure the directory exists. - test -f "$dirpart/$file" && continue - fdir=`$as_dirname -- "$file" || -$as_expr X"$file" : 'X\(.*[^/]\)//*[^/][^/]*/*$' \| \ - X"$file" : 'X\(//\)[^/]' \| \ - X"$file" : 'X\(//\)$' \| \ - X"$file" : 'X\(/\)' \| . 2>/dev/null || -$as_echo X"$file" | - sed '/^X\(.*[^/]\)\/\/*[^/][^/]*\/*$/{ - s//\1/ - q - } - /^X\(\/\/\)[^/].*/{ + am_filepart=`$as_basename -- "$am_mf" || +$as_expr X/"$am_mf" : '.*/\([^/][^/]*\)/*$' \| \ + X"$am_mf" : 'X\(//\)$' \| \ + X"$am_mf" : 'X\(/\)' \| . 2>/dev/null || +$as_echo X/"$am_mf" | + sed '/^.*\/\([^/][^/]*\)\/*$/{ s//\1/ q } - /^X\(\/\/\)$/{ + /^X\/\(\/\/\)$/{ s//\1/ q } - /^X\(\/\).*/{ + /^X\/\(\/\).*/{ s//\1/ q } s/.*/./; q'` - as_dir=$dirpart/$fdir; as_fn_mkdir_p - # echo "creating $dirpart/$file" - echo '# dummy' > "$dirpart/$file" - done + { echo "$as_me:$LINENO: cd "$am_dirpart" \ + && sed -e '/# am--include-marker/d' "$am_filepart" \ + | $MAKE -f - am--depfiles" >&5 + (cd "$am_dirpart" \ + && sed -e '/# am--include-marker/d' "$am_filepart" \ + | $MAKE -f - am--depfiles) >&5 2>&5 + ac_status=$? + echo "$as_me:$LINENO: \$? = $ac_status" >&5 + (exit $ac_status); } || am_rc=$? done + if test $am_rc -ne 0; then + { { $as_echo "$as_me:${as_lineno-$LINENO}: error: in \`$ac_pwd':" >&5 +$as_echo "$as_me: error: in \`$ac_pwd':" >&2;} +as_fn_error $? "Something went wrong bootstrapping makefile fragments + for automatic dependency tracking. Try re-running configure with the + '--disable-dependency-tracking' option to at least be able to build + the package (albeit without support for automatic dependency tracking). +See \`config.log' for more details" "$LINENO" 5; } + fi + { am_dirpart=; unset am_dirpart;} + { am_filepart=; unset am_filepart;} + { am_mf=; unset am_mf;} + { am_rc=; unset am_rc;} + rm -f conftest-deps.mk } ;; diff --git a/configure.ac b/configure.ac index 18dbfdd3..9edcb290 100755 --- a/configure.ac +++ b/configure.ac @@ -36,11 +36,11 @@ dnl NOTE : There is no cross compilation support. ############################################################################### # NOTE: the literal for version (the second argument to AC_INIT should be a literal!) -AC_INIT([PSBLAS],3.5, [https://github.com/sfilippone/psblas3/issues]) +AC_INIT([PSBLAS],3.7.0, [https://github.com/sfilippone/psblas3/issues]) # VERSION is the file containing the PSBLAS version code # FIXME -psblas_cv_version="3.5" +psblas_cv_version="3.7.0" # A sample source file AC_CONFIG_SRCDIR([base/modules/psb_base_mod.f90]) @@ -561,13 +561,26 @@ PAC_FORTRAN_TEST_MOLD( Please get a Fortran compiler that supports it, e.g. GNU Fortran 4.8.])] ) +PAC_FORTRAN_TEST_VOLATILE( + [], + [AC_MSG_ERROR([Sorry, cannot build PSBLAS without support for VOLATILE])] +) + +PAC_FORTRAN_TEST_ISO_FORTRAN_ENV( + [], + [AC_MSG_ERROR([Sorry, cannot build PSBLAS without support for ISO_FORTRAN_ENV])] +) + +PAC_FORTRAN_TEST_FINAL( + [], + [AC_MSG_ERROR([Sorry, cannot build PSBLAS without support for FINAL])] +) + + # # Optional features # -PAC_FORTRAN_TEST_VOLATILE( - [FDEFINES="$psblas_cv_define_prepend-DHAVE_VOLATILE $FDEFINES"], -) PAC_FORTRAN_TEST_GENERICS( [], @@ -578,14 +591,6 @@ PAC_FORTRAN_TEST_FLUSH( [FDEFINES="$psblas_cv_define_prepend-DHAVE_FLUSH_STMT $FDEFINES"], ) -PAC_FORTRAN_TEST_ISO_FORTRAN_ENV( - [FDEFINES="$psblas_cv_define_prepend-DHAVE_ISO_FORTRAN_ENV $FDEFINES"], -) - -PAC_FORTRAN_TEST_FINAL( - [FDEFINES="$psblas_cv_define_prepend-DHAVE_FINAL $FDEFINES"], -) - ############################################################################### # Additional pathname stuff (yes, it is redundant and confusing...) @@ -689,9 +694,12 @@ LIBS="$RSB_LIBS ${LIBS}" dnl AC_CHECK_HEADERS([rsb.h], [ LIBS="${LIBS} $want_rsb_libs"], []) PAC_CHECK_METIS -if test "x$psblas_cv_have_metis" == "xyes" ; then - FDEFINES="$psblas_cv_define_prepend-DHAVE_METIS $FDEFINES" - CDEFINES="-DHAVE_METIS_ $psblas_cv_metis_includes $CDEFINES" +if test "x$pac_cv_metis_idx" == "xunknown" ; then + true ; # do nothing +elif test "x$psblas_cv_have_metis" == "xyes" ; then + FDEFINES="$psblas_cv_define_prepend-DHAVE_METIS $psblas_cv_define_prepend-DMETIS_$pac_cv_metis_idx $FDEFINES" + CDEFINES="-DHAVE_METIS_ $psblas_cv_metis_includes $CDEFINES -DMETIS_$pac_cv_metis_idx" + METISINCFILE=$psblas_cv_metisincfile fi PAC_CHECK_AMD @@ -702,11 +710,6 @@ fi -# Note : We would like to detect PSBLAS, but this is complicated by the -# module symbols mangling rules, which are compiler specific ! -# -# Moreover, the PSBLAS doesn't have an installer, currently. - ############################################################################### # Library target directory and archive files. ############################################################################### @@ -761,7 +764,6 @@ AC_SUBST(LAPACK_LIBS) AC_SUBST(CINCLUDES) AC_SUBST(FINCLUDES) - PSBLASRULES=' PSBLDLIBS=$(LAPACK) $(BLAS) $(METIS_LIB) $(AMD_LIB) $(LIBS) CDEFINES=$(PSBCDEFINES) @@ -784,11 +786,13 @@ AC_SUBST(BASELIBNAME) AC_SUBST(PRECLIBNAME) AC_SUBST(METHDLIBNAME) AC_SUBST(UTILLIBNAME) +AC_SUBST(METISINCFILE) ############################################################################### # the following files will be created by Automake AC_CONFIG_FILES([Make.inc]) +AC_CONFIG_FILES([util/psb_metis_int.h]) AC_OUTPUT() #AC_OUTPUT(Make.inc Makefile) ############################################################################### diff --git a/docs/html/cmsy10-42.png b/docs/html/cmsy10-42.png index 3e3ff34f..bd31f1f8 100644 Binary files a/docs/html/cmsy10-42.png and b/docs/html/cmsy10-42.png differ diff --git a/docs/html/cmsy10-48.png b/docs/html/cmsy10-48.png index b5d0fbf1..df31f928 100644 Binary files a/docs/html/cmsy10-48.png and b/docs/html/cmsy10-48.png differ diff --git a/docs/html/cmsy10-49.png b/docs/html/cmsy10-49.png index b44b91b4..cb292c09 100644 Binary files a/docs/html/cmsy10-49.png and b/docs/html/cmsy10-49.png differ diff --git a/docs/html/userhtml.css b/docs/html/userhtml.css index 201b1c23..86898d97 100644 --- a/docs/html/userhtml.css +++ b/docs/html/userhtml.css @@ -8,13 +8,11 @@ .cmbx-12x-x-144{font-size:172%; font-weight: bold;} .cmbx-12x-x-144{ font-weight: bold;} .cmbx-12x-x-144{ font-weight: bold;} -.cmbx-12x-x-144{ font-weight: bold;} .cmti-10{ font-style: italic;} .cmti-12{font-size:120%; font-style: italic;} .cmbx-10{ font-weight: bold;} .cmbx-10{ font-weight: bold;} .cmbx-10{ font-weight: bold;} -.cmbx-10{ font-weight: bold;} .cmtt-10{font-family: monospace;} .cmtt-10{font-family: monospace;} .cmtt-10{font-family: monospace;} @@ -23,7 +21,6 @@ .cmbx-12{font-size:120%; font-weight: bold;} .cmbx-12{ font-weight: bold;} .cmbx-12{ font-weight: bold;} -.cmbx-12{ font-weight: bold;} .cmtt-8{font-size:80%;font-family: monospace;} .cmtt-8{font-family: monospace;} .cmtt-8{font-family: monospace;} @@ -31,10 +28,6 @@ .cmtt-9{font-family: monospace;} .cmtt-9{font-family: monospace;} .cmmi-8{font-size:80%;font-style: italic;} -.cmbx-9{font-size:90%; font-weight: bold;} -.cmbx-9{ font-weight: bold;} -.cmbx-9{ font-weight: bold;} -.cmbx-9{ font-weight: bold;} p.noindent { text-indent: 0em } td p.noindent { text-indent: 0em; margin-top:0em; } p.nopar { text-indent: 0em; } @@ -104,8 +97,9 @@ span.LATEX span.TEX{ position:relative; left: -0.4em; } div.float, div.figure {margin-left: auto; margin-right: auto;} div.float img {text-align:center;} div.figure img {text-align:center;} -.marginpar {width:20%; float:right; text-align:left; margin-left:auto; margin-top:0.5em; font-size:85%; text-decoration:underline;} -.marginpar p{margin-top:0.4em; margin-bottom:0.4em;} +.marginpar,.reversemarginpar {width:20%; float:right; text-align:left; margin-left:auto; margin-top:0.5em; font-size:85%; text-decoration:underline;} +.marginpar p,.reversemarginpar p{margin-top:0.4em; margin-bottom:0.4em;} +.reversemarginpar{float:left;} table.equation {width:100%;} .equation td{text-align:center; } td.equation { margin-top:1em; margin-bottom:1em; } @@ -148,6 +142,7 @@ div.abstract {width:100%;} .Ovalbox-thick { padding-left:3pt; padding-right:3pt; border:solid thick; } .shadowbox { padding-left:3pt; padding-right:3pt; border:solid thin; border-right:solid thick; border-bottom:solid thick; } .doublebox { padding-left:3pt; padding-right:3pt; border-style:double; border:solid thick; } +.rotatebox{display: inline-block;} .figure img.graphics {margin-left:10%;} .lstlisting .label{margin-right:0.5em; } div.lstlisting{font-family: monospace; white-space: nowrap; margin-top:0.5em; margin-bottom:0.5em; } diff --git a/docs/html/userhtml0x.png b/docs/html/userhtml0x.png index 236023f4..bed5b304 100644 Binary files a/docs/html/userhtml0x.png and b/docs/html/userhtml0x.png differ diff --git a/docs/html/userhtml10x.png b/docs/html/userhtml10x.png index bc555490..a0d4d0d4 100644 Binary files a/docs/html/userhtml10x.png and b/docs/html/userhtml10x.png differ diff --git a/docs/html/userhtml11x.png b/docs/html/userhtml11x.png index 02e3ca83..634594e1 100644 Binary files a/docs/html/userhtml11x.png and b/docs/html/userhtml11x.png differ diff --git a/docs/html/userhtml12x.png b/docs/html/userhtml12x.png index 2fe0d741..7257667a 100644 Binary files a/docs/html/userhtml12x.png and b/docs/html/userhtml12x.png differ diff --git a/docs/html/userhtml13x.png b/docs/html/userhtml13x.png index 2640335f..4c70ecb3 100644 Binary files a/docs/html/userhtml13x.png and b/docs/html/userhtml13x.png differ diff --git a/docs/html/userhtml14x.png b/docs/html/userhtml14x.png index bc55f7d1..711ce998 100644 Binary files a/docs/html/userhtml14x.png and b/docs/html/userhtml14x.png differ diff --git a/docs/html/userhtml15x.png b/docs/html/userhtml15x.png index 0b01ec07..00e0d73d 100644 Binary files a/docs/html/userhtml15x.png and b/docs/html/userhtml15x.png differ diff --git a/docs/html/userhtml16x.png b/docs/html/userhtml16x.png index a18061f0..95794ed0 100644 Binary files a/docs/html/userhtml16x.png and b/docs/html/userhtml16x.png differ diff --git a/docs/html/userhtml17x.png b/docs/html/userhtml17x.png index 7fe28a82..27a7387e 100644 Binary files a/docs/html/userhtml17x.png and b/docs/html/userhtml17x.png differ diff --git a/docs/html/userhtml18x.png b/docs/html/userhtml18x.png index 82edd7f7..ebcb935d 100644 Binary files a/docs/html/userhtml18x.png and b/docs/html/userhtml18x.png differ diff --git a/docs/html/userhtml19x.png b/docs/html/userhtml19x.png index 336cb190..5ac0302a 100644 Binary files a/docs/html/userhtml19x.png and b/docs/html/userhtml19x.png differ diff --git a/docs/html/userhtml1x.png b/docs/html/userhtml1x.png index 00fc2e33..465c255b 100644 Binary files a/docs/html/userhtml1x.png and b/docs/html/userhtml1x.png differ diff --git a/docs/html/userhtml20x.png b/docs/html/userhtml20x.png index b19fd75c..f9f0de44 100644 Binary files a/docs/html/userhtml20x.png and b/docs/html/userhtml20x.png differ diff --git a/docs/html/userhtml21x.png b/docs/html/userhtml21x.png index c7b037cf..85201b8c 100644 Binary files a/docs/html/userhtml21x.png and b/docs/html/userhtml21x.png differ diff --git a/docs/html/userhtml22x.png b/docs/html/userhtml22x.png index 6b413586..8c882bc2 100644 Binary files a/docs/html/userhtml22x.png and b/docs/html/userhtml22x.png differ diff --git a/docs/html/userhtml23x.png b/docs/html/userhtml23x.png index ed4251c5..30c84d3d 100644 Binary files a/docs/html/userhtml23x.png and b/docs/html/userhtml23x.png differ diff --git a/docs/html/userhtml24x.png b/docs/html/userhtml24x.png index 3c6bb3bd..0910fe9c 100644 Binary files a/docs/html/userhtml24x.png and b/docs/html/userhtml24x.png differ diff --git a/docs/html/userhtml25x.png b/docs/html/userhtml25x.png index 37836c81..88630c50 100644 Binary files a/docs/html/userhtml25x.png and b/docs/html/userhtml25x.png differ diff --git a/docs/html/userhtml26x.png b/docs/html/userhtml26x.png index 75ae8f6b..ac0518e5 100644 Binary files a/docs/html/userhtml26x.png and b/docs/html/userhtml26x.png differ diff --git a/docs/html/userhtml27x.png b/docs/html/userhtml27x.png index 51d8b238..57beeb01 100644 Binary files a/docs/html/userhtml27x.png and b/docs/html/userhtml27x.png differ diff --git a/docs/html/userhtml28x.png b/docs/html/userhtml28x.png index 2ea69a48..b2106916 100644 Binary files a/docs/html/userhtml28x.png and b/docs/html/userhtml28x.png differ diff --git a/docs/html/userhtml29x.png b/docs/html/userhtml29x.png index b890fc38..64db17ce 100644 Binary files a/docs/html/userhtml29x.png and b/docs/html/userhtml29x.png differ diff --git a/docs/html/userhtml2x.png b/docs/html/userhtml2x.png index 31c7ec89..49dbf305 100644 Binary files a/docs/html/userhtml2x.png and b/docs/html/userhtml2x.png differ diff --git a/docs/html/userhtml30x.png b/docs/html/userhtml30x.png index 059fe3e0..509285f3 100644 Binary files a/docs/html/userhtml30x.png and b/docs/html/userhtml30x.png differ diff --git a/docs/html/userhtml31x.png b/docs/html/userhtml31x.png index f50412f3..72203987 100644 Binary files a/docs/html/userhtml31x.png and b/docs/html/userhtml31x.png differ diff --git a/docs/html/userhtml32x.png b/docs/html/userhtml32x.png index b4a91874..a92c5455 100644 Binary files a/docs/html/userhtml32x.png and b/docs/html/userhtml32x.png differ diff --git a/docs/html/userhtml3x.png b/docs/html/userhtml3x.png index 55f1e95a..3676aa53 100644 Binary files a/docs/html/userhtml3x.png and b/docs/html/userhtml3x.png differ diff --git a/docs/html/userhtml4x.png b/docs/html/userhtml4x.png index cdca4696..1ce118eb 100644 Binary files a/docs/html/userhtml4x.png and b/docs/html/userhtml4x.png differ diff --git a/docs/html/userhtml5x.png b/docs/html/userhtml5x.png index 403b9248..cabc66d4 100644 Binary files a/docs/html/userhtml5x.png and b/docs/html/userhtml5x.png differ diff --git a/docs/html/userhtml6x.png b/docs/html/userhtml6x.png index 60d3b2c5..a0ef30ce 100644 Binary files a/docs/html/userhtml6x.png and b/docs/html/userhtml6x.png differ diff --git a/docs/html/userhtml7x.png b/docs/html/userhtml7x.png index bc555490..a0d4d0d4 100644 Binary files a/docs/html/userhtml7x.png and b/docs/html/userhtml7x.png differ diff --git a/docs/html/userhtml8x.png b/docs/html/userhtml8x.png index c16766fd..147d6b03 100644 Binary files a/docs/html/userhtml8x.png and b/docs/html/userhtml8x.png differ diff --git a/docs/html/userhtml9x.png b/docs/html/userhtml9x.png index aaff2dea..dcf8438f 100644 Binary files a/docs/html/userhtml9x.png and b/docs/html/userhtml9x.png differ diff --git a/docs/html/userhtmlse8.html b/docs/html/userhtmlse8.html index 504703f8..775422c8 100644 --- a/docs/html/userhtmlse8.html +++ b/docs/html/userhtmlse8.html @@ -70,264 +70,254 @@ explicitly.
subroutine psb_foo(some args, info) 
subroutine psb_foo(some args, info) 
   ... 
   ... 
   if(error detected) then 
   if(error detected) then 
      info=errcode1 
      info=errcode1 
      call psb_errpush(’psb_foo’, errcode1) 
      call psb_errpush(’psb_foo’, errcode1) 
      goto 9999 
      goto 9999 
   end if 
   end if 
   ... 
   ... 
   call psb_bar(some args, info) 
   call psb_bar(some args, info) 
   if(info .ne. zero) then 
   if(info .ne. zero) then 
      info=errcode2 
      info=errcode2 
      call psb_errpush(’psb_foo’, errcode2) 
      call psb_errpush(’psb_foo’, errcode2) 
      goto 9999 
      goto 9999 
   end if 
   end if 
   ... 
   ... 
9999 continue 
9999 continue 
   if (err_act .eq. act_abort) then 
   if (err_act .eq. act_abort) then 
     call psb_error(icontxt) 
     call psb_error(icontxt) 
     return 
     return 
   else 
   else 
     return 
     return 
   end if 
   end if 
 
 
end subroutine psb_foo +class="cmtt-9">end subroutine psb_foo

Notes possible to improve the runtime efficiency by using the following scheme:
    vres(1) = psb_gedot(x1,y1,desc_a,info,global=.false.) 
    vres(2) = psb_gedot(x2,y2,desc_a,info,global=.false.) 
    vres(3) = psb_gedot(x3,y3,desc_a,info,global=.false.) 
    call psb_sum(ictxt,vres(1:3)) + id="x20-57004r1">   vres(1) = psb_gedot(x1,y1,desc_a,info,global=.false.) 
   vres(2) = psb_gedot(x2,y2,desc_a,info,global=.false.) 
   vres(3) = psb_gedot(x3,y3,desc_a,info,global=.false.) 
   call psb_sum(ictxt,vres(1:3))

In this way the global communication, which for small sizes is a latency-bound diff --git a/docs/html/userhtmlsu13.html b/docs/html/userhtmlsu13.html index db09267b..0cb73107 100644 --- a/docs/html/userhtmlsu13.html +++ b/docs/html/userhtmlsu13.html @@ -194,14 +194,84 @@ class="cmbx-12">Notes multiple norms at the same time; in this case, it is possible to improve the runtime efficiency by using the following scheme:

    vres(1) = psb_geamax(x1,desc_a,info,global=.false.) 
    vres(2) = psb_geamax(x2,desc_a,info,global=.false.) 
    vres(3) = psb_geamax(x3,desc_a,info,global=.false.) 
    call psb_amx(ictxt,vres(1:3)) + id="x22-59004r1">   vres(1) = psb_geamax(x1,desc_a,info,global=.false.) 
   vres(2) = psb_geamax(x2,desc_a,info,global=.false.) 
   vres(3) = psb_geamax(x3,desc_a,info,global=.false.) 
   call psb_amx(ictxt,vres(1:3))

In this way the global communication, which for small sizes is a latency-bound diff --git a/docs/html/userhtmlsu15.html b/docs/html/userhtmlsu15.html index 61096aa2..55c54a24 100644 --- a/docs/html/userhtmlsu15.html +++ b/docs/html/userhtmlsu15.html @@ -195,14 +195,84 @@ class="cmbx-12">Notes multiple norms at the same time; in this case, it is possible to improve the runtime efficiency by using the following scheme:

    vres(1) = psb_geasum(x1,desc_a,info,global=.false.) 
    vres(2) = psb_geasum(x2,desc_a,info,global=.false.) 
    vres(3) = psb_geasum(x3,desc_a,info,global=.false.) 
    call psb_sum(ictxt,vres(1:3)) + id="x24-61004r1">   vres(1) = psb_geasum(x1,desc_a,info,global=.false.) 
   vres(2) = psb_geasum(x2,desc_a,info,global=.false.) 
   vres(3) = psb_geasum(x3,desc_a,info,global=.false.) 
   call psb_sum(ictxt,vres(1:3))

In this way the global communication, which for small sizes is a latency-bound diff --git a/docs/html/userhtmlsu17.html b/docs/html/userhtmlsu17.html index f040eafb..44716f08 100644 --- a/docs/html/userhtmlsu17.html +++ b/docs/html/userhtmlsu17.html @@ -198,14 +198,84 @@ class="cmbx-12">Notes multiple norms at the same time; in this case, it is possible to improve the runtime efficiency by using the following scheme:

    vres(1) = psb_genrm2(x1,desc_a,info,global=.false.) 
    vres(2) = psb_genrm2(x2,desc_a,info,global=.false.) 
    vres(3) = psb_genrm2(x3,desc_a,info,global=.false.) 
    call psb_nrm2(ictxt,vres(1:3)) + id="x26-63004r1">   vres(1) = psb_genrm2(x1,desc_a,info,global=.false.) 
   vres(2) = psb_genrm2(x2,desc_a,info,global=.false.) 
   vres(3) = psb_genrm2(x3,desc_a,info,global=.false.) 
   call psb_nrm2(ictxt,vres(1:3))

In this way the global communication, which for small sizes is a latency-bound diff --git a/docs/html/userhtmlsu26.html b/docs/html/userhtmlsu26.html index 2287390e..d5fe2035 100644 --- a/docs/html/userhtmlsu26.html +++ b/docs/html/userhtmlsu26.html @@ -87,10 +87,38 @@ class="content">Data types

call psb_halo(x, desc_a, info) 
call psb_halo(x, desc_a, info) 
call psb_halo(x, desc_a, info, work, data) +class="cmtt-10">call psb_halo(x, desc_a, info, work, data)
diff --git a/docs/html/userhtmlsu27.html b/docs/html/userhtmlsu27.html index 2bb4e849..f0fe399b 100644 --- a/docs/html/userhtmlsu27.html +++ b/docs/html/userhtmlsu27.html @@ -92,9 +92,42 @@ class="content">Data types
call psb_ovrl(x, desc_a, info) 
call psb_ovrl(x, desc_a, info) 
call psb_ovrl(x, desc_a, info, update=update_type, work=work) +class="cmtt-10">call psb_ovrl(x, desc_a, info, update=update_type, work=work)
diff --git a/docs/html/userhtmlsu28.html b/docs/html/userhtmlsu28.html index 45fa580e..513f7e28 100644 --- a/docs/html/userhtmlsu28.html +++ b/docs/html/userhtmlsu28.html @@ -109,9 +109,44 @@ class="content">Data types
call psb_gather(glob_x, loc_x, desc_a, info, root) 
call psb_gather(glob_x, loc_x, desc_a, info, root) 
call psb_gather(glob_x, loc_x, desc_a, info, root) +class="cmtt-10">call psb_gather(glob_x, loc_x, desc_a, info, root)
diff --git a/docs/html/userhtmlsu29.html b/docs/html/userhtmlsu29.html index 7186e307..baeee944 100644 --- a/docs/html/userhtmlsu29.html +++ b/docs/html/userhtmlsu29.html @@ -108,7 +108,27 @@ class="content">Data types
call psb_scatter(glob_x, loc_x, desc_a, info, root, mold) +class="cmtt-10">call psb_scatter(glob_x, loc_x, desc_a, info, root, mold)
diff --git a/docs/html/userhtmlsu6.html b/docs/html/userhtmlsu6.html index 18b94c13..f0ba830f 100644 --- a/docs/html/userhtmlsu6.html +++ b/docs/html/userhtmlsu6.html @@ -45,8 +45,14 @@ class="description">Single precision complex;
Z
Double precision complex.
-

The actual data is contained in the polymorphic component Double precision complex; +

+LS,LD,LC,LZ
Same numeric type as above, but with psb_lpk_ integer + indices.
+

The actual data is contained in the polymorphic component a%a of type psb_spasb routine.

-

+

  type :: psb_Tspmat_type  
    class(psb_T_base_sparse_mat), allocatable  :: a  
  end type  psb_Tspmat_type
-

+


Figure 4: The PSBLAS defined data type that contains a sparse matrix.
-

-

The following very common formats are precompiled in PSBLAS and thus are +

+

The following very common formats are precompiled in PSBLAS and thus are always available:

_coo_sparse_mat
Coordinate storage; + + +
psb_csc_sparse_mat
Compressed storage by columns;
- - - -

The inner sparse matrix has an associated state, which can take the following +

The inner sparse matrix has an associated state, which can take the following values:

State entered after a reinitalization; this is used to handl in which the same sparsity pattern is used multiple times with different coefficients. In this state it is only possible to enter coefficients for already existing nonzero entries.
-

The only storage variant supporting the build state is COO; all other variants are +

The only storage variant supporting the build state is COO; all other variants are obtained by conversion to/from it. -

+

3.2.1 Sparse Matrix Methods
-

+

3.2.2 get_nrows — Get number of rows in a sparse matrix
@@ -135,8 +141,8 @@ obtained by conversion to/from it.
nr = a%get_nrows()
-

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
class="cmbx-10">Function value
The number of rows of sparse matrix a.
-

+

3.2.3 get_ncols — Get number of columns in a sparse matrix
@@ -171,8 +177,8 @@ class="cmtt-10">a.
nc = a%get_ncols()
-

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
class="cmbx-10">Function value
The number of columns of sparse matrix a.
-

+

3.2.4 get_nnzeros — Get number of nonzero elements in a sparse matrix
@@ -207,8 +213,8 @@ class="cmtt-10">a.
nz = a%get_nnzeros()
-

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
class="cmbx-10">Function value
The number of nonzero elements stored in sparse matrix a.
-

Notes

  1. Notes class="cmtt-10">a; some storage formats employ padding, thus the returned value for the same matrix may be different for different storage choices.
-

+

3.2.5 get_size — Get maximum number of nonzero elements in a sparse matrix
@@ -252,8 +258,8 @@ matrix
maxnz = a%get_size()
-

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
Function value
The maximum number of nonzero elements that can be stored in sparse matrix a using its current memory allocation.
-

+

3.2.6 sizeof — Get memory occupation in bytes of a sparse matrix
@@ -289,8 +295,8 @@ class="cmtt-10">a using its current memory allocation.
memory_size = a%sizeof()
-

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
Function value
The memory occupation in bytes.
-

+

3.2.7 get_fmt — Short description of the dynamic type
@@ -324,8 +330,8 @@ class="description">The memory occupation in bytes.
write(*,*) a%get_fmt()
-

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
NULL, COO, CSR and CSC.
-

+

3.2.8 is_bld, is_upd, is_asb — Status check
@@ -366,8 +372,8 @@ if (a%is_bld()) then  
if (a%is_upd()) then  
if (a%is_asb()) then -

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
Function value
A logical value indicating whether the matrix is in the Build, Update or Assembled state, respectively.
-

+

3.2.9 is_lower, is_upper, is_triangle, is_unit — Format check
@@ -406,8 +412,8 @@ if (a%is_triangle()) then  
if (a%is_lower()) then  
if (a%is_unit()) then -

-

+

+

Type:
the sparse matrix
Scope: local
-

+

On Return
logical value indicating whether the matrix class="cmtt-10">is_triangle() returns .true. check also if it is lower, upper and with a unit (i.e. assumed) diagonal.
-

+

3.2.10 cscnv — Convert to a different storage format
@@ -447,8 +453,8 @@ class="cmtt-10">.true. check also if it is lower, upper and call  a%cscnv(b,info [, type, mold, dupl])  
call  a%cscnv(info [, type, mold, dupl]) -

-

+

+

Type:
Type: optional. class="cmbx-10">dupl
an integer value specifing how to handle duplicates (see Named Constants below)
-

+

On Return
psb_Tspmat_type. info
Return code.
-

The

The mold arguments may be employed to interface with special devices, such as GPUs and other accelerators. -

+

3.2.11 csclip — Reduce to a submatrix
@@ -511,8 +517,8 @@ and other accelerators.     call a%csclip(b,info[,&  
       & imin,imax,jmin,jmax,rscale,cscale]) -

-

Returns the submatrix

+

Returns the submatrix A(imin:imax,jmin:jmax), optionally rescaling row/col indices to the range 1:imax-imin+1,1:jmax-jmin+1. @@ -542,7 +548,7 @@ class="newline" />Type: optional. rscale,cscale

Whether to rescale row/column indices. Type: optional.
-

+

On Return
psb_Tspmat_type. info
Return code.
-

+

3.2.12 clean_zeros — Eliminate zero coefficients
@@ -567,11 +573,11 @@ class="description">Return code.
    call a%clean_zeros(info)
-

-

Eliminates zero coefficients in the input matrix. Note that depending on the +

+

Eliminates zero coefficients in the input matrix. Note that depending on the internal storage format, there may still be some amount of zero padding in the output. -

+

Type:
psb_Tspmat_type.
Scope: local.
-

+

On Return
psb_Tspmat_type. info
Return code.
-

+

3.2.13 get_diag — Get main diagonal
@@ -614,8 +620,8 @@ class="description">Return code.
    call a%get_diag(d,info)
-

-

Returns a copy of the main diagonal. +

+

Returns a copy of the main diagonal.

Type:
psb_Tspmat_type.
Scope: local.
-

+

On Return
A one-dimensional array of the appropriate type. info
Return code.
-

+

3.2.14 clip_diag — Cut out main diagonal
@@ -656,8 +662,8 @@ class="description">Return code.
    call a%clip_diag(b,info)
-

-

Returns a copy of

+

Returns a copy of a without the main diagonal.

psb_Tspmat_type.
Scope: local.
-

+

On Return
psb_Tspmat_type. info
Return code.
-

+

3.2.15 tril — Return the lower triangle
@@ -702,8 +708,8 @@ class="description">Return code.     call a%tril(l,info[,&  
       & diag,imin,imax,jmin,jmax,rscale,cscale,u]) -

-

Returns the lower triangular part of submatrix

+

Returns the lower triangular part of submatrix A(imin:imax,jmin:jmax), optionally rescaling row/col indices to the range 1:imax-imin+1,1:jmax-jmin+1 and @@ -741,7 +747,7 @@ class="newline" />Type: optional. rscale,cscale

Whether to rescale row/column indices. Type: optional.
-

+

On Return
psb_Tspmat_type. info
Return code.
-

+

3.2.16 triu — Return the upper triangle
@@ -774,8 +780,8 @@ class="description">Return code.     call a%triu(u,info[,&  
       & diag,imin,imax,jmin,jmax,rscale,cscale,l]) -

-

Returns the upper triangular part of submatrix

+

Returns the upper triangular part of submatrix A(imin:imax,jmin:jmax), optionally rescaling row/col indices to the range 1:imax-imin+1,1:jmax-jmin+1, @@ -813,7 +819,7 @@ class="newline" />Type: optional. rscale,cscale

Whether to rescale row/column indices. Type: optional.
-

+

On Return
psb_Tspmat_type. info
Return code.
-

+

3.2.17 psb_set_mat_default — Set default storage format
@@ -845,8 +851,8 @@ class="description">Return code.
call  psb_set_mat_default(a)
-

-

+

+

Type:
a variable of class(psb_T_base_sparse_mat) requesting a new default storage format.
Type: required.
-

+

3.2.18 clone — Clone current object
@@ -871,8 +877,8 @@ class="newline" />Type: required.
call  a%clone(b,info)
-

-

+

+

Type:
the sparse matrix.
Scope: local.
-

+

On Return
A copy of the input object. info
Return code.
-

+

3.2.19 Named Constants
-

+

psbUpdate strategy based on additional permutation data (see -

[

[next] [prev] [prev-tail] [front] [up]

-

diff --git a/docs/html/userhtmlsu7.html b/docs/html/userhtmlsu7.html index 87e649d8..439e32c9 100644 --- a/docs/html/userhtmlsu7.html +++ b/docs/html/userhtmlsu7.html @@ -10,7 +10,7 @@ -

[

[next] [prev] [prev-tail] [tail] [up]

3.3 Dense Vector Data Structure

-

The

The psb_T6. Among other simple things, we define here an extraction method that can be used to get a full copy of the part of the vector stored on the local process. -

The type declaration is shown in figure 

The type declaration is shown in figure 5 where T is a placeholder for the data type and precision variants @@ -60,7 +60,7 @@ class="description">Single precision complex; Z

Double precision complex.
-

The actual data is contained in the polymorphic component

The actual data is contained in the polymorphic component v%v; the separation between the application and the actual data is essential for cases where it is necessary to link to data storage made available elsewhere outside the direct control of the @@ -70,7 +70,7 @@ compiler/application, e.g. data stored in a graphics accelerator’s private id="x15-460015">

-

+

  type psb_T_base_vect_type  
    TYPE(KIND_), allocatable :: v(:) @@ -81,18 +81,18 @@ compiler/application, e.g. data stored in a graphics accelerator’s private  
  end type  psb_T_vect_type  
-

+


Figure 5: The PSBLAS defined data type that contains a dense vector.
-

-

+

+

3.3.1 Vector Methods
-

+

3.3.2 get_nrows — Get number of rows in a dense vector
@@ -101,8 +101,8 @@ class="content"> The PSBLAS defined data type that contains a dense vector. nr = v%get_nrows() -

-

+

+

Type:
the dense vector
Scope: local
-

+

On Return
class="cmbx-10">Function value
The number of rows of dense vector v.
-

+

3.3.3 sizeof — Get memory occupation in bytes of a dense vector
@@ -137,8 +137,8 @@ class="cmtt-10">v.
memory_size = v%sizeof()
-

-

+

+

Type:
the dense vector
Scope: local
-

+

On Return
Function value
The memory occupation in bytes.
-

+

3.3.4 set — Set contents of the vector
@@ -174,8 +174,8 @@ class="description">The memory occupation in bytes.  
 call  v%set(vect[,first,last])  
 call  v%zero() -

-

+

+

Type:
Intent: in.
Specified as: a number of the data type indicated in Table 1.
-

Note that a call to

Note that a call to v%zero() is provided as a shorthand, but is equivalent to a call to v%set(zero) with the zero constant having the appropriate type and kind. -

+

On Return
-

+

3.3.5 get_vect — Get a copy of the vector contents
@@ -256,8 +256,8 @@ class="newline" />
extv = v%get_vect([n])
-

-

+

+

Type:
Type: optional; default: entire vector.
-

+

On Return
n and the internal size of the vector, or 0 if n is negative; otherwise, the size of the array is the same as the internal size of the vector.
-

+

3.3.6 clone — Clone current object
@@ -308,8 +308,8 @@ class="cmmi-10">n is negative; otherwise, the size of the array is the sa
call  x%clone(y,info)
-

-

+

+

Type:
the dense vector.
Scope: local.
-

+

On Return
Return code.
-

[

[next] [prev] [prev-tail] [front] [up]

-

diff --git a/docs/html/userhtmlsu74.html b/docs/html/userhtmlsu74.html index e9a4bf75..402007ce 100644 --- a/docs/html/userhtmlsu74.html +++ b/docs/html/userhtmlsu74.html @@ -19,7 +19,21 @@ href="userhtmlsu76.html#userhtmlsu77.html" >up]

call psb_errpush(err_c, r_name, i_err, a_err) +class="cmtt-10">call psb_errpush(err_c, r_name, i_err, a_err)

diff --git a/docs/html/userhtmlsu75.html b/docs/html/userhtmlsu75.html index b4726464..1b0b6bcd 100644 --- a/docs/html/userhtmlsu75.html +++ b/docs/html/userhtmlsu75.html @@ -21,7 +21,12 @@ href="userhtmlsu76.html#userhtmlse9.html" >up]

call psb_error(icontxt) +class="cmtt-10">call psb_error(icontxt)

diff --git a/docs/html/userhtmlsu76.html b/docs/html/userhtmlsu76.html index 218404fa..6ae8a04e 100644 --- a/docs/html/userhtmlsu76.html +++ b/docs/html/userhtmlsu76.html @@ -21,7 +21,12 @@ href="userhtmlsu76.html#userhtmlsu78.html" >up]

call psb_set_errverbosity(v) +class="cmtt-10">call psb_set_errverbosity(v)

diff --git a/docs/html/userhtmlsu77.html b/docs/html/userhtmlsu77.html index e1e62e1b..e3002b2c 100644 --- a/docs/html/userhtmlsu77.html +++ b/docs/html/userhtmlsu77.html @@ -21,7 +21,12 @@ condition

call psb_set_erraction(err_act) +class="cmtt-10">call psb_set_erraction(err_act)

diff --git a/docs/html/userhtmlsu78.html b/docs/html/userhtmlsu78.html index d7e6fd9f..d793b344 100644 --- a/docs/html/userhtmlsu78.html +++ b/docs/html/userhtmlsu78.html @@ -21,7 +21,27 @@ format

call hb_read(a, iret, iunit, filename, b, mtitle) +class="cmtt-10">call hb_read(a, iret, iunit, filename, b, mtitle)

diff --git a/docs/html/userhtmlsu79.html b/docs/html/userhtmlsu79.html index 47418d21..55d166d1 100644 --- a/docs/html/userhtmlsu79.html +++ b/docs/html/userhtmlsu79.html @@ -22,7 +22,30 @@ format

call hb_write(a, iret, iunit, filename, key, rhs, mtitle) +class="cmtt-10">call hb_write(a, iret, iunit, filename, key, rhs, mtitle)

diff --git a/docs/html/userhtmlsu8.html b/docs/html/userhtmlsu8.html index 5c0ff71c..5a02be22 100644 --- a/docs/html/userhtmlsu8.html +++ b/docs/html/userhtmlsu8.html @@ -10,7 +10,7 @@ -

[

[next] [prev] [prev-tail] [tail] [up]

3.4 Preconditioner data structure

-

Our base library offers support for simple well known preconditioners like Diagonal +

Our base library offers support for simple well known preconditioners like Diagonal Scaling or Block Jacobi with incomplete factorization ILU(0). -

A preconditioner is held in the

A preconditioner is held in the psb_prec

-

+

 
  type psb_Tprec_type  
    class(psb_T_base_prec_type), allocatable :: prec  
  end type psb_Tprec_type  
-

+


Figure 6: The PSBLAS defined data type that contains a preconditioner.
-

+

-

[

[next] [prev] [prev-tail] [front] [up]

-

diff --git a/docs/html/userhtmlsu80.html b/docs/html/userhtmlsu80.html index ee9c2c16..547efe1d 100644 --- a/docs/html/userhtmlsu80.html +++ b/docs/html/userhtmlsu80.html @@ -23,7 +23,21 @@ format

call mm_mat_read(a, iret, iunit, filename) +class="cmtt-10">call mm_mat_read(a, iret, iunit, filename)

diff --git a/docs/html/userhtmlsu81.html b/docs/html/userhtmlsu81.html index b74cc797..38d5c58f 100644 --- a/docs/html/userhtmlsu81.html +++ b/docs/html/userhtmlsu81.html @@ -23,7 +23,21 @@ format

call mm_array_read(b, iret, iunit, filename) +class="cmtt-10">call mm_array_read(b, iret, iunit, filename)

diff --git a/docs/html/userhtmlsu82.html b/docs/html/userhtmlsu82.html index 6b3014e0..8773da4a 100644 --- a/docs/html/userhtmlsu82.html +++ b/docs/html/userhtmlsu82.html @@ -23,7 +23,24 @@ format

call mm_mat_write(a, mtitle, iret, iunit, filename) +class="cmtt-10">call mm_mat_write(a, mtitle, iret, iunit, filename)

diff --git a/docs/html/userhtmlsu83.html b/docs/html/userhtmlsu83.html index f24fcc26..a82cb59f 100644 --- a/docs/html/userhtmlsu83.html +++ b/docs/html/userhtmlsu83.html @@ -22,7 +22,24 @@ format

call mm_array_write(b, vtitle, iret, iunit, filename) +class="cmtt-10">call mm_array_write(b, vtitle, iret, iunit, filename)

diff --git a/docs/html/userhtmlsu9.html b/docs/html/userhtmlsu9.html index e78edec3..311ae768 100644 --- a/docs/html/userhtmlsu9.html +++ b/docs/html/userhtmlsu9.html @@ -10,14 +10,14 @@ -

[

[prev] [prev-tail] [tail] [up]

3.5 Heap data structure

-

Among the tools routines of sec. 

Among the tools routines of sec. 6, we have a number of sorting utilities; the heap sort is implemented in terms of heaps having the following signatures:

@@ -36,7 +36,7 @@ class="cmtt-10">_idx_heap
: a heap containing elements of type T, as above, together with an integer index.
-

Given a heap object, the following methods are defined on it: +

Given a heap object, the following methods are defined on it:

init
Print on file; free
Release memory.
-

These objects are used in MLD2P4 to implement the factorization algorithms. +

These objects are used in MLD2P4 to implement the factorization algorithms. diff --git a/docs/html/userhtmlsu90.html b/docs/html/userhtmlsu90.html index 751fcd56..b9805544 100644 --- a/docs/html/userhtmlsu90.html +++ b/docs/html/userhtmlsu90.html @@ -63,9 +63,48 @@ iteration.

call psb_krylov(method,a,prec,b,x,eps,desc_a,info,& 
     & itmax,iter,err,itrace,irst,istop,cond) +class="cmtt-10">call psb_krylov(method,a,prec,b,x,eps,desc_a,info,& 
    & itmax,iter,err,itrace,irst,istop,cond)

diff --git a/docs/psblas-3.7.pdf b/docs/psblas-3.7.pdf index 94dab2e2..0da729c5 100644 --- a/docs/psblas-3.7.pdf +++ b/docs/psblas-3.7.pdf @@ -861,7 +861,7 @@ BT 0 g 0 G [-269(.)-500(.)]TJ 0 g 0 G - [-1084(17)]TJ + [-1084(18)]TJ 0 g 0 G 0 0 1 rg 0 0 1 RG -48.751 -12.091 Td [(3.2.4)-1144(ge)-1(t)]TJ @@ -6353,7 +6353,7 @@ endstream endobj 945 0 obj << -/Length 9167 +/Length 7975 >> stream 0 g 0 G @@ -6363,221 +6363,221 @@ BT /F27 9.9626 Tf 150.705 706.129 Td [(On)-383(Return)]TJ 0 g 0 G 0 g 0 G - 0 -19.686 Td [(iprc)]TJ + 0 -21.528 Td [(iprc)]TJ 0 g 0 G -/F8 9.9626 Tf 24.339 0 Td [(the)-333(list)-334(of)-333(pro)-28(cesses)-333(o)28(w)-1(n)1(ing)-334(the)-333(indices)-333(in)]TJ/F30 9.9626 Tf 183.783 0 Td [(idx)]TJ/F8 9.9626 Tf 15.691 0 Td [(.)]TJ -198.907 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.379 0 Td [(lo)-32(cal)]TJ/F8 9.9626 Tf 23.073 0 Td [(.)]TJ -55.452 -11.955 Td [(T)28(yp)-28(e:)]TJ/F27 9.9626 Tf 29.612 0 Td [(required)]TJ/F8 9.9626 Tf 41.898 0 Td [(.)]TJ -71.51 -11.955 Td [(In)28(ten)28(t:)]TJ/F27 9.9626 Tf 33.486 0 Td [(in)]TJ/F8 9.9626 Tf 9.547 0 Td [(.)]TJ -43.033 -11.955 Td [(Sp)-28(eci\014ed)-525(as:)-828(an)-526(allo)-27(catable)-526(one-dimensional)-525(arra)28(y)-525(of)-526(in)28(tegers)-525(of)-525(kind)]TJ/F30 9.9626 Tf 0 -11.955 Td [(psb_ipk_)]TJ/F8 9.9626 Tf 41.843 0 Td [(.)]TJ -66.749 -21.318 Td [(Note:)-417(th)1(is)-278(metho)-28(d)-277(ma)28(y)-278(or)-277(ma)27(y)-277(not)-278(actual)1(ly)-278(require)-277(comm)27(unications,)-288(dep)-28(ending)]TJ 0 -11.956 Td [(on)-466(the)-465(exact)-466(in)28(ternal)-466(data)-465(storage;)-532(giv)28(en)-466(that)-465(the)-466(c)28(hoice)-466(of)-466(storage)-465(ma)27(y)-465(b)-28(e)]TJ 0 -11.955 Td [(altered)-440(b)28(y)-440(run)28(time)-440(parameters,)-467(it)-439(is)-440(nece)-1(ssary)-439(for)-440(safet)27(y)-439(that)-440(this)-440(metho)-28(d)-440(is)]TJ 0 -11.955 Td [(called)-333(b)27(y)-333(all)-333(pro)-28(cesses.)]TJ/F27 9.9626 Tf 0 -25.793 Td [(3.1.15)-1150(Named)-383(Constan)32(ts)]TJ +/F8 9.9626 Tf 24.339 0 Td [(the)-333(list)-334(of)-333(pro)-28(cesses)-333(o)28(w)-1(n)1(ing)-334(the)-333(indices)-333(in)]TJ/F30 9.9626 Tf 183.783 0 Td [(idx)]TJ/F8 9.9626 Tf 15.691 0 Td [(.)]TJ -198.907 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.379 0 Td [(lo)-32(cal)]TJ/F8 9.9626 Tf 23.073 0 Td [(.)]TJ -55.452 -11.955 Td [(T)28(yp)-28(e:)]TJ/F27 9.9626 Tf 29.612 0 Td [(required)]TJ/F8 9.9626 Tf 41.898 0 Td [(.)]TJ -71.51 -11.955 Td [(In)28(ten)28(t:)]TJ/F27 9.9626 Tf 33.486 0 Td [(in)]TJ/F8 9.9626 Tf 9.547 0 Td [(.)]TJ -43.033 -11.956 Td [(Sp)-28(eci\014ed)-525(as:)-828(an)-526(allo)-27(catable)-526(one-dimensional)-525(arra)28(y)-525(of)-526(in)28(tegers)-525(of)-525(kind)]TJ/F30 9.9626 Tf 0 -11.955 Td [(psb_ipk_)]TJ/F8 9.9626 Tf 41.843 0 Td [(.)]TJ -66.749 -23.52 Td [(Note:)-417(th)1(is)-278(metho)-28(d)-277(ma)28(y)-278(or)-277(ma)27(y)-277(not)-278(actual)1(ly)-278(require)-277(comm)27(unications,)-288(dep)-28(ending)]TJ 0 -11.955 Td [(on)-466(the)-465(exact)-466(in)28(ternal)-466(data)-465(storage;)-532(giv)28(en)-466(that)-465(the)-466(c)28(hoice)-466(of)-466(storage)-465(ma)27(y)-465(b)-28(e)]TJ 0 -11.955 Td [(altered)-440(b)28(y)-440(run)28(time)-440(parameters,)-467(it)-439(is)-440(nece)-1(ssary)-439(for)-440(safet)27(y)-439(that)-440(this)-440(metho)-28(d)-440(is)]TJ 0 -11.956 Td [(called)-333(b)27(y)-333(all)-333(pro)-28(cesses.)]TJ/F27 9.9626 Tf 0 -28.021 Td [(3.1.15)-1150(Named)-383(Constan)32(ts)]TJ 0 g 0 G - 0 -18.389 Td [(psb)]TJ + 0 -19.135 Td [(psb)]TJ ET q -1 0 0 1 168.641 525.501 cm +1 0 0 1 168.641 518.482 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 525.302 Td [(none)]TJ +/F27 9.9626 Tf 172.078 518.283 Td [(none)]TJ ET q -1 0 0 1 196.475 525.501 cm +1 0 0 1 196.475 518.482 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 204.894 525.302 Td [(Generic)-333(no-op;)]TJ +/F8 9.9626 Tf 204.894 518.283 Td [(Generic)-333(no-op;)]TJ 0 g 0 G -/F27 9.9626 Tf -54.189 -19.685 Td [(psb)]TJ +/F27 9.9626 Tf -54.189 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 505.816 cm +1 0 0 1 168.641 496.955 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 505.617 Td [(ro)-32(ot)]TJ +/F27 9.9626 Tf 172.078 496.755 Td [(ro)-32(ot)]TJ ET q -1 0 0 1 193.715 505.816 cm +1 0 0 1 193.715 496.955 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 202.133 505.617 Td [(Default)-333(ro)-28(ot)-333(pro)-28(cess)-334(for)-333(broadcast)-333(and)-333(sc)-1(atter)-333(op)-28(erations;)]TJ +/F8 9.9626 Tf 202.133 496.755 Td [(Default)-333(ro)-28(ot)-333(pro)-28(cess)-334(for)-333(broadcast)-333(and)-333(sc)-1(atter)-333(op)-28(erations;)]TJ 0 g 0 G -/F27 9.9626 Tf -51.428 -19.686 Td [(psb)]TJ +/F27 9.9626 Tf -51.428 -21.527 Td [(psb)]TJ ET q -1 0 0 1 168.641 486.131 cm +1 0 0 1 168.641 475.427 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 485.931 Td [(nohalo)]TJ +/F27 9.9626 Tf 172.078 475.228 Td [(nohalo)]TJ ET q -1 0 0 1 205.705 486.131 cm +1 0 0 1 205.705 475.427 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 214.123 485.931 Td [(Do)-333(not)-334(fetc)28(h)-333(halo)-333(e)-1(l)1(e)-1(men)28(ts;)]TJ +/F8 9.9626 Tf 214.123 475.228 Td [(Do)-333(not)-334(fetc)28(h)-333(halo)-333(e)-1(l)1(e)-1(men)28(ts;)]TJ 0 g 0 G -/F27 9.9626 Tf -63.418 -19.685 Td [(psb)]TJ +/F27 9.9626 Tf -63.418 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 466.445 cm +1 0 0 1 168.641 453.899 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 466.246 Td [(halo)]TJ +/F27 9.9626 Tf 172.078 453.7 Td [(halo)]TJ ET q -1 0 0 1 193.611 466.445 cm +1 0 0 1 193.611 453.899 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 202.03 466.246 Td [(F)83(etc)28(h)-333(halo)-333(e)-1(lemen)28(ts)-333(from)-333(neigh)27(b)-27(ouring)-334(pro)-27(cesse)-1(s;)]TJ +/F8 9.9626 Tf 202.03 453.7 Td [(F)83(etc)28(h)-333(halo)-333(e)-1(lemen)28(ts)-333(from)-333(neigh)27(b)-27(ouring)-334(pro)-27(cesse)-1(s;)]TJ 0 g 0 G -/F27 9.9626 Tf -51.325 -19.686 Td [(psb)]TJ +/F27 9.9626 Tf -51.325 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 446.76 cm +1 0 0 1 168.641 432.371 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 446.56 Td [(sum)]TJ +/F27 9.9626 Tf 172.078 432.172 Td [(sum)]TJ ET q -1 0 0 1 193.197 446.76 cm +1 0 0 1 193.197 432.371 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 201.616 446.56 Td [(Sum)-333(o)27(v)28(erlapp)-27(ed)-334(elemen)28(ts)]TJ +/F8 9.9626 Tf 201.616 432.172 Td [(Sum)-333(o)27(v)28(erlapp)-27(ed)-334(elemen)28(ts)]TJ 0 g 0 G -/F27 9.9626 Tf -50.911 -19.685 Td [(psb)]TJ +/F27 9.9626 Tf -50.911 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 427.074 cm +1 0 0 1 168.641 410.844 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 426.875 Td [(a)32(vg)]TJ +/F27 9.9626 Tf 172.078 410.644 Td [(a)32(vg)]TJ ET q -1 0 0 1 189.792 427.074 cm +1 0 0 1 189.792 410.844 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 198.211 426.875 Td [(Av)28(erage)-334(o)28(v)28(erlapp)-28(ed)-333(elemen)28(ts)]TJ +/F8 9.9626 Tf 198.211 410.644 Td [(Av)28(erage)-334(o)28(v)28(erlapp)-28(ed)-333(elemen)28(ts)]TJ 0 g 0 G -/F27 9.9626 Tf -47.506 -19.685 Td [(psb)]TJ +/F27 9.9626 Tf -47.506 -21.527 Td [(psb)]TJ ET q -1 0 0 1 168.641 407.389 cm +1 0 0 1 168.641 389.316 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 407.19 Td [(comm)]TJ +/F27 9.9626 Tf 172.078 389.117 Td [(comm)]TJ ET q -1 0 0 1 202.681 407.389 cm +1 0 0 1 202.681 389.316 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 206.118 407.19 Td [(halo)]TJ +/F27 9.9626 Tf 206.118 389.117 Td [(halo)]TJ ET q -1 0 0 1 227.651 407.389 cm +1 0 0 1 227.651 389.316 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 236.07 407.19 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(halo_index)]TJ/F8 9.9626 Tf 55.624 0 Td [(list;)]TJ +/F8 9.9626 Tf 236.07 389.117 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(halo_index)]TJ/F8 9.9626 Tf 55.624 0 Td [(list;)]TJ 0 g 0 G -/F27 9.9626 Tf -267.376 -19.686 Td [(psb)]TJ +/F27 9.9626 Tf -267.376 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 387.704 cm +1 0 0 1 168.641 367.788 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 387.504 Td [(comm)]TJ +/F27 9.9626 Tf 172.078 367.589 Td [(comm)]TJ ET q -1 0 0 1 202.681 387.704 cm +1 0 0 1 202.681 367.788 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 206.118 387.504 Td [(ext)]TJ +/F27 9.9626 Tf 206.118 367.589 Td [(ext)]TJ ET q -1 0 0 1 222.559 387.704 cm +1 0 0 1 222.559 367.788 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 230.978 387.504 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(ext_index)]TJ/F8 9.9626 Tf 50.394 0 Td [(list;)]TJ +/F8 9.9626 Tf 230.978 367.589 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(ext_index)]TJ/F8 9.9626 Tf 50.394 0 Td [(list;)]TJ 0 g 0 G -/F27 9.9626 Tf -257.054 -19.685 Td [(psb)]TJ +/F27 9.9626 Tf -257.054 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 368.018 cm +1 0 0 1 168.641 346.26 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 367.819 Td [(comm)]TJ +/F27 9.9626 Tf 172.078 346.061 Td [(comm)]TJ ET q -1 0 0 1 202.681 368.018 cm +1 0 0 1 202.681 346.26 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 206.118 367.819 Td [(o)32(vr)]TJ +/F27 9.9626 Tf 206.118 346.061 Td [(o)32(vr)]TJ ET q -1 0 0 1 222.981 368.018 cm +1 0 0 1 222.981 346.26 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 231.4 367.819 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(ovrlap_index)]TJ/F8 9.9626 Tf 66.085 0 Td [(list;)]TJ +/F8 9.9626 Tf 231.4 346.061 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(ovrlap_index)]TJ/F8 9.9626 Tf 66.085 0 Td [(list;)]TJ 0 g 0 G -/F27 9.9626 Tf -273.167 -19.685 Td [(psb)]TJ +/F27 9.9626 Tf -273.167 -21.528 Td [(psb)]TJ ET q -1 0 0 1 168.641 348.333 cm +1 0 0 1 168.641 324.733 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 172.078 348.134 Td [(comm)]TJ +/F27 9.9626 Tf 172.078 324.533 Td [(comm)]TJ ET q -1 0 0 1 202.681 348.333 cm +1 0 0 1 202.681 324.733 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 206.118 348.134 Td [(mo)32(v)]TJ +/F27 9.9626 Tf 206.118 324.533 Td [(mo)32(v)]TJ ET q -1 0 0 1 227.81 348.333 cm +1 0 0 1 227.81 324.733 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q 0 g 0 G BT -/F8 9.9626 Tf 236.229 348.134 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(ovr_mst_idx)]TJ/F8 9.9626 Tf 60.855 0 Td [(list;)]TJ/F16 11.9552 Tf -272.766 -27.786 Td [(3.2)-1125(Sparse)-375(Matrix)-375(class)]TJ/F8 9.9626 Tf 0 -18.389 Td [(The)]TJ/F30 9.9626 Tf 20.653 0 Td [(psb)]TJ +/F8 9.9626 Tf 236.229 324.533 Td [(Exc)28(hange)-333(data)-334(based)-333(on)-333(the)]TJ/F30 9.9626 Tf 126.387 0 Td [(ovr_mst_idx)]TJ/F8 9.9626 Tf 60.855 0 Td [(list;)]TJ/F16 11.9552 Tf -272.766 -30.014 Td [(3.2)-1125(Sparse)-375(Matrix)-375(class)]TJ/F8 9.9626 Tf 0 -19.135 Td [(The)]TJ/F30 9.9626 Tf 20.653 0 Td [(psb)]TJ ET q -1 0 0 1 187.676 302.158 cm +1 0 0 1 187.676 275.584 cm []0 d 0 J 0.398 w 0 0 m 3.138 0 l S Q BT -/F30 9.9626 Tf 190.814 301.959 Td [(Tspmat)]TJ +/F30 9.9626 Tf 190.814 275.384 Td [(Tspmat)]TJ ET q -1 0 0 1 222.824 302.158 cm +1 0 0 1 222.824 275.584 cm []0 d 0 J 0.398 w 0 0 m 3.138 0 l S Q BT -/F30 9.9626 Tf 225.962 301.959 Td [(type)]TJ/F8 9.9626 Tf 24.416 0 Td [(class)-351(con)28(tains)-351(all)-351(information)-350(ab)-28(out)-351(the)-351(lo)-27(cal)-351(p)-28(ortion)-351(of)]TJ -99.673 -11.955 Td [(the)-249(sparse)-249(matrix)-248(and)-249(its)-249(storage)-249(mo)-27(de.)-417(Its)-248(design)-249(is)-249(based)-249(on)-248(the)-249(ST)83(A)84(TE)-249(design)]TJ 0 -11.955 Td [(pattern)-347([)]TJ +/F30 9.9626 Tf 225.962 275.384 Td [(type)]TJ/F8 9.9626 Tf 24.416 0 Td [(class)-351(con)28(tains)-351(all)-351(information)-350(ab)-28(out)-351(the)-351(lo)-27(cal)-351(p)-28(ortion)-351(of)]TJ -99.673 -11.955 Td [(the)-249(sparse)-249(matrix)-248(and)-249(its)-249(storage)-249(mo)-27(de.)-417(Its)-248(design)-249(is)-249(based)-249(on)-248(the)-249(ST)83(A)84(TE)-249(design)]TJ 0 -11.955 Td [(pattern)-347([)]TJ 1 0 0 rg 1 0 0 RG [(13)]TJ 0 g 0 G @@ -6589,61 +6589,37 @@ BT 0 0 1 rg 0 0 1 RG [-347(4)]TJ 0 g 0 G - [-347(where)]TJ/F30 9.9626 Tf 0 -11.956 Td [(T)]TJ/F8 9.9626 Tf 8.551 0 Td [(is)-333(a)-334(placeholder)-333(for)-333(the)-334(data)-333(t)28(yp)-28(e)-333(and)-333(precision)-334(v)56(arian)28(ts)]TJ + [-347(where)]TJ/F30 9.9626 Tf 0 -11.955 Td [(T)]TJ/F8 9.9626 Tf 8.551 0 Td [(is)-333(a)-334(placeholder)-333(for)-333(the)-334(data)-333(t)28(yp)-28(e)-333(and)-333(precision)-334(v)56(arian)28(ts)]TJ 0 g 0 G -/F27 9.9626 Tf -8.551 -19.445 Td [(S)]TJ +/F27 9.9626 Tf -8.551 -21.127 Td [(S)]TJ 0 g 0 G /F8 9.9626 Tf 11.346 0 Td [(Single)-333(precision)-334(real;)]TJ 0 g 0 G -/F27 9.9626 Tf -11.346 -19.685 Td [(D)]TJ +/F27 9.9626 Tf -11.346 -21.528 Td [(D)]TJ 0 g 0 G /F8 9.9626 Tf 13.768 0 Td [(Double)-333(precision)-334(real;)]TJ 0 g 0 G -/F27 9.9626 Tf -13.768 -19.686 Td [(C)]TJ +/F27 9.9626 Tf -13.768 -21.528 Td [(C)]TJ 0 g 0 G /F8 9.9626 Tf 13.256 0 Td [(Single)-333(precision)-334(complex;)]TJ 0 g 0 G -/F27 9.9626 Tf -13.256 -19.685 Td [(Z)]TJ +/F27 9.9626 Tf -13.256 -21.527 Td [(Z)]TJ 0 g 0 G -/F8 9.9626 Tf 11.983 0 Td [(Double)-333(precision)-334(complex.)]TJ -11.983 -19.446 Td [(The)-222(actual)-222(data)-223(is)-222(con)28(tained)-222(in)-222(the)-223(p)-27(olymorphic)-223(comp)-27(onen)27(t)]TJ/F30 9.9626 Tf 255.514 0 Td [(a%a)]TJ/F8 9.9626 Tf 17.905 0 Td [(of)-222(t)28(yp)-28(e)]TJ/F30 9.9626 Tf 31.549 0 Td [(psb)]TJ -ET -q -1 0 0 1 471.991 168.346 cm -[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S -Q -BT -/F30 9.9626 Tf 475.13 168.146 Td [(T)]TJ -ET -q -1 0 0 1 480.988 168.346 cm -[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S -Q -BT -/F30 9.9626 Tf 484.126 168.146 Td [(base)]TJ -ET -q -1 0 0 1 505.675 168.346 cm -[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S -Q -BT -/F30 9.9626 Tf 508.813 168.146 Td [(sparse)]TJ -ET -q -1 0 0 1 540.823 168.346 cm -[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S -Q -BT -/F30 9.9626 Tf 543.961 168.146 Td [(mat)]TJ/F8 9.9626 Tf 15.691 0 Td [(;)]TJ -408.947 -11.955 Td [(its)-300(sp)-28(eci\014c)-301(la)28(y)28(out)-300(can)-301(b)-27(e)-301(c)28(hosen)-301(d)1(ynamically)-301(among)-300(the)-301(pr)1(e)-1(d)1(e)-1(\014n)1(e)-1(d)-300(t)28(yp)-28(es,)-307(or)-300(an)]TJ 0 -11.955 Td [(en)28(tirely)-419(new)-419(storage)-419(la)28(y)27(out)-419(can)-419(b)-27(e)-419(implemen)27(ted)-419(an)1(d)-419(passe)-1(d)-418(to)-419(the)-419(library)-419(at)]TJ 0 -11.955 Td [(run)28(time)-420(via)-419(the)]TJ/F30 9.9626 Tf 73.447 0 Td [(psb_spasb)]TJ/F8 9.9626 Tf 51.252 0 Td [(routine.)-703(The)-419(follo)28(wing)-420(v)28(ery)-419(common)-420(formats)-419(are)]TJ -124.699 -11.955 Td [(precompiled)-333(in)-334(PSBLAS)-333(and)-333(th)28(us)-334(are)-333(alw)28(a)28(ys)-334(a)28(v)56(ailable:)]TJ +/F8 9.9626 Tf 11.983 0 Td [(Double)-333(precision)-334(complex;)]TJ +0 g 0 G +/F27 9.9626 Tf -11.983 -21.528 Td [(LS,LD,LC,LZ)]TJ +0 g 0 G +/F8 9.9626 Tf 72.519 0 Td [(Same)-446(n)28(umeric)-445(t)27(y)1(p)-28(e)-446(as)-445(ab)-28(o)28(v)28(e)-1(,)-473(but)-445(with)]TJ/F30 9.9626 Tf 179.211 0 Td [(psb_lpk_)]TJ/F8 9.9626 Tf 46.281 0 Td [(in)28(teger)-446(in-)]TJ -273.104 -11.955 Td [(dices.)]TJ 0 g 0 G - 166.874 -29.888 Td [(16)]TJ + 141.967 -29.888 Td [(16)]TJ 0 g 0 G ET endstream endobj -953 0 obj +952 0 obj << -/Length 5963 +/Length 6883 >> stream 0 g 0 G @@ -6659,161 +6635,173 @@ BT /F8 9.9626 Tf -24.739 -30.054 Td [(Figure)-333(4:)-889(The)-333(PSBLAS)-334(de\014ned)-333(data)-333(t)28(yp)-28(e)-334(that)-333(con)28(tains)-333(a)-334(sparse)-333(matrix.)]TJ 0 g 0 G 0 g 0 G + -11.499 -32.624 Td [(The)-222(actual)-222(data)-223(is)-222(con)28(tained)-222(in)-222(the)-223(p)-27(olymorphic)-223(comp)-27(onen)27(t)]TJ/F30 9.9626 Tf 255.515 0 Td [(a%a)]TJ/F8 9.9626 Tf 17.905 0 Td [(of)-222(t)28(yp)-28(e)]TJ/F30 9.9626 Tf 31.548 0 Td [(psb)]TJ +ET +q +1 0 0 1 421.182 623.614 cm +[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S +Q +BT +/F30 9.9626 Tf 424.32 623.415 Td [(T)]TJ +ET +q +1 0 0 1 430.178 623.614 cm +[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S +Q +BT +/F30 9.9626 Tf 433.316 623.415 Td [(base)]TJ +ET +q +1 0 0 1 454.865 623.614 cm +[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S +Q +BT +/F30 9.9626 Tf 458.003 623.415 Td [(sparse)]TJ +ET +q +1 0 0 1 490.013 623.614 cm +[]0 d 0 J 0.398 w 0 0 m 3.138 0 l S +Q +BT +/F30 9.9626 Tf 493.151 623.415 Td [(mat)]TJ/F8 9.9626 Tf 15.691 0 Td [(;)]TJ -408.947 -11.955 Td [(its)-300(sp)-28(eci\014c)-301(la)28(y)28(out)-300(can)-301(b)-28(e)-300(c)28(hosen)-301(dyn)1(am)-1(ically)-300(among)-300(the)-301(prede\014ned)-300(t)28(yp)-28(es,)-307(or)-300(an)]TJ 0 -11.956 Td [(en)28(tirely)-419(new)-419(storage)-419(la)28(y)27(out)-419(can)-419(b)-27(e)-419(implemen)27(ted)-419(and)-418(pass)-1(ed)-418(to)-419(the)-419(library)-419(at)]TJ 0 -11.955 Td [(run)28(time)-420(via)-419(the)]TJ/F30 9.9626 Tf 73.447 0 Td [(psb_spasb)]TJ/F8 9.9626 Tf 51.252 0 Td [(routine.)-703(The)-419(follo)28(wing)-420(v)28(ery)-419(common)-420(formats)-419(are)]TJ -124.699 -11.955 Td [(precompiled)-333(in)-334(PSBLAS)-333(and)-333(th)28(us)-334(are)-333(alw)28(a)28(ys)-334(a)28(v)56(ailable:)]TJ 0 g 0 G -/F27 9.9626 Tf -11.499 -31.954 Td [(psb)]TJ +/F27 9.9626 Tf 0 -21.124 Td [(psb)]TJ ET q -1 0 0 1 117.832 624.284 cm +1 0 0 1 117.832 554.67 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 121.269 624.085 Td [(T)]TJ +/F27 9.9626 Tf 121.269 554.47 Td [(T)]TJ ET q -1 0 0 1 129.926 624.284 cm +1 0 0 1 129.926 554.67 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 133.364 624.085 Td [(co)-32(o)]TJ +/F27 9.9626 Tf 133.364 554.47 Td [(co)-32(o)]TJ ET q -1 0 0 1 150.918 624.284 cm +1 0 0 1 150.918 554.67 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 154.355 624.085 Td [(sparse)]TJ +/F27 9.9626 Tf 154.355 554.47 Td [(sparse)]TJ ET q -1 0 0 1 185.985 624.284 cm +1 0 0 1 185.985 554.67 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 189.422 624.085 Td [(mat)]TJ +/F27 9.9626 Tf 189.422 554.47 Td [(mat)]TJ 0 g 0 G /F8 9.9626 Tf 24.554 0 Td [(Co)-28(ordinate)-333(storage;)]TJ 0 g 0 G -/F27 9.9626 Tf -114.081 -20.184 Td [(psb)]TJ +/F27 9.9626 Tf -114.081 -21.523 Td [(psb)]TJ ET q -1 0 0 1 117.832 604.101 cm +1 0 0 1 117.832 533.146 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 121.269 603.901 Td [(T)]TJ +/F27 9.9626 Tf 121.269 532.947 Td [(T)]TJ ET q -1 0 0 1 129.926 604.101 cm +1 0 0 1 129.926 533.146 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 133.364 603.901 Td [(csr)]TJ +/F27 9.9626 Tf 133.364 532.947 Td [(csr)]TJ ET q -1 0 0 1 148.38 604.101 cm +1 0 0 1 148.38 533.146 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 151.818 603.901 Td [(sparse)]TJ +/F27 9.9626 Tf 151.818 532.947 Td [(sparse)]TJ ET q -1 0 0 1 183.447 604.101 cm +1 0 0 1 183.447 533.146 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 186.884 603.901 Td [(mat)]TJ +/F27 9.9626 Tf 186.884 532.947 Td [(mat)]TJ 0 g 0 G /F8 9.9626 Tf 24.554 0 Td [(Compressed)-333(storage)-334(b)28(y)-333(ro)27(ws;)]TJ 0 g 0 G -/F27 9.9626 Tf -111.543 -20.183 Td [(psb)]TJ +/F27 9.9626 Tf -111.543 -21.523 Td [(psb)]TJ ET q -1 0 0 1 117.832 583.917 cm +1 0 0 1 117.832 511.623 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 121.269 583.718 Td [(T)]TJ +/F27 9.9626 Tf 121.269 511.424 Td [(T)]TJ ET q -1 0 0 1 129.926 583.917 cm +1 0 0 1 129.926 511.623 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 133.364 583.718 Td [(csc)]TJ +/F27 9.9626 Tf 133.364 511.424 Td [(csc)]TJ ET q -1 0 0 1 148.754 583.917 cm +1 0 0 1 148.754 511.623 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 152.191 583.718 Td [(sparse)]TJ +/F27 9.9626 Tf 152.191 511.424 Td [(sparse)]TJ ET q -1 0 0 1 183.821 583.917 cm +1 0 0 1 183.821 511.623 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 187.258 583.718 Td [(mat)]TJ +/F27 9.9626 Tf 187.258 511.424 Td [(mat)]TJ 0 g 0 G -/F8 9.9626 Tf 24.553 0 Td [(Compressed)-334(storage)-333(b)28(y)-333(columns;)]TJ -111.916 -20.119 Td [(The)-373(inner)-373(sparse)-373(matrix)-373(has)-373(an)-373(asso)-28(ciated)-373(state,)-383(whic)28(h)-373(can)-373(tak)28(e)-373(the)-373(follo)27(wing)]TJ 0 -11.955 Td [(v)56(alues:)]TJ +/F8 9.9626 Tf 24.553 0 Td [(Compressed)-334(storage)-333(b)28(y)-333(columns;)]TJ -111.916 -21.124 Td [(The)-373(inner)-373(sparse)-373(matrix)-373(has)-373(an)-373(asso)-28(ciated)-373(state,)-383(whic)28(h)-373(can)-373(tak)28(e)-373(the)-373(follo)27(wing)]TJ 0 -11.955 Td [(v)56(alues:)]TJ 0 g 0 G -/F27 9.9626 Tf 0 -20.119 Td [(Build:)]TJ +/F27 9.9626 Tf 0 -21.124 Td [(Build:)]TJ 0 g 0 G /F8 9.9626 Tf 35.409 0 Td [(State)-306(en)28(tered)-306(after)-307(the)-306(\014rst)-306(allo)-28(cation,)-311(and)-306(b)-28(efore)-306(the)-306(\014rst)-306(assem)27(bly;)-315(in)]TJ -10.502 -11.955 Td [(this)-333(state)-334(it)-333(is)-333(p)-28(ossible)-334(to)-333(add)-333(nonzero)-333(en)27(tries.)]TJ 0 g 0 G -/F27 9.9626 Tf -24.907 -20.183 Td [(Assem)32(bled:)]TJ +/F27 9.9626 Tf -24.907 -21.523 Td [(Assem)32(bled:)]TJ 0 g 0 G /F8 9.9626 Tf 61.508 0 Td [(State)-373(en)27(tered)-373(after)-373(the)-374(assem)28(bly;)-393(computations)-373(us)-1(i)1(ng)-374(the)-373(sparse)]TJ -36.601 -11.956 Td [(matrix,)-333(suc)27(h)-333(as)-333(matrix-v)28(ec)-1(tor)-333(pro)-28(d)1(ucts)-1(,)-333(are)-333(only)-333(p)-28(ossible)-334(in)-333(this)-333(state;)]TJ 0 g 0 G -/F27 9.9626 Tf -24.907 -20.183 Td [(Up)-32(date:)]TJ +/F27 9.9626 Tf -24.907 -21.523 Td [(Up)-32(date:)]TJ 0 g 0 G -/F8 9.9626 Tf 45.302 0 Td [(State)-233(en)27(tered)-233(after)-233(a)-234(reinit)1(aliz)-1(at)1(ion;)-267(this)-233(is)-234(used)-233(to)-233(handle)-234(applications)]TJ -20.395 -11.955 Td [(in)-395(whic)28(h)-396(th)1(e)-396(same)-395(sparsit)28(y)-395(pattern)-396(is)-395(used)-395(m)28(ultiple)-395(times)-396(with)-395(di\013eren)28(t)]TJ 0 -11.955 Td [(co)-28(e\016cien)28(ts.)-427(In)-280(this)-280(state)-280(it)-280(is)-281(only)-280(p)-27(oss)-1(ib)1(le)-281(to)-280(en)28(ter)-280(co)-28(e\016cien)28(ts)-280(for)-281(already)]TJ 0 -11.956 Td [(existing)-333(nonzero)-334(en)28(tries.)]TJ -24.907 -20.118 Td [(The)-358(only)-357(storage)-358(v)56(arian)27(t)-357(supp)-28(orting)-357(the)-358(build)-357(s)-1(tate)-357(is)-358(COO;)-357(all)-358(other)-358(v)56(arian)28(ts)]TJ 0 -11.956 Td [(are)-333(obtained)-334(b)28(y)-333(con)28(v)27(ersion)-333(to/from)-333(it.)]TJ/F27 9.9626 Tf 0 -26.238 Td [(3.2.1)-1150(Sparse)-383(Matri)-1(x)-383(Metho)-32(ds)]TJ 0 -18.509 Td [(3.2.2)-1150(get)]TJ +/F8 9.9626 Tf 45.302 0 Td [(State)-233(en)27(tered)-233(after)-233(a)-234(reinit)1(aliz)-1(at)1(ion;)-267(this)-233(is)-234(used)-233(to)-233(handle)-234(applications)]TJ -20.395 -11.955 Td [(in)-395(whic)28(h)-396(th)1(e)-396(same)-395(sparsit)28(y)-395(pattern)-396(is)-395(used)-395(m)28(ultiple)-395(times)-396(with)-395(di\013eren)28(t)]TJ 0 -11.955 Td [(co)-28(e\016cien)28(ts.)-427(In)-280(this)-280(state)-280(it)-280(is)-281(only)-280(p)-27(oss)-1(ib)1(le)-281(to)-280(en)28(ter)-280(co)-28(e\016cien)28(ts)-280(for)-281(already)]TJ 0 -11.955 Td [(existing)-333(nonzero)-334(en)28(tries.)]TJ -24.907 -21.124 Td [(The)-358(only)-357(storage)-358(v)56(arian)27(t)-357(supp)-28(orting)-357(the)-358(build)-357(s)-1(tate)-357(is)-358(COO;)-357(all)-358(other)-358(v)56(arian)28(ts)]TJ 0 -11.955 Td [(are)-333(obtained)-334(b)28(y)-333(con)28(v)27(ersion)-333(to/from)-333(it.)]TJ/F27 9.9626 Tf 0 -28.016 Td [(3.2.1)-1150(Sparse)-383(Matri)-1(x)-383(Metho)-32(ds)]TJ 0 -19.133 Td [(3.2.2)-1150(get)]TJ ET q -1 0 0 1 151.025 354.76 cm +1 0 0 1 151.025 274.371 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 154.462 354.561 Td [(nro)32(ws)-383(|)-384(Get)-383(n)32(um)32(b)-32(er)-383(of)-384(ro)32(ws)-383(in)-383(a)-384(sparse)-383(matrix)]TJ +/F27 9.9626 Tf 154.462 274.171 Td [(nro)32(ws)-383(|)-384(Get)-383(n)32(um)32(b)-32(er)-383(of)-384(ro)32(ws)-383(in)-383(a)-384(sparse)-383(matrix)]TJ 0 g 0 G 0 g 0 G -/F30 9.9626 Tf -54.567 -18.51 Td [(nr)-525(=)-525(a%get_nrows\050\051)]TJ +/F30 9.9626 Tf -54.567 -19.133 Td [(nr)-525(=)-525(a%get_nrows\050\051)]TJ 0 g 0 G -/F27 9.9626 Tf 0 -22.111 Td [(T)32(yp)-32(e:)]TJ +/F27 9.9626 Tf 0 -23.116 Td [(T)32(yp)-32(e:)]TJ 0 g 0 G /F8 9.9626 Tf 33.797 0 Td [(Async)28(hronous.)]TJ 0 g 0 G -/F27 9.9626 Tf -33.797 -20.183 Td [(On)-383(En)32(try)]TJ -0 g 0 G -0 g 0 G - 0 -20.184 Td [(a)]TJ -0 g 0 G -/F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.356 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.379 0 Td [(lo)-32(cal)]TJ -0 g 0 G - -57.286 -34.067 Td [(On)-383(Return)]TJ +/F27 9.9626 Tf -33.797 -21.523 Td [(On)-383(En)32(try)]TJ 0 g 0 G 0 g 0 G - 0 -20.183 Td [(F)96(unction)-384(v)64(alue)]TJ + 0 -21.523 Td [(a)]TJ 0 g 0 G -/F8 9.9626 Tf 78.387 0 Td [(The)-333(n)28(um)27(b)-27(e)-1(r)-333(of)-333(ro)28(ws)-334(of)-333(sparse)-333(matrix)]TJ/F30 9.9626 Tf 164.937 0 Td [(a)]TJ/F8 9.9626 Tf 5.23 0 Td [(.)]TJ/F27 9.9626 Tf -248.554 -26.238 Td [(3.2.3)-1150(get)]TJ -ET -q -1 0 0 1 151.025 181.329 cm -[]0 d 0 J 0.398 w 0 0 m 3.437 0 l S -Q -BT -/F27 9.9626 Tf 154.462 181.13 Td [(ncols)-383(|)-384(Get)-383(n)32(um)32(b)-32(er)-383(of)-384(columns)-383(in)-384(a)-383(sparse)-383(matrix)]TJ -0 g 0 G -0 g 0 G -/F30 9.9626 Tf -54.567 -18.51 Td [(nc)-525(=)-525(a%get_ncols\050\051)]TJ +/F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.356 -11.956 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.379 0 Td [(lo)-32(cal)]TJ 0 g 0 G -/F27 9.9626 Tf 0 -22.111 Td [(T)32(yp)-32(e:)]TJ + -57.286 -35.071 Td [(On)-383(Return)]TJ 0 g 0 G -/F8 9.9626 Tf 33.797 0 Td [(Async)28(hronous.)]TJ 0 g 0 G -/F27 9.9626 Tf -33.797 -20.183 Td [(On)-383(En)32(try)]TJ + 0 -21.523 Td [(F)96(unction)-384(v)64(alue)]TJ 0 g 0 G +/F8 9.9626 Tf 78.387 0 Td [(The)-333(n)28(um)27(b)-27(e)-1(r)-333(of)-333(ro)28(ws)-334(of)-333(sparse)-333(matrix)]TJ/F30 9.9626 Tf 164.937 0 Td [(a)]TJ/F8 9.9626 Tf 5.23 0 Td [(.)]TJ 0 g 0 G -/F8 9.9626 Tf 166.875 -29.888 Td [(17)]TJ + -81.679 -29.888 Td [(17)]TJ 0 g 0 G ET @@ -6821,81 +6809,99 @@ endstream endobj 957 0 obj << -/Length 3543 +/Length 4079 >> stream 0 g 0 G 0 g 0 G -0 g 0 G BT -/F27 9.9626 Tf 150.705 706.129 Td [(a)]TJ +/F27 9.9626 Tf 150.705 706.129 Td [(3.2.3)-1150(get)]TJ +ET +q +1 0 0 1 201.835 706.328 cm +[]0 d 0 J 0.398 w 0 0 m 3.437 0 l S +Q +BT +/F27 9.9626 Tf 205.272 706.129 Td [(ncols)-383(|)-384(Get)-383(n)32(um)32(b)-32(er)-383(of)-384(columns)-383(in)-383(a)-384(sparse)-383(matrix)]TJ 0 g 0 G -/F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.355 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.379 0 Td [(lo)-32(cal)]TJ 0 g 0 G - -57.285 -37.362 Td [(On)-383(Return)]TJ +/F30 9.9626 Tf -54.567 -18.389 Td [(nc)-525(=)-525(a%get_ncols\050\051)]TJ +0 g 0 G +/F27 9.9626 Tf 0 -21.464 Td [(T)32(yp)-32(e:)]TJ +0 g 0 G +/F8 9.9626 Tf 33.797 0 Td [(Async)28(hronous.)]TJ +0 g 0 G +/F27 9.9626 Tf -33.797 -19.743 Td [(On)-383(En)32(try)]TJ +0 g 0 G +0 g 0 G + 0 -19.743 Td [(a)]TJ +0 g 0 G +/F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.356 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.378 0 Td [(lo)-32(cal)]TJ +0 g 0 G + -57.285 -33.419 Td [(On)-383(Return)]TJ 0 g 0 G 0 g 0 G - 0 -24.577 Td [(F)96(unction)-384(v)64(alue)]TJ + 0 -19.743 Td [(F)96(unction)-384(v)64(alue)]TJ 0 g 0 G -/F8 9.9626 Tf 78.386 0 Td [(The)-333(n)27(um)28(b)-28(er)-333(of)-333(columns)-334(of)-333(sparse)-333(matrix)]TJ/F30 9.9626 Tf 180.684 0 Td [(a)]TJ/F8 9.9626 Tf 5.23 0 Td [(.)]TJ/F27 9.9626 Tf -264.3 -32.066 Td [(3.2.4)-1150(get)]TJ +/F8 9.9626 Tf 78.386 0 Td [(The)-333(n)27(um)28(b)-28(er)-333(of)-333(columns)-334(of)-333(sparse)-333(matrix)]TJ/F30 9.9626 Tf 180.684 0 Td [(a)]TJ/F8 9.9626 Tf 5.23 0 Td [(.)]TJ/F27 9.9626 Tf -264.3 -25.818 Td [(3.2.4)-1150(get)]TJ ET q -1 0 0 1 201.835 600.368 cm +1 0 0 1 201.835 536.055 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 205.272 600.169 Td [(nnzeros)-497(|)-498(Get)-497(n)32(um)32(b)-32(er)-498(of)-497(nonzero)-497(elemen)32(ts)-498(in)-497(a)-498(sparse)]TJ -19.56 -11.955 Td [(matrix)]TJ +/F27 9.9626 Tf 205.272 535.855 Td [(nnzeros)-497(|)-498(Get)-497(n)32(um)32(b)-32(er)-498(of)-497(nonzero)-497(elemen)32(ts)-498(in)-497(a)-498(sparse)]TJ -19.56 -11.955 Td [(matrix)]TJ 0 g 0 G 0 g 0 G -/F30 9.9626 Tf -35.007 -20.554 Td [(nz)-525(=)-525(a%get_nnzeros\050\051)]TJ +/F30 9.9626 Tf -35.007 -18.389 Td [(nz)-525(=)-525(a%get_nnzeros\050\051)]TJ 0 g 0 G -/F27 9.9626 Tf 0 -25.407 Td [(T)32(yp)-32(e:)]TJ +/F27 9.9626 Tf 0 -21.463 Td [(T)32(yp)-32(e:)]TJ 0 g 0 G /F8 9.9626 Tf 33.797 0 Td [(Async)28(hronous.)]TJ 0 g 0 G -/F27 9.9626 Tf -33.797 -24.577 Td [(On)-383(En)32(try)]TJ +/F27 9.9626 Tf -33.797 -19.744 Td [(On)-383(En)32(try)]TJ 0 g 0 G 0 g 0 G - 0 -24.577 Td [(a)]TJ + 0 -19.743 Td [(a)]TJ 0 g 0 G /F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.356 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.378 0 Td [(lo)-32(cal)]TJ 0 g 0 G - -57.285 -37.362 Td [(On)-383(Return)]TJ + -57.285 -33.418 Td [(On)-383(Return)]TJ 0 g 0 G 0 g 0 G - 0 -24.577 Td [(F)96(unction)-384(v)64(alue)]TJ + 0 -19.744 Td [(F)96(unction)-384(v)64(alue)]TJ 0 g 0 G -/F8 9.9626 Tf 78.386 0 Td [(The)-333(n)27(um)28(b)-28(er)-333(of)-333(nonzero)-333(e)-1(l)1(e)-1(men)28(ts)-333(stored)-334(in)-333(sparse)-333(matrix)]TJ/F30 9.9626 Tf 249.98 0 Td [(a)]TJ/F8 9.9626 Tf 5.23 0 Td [(.)]TJ/F27 9.9626 Tf -333.596 -26.57 Td [(Notes)]TJ +/F8 9.9626 Tf 78.386 0 Td [(The)-333(n)27(um)28(b)-28(er)-333(of)-333(nonzero)-333(e)-1(l)1(e)-1(men)28(ts)-333(stored)-334(in)-333(sparse)-333(matrix)]TJ/F30 9.9626 Tf 249.98 0 Td [(a)]TJ/F8 9.9626 Tf 5.23 0 Td [(.)]TJ/F27 9.9626 Tf -333.596 -21.463 Td [(Notes)]TJ 0 g 0 G -/F8 9.9626 Tf 12.176 -23.414 Td [(1.)]TJ +/F8 9.9626 Tf 12.176 -19.561 Td [(1.)]TJ 0 g 0 G - [-500(The)-462(function)-462(v)55(alue)-462(is)-462(sp)-28(eci\014c)-462(to)-463(th)1(e)-463(storage)-462(format)-462(of)-462(matrix)]TJ/F30 9.9626 Tf 296.649 0 Td [(a)]TJ/F8 9.9626 Tf 5.231 0 Td [(;)-527(some)]TJ -289.149 -11.955 Td [(storage)-465(formats)-466(emplo)28(y)-465(padding,)-498(th)28(us)-466(the)-465(returned)-465(v)55(alue)-465(for)-465(the)-466(same)]TJ 0 -11.955 Td [(matrix)-333(ma)27(y)-333(b)-28(e)-333(di\013eren)28(t)-333(for)-334(di\013eren)28(t)-333(storage)-334(c)28(hoices.)]TJ/F27 9.9626 Tf -24.907 -32.066 Td [(3.2.5)-1150(get)]TJ + [-500(The)-462(function)-462(v)55(alue)-462(is)-462(sp)-28(eci\014c)-462(to)-463(th)1(e)-463(storage)-462(format)-462(of)-462(matrix)]TJ/F30 9.9626 Tf 296.649 0 Td [(a)]TJ/F8 9.9626 Tf 5.231 0 Td [(;)-527(some)]TJ -289.149 -11.955 Td [(storage)-465(formats)-466(emplo)28(y)-465(padding,)-498(th)28(us)-466(the)-465(returned)-465(v)55(alue)-465(for)-465(the)-466(same)]TJ 0 -11.956 Td [(matrix)-333(ma)27(y)-333(b)-28(e)-333(di\013eren)28(t)-333(for)-334(di\013eren)28(t)-333(storage)-334(c)28(hoices.)]TJ/F27 9.9626 Tf -24.907 -25.817 Td [(3.2.5)-1150(get)]TJ ET q -1 0 0 1 201.835 313.444 cm +1 0 0 1 201.835 288.891 cm []0 d 0 J 0.398 w 0 0 m 3.437 0 l S Q BT -/F27 9.9626 Tf 205.272 313.245 Td [(size)-514(|)-513(Get)-514(maxim)32(um)-513(n)31(u)1(m)31(b)-32(er)-513(of)-514(nonzero)-513(elemen)32(ts)-514(in)-513(a)]TJ -19.56 -11.955 Td [(sparse)-383(matri)-1(x)]TJ +/F27 9.9626 Tf 205.272 288.692 Td [(size)-514(|)-513(Get)-514(maxim)32(um)-513(n)31(u)1(m)31(b)-32(er)-513(of)-514(nonzero)-513(elemen)32(ts)-514(in)-513(a)]TJ -19.56 -11.955 Td [(sparse)-383(matri)-1(x)]TJ 0 g 0 G 0 g 0 G -/F30 9.9626 Tf -35.007 -20.554 Td [(maxnz)-525(=)-525(a%get_size\050\051)]TJ +/F30 9.9626 Tf -35.007 -18.39 Td [(maxnz)-525(=)-525(a%get_size\050\051)]TJ 0 g 0 G -/F27 9.9626 Tf 0 -25.407 Td [(T)32(yp)-32(e:)]TJ +/F27 9.9626 Tf 0 -21.463 Td [(T)32(yp)-32(e:)]TJ 0 g 0 G /F8 9.9626 Tf 33.797 0 Td [(Async)28(hronous.)]TJ 0 g 0 G -/F27 9.9626 Tf -33.797 -24.577 Td [(On)-383(En)32(try)]TJ +/F27 9.9626 Tf -33.797 -19.743 Td [(On)-383(En)32(try)]TJ 0 g 0 G 0 g 0 G - 0 -24.577 Td [(a)]TJ + 0 -19.743 Td [(a)]TJ 0 g 0 G -/F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.356 -11.955 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.378 0 Td [(lo)-32(cal)]TJ +/F8 9.9626 Tf 10.551 0 Td [(the)-333(sparse)-334(matrix)]TJ 14.356 -11.956 Td [(Scop)-28(e:)]TJ/F27 9.9626 Tf 32.378 0 Td [(lo)-32(cal)]TJ 0 g 0 G - -57.285 -37.362 Td [(On)-383(Return)]TJ + -57.285 -33.418 Td [(On)-383(Return)]TJ 0 g 0 G 0 g 0 G - 0 -24.577 Td [(F)96(unction)-384(v)64(alue)]TJ + 0 -19.743 Td [(F)96(unction)-384(v)64(alue)]TJ 0 g 0 G /F8 9.9626 Tf 78.386 0 Td [(The)-253(maxim)27(um)-253(n)28(um)28(b)-28(er)-253(of)-253(nonzero)-254(elemen)28(ts)-253(that)-253(can)-254(b)-27(e)-254(stored)]TJ -53.479 -11.955 Td [(in)-333(sparse)-334(matrix)]TJ/F30 9.9626 Tf 74.055 0 Td [(a)]TJ/F8 9.9626 Tf 8.551 0 Td [(using)-333(its)-334(curren)28(t)-333(memory)-334(allo)-27(cation.)]TJ 0 g 0 G @@ -8000,19 +8006,19 @@ endobj /Type /ObjStm /N 100 /First 885 -/Length 8978 +/Length 8979 >> stream 931 0 933 114 71 171 75 227 934 283 79 339 83 395 930 451 937 543 939 657 87 713 91 768 95 823 936 878 944 970 940 1120 941 1276 942 1426 946 1571 99 1628 -103 1684 947 1741 948 1798 943 1855 952 1960 954 2074 950 2130 107 2187 111 2243 115 2299 -951 2355 956 2447 958 2561 119 2618 959 2674 123 2731 955 2788 961 2880 963 2994 127 3050 -131 3106 135 3162 960 3218 966 3310 968 3424 139 3481 143 3538 965 3595 970 3687 972 3801 -147 3857 151 3913 969 3969 974 4061 976 4175 155 4232 159 4289 163 4346 973 4403 978 4495 -980 4609 167 4665 977 4721 984 4813 981 4955 982 5102 986 5248 171 5305 175 5362 179 5419 -183 5476 987 5533 983 5590 990 5695 992 5809 988 5865 187 5922 191 5978 195 6034 989 6090 -997 6182 994 6324 995 6469 999 6614 199 6671 996 6728 1004 6820 1001 6958 1006 7105 204 7163 -208 7220 212 7276 1007 7333 1003 7392 1010 7511 1002 7649 1012 7797 1008 7856 216 7915 1009 7973 +103 1684 947 1741 943 1798 951 1903 953 2017 949 2073 954 2130 107 2187 111 2243 950 2299 +956 2391 958 2505 115 2562 119 2619 959 2676 123 2733 955 2789 961 2881 963 2995 127 3051 +131 3107 135 3163 960 3219 966 3311 968 3425 139 3482 143 3539 965 3596 970 3688 972 3802 +147 3858 151 3914 969 3970 974 4062 976 4176 155 4233 159 4290 163 4347 973 4404 978 4496 +980 4610 167 4666 977 4722 984 4814 981 4956 982 5103 986 5249 171 5306 175 5363 179 5420 +183 5477 987 5534 983 5591 990 5696 992 5810 988 5866 187 5923 191 5979 195 6035 989 6091 +997 6183 994 6325 995 6470 999 6615 199 6672 996 6729 1004 6821 1001 6959 1006 7106 204 7164 +208 7221 212 7277 1007 7334 1003 7393 1010 7512 1002 7650 1012 7798 1008 7857 216 7916 1009 7974 % 931 0 obj << /Type /Page @@ -8093,7 +8099,7 @@ stream /Type /Annot /Subtype /Link /Border[0 0 0]/H/I/C[0 1 0] -/Rect [188.06 277.052 200.015 285.465] +/Rect [188.06 250.478 200.015 258.891] /A << /S /GoTo /D (cite.DesignPatterns) >> >> % 941 0 obj @@ -8101,7 +8107,7 @@ stream /Type /Annot /Subtype /Link /Border[0 0 0]/H/I/C[0 1 0] -/Rect [268.904 277.052 280.86 285.465] +/Rect [268.904 250.478 280.86 258.891] /A << /S /GoTo /D (cite.Sparse03) >> >> % 942 0 obj @@ -8109,7 +8115,7 @@ stream /Type /Annot /Subtype /Link /Border[0 0 0]/H/I/C[1 0 0] -/Rect [459.496 274.562 466.47 286.517] +/Rect [459.496 247.987 466.47 259.942] /A << /S /GoTo /D (figure.4) >> >> % 946 0 obj @@ -8118,54 +8124,50 @@ stream >> % 99 0 obj << -/D [944 0 R /XYZ 150.705 553.709 null] +/D [944 0 R /XYZ 150.705 547.837 null] >> % 103 0 obj << -/D [944 0 R /XYZ 150.705 332.359 null] +/D [944 0 R /XYZ 150.705 306.931 null] >> % 947 0 obj << -/D [944 0 R /XYZ 171.358 301.959 null] ->> -% 948 0 obj -<< -/D [944 0 R /XYZ 455.673 168.146 null] +/D [944 0 R /XYZ 171.358 275.384 null] >> % 943 0 obj << /Font << /F27 584 0 R /F8 585 0 R /F30 811 0 R /F16 582 0 R >> /ProcSet [ /PDF /Text ] >> -% 952 0 obj +% 951 0 obj << /Type /Page -/Contents 953 0 R -/Resources 951 0 R +/Contents 952 0 R +/Resources 950 0 R /MediaBox [0 0 595.276 841.89] /Parent 929 0 R >> -% 954 0 obj +% 953 0 obj << -/D [952 0 R /XYZ 98.895 753.953 null] +/D [951 0 R /XYZ 98.895 753.953 null] >> -% 950 0 obj +% 949 0 obj << -/D [952 0 R /XYZ 155.008 667.994 null] +/D [951 0 R /XYZ 155.008 667.994 null] >> -% 107 0 obj +% 954 0 obj << -/D [952 0 R /XYZ 99.895 382.599 null] +/D [951 0 R /XYZ 404.863 623.415 null] >> -% 111 0 obj +% 107 0 obj << -/D [952 0 R /XYZ 99.895 364.643 null] +/D [951 0 R /XYZ 99.895 303.168 null] >> -% 115 0 obj +% 111 0 obj << -/D [952 0 R /XYZ 99.895 191.212 null] +/D [951 0 R /XYZ 99.895 284.589 null] >> -% 951 0 obj +% 950 0 obj << /Font << /F30 811 0 R /F8 585 0 R /F27 584 0 R >> /ProcSet [ /PDF /Text ] @@ -8182,21 +8184,25 @@ stream << /D [956 0 R /XYZ 149.705 753.953 null] >> +% 115 0 obj +<< +/D [956 0 R /XYZ 150.705 716.092 null] +>> % 119 0 obj << -/D [956 0 R /XYZ 150.705 611.35 null] +/D [956 0 R /XYZ 150.705 545.873 null] >> % 959 0 obj << -/D [956 0 R /XYZ 150.705 387.487 null] +/D [956 0 R /XYZ 150.705 354.269 null] >> % 123 0 obj << -/D [956 0 R /XYZ 150.705 324.426 null] +/D [956 0 R /XYZ 150.705 298.71 null] >> % 955 0 obj << -/Font << /F27 584 0 R /F8 585 0 R /F30 811 0 R >> +/Font << /F27 584 0 R /F30 811 0 R /F8 585 0 R >> /ProcSet [ /PDF /Text ] >> % 961 0 obj @@ -26403,7 +26409,7 @@ stream 1914 3023 1915 3165 1916 3309 1920 3454 1917 3513 1923 3632 1921 3771 1925 3915 1922 3973 1928 4079 1930 4197 1927 4256 1932 4323 1934 4441 1935 4499 781 4557 1936 4614 833 4671 831 4728 787 4785 788 4842 804 4899 784 4956 785 5013 1937 5070 780 5128 1938 5185 1931 5243 1940 5336 1942 5454 -949 5513 819 5571 786 5629 783 5687 779 5745 782 5803 1943 5861 1939 5920 1944 6013 1945 6058 +948 5513 819 5571 786 5629 783 5687 779 5745 782 5803 1943 5861 1939 5920 1944 6013 1945 6058 1946 6197 1947 6384 1948 6878 1949 7207 1950 7550 1951 7679 1952 7700 1953 8206 1954 8251 1955 8941 1956 9269 1957 9350 1958 9725 1959 10362 1960 11021 1961 11649 1962 12293 1963 12578 1965 13226 1967 13656 1969 14083 1971 14332 1973 14574 1975 14908 1977 15125 1979 15364 1981 15586 1983 16122 1985 16359 1987 16607 @@ -26736,7 +26742,7 @@ stream << /D [1940 0 R /XYZ 149.705 753.953 null] >> -% 949 0 obj +% 948 0 obj << /D [1940 0 R /XYZ 150.705 716.092 null] >> @@ -28308,7 +28314,7 @@ stream /Type /Pages /Count 6 /Parent 2005 0 R -/Kids [926 0 R 931 0 R 937 0 R 944 0 R 952 0 R 956 0 R] +/Kids [926 0 R 931 0 R 937 0 R 944 0 R 951 0 R 956 0 R] >> % 964 0 obj << @@ -29835,7 +29841,7 @@ stream >> % 2037 0 obj << -/Names [(cite.BLAS2) 787 0 R (cite.BLAS3) 788 0 R (cite.DesPat:11) 781 0 R (cite.DesignPatterns) 949 0 R (cite.KIVA3PSBLAS) 1938 0 R (cite.METIS) 819 0 R] +/Names [(cite.BLAS2) 787 0 R (cite.BLAS3) 788 0 R (cite.DesPat:11) 781 0 R (cite.DesignPatterns) 948 0 R (cite.KIVA3PSBLAS) 1938 0 R (cite.METIS) 819 0 R] /Limits [(cite.BLAS2) (cite.METIS)] >> % 2038 0 obj @@ -29850,7 +29856,7 @@ stream >> % 2040 0 obj << -/Names [(equation.4.3) 1184 0 R (figure.1) 798 0 R (figure.10) 1746 0 R (figure.2) 827 0 R (figure.3) 918 0 R (figure.4) 950 0 R] +/Names [(equation.4.3) 1184 0 R (figure.1) 798 0 R (figure.10) 1746 0 R (figure.2) 827 0 R (figure.3) 918 0 R (figure.4) 949 0 R] /Limits [(equation.4.3) (figure.4)] >> % 2041 0 obj @@ -29975,7 +29981,7 @@ stream >> % 2065 0 obj << -/Names [(page.153) 1934 0 R (page.154) 1942 0 R (page.16) 946 0 R (page.17) 954 0 R (page.18) 958 0 R (page.19) 963 0 R] +/Names [(page.153) 1934 0 R (page.154) 1942 0 R (page.16) 946 0 R (page.17) 953 0 R (page.18) 958 0 R (page.19) 963 0 R] /Limits [(page.153) (page.19)] >> % 2066 0 obj @@ -30065,7 +30071,7 @@ stream >> % 2083 0 obj << -/Names [(section.6) 308 0 R (section.7) 421 0 R (section.8) 493 0 R (section.9) 513 0 R (spbasedata) 948 0 R (spdata) 947 0 R] +/Names [(section.6) 308 0 R (section.7) 421 0 R (section.8) 493 0 R (section.9) 513 0 R (spbasedata) 954 0 R (spdata) 947 0 R] /Limits [(section.6) (spdata)] >> % 2084 0 obj @@ -30198,11 +30204,11 @@ endstream endobj 2133 0 obj << - /Title (Parallel Sparse BLAS V. 3.7.0) /Subject (Parallel Sparse Basic Linear Algebra Subroutines) /Keywords (Computer Science Linear Algebra Fluid Dynamics Parallel Linux MPI PSBLAS Iterative Solvers Preconditioners) /Creator (pdfLaTeX) /Producer ($Id$) /Author()/Title()/Subject()/Creator(LaTeX with hyperref package)/Producer(pdfTeX-1.40.19)/Keywords() -/CreationDate (D:20200623125345+02'00') -/ModDate (D:20200623125345+02'00') + /Title (Parallel Sparse BLAS V. 3.7.0) /Subject (Parallel Sparse Basic Linear Algebra Subroutines) /Keywords (Computer Science Linear Algebra Fluid Dynamics Parallel Linux MPI PSBLAS Iterative Solvers Preconditioners) /Creator (pdfLaTeX) /Producer ($Id$) /Author()/Title()/Subject()/Creator(LaTeX with hyperref)/Keywords() +/CreationDate (D:20200723115613+02'00') +/ModDate (D:20200723115613+02'00') /Trapped /False -/PTEX.Fullbanner (This is pdfTeX, Version 3.14159265-2.6-1.40.19 (TeX Live 2018) kpathsea version 6.3.0) +/PTEX.Fullbanner (This is pdfTeX, Version 3.14159265-2.6-1.40.20 (TeX Live 2019) kpathsea version 6.3.1) >> endobj 2110 0 obj @@ -30345,13 +30351,13 @@ endobj /W [1 3 1] /Root 2132 0 R /Info 2133 0 R -/ID [ ] +/ID [<21388786B3B1D6C580F616CB7C8CC241> <21388786B3B1D6C580F616CB7C8CC241>] /Length 10675 >> stream ZIS[@  @ c  @,b @3a@8`@9_@Y^@Z]@[\@\[@`Z@aY @bX!"W#$V%&U'(T)* -S+, R-. Q/0P12O34N56M78L9:!K;<#J=>'I?@(HAB)GCD-FEF.EGH2DIJ3CKL7BMN8AOP9@QR=?STC>UVD=WXE<YZF;[\L:]^M9_`N8abT7cRY6Z5[4b32  +S+, R-. Q/0P12O34N56M78 L9:!K;<#J=>'I?@(HAB)GCD-FEF.EGH2DIJ3CKL7BMN8AOP9@QR=?STC>UVD=WXE<YZF;[\L:]^M9_`N8abT7cRY6Z5[4b32   1  0 1/<.M-U,f+f *f)f%(f,' f5&!"fI%#$$%& #'(")*"!+,) -.C/0I12I 34I56I78I19:I>;<IE=>IJ?@IPABI[CDIbEFGH'IJ-KL4MN:OPFQRL STP UVU WX_ YZ! [\!]^! _`!ab!cs!!&!3!7!G  !Kc  !Qb !Wa!^`!b_^] @@ -30359,35 +30365,37 @@ stream     !"#$%&'().,_5*+/0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\`^-]abc     !#$&(-+_* "%')/13456789:;=,W<.02?@ABCDEFGHJ>h_@_7_-_A_?_4_5_>_1_2PQR WUK TLMNO_3XY^VcZ\]_=c@S@_5@@`ab_0_/@@@ @*@@@@ -@z@  @ @@@@@@@@)@'@b@@@@@@@@@ @!@"@#@$@%@&@*@4@(e@+@-@.@/@0@1@2@6@;@5@7@:@=@>@?@F@<@@@A@B@C@D@E@H@I@T@G@J@K@L@M@N@O@P@Q@R@S@V@W@]@U{@X@c@^+@_9b MM d_<$9 "*%L& /+-,40ρ1:56>;<@AH?BGKOI$J!QRUP2S W_]VBX\ac^Tm`k  -E %" (!"#$%&'*+,-./3)a0264589:>75;=G? M@ABCDEF#IJKOH2OLNQRSWPK<TV`XgyYZ[\]^_bcffavfftff8ff f -f fff f$ffpfffffffffff ffff"f#f'f! f$f&f)f*f.f(f+f-f0f1f2f3f<f:f/3f4f6f7f8f9f=f>f?f@fAfBfCfEf;PfD%fGfMfKfFl9fHfJfNfOfPfQfRfSfUfLfTfWfXfZfVfYf\f]f^f_f`fafbf[#fcΌg  -   &&Z,d #I!%&'6.$L(*+,-v7:/iV89012345j=;q<'?@AH>BDEFGJ: -KLQIMNOP]TR*(S^U8\VWXYZ[H|`abI_cIIII -:II I(I -III  I I IIIIII $III 5CIII I OIII"I%I! iAI#I$I-I& I'I(I)I*I+I,)I/I2I. `I0I4I9I3 I5I6I7I8I;I<I@I: ;I=I?IBICIFIA IDIHIKIG IIIMINIQIL IO*IVIR - -NISITIUIXIYI_IW -IZI\I]I^IcI` -$Ia +@z@  @ @@@@@@@@)@'@b@@@@@@@@@ @!@"@#@$@%@&@*@4@(e@+@-@.@/@0@1@2@6@;@5@7@:@=@>@?@F@<@@@A@B@C@D@E@H@I@T@G@J@K@L@M@N@O@P@Q@R@S@V@W@]@U{@X@c@^+@_9d MM d_< +$)"*%T& /+5,40Љ1:56>;<@AH?BGKOI% J!QRUP3S W_]VC +X\ac^Uu`t  +N ." ((!"#$%&'*+,-./3)j0264589:>7>;=G?!V@ABCDEF#IJKOH3XLNQRSWPLETV`XhYZ[\]^_bcffawff}ffAff f +f fff f$ffyfffffffffff ffff"f#f'f! f$f&f)f*f.f(f+f-f0f1f2f3f<f:f/4f4f6f7f8f9f=f>f?f@fAfBfCfEf;QfD%fGfMfKfFmBfHfJfNfOfPfQfRfSfUfLfTfWfXfZfVfYf\f]f^f_f`fafbf[,fcϕp  +  + &'c-m #J!%&'6.$M(*+,-w7:/j_89012345s=;z<'?@AH>BDEFGJ;KLQI +MNOP]TR+1S^U9\VWXYZ[I} `abI_cIIII +;II  I(I +III  I I IIIIII %III 6LIII I PIII"I%I! jJI#I$I-I& I'I(I)I*I+I,)I/I2I. iI0I4I9I3 I5I6I7I8I;I<I@I: DI=I?IBICIFIA IDIHIKIG IIIMINIQIL IO*IVIR + WISITIUIXIYI_IW +IZI\I]I^IcI` +%%Ia  -d*  t  -9   -X+" - !$%(# -&*+/) -,.1250 -378;6 -9@< -=>?,BCDGA EIJMH 7KQN %OSVR 3xT[W MMXYZ]`\ PK^-b!a fc! !!! !!! -!!  ! ! !!! !!!!! _!!!!!"! !! !!.!$!'!# 1!%!-!( 2!)!*!+!,!/!0!1!4!. !2!8!5 &L!6!D!9 =!:!;!<!=!>!?!@!A!B!C!H!E X!F/!N!I Y!J!L!M!T!O o!P!R!S![!U !V!X!Y!Z!_!\ e!]!c!` !a  0         ! | 1'" 1#%&-( F)+,4. [/12 395 tH68>: r;=@AD? B2`abE FGH -IJKLMNOPQRSTUVWXYZ[\]^_c <5 -    W3!) #("0$&'*/)D+-.160Sr2458=7d9;<?D>u@@BC4HEFJKOI LNQRSTWP+UZXY\]^_b[%`__cH_5__?__ _ -_ _b_ __O______________ _"_@_!_$_&_#\_%6_)_'oO_(_9_*o_+_,_._6_8_C_:_;_B_D_E_F_G_H_I_J_K_L_M_N_O_P_Q_R_S_T_U_V_W_X0_Y>_Zf<_[_\Ń_]_^ __)Y_`O_a<_b_cXSe-(\789:;<=>*  -    !"#$%&'()*+,-./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abc>>>>>>>>>> > -> > > >>>>>>>>>ZL +e3  u  +B   +a+" + !$%(# +&*+/) +,.1250 +378;6 +9@< =>?,BCDGA EIJMH @KQN &OSVR 4T[W NVXYZ]`\ QT^-b!a gc! !!! !!! +!!  ! ! !!! !!!!! h!!!!!"! !! !!.!$!'!# :!%!-!( ;!)!*!+!,!/!0!1!4!. !2!8!5 'U!6!D!9 >!:!;!<!=!>!?!@!A!B!C!H!E Y!F/!N!I Z!J!L!M!T!O p!P!R!S![!U $!V!X!Y!Z!_!\ n!]!c!` !a  0         !  1'" 2#%&-( G)+,4. \/12 395 uQ68>: {;=@AD? B2`abE FGH +IJKLMNOPQRSTUVWXYZ[\]^_c E> +    `3!2 #("1$&'*/)E+-.160T{2458=7e9;<?D>vI@BC4HEFJKOILNQRSTWP4UZXY\]^_b[.`__cQ_5__H__ _ +_ _k_ __X______ ________ _"_A_!_$_&_#]_%6_)_'pX_(_9_*p_+_,_._6_8_C_:_;_B_D_E_F_G_H_I_J_K_L_M_N_O_P_Q_R_S_T_U_V_W_X9_Y?_ZgE_[_\ƌ_]_^ +__*b_`X_aE_b _cY\"n.1]789:;<=?3  +    !"#$%&'()*+,-./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abc>>>>>>>>>> > +> > > >>>>>>>>>c4 endstream endobj startxref -1350476 +1350708 %%EOF diff --git a/docs/src/datastruct.tex b/docs/src/datastruct.tex index 72ee814b..613b54ee 100644 --- a/docs/src/datastruct.tex +++ b/docs/src/datastruct.tex @@ -530,7 +530,9 @@ data type and precision variants \item[S] Single precision real; \item[D] Double precision real; \item[C] Single precision complex; -\item[Z] Double precision complex. +\item[Z] Double precision complex; +\item[LS,LD,LC,LZ] Same numeric type as above, but with + \verb|psb_lpk_| integer indices. \end{description} The actual data is contained in the polymorphic component \verb|a%a| of type \hypertarget{spbasedata}{{\tt psb\_T\_base\_sparse\_mat}}; its diff --git a/test/fileread/psb_cf_sample.f90 b/test/fileread/psb_cf_sample.f90 index 8170a446..699c4347 100644 --- a/test/fileread/psb_cf_sample.f90 +++ b/test/fileread/psb_cf_sample.f90 @@ -38,10 +38,11 @@ program psb_cf_sample implicit none ! input parameters - character(len=40) :: kmethd, ptype, mtrx_file, rhs_file,renum + character(len=40) :: kmethd, ptype, mtrx_file, rhs_file ! sparse matrices - type(psb_cspmat_type) :: a, aux_a + type(psb_cspmat_type) :: a + type(psb_lcspmat_type) :: aux_a ! preconditioner data type(psb_cprec_type) :: prec @@ -56,6 +57,7 @@ program psb_cf_sample type(psb_desc_type):: desc_a integer(psb_ipk_) :: ictxt, iam, np + integer(psb_lpk_) :: lnp ! solver paramters integer(psb_ipk_) :: iter, itmax, ierr, itrace, ircode,& @@ -140,7 +142,6 @@ program psb_cf_sample m_problem = aux_a%get_nrows() call psb_bcast(ictxt,m_problem) - call psb_mat_renum(psb_mat_renum_identity_,aux_a,info,perm) ! At this point aux_b may still be unallocated if (size(aux_b,dim=1) == m_problem) then @@ -180,7 +181,9 @@ program psb_cf_sample write(psb_out_unit,'("Partition type: graph vector")') write(psb_out_unit,'(" ")') ! write(psb_err_unit,'("Build type: graph")') - call build_mtpart(aux_a,np) + call aux_a%cscnv(info,type='csr') + lnp = np + call build_mtpart(aux_a,lnp) endif call psb_barrier(ictxt) diff --git a/test/fileread/psb_df_sample.f90 b/test/fileread/psb_df_sample.f90 index 17f821b3..8fbf8bc4 100644 --- a/test/fileread/psb_df_sample.f90 +++ b/test/fileread/psb_df_sample.f90 @@ -38,10 +38,11 @@ program psb_df_sample implicit none ! input parameters - character(len=40) :: kmethd, ptype, mtrx_file, rhs_file,renum + character(len=40) :: kmethd, ptype, mtrx_file, rhs_file ! sparse matrices - type(psb_dspmat_type) :: a, aux_a + type(psb_dspmat_type) :: a + type(psb_ldspmat_type) :: aux_a ! preconditioner data type(psb_dprec_type) :: prec @@ -56,6 +57,7 @@ program psb_df_sample type(psb_desc_type):: desc_a integer(psb_ipk_) :: ictxt, iam, np + integer(psb_lpk_) :: lnp ! solver paramters integer(psb_ipk_) :: iter, itmax, ierr, itrace, ircode,& @@ -140,7 +142,6 @@ program psb_df_sample m_problem = aux_a%get_nrows() call psb_bcast(ictxt,m_problem) - call psb_mat_renum(psb_mat_renum_identity_,aux_a,info,perm) ! At this point aux_b may still be unallocated if (size(aux_b,dim=1) == m_problem) then @@ -180,7 +181,9 @@ program psb_df_sample write(psb_out_unit,'("Partition type: graph vector")') write(psb_out_unit,'(" ")') ! write(psb_err_unit,'("Build type: graph")') - call build_mtpart(aux_a,np) + call aux_a%cscnv(info,type='csr') + lnp = np + call build_mtpart(aux_a,lnp) endif call psb_barrier(ictxt) diff --git a/test/fileread/psb_sf_sample.f90 b/test/fileread/psb_sf_sample.f90 index 05127088..71dd646b 100644 --- a/test/fileread/psb_sf_sample.f90 +++ b/test/fileread/psb_sf_sample.f90 @@ -38,10 +38,11 @@ program psb_sf_sample implicit none ! input parameters - character(len=40) :: kmethd, ptype, mtrx_file, rhs_file,renum + character(len=40) :: kmethd, ptype, mtrx_file, rhs_file ! sparse matrices - type(psb_sspmat_type) :: a, aux_a + type(psb_sspmat_type) :: a + type(psb_lsspmat_type) :: aux_a ! preconditioner data type(psb_sprec_type) :: prec @@ -56,6 +57,7 @@ program psb_sf_sample type(psb_desc_type):: desc_a integer(psb_ipk_) :: ictxt, iam, np + integer(psb_lpk_) :: lnp ! solver paramters integer(psb_ipk_) :: iter, itmax, ierr, itrace, ircode,& @@ -140,7 +142,6 @@ program psb_sf_sample m_problem = aux_a%get_nrows() call psb_bcast(ictxt,m_problem) - call psb_mat_renum(psb_mat_renum_identity_,aux_a,info,perm) ! At this point aux_b may still be unallocated if (size(aux_b,dim=1) == m_problem) then @@ -180,7 +181,9 @@ program psb_sf_sample write(psb_out_unit,'("Partition type: graph vector")') write(psb_out_unit,'(" ")') ! write(psb_err_unit,'("Build type: graph")') - call build_mtpart(aux_a,np) + call aux_a%cscnv(info,type='csr') + lnp = np + call build_mtpart(aux_a,lnp) endif call psb_barrier(ictxt) diff --git a/test/fileread/psb_zf_sample.f90 b/test/fileread/psb_zf_sample.f90 index a2d54e6a..860405ba 100644 --- a/test/fileread/psb_zf_sample.f90 +++ b/test/fileread/psb_zf_sample.f90 @@ -38,10 +38,11 @@ program psb_zf_sample implicit none ! input parameters - character(len=40) :: kmethd, ptype, mtrx_file, rhs_file,renum + character(len=40) :: kmethd, ptype, mtrx_file, rhs_file ! sparse matrices - type(psb_zspmat_type) :: a, aux_a + type(psb_zspmat_type) :: a + type(psb_lzspmat_type) :: aux_a ! preconditioner data type(psb_zprec_type) :: prec @@ -56,6 +57,7 @@ program psb_zf_sample type(psb_desc_type):: desc_a integer(psb_ipk_) :: ictxt, iam, np + integer(psb_lpk_) :: lnp ! solver paramters integer(psb_ipk_) :: iter, itmax, ierr, itrace, ircode,& @@ -140,7 +142,6 @@ program psb_zf_sample m_problem = aux_a%get_nrows() call psb_bcast(ictxt,m_problem) - call psb_mat_renum(psb_mat_renum_identity_,aux_a,info,perm) ! At this point aux_b may still be unallocated if (size(aux_b,dim=1) == m_problem) then @@ -180,7 +181,9 @@ program psb_zf_sample write(psb_out_unit,'("Partition type: graph vector")') write(psb_out_unit,'(" ")') ! write(psb_err_unit,'("Build type: graph")') - call build_mtpart(aux_a,np) + call aux_a%cscnv(info,type='csr') + lnp = np + call build_mtpart(aux_a,lnp) endif call psb_barrier(ictxt) diff --git a/util/Makefile b/util/Makefile index 0f572e9c..9809a34f 100644 --- a/util/Makefile +++ b/util/Makefile @@ -21,7 +21,7 @@ IMPLOBJS= psb_s_hbio_impl.o psb_d_hbio_impl.o \ psb_c_renum_impl.o psb_z_renum_impl.o psi_build_mtpart.o MODOBJS=psb_util_mod.o $(BASEOBJS) -COBJS=metis_int.o psb_amd_order.o +COBJS=psb_metis_int.o psb_amd_order.o OBJS=$(COBJS) $(MODOBJS) $(IMPLOBJS) LOCAL_MODS=$(MODOBJS:.o=$(.mod)) LIBNAME=$(UTILLIBNAME) @@ -38,7 +38,7 @@ $(HERE)/$(LIBNAME): $(OBJS) $(OBJS): $(MODDIR)/$(BASEMODNAME)$(.mod) psb_util_mod.o: $(BASEOBJS) -psb_metispart_mod.o: metis_int.o +psb_metispart_mod.o: psb_metis_int.o psb_mat_dist_mod.o: psb_s_mat_dist_mod.o psb_d_mat_dist_mod.o psb_c_mat_dist_mod.o psb_z_mat_dist_mod.o $(IMPLOBJS): $(BASEOBJS) diff --git a/util/metis_int.c b/util/psb_metis_int.c similarity index 73% rename from util/metis_int.c rename to util/psb_metis_int.c index b9d258b4..0b76c2a6 100644 --- a/util/metis_int.c +++ b/util/psb_metis_int.c @@ -1,13 +1,10 @@ #include #if defined(HAVE_METIS_) -#include "metis.h" +#include "psb_metis_int.h" -/* extern int METIS_PartGraphKway(int *, int *, int *, int *, int *, int *, int *, int *, float *, float, int *, int *, int *); */ - - -int metis_PartGraphKway_C(int *n, int *ixadj, int *iadj, int *ivwg, - int *iajw, int *nparts, float *weights, - int *graphpart) +int metis_PartGraphKway_C(idx_t *n, idx_t *ixadj, idx_t *iadj, idx_t *ivwg, + idx_t *iajw, idx_t *nparts, float *weights, + idx_t *graphpart) { int res = -1; idx_t objval = 0; @@ -40,9 +37,9 @@ int metis_PartGraphKway_C(int *n, int *ixadj, int *iadj, int *ivwg, #else -int metis_PartGraphKway_C(int *n, int *ixadj, int *iadj, int *ivwg, - int *iajw, int *nparts, float *weights, - int *graphpart) +int metis_PartGraphKway_C(idx_t *n, idx_t *ixadj, idx_t *iadj, idx_t *ivwg, + idx_t *iajw, idx_t *nparts, float *weights, + idx_t *graphpart) { return(-1); } diff --git a/util/psb_metispart_mod.F90 b/util/psb_metispart_mod.F90 index 425d8106..fa6581e9 100644 --- a/util/psb_metispart_mod.F90 +++ b/util/psb_metispart_mod.F90 @@ -54,31 +54,32 @@ ! uses information prepared by the previous two subroutines. ! module psb_metispart_mod - use psb_base_mod, only : psb_sspmat_type, psb_cspmat_type,& - & psb_dspmat_type, psb_zspmat_type, psb_err_unit, & - & psb_ipk_, psb_lpk_, psb_mpk_, psb_epk_, psb_spk_,& - & psb_s_csr_sparse_mat, psb_d_csr_sparse_mat, & - & psb_c_csr_sparse_mat, psb_z_csr_sparse_mat + use psb_base_mod, only : psb_ipk_, psb_lpk_, psb_mpk_, psb_epk_, & + & psb_err_unit, psb_spk_,& + & psb_lsspmat_type, psb_lcspmat_type,& + & psb_ldspmat_type, psb_lzspmat_type, & + & psb_ls_csr_sparse_mat, psb_ld_csr_sparse_mat, & + & psb_lc_csr_sparse_mat, psb_lz_csr_sparse_mat public part_graph, build_mtpart, distr_mtpart,& & getv_mtpart, free_part private - integer(psb_ipk_), allocatable, save :: graph_vect(:) + integer(psb_lpk_), allocatable, save :: graph_vect(:) interface build_mtpart - module procedure d_mat_build_mtpart, s_mat_build_mtpart,& - & z_mat_build_mtpart, c_mat_build_mtpart + module procedure ld_mat_build_mtpart, ls_mat_build_mtpart,& + & lz_mat_build_mtpart, lc_mat_build_mtpart end interface - interface - subroutine psi_build_mtpart(n,ja,irp,nparts,vect, weights) - import :: psb_ipk_, psb_spk_ + interface psi_build_mtpart + subroutine psi_l_build_mtpart(n,ja,irp,nparts,vect, weights) + import :: psb_lpk_, psb_spk_ implicit none - integer(psb_ipk_), intent(in) :: n, nparts - integer(psb_ipk_), intent(in) :: ja(:), irp(:) - integer(psb_ipk_), allocatable, intent(inout) :: vect(:) + integer(psb_lpk_), intent(in) :: n, nparts + integer(psb_lpk_), intent(in) :: ja(:), irp(:) + integer(psb_lpk_), allocatable, intent(inout) :: vect(:) real(psb_spk_),optional, intent(in) :: weights(:) - end subroutine psi_build_mtpart + end subroutine psi_l_build_mtpart end interface contains @@ -155,29 +156,30 @@ contains ivg(:) = graph_vect(:) end if end subroutine getv_mtpart - - subroutine d_mat_build_mtpart(a,nparts,weights) + + + subroutine ld_mat_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_dspmat_type), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_ldspmat_type), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_dpk_), optional :: weights(:) select type (aa=>a%a) - type is (psb_d_csr_sparse_mat) - call d_csr_build_mtpart(aa,nparts,weights) + type is (psb_ld_csr_sparse_mat) + call ld_csr_build_mtpart(aa,nparts,weights) class default write(psb_err_unit,*) 'Sorry, right now we only take CSR input!' end select - end subroutine d_mat_build_mtpart + end subroutine ld_mat_build_mtpart - subroutine d_csr_build_mtpart(a,nparts,weights) + subroutine ld_csr_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_d_csr_sparse_mat), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_ld_csr_sparse_mat), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_dpk_), optional :: weights(:) real(psb_spk_), allocatable :: wgh_(:) @@ -193,30 +195,30 @@ contains call psi_build_mtpart(a%get_nrows(),a%ja,a%irp,nparts,graph_vect) end if - end subroutine d_csr_build_mtpart + end subroutine ld_csr_build_mtpart - subroutine z_mat_build_mtpart(a,nparts,weights) + subroutine lz_mat_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_zspmat_type), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_lzspmat_type), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_dpk_), optional :: weights(:) select type (aa=>a%a) - type is (psb_z_csr_sparse_mat) - call z_csr_build_mtpart(aa,nparts,weights) + type is (psb_lz_csr_sparse_mat) + call lz_csr_build_mtpart(aa,nparts,weights) class default write(psb_err_unit,*) 'Sorry, right now we only take CSR input!' end select - end subroutine z_mat_build_mtpart + end subroutine lz_mat_build_mtpart - subroutine z_csr_build_mtpart(a,nparts,weights) + subroutine lz_csr_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_z_csr_sparse_mat), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_lz_csr_sparse_mat), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_dpk_), optional :: weights(:) real(psb_spk_), allocatable :: wgh_(:) @@ -232,66 +234,66 @@ contains call psi_build_mtpart(a%get_nrows(),a%ja,a%irp,nparts,graph_vect) end if - end subroutine z_csr_build_mtpart + end subroutine lz_csr_build_mtpart - subroutine s_mat_build_mtpart(a,nparts,weights) + subroutine ls_mat_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_sspmat_type), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_lsspmat_type), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_spk_), optional :: weights(:) select type (aa=>a%a) - type is (psb_s_csr_sparse_mat) - call s_csr_build_mtpart(aa,nparts,weights) + type is (psb_ls_csr_sparse_mat) + call ls_csr_build_mtpart(aa,nparts,weights) class default write(psb_err_unit,*) 'Sorry, right now we only take CSR input!' end select - end subroutine s_mat_build_mtpart + end subroutine ls_mat_build_mtpart - subroutine c_mat_build_mtpart(a,nparts,weights) + subroutine lc_mat_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_cspmat_type), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_lcspmat_type), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_spk_), optional :: weights(:) select type (aa=>a%a) - type is (psb_c_csr_sparse_mat) - call c_csr_build_mtpart(aa,nparts,weights) + type is (psb_lc_csr_sparse_mat) + call lc_csr_build_mtpart(aa,nparts,weights) class default write(psb_err_unit,*) 'Sorry, right now we only take CSR input!' end select - end subroutine c_mat_build_mtpart + end subroutine lc_mat_build_mtpart - subroutine c_csr_build_mtpart(a,nparts,weights) + subroutine lc_csr_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_c_csr_sparse_mat), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_lc_csr_sparse_mat), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_spk_), optional :: weights(:) call psi_build_mtpart(a%get_nrows(),a%ja,a%irp,nparts,graph_vect,weights) - end subroutine c_csr_build_mtpart + end subroutine lc_csr_build_mtpart - subroutine s_csr_build_mtpart(a,nparts,weights) + subroutine ls_csr_build_mtpart(a,nparts,weights) use psb_base_mod implicit none - type(psb_s_csr_sparse_mat), intent(in) :: a - integer(psb_ipk_) :: nparts + type(psb_ls_csr_sparse_mat), intent(in) :: a + integer(psb_lpk_) :: nparts real(psb_spk_), optional :: weights(:) - + call psi_build_mtpart(a%get_nrows(),a%ja,a%irp,nparts,graph_vect,weights) - end subroutine s_csr_build_mtpart - + end subroutine ls_csr_build_mtpart + ! ! WARNING: called IRET otherwise Intel compiler complains, ! methinks it's a compiler bug, will need to report. diff --git a/util/psi_build_mtpart.F90 b/util/psi_build_mtpart.F90 index b8974185..d01c3d85 100644 --- a/util/psi_build_mtpart.F90 +++ b/util/psi_build_mtpart.F90 @@ -1,31 +1,21 @@ - -subroutine psi_build_mtpart(n,ja,irp,nparts,graph_vect,weights) +subroutine psi_l_build_mtpart(n,ja,irp,nparts,graph_vect,weights) use psb_base_mod use iso_c_binding implicit none - integer(psb_ipk_), intent(in) :: n, nparts - integer(psb_ipk_), intent(in) :: ja(:), irp(:) - integer(psb_ipk_), allocatable, intent(inout) :: graph_vect(:) + integer(psb_lpk_), intent(in) :: n, nparts + integer(psb_lpk_), intent(in) :: ja(:), irp(:) + integer(psb_lpk_), allocatable, intent(inout) :: graph_vect(:) real(psb_spk_),optional, intent(in) :: weights(:) ! local variables - integer(psb_ipk_) :: i,numflag, nedc,wgflag - integer(psb_ipk_) :: iopt(10),idummy(2),jdummy(2), info - integer(psb_ipk_) :: nl,nptl - integer(psb_ipk_), allocatable :: irpl(:),jal(:),gvl(:) + integer(psb_lpk_) :: i,numflag, nedc,wgflag + integer(psb_lpk_) :: iopt(10),idummy(2),jdummy(2) + integer(psb_ipk_) :: info + integer(psb_lpk_) :: nl,nptl + integer(psb_lpk_), allocatable :: irpl(:),jal(:),gvl(:) real(psb_spk_),allocatable :: wgh_(:) -#if defined(HAVE_METIS) && defined(IPK4) +#if defined(HAVE_METIS) && defined(LPK4) && defined(METIS_32) interface - ! subroutine METIS_PartGraphKway(n,ixadj,iadj,ivwg,iajw,& - ! & wgflag,numflag,nparts,weights,iopt,nedc,part) bind(c) - ! use iso_c_binding - ! integer(c_int) :: n,wgflag,numflag,nparts,nedc - ! integer(c_int) :: ixadj(*),iadj(*),ivwg(*),iajw(*),iopt(*),part(*) - ! real(c_float) :: weights(*) - ! !integer(psb_ipk_) :: n,wgflag,numflag,nparts,nedc - ! !integer(psb_ipk_) :: ixadj(*),iadj(*),ivwg(*),iajw(*),iopt(*),part(*) - ! end subroutine METIS_PartGraphKway - function METIS_PartGraphKway(n,ixadj,iadj,ivwg,iajw,& & nparts,weights,part) bind(c,name="metis_PartGraphKway_C") result(res) use iso_c_binding @@ -86,10 +76,76 @@ subroutine psi_build_mtpart(n,ja,irp,nparts,graph_vect,weights) graph_vect(i) = 0 enddo endif + +#elif defined(HAVE_METIS) && defined(LPK8) && defined(METIS_64) + + interface + function METIS_PartGraphKway(n,ixadj,iadj,ivwg,iajw,& + & nparts,weights,part) bind(c,name="metis_PartGraphKway_C") result(res) + use iso_c_binding + integer(c_long_long) :: res + integer(c_long_long) :: n,nparts + integer(c_long_long) :: ixadj(*),iadj(*),ivwg(*),iajw(*),part(*) + real(c_float) :: weights(*) + !integer(psb_ipk_) :: n,wgflag,numflag,nparts,nedc + !integer(psb_ipk_) :: ixadj(*),iadj(*),ivwg(*),iajw(*),iopt(*),part(*) + end function METIS_PartGraphKway + end interface + + call psb_realloc(n,graph_vect,info) + if (info == psb_success_) allocate(gvl(n),wgh_(nparts),stat=info) + + if (info /= psb_success_) then + write(psb_err_unit,*) 'Fatal error in BUILD_MTPART: memory allocation ',& + & ' failure.' + return + endif + if (nparts > 1) then + iopt(1) = 0 + numflag = 1 + wgflag = 0 + +!!$ write(*,*) 'Before allocation',nparts + + irpl=irp + jal = ja + nl = n + nptl = nparts + wgh_ = -1.0 + if(present(weights)) then + if (size(weights) == nptl) then +!!$ write(*,*) 'weights present',weights + ! call METIS_PartGraphKway(n,irp,ja,idummy,jdummy,& + ! & wgflag,numflag,nparts,weights,iopt,nedc,graph_vect) + info = METIS_PartGraphKway(nl,irpl,jal,idummy,jdummy,& + & nptl,weights,gvl) + + else +!!$ write(*,*) 'weights absent',wgh_ + info = METIS_PartGraphKway(nl,irpl,jal,idummy,jdummy,& + & nptl,wgh_,gvl) + end if + else +!!$ write(*,*) 'weights absent',wgh_ + info = METIS_PartGraphKway(nl,irpl,jal,idummy,jdummy,& + & nptl,wgh_,gvl) + endif +!!$ write(*,*) 'after allocation',info + + do i=1, n + graph_vect(i) = gvl(i) - 1 + enddo + else + do i=1, n + graph_vect(i) = 0 + enddo + endif + #else - write(psb_err_unit,*) 'Warning: METIS was not configured at PSBLAS compile time !' + + write(psb_err_unit,*) 'Warning: no suitable METIS interface for LPK indices' #endif return -end subroutine psi_build_mtpart +end subroutine psi_l_build_mtpart