[petsc-dev] [petsc-maint #40976] convert a matrix from CSR format to BCSR format
Satish Balay
balay at mcs.anl.gov
Tue Feb 16 21:17:58 CST 2010
On Tue, 16 Feb 2010, Barry Smith wrote:
>
> --with-blas-lapack=0 doesn't work.
I guess we should add that check - to be consistant? [blaslapack.py
doesn't use package.py - so is not completely consistant to all other
packages]
>
> --with-mpi=0 is actually a lie because there is an MPI package that is used;
Its not really a separate package because its part of PETSc
sources... Its tied into PETSc sources - and need petsconf.h. And it
doesn't have separate makefiles to build it - like other pacakges do.
But the primary user interface it provides is 'sequential PETSc' -
which most users relate to. [download-mpiuni is more confusing than
with-mpi=0. For one there is no download - and then one has to know
implicitly mpiuni - is for mpi substitute for sequential mpi]
Satish
> yes it is limited and doesn't work in parallel, it is like the old
> --download-c-blas-lapack
>
> Barry
>
> On Feb 16, 2010, at 9:04 PM, Satish Balay wrote:
>
> > Added note:
> >
> > we do have --with-blacs=0 --with-hypre=0 etc - so --with-mpi=0 is not
> > one of a kind option that need eliminating.
> >
> > We default to --with-mpi=1 --with-blas-lapack=1 [as there are the
> > basic dependencies for PETSc ], and --with-PACKAGENAME=0 for every
> > other PACKAGENAME.
> >
> > The user has to explicitly enable other packages. [either by
> > --download-PACAKGE or --with-PACKAGE=1 option - if its installed in
> > standard search PATHs]
> >
> > --with-x is the funky one thats automatically enabled - if configure
> > is able to find it.
> >
> > satish
> >
> > On Tue, 16 Feb 2010, Satish Balay wrote:
> >
> > > For one - the interface most users are looking for is
> > > --with-sequential - or --without-mpi - which '-with-mpi=0'
> > > represents.. They don't know are care about mpiuni [so download-mpiuni
> > > is not really better]
> > >
> > > Also mpiuni is not distributed or built separately from petsc
> > > sources/makefiles - and is akin to --with-sieve=1
> > >
> > > Blas/lapack is not optional - hence no --with-blas-lapack=0 support.
> > >
> > > Satish
> > >
> > > On Tue, 16 Feb 2010, Barry Smith wrote:
> > >
> > > >
> > > > --download-c-blas-lapack puts "bad" Blas/LAPACK into its own
> > > > library. There is no --with-blaslapack=0
> > > >
> > > > Shouldn't we eliminate -with-mpi=0 and instead have --download-
> > > > mpiuni
> > > >
> > > > Barry
> > > >
> > > > Removes one more special case, and special cases are BAD! So
> > > > obviously it is a good idea, let's see if anyone has a legitimate
> > > > reason for rejecting it instead of knee-jerk conservatism.
> > > >
> > > >
> > > >
> > > > On Feb 16, 2010, at 8:31 PM, Satish Balay wrote:
> > > >
> > > > > This is a conflict with --with-mpi=0 and the curent default
> > > > > --with-single-library=1
> > > > >
> > > > > Ideally we should create libmpiuni.a [and a .so] even with
> > > > > --with-single-library=1 - but I'm not sure how to do that cleanly. So
> > > > > my current fix is to remove libmpiuni.a for
> > > > > --with-single-library=1. [i.e libpetsc.a will have the mpiuni stuff as
> > > > > well]
> > > > >
> > > > > Testing this change - and will push once I verify.
> > > > >
> > > > > Satish
> > > > >
> > > > > On Tue, 16 Feb 2010, Matthew Knepley wrote:
> > > > >
> > > > > > On Tue, Feb 16, 2010 at 8:16 PM, Dahai Guo <dhguo at ncsa.uiuc.edu>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Thanks, Barry.
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > I tried to configure and built petsc-dev on AIX with the following
> > > > > > > options,
> > > > > > > but failed when "make test" as follow "ld: 0706-006 Cannot find or
> > > > > > > open
> > > > > > > library file: -l mpiuni". I used the same options for petsc-3.0,
> > > > > > > everything
> > > > > > > is fine. Would you tell me how to fix it?
> > > > > > >
> > > > > >
> > > > > > Always send configure.log
> > > > > >
> > > > > > Matt
> > > > > >
> > > > > >
> > > > > > > Dahai
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > export OBJECT_MODE=64 # for 64/32-bit ar
> > > > > > > export PETSC_DIR=$PWD
> > > > > > >
> > > > > > > config/configure.py --with-mpi=0 --with-64-bit-pointers=1
> > > > > > > --with-
> > > > > > > cc=xlc_r
> > > > > > > --with-cxx=xlC_r --with-fc=xlf_r COPTFLAGS='-q64 -O3 -qstrict -
> > > > > > > qmaxmem=-1
> > > > > > > -qarch=pwr6 -qtune=pwr6' FOPTFLAGS='-q64 -O3 -qstrict -qmaxmem=-1
> > > > > > > -qarch=pwr6 -qtune=pwr6' CXXOPTFLAGS='-q64 -O3 -qstrict
> > > > > > > -qmaxmem=-1
> > > > > > > -qarch=pwr6 -qtune=pwr6'
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > bash-3.1$ make test
> > > > > > > Running test examples to verify correct installation
> > > > > > > --------------Error detected during compile or
> > > > > > > link!-----------------------
> > > > > > > See
> > > > > > > http://www.mcs.anl.gov/petsc/petsc-2/documentation/troubleshooting.html
> > > > > > > xlc_r -o ex19.o -c -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6
> > > > > > > -I/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/include
> > > > > > > -I/blhome/dhguo/petsc-dev/include -I/blhome/dhguo/petsc-dev/
> > > > > > > include/mpiuni
> > > > > > > -D__SDIR__="src/snes/examples/tutorials/" ex19.c
> > > > > > > xlc_r -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6 -o
> > > > > > > ex19 ex19.o -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lpetsc -lX11 -
> > > > > > > lessl
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lmpiuni -lm -ldl
> > > > > > > -L/usr/vac/lib -lxlopt -lxl -lpthreads -lxlf90 -L/usr/lpp/xlf/lib
> > > > > > > -
> > > > > > > lxlf
> > > > > > > -lxlomp_ser -lm -ldl -lxlopt -lxl -lpthreads -ldl
> > > > > > > ld: 0706-006 Cannot find or open library file: -l mpiuni
> > > > > > > ld:open(): No such file or directory
> > > > > > > make: The error code from the last command is 255.
> > > > > > > make: Ignored error code 255 from last command.
> > > > > > > /usr/linux/bin/rm -f ex19.o
> > > > > > > --------------Error detected during compile or
> > > > > > > link!-----------------------
> > > > > > > See
> > > > > > > http://www.mcs.anl.gov/petsc/petsc-2/documentation/troubleshooting.html
> > > > > > > xlf_r -c -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6
> > > > > > > -I/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/include
> > > > > > > -I/blhome/dhguo/petsc-dev/include -I/blhome/dhguo/petsc-dev/
> > > > > > > include/mpiuni
> > > > > > > -I/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/include
> > > > > > > -I/blhome/dhguo/petsc-dev/include/mpiuni -o ex5f.o ex5f.F
> > > > > > > ** main === End of Compilation 1 ===
> > > > > > > ** forminitialguess === End of Compilation 2 ===
> > > > > > > ** initialguesslocal === End of Compilation 3 ===
> > > > > > > ** formfunctionlocal === End of Compilation 4 ===
> > > > > > > ** formjacobianlocal === End of Compilation 5 ===
> > > > > > > 1501-510 Compilation successful for file ex5f.F.
> > > > > > > xlf_r -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6 -o
> > > > > > > ex5f ex5f.o -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lpetsc -lX11 -
> > > > > > > lessl
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lmpiuni -lm -ldl
> > > > > > > -L/usr/vac/lib -lxlopt -lxl -lpthreads -lxlf90 -L/usr/lpp/xlf/lib
> > > > > > > -
> > > > > > > lxlf
> > > > > > > -lxlomp_ser -lm -ldl -lxlopt -lxl -lpthreads -ldl
> > > > > > > ld: 0706-006 Cannot find or open library file: -l mpiuni
> > > > > > > ld:open(): No such file or directory
> > > > > > > make: The error code from the last command is 255.
> > > > > > > make: Ignored error code 255 from last command.
> > > > > > > /usr/linux/bin/rm -f ex5f.o
> > > > > > > Completed test examples
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > ----- Original Message -----
> > > > > > > From: "Barry Smith" <bsmith at mcs.anl.gov>
> > > > > > > To: petsc-maint at mcs.anl.gov, "Dahai Guo" <dhguo at ncsa.uiuc.edu>
> > > > > > > Sent: Tuesday, February 16, 2010 1:40:53 PM GMT -06:00 Central
> > > > > > > America
> > > > > > > Subject: Re: [petsc-maint #40976] convert a matrix from CSR format
> > > > > > > to BCSR
> > > > > > > format
> > > > > > >
> > > > > > >
> > > > > > > Sorry, I forgot to check the code last weekend.
> > > > > > >
> > > > > > > I cannot reproduce your problem.
> > > > > > >
> > > > > > > Looks like you are using petsc 3.0.0. That is completely insane!
> > > > > > > This problem has been fixed for a year. You should be using petsc-
> > > > > > > dev
> > > > > > > http://www.mcs.anl.gov/petsc/petsc-as/developers/index.html
> > > > > > > it will make your live much easier.
> > > > > > >
> > > > > > >
> > > > > > > Barry
> > > > > > >
> > > > > > > On Feb 15, 2010, at 10:38 AM, Dahai Guo wrote:
> > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Barry:
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > In the function “MatGetFactor” in src/mat/interface/matrix.c, it
> > > > > > > > only checks if the matsolver is available for subclass, but NOT
> > > > > > > > further check for the AIJ/BAIJ. I think it is one of the
> > > > > > > > reasons
> > > > > > > > that a subclass(e.g. crl ) inherited from AIJ/BAIJ can not
> > > > > > > > inherit
> > > > > > > > "MatGetFactor" from AIJ/BAIJ.
> > > > > > > >
> > > > > > > > Is there some special reason for designing it in this way? Or it
> > > > > > > > is
> > > > > > > > just a bug? How to improve it? Please advise.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > >
> > > > > > > > Dahai
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > ----- Original Message -----
> > > > > > > > From: "Barry Smith" <bsmith at mcs.anl.gov>
> > > > > > > > To: petsc-maint at mcs.anl.gov, "Dahai Guo" <dhguo at ncsa.uiuc.edu>
> > > > > > > > Sent: Thursday, February 11, 2010 9:45:24 PM GMT -06:00
> > > > > > > > US/Canada
> > > > > > > > Central
> > > > > > > > Subject: Re: [petsc-maint #40976] convert a matrix from CSR
> > > > > > > > format
> > > > > > > > to BCSR format
> > > > > > > >
> > > > > > > >
> > > > > > > > Dahai,
> > > > > > > >
> > > > > > > > You have to switch to a decent email client. Stuff you send is
> > > > > > > > never in attachments, it is always crammed into the message as
> > > > > > > > one
> > > > > > > > long piece of text.
> > > > > > > >
> > > > > > > > I'll try to reproduce the problem tomorrow,
> > > > > > > >
> > > > > > > > Barry
> > > > > > > >
> > > > > > > > On Feb 11, 2010, at 8:01 PM, Dahai Guo wrote:
> > > > > > > >
> > > > > > > > > tried CSR format, it showed the similar problem as follow.???
> > > > > > > > > I
> > > > > > > > > attached the test code and datafile.
> > > > > > > >
> > > > > > > >
> > > > > > > > Barry:
> > > > > > > >
> > > > > > > > In the function “MatGetFactor” in src/mat/interface/matrix.c, it
> > > > > > > > only checks if the matsolver is available for subclass, but NOT
> > > > > > > > further check for the AIJ/BAIJ. I think it is one of the
> > > > > > > > reasons
> > > > > > > > that a subclass(e.g. crl ) inherited from AIJ/BAIJ can not
> > > > > > > > inherit
> > > > > > > > "MatGetFactor" from AIJ/BAIJ.
> > > > > > > >
> > > > > > > > Is there some special reason for designing it in this way? Or it
> > > > > > > > is
> > > > > > > > just a bug? How to improve it? Please advise.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > >
> > > > > > > > Dahai
> > > > > > > >
> > > > > > > >
> > > > > > > > ----- Original Message -----
> > > > > > > > From: "Barry Smith" <bsmith at mcs.anl.gov>
> > > > > > > > To: petsc-maint at mcs.anl.gov, "Dahai Guo" <dhguo at ncsa.uiuc.edu>
> > > > > > > > Sent: Thursday, February 11, 2010 9:45:24 PM GMT -06:00
> > > > > > > > US/Canada
> > > > > > > > Central
> > > > > > > > Subject: Re: [petsc-maint #40976] convert a matrix from CSR
> > > > > > > > format
> > > > > > > > to BCSR format
> > > > > > > >
> > > > > > > >
> > > > > > > > Dahai,
> > > > > > > >
> > > > > > > > You have to switch to a decent email client. Stuff you send is
> > > > > > > > never in attachments, it is always crammed into the message as
> > > > > > > > one
> > > > > > > > long piece of text.
> > > > > > > >
> > > > > > > > I'll try to reproduce the problem tomorrow,
> > > > > > > >
> > > > > > > > Barry
> > > > > > > >
> > > > > > > > On Feb 11, 2010, at 8:01 PM, Dahai Guo wrote:
> > > > > > > >
> > > > > > > > > tried CSR format, it showed the similar problem as follow.???
> > > > > > > > > I
> > > > > > > > > attached the test code and datafile.
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Thanks, Barry.
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > I tried to configure and built petsc-dev on AIX with the following
> > > > > > > options,
> > > > > > > but failed when "make test" as follow "ld: 0706-006 Cannot find or
> > > > > > > open
> > > > > > > library file: -l mpiuni". I used the same options for petsc-3.0,
> > > > > > > everything
> > > > > > > is fine. Would you tell me how to fix it?
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Dahai
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > export OBJECT_MODE=64 # for 64/32-bit ar
> > > > > > > export PETSC_DIR=$PWD
> > > > > > >
> > > > > > > config/configure.py --with-mpi=0 --with-64-bit-pointers=1
> > > > > > > --with-
> > > > > > > cc=xlc_r
> > > > > > > --with-cxx=xlC_r --with-fc=xlf_r COPTFLAGS='-q64 -O3 -qstrict -
> > > > > > > qmaxmem=-1
> > > > > > > -qarch=pwr6 -qtune=pwr6' FOPTFLAGS='-q64 -O3 -qstrict -qmaxmem=-1
> > > > > > > -qarch=pwr6 -qtune=pwr6' CXXOPTFLAGS='-q64 -O3 -qstrict
> > > > > > > -qmaxmem=-1
> > > > > > > -qarch=pwr6 -qtune=pwr6'
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > bash-3.1$ make test
> > > > > > > Running test examples to verify correct installation
> > > > > > > --------------Error detected during compile or
> > > > > > > link!-----------------------
> > > > > > > See
> > > > > > > http://www.mcs.anl.gov/petsc/petsc-2/documentation/troubleshooting.html
> > > > > > > xlc_r -o ex19.o -c -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6
> > > > > > > -I/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/include
> > > > > > > -I/blhome/dhguo/petsc-dev/include -I/blhome/dhguo/petsc-dev/
> > > > > > > include/mpiuni
> > > > > > > -D__SDIR__="src/snes/examples/tutorials/" ex19.c
> > > > > > > xlc_r -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6 -o
> > > > > > > ex19 ex19.o -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lpetsc -lX11 -
> > > > > > > lessl
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lmpiuni -lm -ldl
> > > > > > > -L/usr/vac/lib -lxlopt -lxl -lpthreads -lxlf90 -L/usr/lpp/xlf/lib
> > > > > > > -
> > > > > > > lxlf
> > > > > > > -lxlomp_ser -lm -ldl -lxlopt -lxl -lpthreads -ldl
> > > > > > > ld: 0706-006 Cannot find or open library file: -l mpiuni
> > > > > > > ld:open(): No such file or directory
> > > > > > > make: The error code from the last command is 255.
> > > > > > > make: Ignored error code 255 from last command.
> > > > > > > /usr/linux/bin/rm -f ex19.o
> > > > > > > --------------Error detected during compile or
> > > > > > > link!-----------------------
> > > > > > > See
> > > > > > > http://www.mcs.anl.gov/petsc/petsc-2/documentation/troubleshooting.html
> > > > > > > xlf_r -c -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6
> > > > > > > -I/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/include
> > > > > > > -I/blhome/dhguo/petsc-dev/include -I/blhome/dhguo/petsc-dev/
> > > > > > > include/mpiuni
> > > > > > > -I/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/include
> > > > > > > -I/blhome/dhguo/petsc-dev/include/mpiuni -o ex5f.o ex5f.F
> > > > > > > ** main === End of Compilation 1 ===
> > > > > > > ** forminitialguess === End of Compilation 2 ===
> > > > > > > ** initialguesslocal === End of Compilation 3 ===
> > > > > > > ** formfunctionlocal === End of Compilation 4 ===
> > > > > > > ** formjacobianlocal === End of Compilation 5 ===
> > > > > > > 1501-510 Compilation successful for file ex5f.F.
> > > > > > > xlf_r -q64 -O3 -qstrict -qmaxmem=-1 -qarch=pwr6 -
> > > > > > > qtune=pwr6 -o
> > > > > > > ex5f ex5f.o -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lpetsc -lX11 -
> > > > > > > lessl
> > > > > > > -L/blhome/dhguo/petsc-dev/aix5.3.0.0-c-debug/lib -lmpiuni -lm -ldl
> > > > > > > -L/usr/vac/lib -lxlopt -lxl -lpthreads -lxlf90 -L/usr/lpp/xlf/lib
> > > > > > > -
> > > > > > > lxlf
> > > > > > > -lxlomp_ser -lm -ldl -lxlopt -lxl -lpthreads -ldl
> > > > > > > ld: 0706-006 Cannot find or open library file: -l mpiuni
> > > > > > > ld:open(): No such file or directory
> > > > > > > make: The error code from the last command is 255.
> > > > > > > make: Ignored error code 255 from last command.
> > > > > > > /usr/linux/bin/rm -f ex5f.o
> > > > > > > Completed test examples
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > ----- Original Message -----
> > > > > > > From: "Barry Smith" <bsmith at mcs.anl.gov>
> > > > > > > To: petsc-maint at mcs.anl.gov, "Dahai Guo" <dhguo at ncsa.uiuc.edu>
> > > > > > > Sent: Tuesday, February 16, 2010 1:40:53 PM GMT -06:00 Central
> > > > > > > America
> > > > > > > Subject: Re: [petsc-maint #40976] convert a matrix from CSR format
> > > > > > > to BCSR
> > > > > > > format
> > > > > > >
> > > > > > >
> > > > > > > Sorry, I forgot to check the code last weekend.
> > > > > > >
> > > > > > > I cannot reproduce your problem.
> > > > > > >
> > > > > > > Looks like you are using petsc 3.0.0. That is completely insane!
> > > > > > > This problem has been fixed for a year. You should be using petsc-
> > > > > > > dev
> > > > > > > http://www.mcs.anl.gov/petsc/petsc-as/developers/index.html
> > > > > > > it will make your live much easier.
> > > > > > >
> > > > > > >
> > > > > > > Barry
> > > > > > >
> > > > > > > On Feb 15, 2010, at 10:38 AM, Dahai Guo wrote:
> > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Barry:
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > In the function “MatGetFactor” in src/mat/interface/matrix.c, it
> > > > > > > > only checks if the matsolver is available for subclass, but NOT
> > > > > > > > further check for the AIJ/BAIJ. I think it is one of the
> > > > > > > > reasons
> > > > > > > > that a subclass(e.g. crl ) inherited from AIJ/BAIJ can not
> > > > > > > > inherit
> > > > > > > > "MatGetFactor" from AIJ/BAIJ.
> > > > > > > >
> > > > > > > > Is there some special reason for designing it in this way? Or it
> > > > > > > > is
> > > > > > > > just a bug? How to improve it? Please advise.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > >
> > > > > > > > Dahai
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > ----- Original Message -----
> > > > > > > > From: "Barry Smith" <bsmith at mcs.anl.gov>
> > > > > > > > To: petsc-maint at mcs.anl.gov, "Dahai Guo" <dhguo at ncsa.uiuc.edu>
> > > > > > > > Sent: Thursday, February 11, 2010 9:45:24 PM GMT -06:00
> > > > > > > > US/Canada
> > > > > > > > Central
> > > > > > > > Subject: Re: [petsc-maint #40976] convert a matrix from CSR
> > > > > > > > format
> > > > > > > > to BCSR format
> > > > > > > >
> > > > > > > >
> > > > > > > > Dahai,
> > > > > > > >
> > > > > > > > You have to switch to a decent email client. Stuff you send is
> > > > > > > > never in attachments, it is always crammed into the message as
> > > > > > > > one
> > > > > > > > long piece of text.
> > > > > > > >
> > > > > > > > I'll try to reproduce the problem tomorrow,
> > > > > > > >
> > > > > > > > Barry
> > > > > > > >
> > > > > > > > On Feb 11, 2010, at 8:01 PM, Dahai Guo wrote:
> > > > > > > >
> > > > > > > > > tried CSR format, it showed the similar problem as follow.???
> > > > > > > > > I
> > > > > > > > > attached the test code and datafile.
> > > > > > > >
> > > > > > > >
> > > > > > > > Barry:
> > > > > > > >
> > > > > > > > In the function “MatGetFactor” in src/mat/interface/matrix.c, it
> > > > > > > > only checks if the matsolver is available for subclass, but NOT
> > > > > > > > further check for the AIJ/BAIJ. I think it is one of the
> > > > > > > > reasons
> > > > > > > > that a subclass(e.g. crl ) inherited from AIJ/BAIJ can not
> > > > > > > > inherit
> > > > > > > > "MatGetFactor" from AIJ/BAIJ.
> > > > > > > >
> > > > > > > > Is there some special reason for designing it in this way? Or it
> > > > > > > > is
> > > > > > > > just a bug? How to improve it? Please advise.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > >
> > > > > > > > Dahai
> > > > > > > >
> > > > > > > >
> > > > > > > > ----- Original Message -----
> > > > > > > > From: "Barry Smith" <bsmith at mcs.anl.gov>
> > > > > > > > To: petsc-maint at mcs.anl.gov, "Dahai Guo" <dhguo at ncsa.uiuc.edu>
> > > > > > > > Sent: Thursday, February 11, 2010 9:45:24 PM GMT -06:00
> > > > > > > > US/Canada
> > > > > > > > Central
> > > > > > > > Subject: Re: [petsc-maint #40976] convert a matrix from CSR
> > > > > > > > format
> > > > > > > > to BCSR format
> > > > > > > >
> > > > > > > >
> > > > > > > > Dahai,
> > > > > > > >
> > > > > > > > You have to switch to a decent email client. Stuff you send is
> > > > > > > > never in attachments, it is always crammed into the message as
> > > > > > > > one
> > > > > > > > long piece of text.
> > > > > > > >
> > > > > > > > I'll try to reproduce the problem tomorrow,
> > > > > > > >
> > > > > > > > Barry
> > > > > > > >
> > > > > > > > On Feb 11, 2010, at 8:01 PM, Dahai Guo wrote:
> > > > > > > >
> > > > > > > > > tried CSR format, it showed the similar problem as follow.???
> > > > > > > > > I
> > > > > > > > > attached the test code and datafile.
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > >
> > > >
> > > >
> > >
> > >
>
More information about the petsc-dev
mailing list