April 2012 - petsc-dev - mailman3.cels.anl.gov

Re: [petsc-dev] make reconfigure
by Satish Balay 10 Apr '12

10 Apr '12

On Tue, 10 Apr 2012, Lisandro Dalcin wrote: > On 10 April 2012 21:12, Satish Balay <balay(a)mcs.anl.gov> wrote: > > Lisandro has commit access - so can push directly. [Lisandro - can you > > commit/push?] > > > > Yes, but I'm visiting KAUST, and while being at home, sometimes my ssh > connections randomly stop working for a random period of time. I can push your diffs - but perhaps you might want to test if ssh on port 443 [https port] works arround the difficulties you are having? hg pull -u ssh://[email protected]:443//hg/petsc/petsc-dev Satish

1 0

Re: [petsc-dev] make reconfigure
by Lisandro Dalcin 10 Apr '12

10 Apr '12

On 10 April 2012 21:12, Satish Balay <balay(a)mcs.anl.gov> wrote: > Lisandro has commit access - so can push directly. [Lisandro - can you > commit/push?] > Yes, but I'm visiting KAUST, and while being at home, sometimes my ssh connections randomly stop working for a random period of time. -- Lisandro Dalcin --------------- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa Fe Colectora RN 168 Km 472, Paraje El Pozo 3000 Santa Fe, Argentina Tel: +54-342-4511594 (ext 1011) Tel/Fax: +54-342-4511169

1 0

Re: [petsc-dev] make reconfigure
by Satish Balay 10 Apr '12

10 Apr '12

Lisandro has commit access - so can push directly. [Lisandro - can you commit/push?] Satish On Tue, 10 Apr 2012, Barry Smith wrote: > > Satish, > > Please add this rule. > > Barry > > > On Apr 10, 2012, at 12:53 PM, Lisandro Dalcin wrote: > > > Can I add the following target to the top level makefile:? > > > > reconfigure: > > @${PYTHON} ${PETSC_ARCH}/conf/reconfigure-${PETSC_ARCH}.py > > > > -- > > Lisandro Dalcin > > --------------- > > CIMEC (INTEC/CONICET-UNL) > > Predio CONICET-Santa Fe > > Colectora RN 168 Km 472, Paraje El Pozo > > 3000 Santa Fe, Argentina > > Tel: +54-342-4511594 (ext 1011) > > Tel/Fax: +54-342-4511169 > >

1 0

Fwd: A few remaining issues with Mat block sizes
by Barry Smith 10 Apr '12

10 Apr '12

Satish, Please add this patch. Thanks Barry Begin forwarded message: > From: Lisandro Dalcin <dalcinl(a)gmail.com> > Date: April 10, 2012 12:43:39 PM CDT > To: For users of the development version of PETSc <petsc-dev(a)mcs.anl.gov> > Subject: Re: [petsc-dev] A few remaining issues with Mat block sizes > Reply-To: For users of the development version of PETSc <petsc-dev(a)mcs.anl.gov> > > On 10 April 2012 17:58, Barry Smith <bsmith(a)mcs.anl.gov> wrote: >> >> Lisandro, >> >> Thanks, I have incorporated these and a few other blocksize fixes. >> > > Many thanks, Barry! However, I still think you are missing this one (I > found the issue after running petsc4py testsuite). Pasted below and > also attached. Please review and push. > > > diff -r 2d1c2b03ef3a src/mat/impls/aij/seq/aij.c > --- a/src/mat/impls/aij/seq/aij.c Tue Apr 10 09:57:31 2012 -0500 > +++ b/src/mat/impls/aij/seq/aij.c Tue Apr 10 20:36:48 2012 +0300 > @@ -3477,9 +3477,12 @@ > PetscErrorCode ierr; > > PetscFunctionBegin; > + if (Ii[0]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE, > "Ii[0] must be 0 it is %D", Ii[0]); > + > + ierr = PetscLayoutSetUp(B->rmap);CHKERRQ(ierr); > + ierr = PetscLayoutSetUp(B->cmap);CHKERRQ(ierr); > + > ierr = MatGetSize(B, &m, &n);CHKERRQ(ierr); > - > - if (Ii[0]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE, > "Ii[0] must be 0 it is %D", Ii[0]); > ierr = PetscMalloc((m+1) * sizeof(PetscInt), &nnz);CHKERRQ(ierr); > for(i = 0; i < m; i++) { > nz = Ii[i+1]- Ii[i]; > > > -- > Lisandro Dalcin > --------------- > CIMEC (INTEC/CONICET-UNL) > Predio CONICET-Santa Fe > Colectora RN 168 Km 472, Paraje El Pozo > 3000 Santa Fe, Argentina > Tel: +54-342-4511594 (ext 1011) > Tel/Fax: +54-342-4511169 [see attached file: matblocksize-seqaij-csr.diff]

1 0

Re: [petsc-dev] make reconfigure
by Barry Smith 10 Apr '12

10 Apr '12

Satish, Please add this rule. Barry On Apr 10, 2012, at 12:53 PM, Lisandro Dalcin wrote: > Can I add the following target to the top level makefile:? > > reconfigure: > @${PYTHON} ${PETSC_ARCH}/conf/reconfigure-${PETSC_ARCH}.py > > -- > Lisandro Dalcin > --------------- > CIMEC (INTEC/CONICET-UNL) > Predio CONICET-Santa Fe > Colectora RN 168 Km 472, Paraje El Pozo > 3000 Santa Fe, Argentina > Tel: +54-342-4511594 (ext 1011) > Tel/Fax: +54-342-4511169

1 0

make reconfigure
by Lisandro Dalcin 10 Apr '12

10 Apr '12

Can I add the following target to the top level makefile:? reconfigure: @${PYTHON} ${PETSC_ARCH}/conf/reconfigure-${PETSC_ARCH}.py -- Lisandro Dalcin --------------- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa Fe Colectora RN 168 Km 472, Paraje El Pozo 3000 Santa Fe, Argentina Tel: +54-342-4511594 (ext 1011) Tel/Fax: +54-342-4511169

1 0

Re: [petsc-dev] A few remaining issues with Mat block sizes
by Lisandro Dalcin 10 Apr '12

10 Apr '12

On 10 April 2012 17:58, Barry Smith <bsmith(a)mcs.anl.gov> wrote: > > Lisandro, > > Thanks, I have incorporated these and a few other blocksize fixes. > Many thanks, Barry! However, I still think you are missing this one (I found the issue after running petsc4py testsuite). Pasted below and also attached. Please review and push. diff -r 2d1c2b03ef3a src/mat/impls/aij/seq/aij.c --- a/src/mat/impls/aij/seq/aij.c Tue Apr 10 09:57:31 2012 -0500 +++ b/src/mat/impls/aij/seq/aij.c Tue Apr 10 20:36:48 2012 +0300 @@ -3477,9 +3477,12 @@ PetscErrorCode ierr; PetscFunctionBegin; + if (Ii[0]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE, "Ii[0] must be 0 it is %D", Ii[0]); + + ierr = PetscLayoutSetUp(B->rmap);CHKERRQ(ierr); + ierr = PetscLayoutSetUp(B->cmap);CHKERRQ(ierr); + ierr = MatGetSize(B, &m, &n);CHKERRQ(ierr); - - if (Ii[0]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE, "Ii[0] must be 0 it is %D", Ii[0]); ierr = PetscMalloc((m+1) * sizeof(PetscInt), &nnz);CHKERRQ(ierr); for(i = 0; i < m; i++) { nz = Ii[i+1]- Ii[i]; -- Lisandro Dalcin --------------- CIMEC (INTEC/CONICET-UNL) Predio CONICET-Santa Fe Colectora RN 168 Km 472, Paraje El Pozo 3000 Santa Fe, Argentina Tel: +54-342-4511594 (ext 1011) Tel/Fax: +54-342-4511169

1 0

Re: [petsc-dev] A few remaining issues with Mat block sizes
by Barry Smith 10 Apr '12

10 Apr '12

Lisandro, Thanks, I have incorporated these and a few other blocksize fixes. Barry On Apr 10, 2012, at 4:36 AM, Lisandro Dalcin wrote: > Barry, could you review the attached patch? Basically, many Mat > preallocation routines should not set the block size to 1, and instead > honor what the user could have set calling MatSetBlockSize(). Please > note that in this patch I'm not addressing all instances of this issue > in the codebase, I've just worked on some of the matrix types > (seq/mpidense, seq/mpiaij, shell) that I know for sure the fix seems > to make sense. > > > > -- > Lisandro Dalcin > --------------- > CIMEC (INTEC/CONICET-UNL) > Predio CONICET-Santa Fe > Colectora RN 168 Km 472, Paraje El Pozo > 3000 Santa Fe, Argentina > Tel: +54-342-4511594 (ext 1011) > Tel/Fax: +54-342-4511169 > <matblocksize.diff>

1 0

Re: [petsc-dev] PCApplyTranspose in BiCG
by Jed Brown 10 Apr '12

10 Apr '12

On Tue, Apr 10, 2012 at 04:47, Alexander Grayver <agrayver(a)gfz-potsdam.de>wrote: > Barry, > > Well, this is not easy to apply bisection here since "good" petsc-dev had > bugs that I patched manually (including Intel MKL BLAS incompatibility) and > only after that I was able to run solver successfully. "Bad" petsc-dev has > all these issues patched and names of some routines have been changed as > well (MatCreateXXX)... However TFQMR doesn't work at all for any system I > tried. > > Unfortunately, the shortest time gap between "good" and "bad" I can give > so far is: > > HG revision: a4dae06b4d05ef5b78a9a9aa2cbc31**b2b158834f HG Date: Wed Feb > 08 11:29:52 2012 -0600 --- GOOD > HG revision: 59ca71b85d784abbe22f12c34ffcb2**ccc6dc68b0 HG Date: Tue Mar > 20 01:05:13 2012 -0500 --- BAD > > I tested them like that: > ./ex10 -f A.dat -rhs RHS.dat -ksp_type tfqmr -pc_type sor > -ksp_monitor_true_residual -ksp_max_it 500 > > 95 KSP preconditioned resid norm 2.117351054162e-18 true resid norm > 1.981564630346e-15 ||r(i)||/||b|| 3.893342156134e-05 ---- GOOD > > 500 KSP preconditioned resid norm 3.000274090666e-13 true resid norm > 4.577554424627e-11 ||r(i)||/||b|| 8.993895702653e-01 ---- BAD > $ ./ex10 -f ~/dl/A.dat -rhs ~/dl/RHS.dat -ksp_type tfqmr -pc_type sor -ksp_monitor_true_residual [0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: Arguments are incompatible! [0]PETSC ERROR: Zero diagonal on row 44! [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Petsc Development HG revision: 21fb71e3b67102b07ab64b3779f68e8719c17739 HG Date: Mon Mar 19 08:24:25 2012 -0500 [0]PETSC ERROR: See docs/changes/index.html for recent updates. [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. [0]PETSC ERROR: See docs/index.html for manual pages. [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: ./ex10 on a mpich named batura by jed Tue Apr 10 07:09:45 2012 [0]PETSC ERROR: Libraries linked from /home/jed/petsc/mpich/lib [0]PETSC ERROR: Configure run at Mon Mar 19 09:50:34 2012 [0]PETSC ERROR: Configure options --download-ams --download-blacs --download-chaco --download-hypre --download-metis --download-ml --download-parmetis --download-spai --download-spooles --download-sundials --download-superlu --download-superlu_dist --download-triangle --with-blas-lapack=/usr --with-c2html --with-cholmod-dir=/usr --with-lgrind --with-mpi-dir=/opt/mpich2 --with-pcbddc --with-pthreadclasses --with-shared-libraries --with-single-library=0 --with-sowing --with-umfpack-dir=/usr --with-x -PETSC_ARCH=mpich [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: MatInvertDiagonal_SeqAIJ() line 1315 in /home/jed/petsc/src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: MatSOR_SeqAIJ() line 1348 in /home/jed/petsc/src/mat/impls/aij/seq/aij.c [0]PETSC ERROR: MatSOR() line 3671 in /home/jed/petsc/src/mat/interface/matrix.c [0]PETSC ERROR: PCApply_SOR() line 35 in /home/jed/petsc/src/ksp/pc/impls/sor/sor.c [0]PETSC ERROR: PCApply() line 384 in /home/jed/petsc/src/ksp/pc/interface/precon.c [0]PETSC ERROR: KSPInitialResidual() line 64 in /home/jed/petsc/src/ksp/ksp/interface/itres.c [0]PETSC ERROR: KSPSolve_TFQMR() line 41 in /home/jed/petsc/src/ksp/ksp/impls/tfqmr/tfqmr.c [0]PETSC ERROR: KSPSolve() line 434 in /home/jed/petsc/src/ksp/ksp/interface/itfunc.c [0]PETSC ERROR: main() line 327 in src/ksp/ksp/examples/tutorials/ex10.c application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0 There is no space allocated at (44,44); note that counting is 0-based. This consistency check has been in PETSc for many years, so I don't know why you didn't trip it in earlier runs. Do you have some options in PETSC_OPTIONS or .petscrc? > > If it is of any help files are here: > http://dl.dropbox.com/u/**60982984/A.dat<http://dl.dropbox.com/u/60982984/A.dat> > http://dl.dropbox.com/u/**60982984/RHS.dat<http://dl.dropbox.com/u/60982984/RHS.dat> > > Is this possible to look up changes in one file only? $ hg log src/ksp/ksp/impls/tfqmr/ But the problem was clearly not in that file.

1 0

Re: [petsc-dev] PCApplyTranspose in BiCG
by Alexander Grayver 10 Apr '12

10 Apr '12

Barry, Well, this is not easy to apply bisection here since "good" petsc-dev had bugs that I patched manually (including Intel MKL BLAS incompatibility) and only after that I was able to run solver successfully. "Bad" petsc-dev has all these issues patched and names of some routines have been changed as well (MatCreateXXX)... However TFQMR doesn't work at all for any system I tried. Unfortunately, the shortest time gap between "good" and "bad" I can give so far is: HG revision: a4dae06b4d05ef5b78a9a9aa2cbc31b2b158834f HG Date: Wed Feb 08 11:29:52 2012 -0600 --- GOOD HG revision: 59ca71b85d784abbe22f12c34ffcb2ccc6dc68b0 HG Date: Tue Mar 20 01:05:13 2012 -0500 --- BAD I tested them like that: ./ex10 -f A.dat -rhs RHS.dat -ksp_type tfqmr -pc_type sor -ksp_monitor_true_residual -ksp_max_it 500 95 KSP preconditioned resid norm 2.117351054162e-18 true resid norm 1.981564630346e-15 ||r(i)||/||b|| 3.893342156134e-05 ---- GOOD 500 KSP preconditioned resid norm 3.000274090666e-13 true resid norm 4.577554424627e-11 ||r(i)||/||b|| 8.993895702653e-01 ---- BAD If it is of any help files are here: http://dl.dropbox.com/u/60982984/A.dat http://dl.dropbox.com/u/60982984/RHS.dat Is this possible to look up changes in one file only? On 07.04.2012 02:46, Barry Smith wrote: > A change in tfqmr that broke it? A change in some underlying thing that tfqmr that broke it. > > The "correct" way to track down the "bad" change is to us bisection with mecurial http://www.selenic.com/mercurial/hg.1.html#bisect it lets you quickly zoom into the change set that causes the change. > > Barry > > On Mar 23, 2012, at 3:35 AM, Alexander Grayver wrote: > >> Barry, >> >> gmres produces the same result. >> What does that mean? >> >> On 22.03.2012 21:10, Barry Smith wrote: >>> Ok, first thing to confirm is if other KSP also produce different results. Run both exactly the same but use -ksp_type gmres on both; does it produce different results or the same? >>> >>> Barry >>> >>> On Mar 22, 2012, at 12:30 PM, Alexander Grayver wrote: >>> >>>> Jed, >>>> >>>> Ok, I see. Initially I've changed solver because I observed one thing. I have two petsc-dev (one is older than the other) versions compiled. >>>> Running the same program with both gives different results: >>>> >>>> ./solveTest -ksp_view -ksp_converged_reason -ksp_monitor_true_residual -log_summary -mat_type aij -ksp_rtol 1.0e-7 -ksp_type tfqmr -pc_type sor -ksp_max_it 750 >>>> 0 KSP preconditioned resid norm 9.760858225033e-10 true resid norm 4.108026443997e-08 ||r(i)||/||b|| 1.000000000000e+00 >>>> 1 KSP preconditioned resid norm 9.497939702484e-10 true resid norm 3.591828545356e-08 ||r(i)||/||b|| 8.743440662618e-01 >>>> 1 KSP preconditioned resid norm 1.341212527158e-09 true resid norm 3.754099340764e-08 ||r(i)||/||b|| 9.138449793207e-01 >>>> 2 KSP preconditioned resid norm 9.481744679086e-10 true resid norm 3.710424003325e-08 ||r(i)||/||b|| 9.032132713620e-01 >>>> 2 KSP preconditioned resid norm 1.340878639500e-09 true resid norm 3.730261027312e-08 ||r(i)||/||b|| 9.080421166135e-01 >>>> 3 KSP preconditioned resid norm 9.481415577317e-10 true resid norm 3.725585425324e-08 ||r(i)||/||b|| 9.069039540306e-01 >>>> ... >>>> 750 KSP preconditioned resid norm 1.340874209607e-09 true resid norm 3.726506102951e-08 ||r(i)||/||b|| 9.071280708031e-01 >>>> Linear solve did not converge due to DIVERGED_ITS iterations 750 >>>> KSP Object: 1 MPI processes >>>> type: tfqmr >>>> maximum iterations=750, initial guess is zero >>>> tolerances: relative=1e-07, absolute=1e-50, divergence=10000 >>>> left preconditioning >>>> using PRECONDITIONED norm type for convergence test >>>> PC Object: 1 MPI processes >>>> type: sor >>>> SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 >>>> linear system matrix = precond matrix: >>>> Matrix Object: 1 MPI processes >>>> type: seqaij >>>> rows=1375920, cols=1375920 >>>> total: nonzeros=17671714, allocated nonzeros=17671714 >>>> total number of mallocs used during MatSetValues calls =0 >>>> not using I-node routines >>>> >>>> Using Petsc Development HG revision: 59ca71b85d784abbe22f12c34ffcb2ccc6dc68b0 HG Date: Tue Mar 20 01:05:13 2012 -0500 >>>> Configure options: --with-petsc-arch=openmpi-intel-complex-release-c --with-mpi-dir=/opt/mpi/intel/openmpi-1.4.2 --with-scalar-type=complex --with-precision=double --with-blas-lapack-dir=/usr/local/lib --with-x=0 --with-debugging=0 >>>> >>>> ####################################### >>>> And now second one: >>>> ./solveTest -ksp_view -ksp_converged_reason -ksp_monitor_true_residual -log_summary -mat_type aij -ksp_rtol 1.0e-7 -ksp_type tfqmr -pc_type sor -ksp_max_it 750 >>>> 0 KSP preconditioned resid norm 9.760858225034e-10 true resid norm 4.108026443997e-08 ||r(i)||/||b|| 1.000000000000e+00 >>>> 1 KSP preconditioned resid norm 9.497930576113e-10 true resid norm 3.591633653854e-08 ||r(i)||/||b|| 8.742966246242e-01 >>>> 1 KSP preconditioned resid norm 1.341211098002e-09 true resid norm 3.754049143693e-08 ||r(i)||/||b|| 9.138327600539e-01 >>>> 2 KSP preconditioned resid norm 9.481662292281e-10 true resid norm 3.709648505700e-08 ||r(i)||/||b|| 9.030244951612e-01 >>>> 2 KSP preconditioned resid norm 1.340863953080e-09 true resid norm 3.730358317575e-08 ||r(i)||/||b|| 9.080657995827e-01 >>>> 3 KSP preconditioned resid norm 9.481308142928e-10 true resid norm 3.725036905334e-08 ||r(i)||/||b|| 9.067704300630e-01 >>>> ... >>>> 230 KSP preconditioned resid norm 8.875661094394e-17 true resid norm 1.062030983067e-14 ||r(i)||/||b|| 2.585258389996e-07 >>>> Linear solve converged due to CONVERGED_RTOL iterations 230 >>>> KSP Object: 1 MPI processes >>>> type: tfqmr >>>> maximum iterations=750, initial guess is zero >>>> tolerances: relative=1e-07, absolute=1e-50, divergence=10000 >>>> left preconditioning >>>> using PRECONDITIONED norm type for convergence test >>>> PC Object: 1 MPI processes >>>> type: sor >>>> SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 >>>> linear system matrix = precond matrix: >>>> Matrix Object: 1 MPI processes >>>> type: seqaij >>>> rows=1375920, cols=1375920 >>>> total: nonzeros=17671714, allocated nonzeros=17671714 >>>> total number of mallocs used during MatSetValues calls =0 >>>> not using I-node routines >>>> >>>> Using Petsc Development HG revision: a4dae06b4d05ef5b78a9a9aa2cbc31b2b158834f HG Date: Wed Feb 08 11:29:52 2012 -0600 >>>> Configure options: --with-petsc-arch=openmpi-intel-complex-release-f-mkl --with-fortran-interfaces=1 --download-superlu --download-superlu_dist --download-mumps --download-parmetis --download-ptscotch --download-metis --with-scalapack-lib=/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_scalapack_lp64.a --with-scalapack-include=/opt/intel/Compiler/11.1/072/mkl/include --with-blacs-lib=/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_blacs_openmpi_lp64.a --with-blacs-include=/opt/intel/Compiler/11.1/072/mkl/include --with-mpi-dir=/opt/mpi/intel/openmpi-1.4.2 --with-scalar-type=complex --with-blas-lapack-lib="[/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_intel_lp64.a,/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_intel_thread.a,/opt/intel/Compiler/11.1/072/mkl/lib/em64t/libmkl_core.a,/opt/intel/Compiler/11.1/072/lib/intel64/libiomp5.a]" --with-precision=double --with-x=0 --with-debugging=0 >>>> >>>> >>>> Well, two version compiled with different lines, but still I don't think this is the reason... >>>> Where actually to look at? >>>> >>>> Thanks. >>>> >>>> On 22.03.2012 15:42, Jed Brown wrote: >>>>> On Thu, Mar 22, 2012 at 04:23, Alexander Grayver<agrayver(a)gfz-potsdam.de> wrote: >>>>> I've got the following error when applying bicg+sor to my problem: >>>>> >>>>> mt-soft/solveTest> ./solveTest -ksp_view -ksp_converged_reason -ksp_monitor_true_residual -log_summary -mat_type aij -ksp_rtol 1.0e-7 -ksp_type bicg -pc_type sor -ksp_max_it 750 >>>>> [0]PETSC ERROR: --------------------- Error Message ------------------------------------ >>>>> [0]PETSC ERROR: No support for this operation for this object type! >>>>> [0]PETSC ERROR: PC does not have apply transpose! >>>>> [0]PETSC ERROR: ------------------------------------------------------------------------ >>>>> [0]PETSC ERROR: Petsc Development HG revision: 59ca71b85d784abbe22f12c34ffcb2ccc6dc68b0 HG Date: Tue Mar 20 01:05:13 2012 -0500 >>>>> [0]PETSC ERROR: See docs/changes/index.html for recent updates. >>>>> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. >>>>> [0]PETSC ERROR: See docs/index.html for manual pages. >>>>> [0]PETSC ERROR: ------------------------------------------------------------------------ >>>>> [0]PETSC ERROR: Configure options --with-petsc-arch=openmpi-intel-complex-release-c --with-mpi-dir=/opt/mpi/intel/openmpi-1.4.2 --with-scalar-type=complex --with-precision=double --with-blas-lapack-dir=/usr/local/lib --with-x=0 --with-debugging=0 >>>>> [0]PETSC ERROR: ------------------------------------------------------------------------ >>>>> [0]PETSC ERROR: PCApplyTranspose() line 515 in /home/mt/agrayver/lib/petsc-dev/src/ksp/pc/interface/precon.c >>>>> [0]PETSC ERROR: KSPSolve_BiCG() line 55 in /home/mt/agrayver/lib/petsc-dev/src/ksp/ksp/impls/bicg/bicg.c >>>>> [0]PETSC ERROR: KSPSolve() line 434 in /home/mt/agrayver/lib/petsc-dev/src/ksp/ksp/interface/itfunc.c >>>>> [0]PETSC ERROR: main() line 114 in /home/solveTest.c >>>>> >>>>> Any ideas? >>>>> >>>>> PCApplyTranspose_SOR is not implemented. It is unfortunately not an operation that is efficient to implement in-place like SOR. >>>> -- >>>> Regards, >>>> Alexander >>>> >> >> -- >> Regards, >> Alexander >> >> -- Regards, Alexander

1 0