[petsc-dev] PETSc VecView with MPIIO and complex numbers broken?

Dave May dave.mayhem23 at gmail.com
Mon Mar 2 22:39:44 CST 2015


Hi Barry,

I fixed all the PetscPrintf's to have the (double) cast.
I am wondering if this could this be an uninitialized variable issue?

In that code, I have the following declaration

  const PetscScalar test_values[] = { 0.311256, 88.068, 11.077444,
9953.62, 7.345, 64.8943, 3.1458, 6699.95, 0.00084, 0.0647 };

which is then used to build a Vec object.
I am not explicitly setting the imaginary part of the vector to zero when
test_values is declared.

Apologize for the crappy code -  I really never ever use complex numbers in
petsc.

If this is in fact the issue, please note that the same time of declaration
appears at the top of
  petsc / src / dm / examples / tutorials / ex15.c


Cheers
   Dave




On 2 March 2015 at 20:31, Barry Smith <bsmith at mcs.anl.gov> wrote:

>
>
> http://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2015/03/01/examples_master_arch-opensolaris-cmplx-pkgs-dbg_n-gage.log
>
> ******* Testing: testexamples_C *******
> 4,7c4,7
> <   min(a)   = +8.40e-04 [loc 8]
> <   max(a)   = +9.95e+03 [loc 3]
> <   min(b)   = +8.40e-04 [loc 8]
> <   max(b)   = +9.95e+03 [loc 3]
> ---
> >   min(a)   = -4.21e+207 [loc 6]
> >   max(a)   = +8.10e-320 [loc 5]
> >   min(b)   = -4.21e+207 [loc 6]
> >   max(b)   = +8.10e-320 [loc 5]
> /export/home/petsc/petsc.clone-3/src/vec/vec/examples/tests
> Possible problem with with ex46, diffs above
> =========================================
> 4,7c4,7
> <   min(a)   = +8.40e-04 [loc 8]
> <   max(a)   = +9.95e+03 [loc 3]
> <   min(b)   = +8.40e-04 [loc 8]
> <   max(b)   = +9.95e+03 [loc 3]
> ---
> >   min(a)   = -4.21e+207 [loc 6]
> >   max(a)   = +8.10e-320 [loc 5]
> >   min(b)   = -4.21e+207 [loc 6]
> >   max(b)   = +8.10e-320 [loc 5]
> /export/home/petsc/petsc.clone-3/src/vec/vec/examples/tests
> Possible problem with with ex46_2, diffs above
> =========================================
> 4,7c4,7
> <   min(a)   = +8.40e-04 [loc 8]
> <   max(a)   = +9.95e+03 [loc 3]
> <   min(b)   = +8.40e-04 [loc 8]
> <   max(b)   = +9.95e+03 [loc 3]
> ---
> >   min(a)   = -4.21e+207 [loc 6]
> >   max(a)   = +8.10e-320 [loc 5]
> >   min(b)   = -4.21e+207 [loc 6]
> >   max(b)   = +8.10e-320 [loc 5]
> /export/home/petsc/petsc.clone-3/src/vec/vec/examples/tests
> Possible problem with with ex46_3, diffs above
> =========================================
> 6,9c6,9
> <   min(a)   = +8.40e-04 [loc 8]
> <   max(a)   = +9.95e+03 [loc 3]
> <   min(b)   = +8.40e-04 [loc 8]
> <   max(b)   = +9.95e+03 [loc 3]
> ---
> >   min(a)   = -4.21e+207 [loc 6]
> >   max(a)   = +8.10e-320 [loc 5]
> >   min(b)   = -4.21e+207 [loc 6]
> >   max(b)   = +8.10e-320 [loc 5]
> /export/home/petsc/petsc.clone-3/src/vec/vec/examples/tests
> Possible problem with with ex46_mpiio, diffs above
> =========================================
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150303/98ffcefe/attachment.html>


More information about the petsc-dev mailing list