[petsc-users] [petsc-maint #61421] ksp/examples/tutorials/Ex2.c: good with with-debugging but error without-debugging

Mon Jan 17 14:27:25 CST 2011

  Compiler bug. Immediately before the call to MatMult() in the code add the two lines

ierr = VecView(u,0);
ierr = MatView(A,0);

how large are the two objects? Given the code it is inconceivable that suddenly the vector length becomes 57.

   Barry

On Jan 17, 2011, at 2:43 AM, Mengda Wu wrote:

> Hi all,
> 
>    I just compiled the debugged and optimized versions of petsc-3.1-p7.
> Both are successful. I am running on Windows Vista 64bit machine.
> The C/C++ compiler is cl.exe from visual studio 2005 (64 bit) and no FORTRAN
> compiler is used. BLAS/LAPACK
> support comes from Intel MKL-10.1.3.028 (under em64t\lib). No MPI is used.
> 
>   The debugged petsc was configured with:
>   $ ./config/configure.py --with-cc='cl' --with-fc=0  --with-cxx='cl'
> --with-mpi=
> 0
> --with-blas-lapack-lib=[mkl_intel_lp64.lib,mkl_intel_thread.lib,mkl_core.lib
> ,libiomp5mt.lib] -CFLAGS='-MDd -W3 -Z7' -CXXFLAGS='-MDd -W3 -Z7'
> 
>   The optimized petsc was configured with:
>   $ ./config/configure.py --with-cc='cl' --with-fc=0  --with-cxx='cl'
> --with-mpi=
> 0 --with-debugging=0
> --with-blas-lapack-lib=[mkl_intel_lp64.lib,mkl_intel_threa
> d.lib,mkl_core.lib,libiomp5mt.lib]  -CFLAGS='-MD -wd4996 -O2' -CXXFLAGS='-MD
> -w
> d4996 -O2'
> 
>   When I run ksp/examples/tutorials/Ex2.c. The result with debugged petsc
> is
> =================================================================
>   Norm of error 0.000156044 iterations 6
> =================================================================
> 
>   However, there are errors with the optimized petsc with the output as
> follows:
> =================================================================
> [0]PETSC ERROR: --------------------- Error Message
> ----------------------------
> --------
> [0]PETSC ERROR: Nonconforming object sizes!
> [0]PETSC ERROR: Mat mat,Vec y: global dim 56 57!
> [0]PETSC ERROR:
> ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 7, Mon Dec 20 14:26:37
> CST 20
> 10
> [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:
> D:\Develop\Test\PETSc\petsc-3.1-p7\src\ksp\ksp\examples\tutorial
> s\ex2.exe on a cygwin-c- named CVBRL-38 by mengda Mon Jan 17 00:35:15 2011
> [0]PETSC ERROR: Libraries linked from
> /cygdrive/d/Develop/Test/PETSc/petsc-3.1-p
> 7/cygwin-c-opt/lib
> [0]PETSC ERROR: Configure run at Sun Jan 16 23:34:25 2011
> [0]PETSC ERROR: Configure options --with-cc=cl --with-fc=0 --with-cxx=cl
> --with-
> mpi=0 --with-debugging=0
> --with-blas-lapack-lib="[mkl_intel_lp64.lib,mkl_intel_t
> hread.lib,mkl_core.lib,libiomp5mt.lib]" -CFLAGS="-MD -wd4996 -O2"
> -CXXFLAGS="-MD
> -wd4996 -O2" --useThreads=0
> [0]PETSC ERROR:
> ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: MatMult() line 1888 in
> src/mat/interface/D:\Develop\Test\PETSc\P
> ETSC-~1.1-P\src\mat\INTERF~1\matrix.c
> [0]PETSC ERROR: main() line 146 in
> src/ksp/ksp/examples/tutorials/D:\Develop\Tes
> t\PETSc\PETSC-~1.1-P\src\ksp\ksp\examples\TUTORI~1\ex2.c
> 
> This application has requested the Runtime to terminate it in an unusual
> way.
> Please contact the application's support team for more information.
> =================================================================
> 
>   I am wondering what problems may lead to the errors. Please let me know
> if you need more
> information.
> 
> Thanks,
> Mengda
> 
> Hi all,
> 
>     I just compiled the debugged and optimized versions of petsc-3.1-p7. Both are successful. I am running on Windows Vista 64bit machine.
> The C/C++ compiler is cl.exe from visual studio 2005 (64 bit) and no FORTRAN compiler is used. BLAS/LAPACK
> support comes from Intel MKL-10.1.3.028 (under em64t\lib). No MPI is used.
> 
>    The debugged petsc was configured with:
>    $ ./config/configure.py --with-cc='cl' --with-fc=0  --with-cxx='cl' --with-mpi=
> 0  --with-blas-lapack-lib=[mkl_intel_lp64.lib,mkl_intel_thread.lib,mkl_core.lib
> ,libiomp5mt.lib] -CFLAGS='-MDd -W3 -Z7' -CXXFLAGS='-MDd -W3 -Z7'
> 
>    The optimized petsc was configured with:
>    $ ./config/configure.py --with-cc='cl' --with-fc=0  --with-cxx='cl' --with-mpi=
> 0 --with-debugging=0 --with-blas-lapack-lib=[mkl_intel_lp64.lib,mkl_intel_threa
> d.lib,mkl_core.lib,libiomp5mt.lib]  -CFLAGS='-MD -wd4996 -O2' -CXXFLAGS='-MD -w
> d4996 -O2'
>   
>    When I run ksp/examples/tutorials/Ex2.c. The result with debugged petsc is
> =================================================================
>    Norm of error 0.000156044 iterations 6
> =================================================================
>   
>    However, there are errors with the optimized petsc with the output as follows:
> =================================================================
> [0]PETSC ERROR: --------------------- Error Message ----------------------------
> --------
> [0]PETSC ERROR: Nonconforming object sizes!
> [0]PETSC ERROR: Mat mat,Vec y: global dim 56 57!
> [0]PETSC ERROR: ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 7, Mon Dec 20 14:26:37 CST 20
> 10
> [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: D:\Develop\Test\PETSc\petsc-3.1-p7\src\ksp\ksp\examples\tutorial
> s\ex2.exe on a cygwin-c- named CVBRL-38 by mengda Mon Jan 17 00:35:15 2011
> [0]PETSC ERROR: Libraries linked from /cygdrive/d/Develop/Test/PETSc/petsc-3.1-p
> 7/cygwin-c-opt/lib
> [0]PETSC ERROR: Configure run at Sun Jan 16 23:34:25 2011
> [0]PETSC ERROR: Configure options --with-cc=cl --with-fc=0 --with-cxx=cl --with-
> mpi=0 --with-debugging=0 --with-blas-lapack-lib="[mkl_intel_lp64.lib,mkl_intel_t
> hread.lib,mkl_core.lib,libiomp5mt.lib]" -CFLAGS="-MD -wd4996 -O2" -CXXFLAGS="-MD
>  -wd4996 -O2" --useThreads=0
> [0]PETSC ERROR: ----------------------------------------------------------------
> --------
> [0]PETSC ERROR: MatMult() line 1888 in src/mat/interface/D:\Develop\Test\PETSc\P
> ETSC-~1.1-P\src\mat\INTERF~1\matrix.c
> [0]PETSC ERROR: main() line 146 in src/ksp/ksp/examples/tutorials/D:\Develop\Tes
> t\PETSc\PETSC-~1.1-P\src\ksp\ksp\examples\TUTORI~1\ex2.c
> 
> This application has requested the Runtime to terminate it in an unusual way.
> Please contact the application's support team for more information.
> =================================================================
> 
>    I am wondering what problems may lead to the errors. Please let me know if you need more
> information.
> 
> Thanks,
> Mengda