[petsc-dev] bullshit in Plex example under complex

[Barry Smith]

> On Aug 29, 2015, at 3:04 PM, Tobin Isaac <tisaac at ices.utexas.edu> wrote:
> 
> On Fri, Aug 28, 2015 at 02:02:49PM -0500, Barry Smith wrote:
>> 
>>   This test should never be run in complex. if it is not coded for it. Please update the makefile rule for running it only without complex
> 
> How am I supposed to add a test to the intersection of a package
> (testexamples_TRIANGLE) and a scalar type (testexamples_C_NoComplex).
> It's pretty frustrating,

    There is no mechanism to handle the intersection of multiple packages or scalar types.

  PetscFVLeastSquaresPseudoInverseSVD_Static

 SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "I don't think this makes sense for complex numbers");

   Is it just this routine where this does not make sense or is that the entire concept of limiters makes no sense for complex numbers?

PetscFVLeastSquaresPseudoInverse_Static() does seem to make sense for complex numbers. How can that be?

 If possible just fix PetscFVLeastSquaresPseudoInverseSVD_Static() to work for complex. Should be just a couple of lines of code?

  Barry





> given how obvious this would be with
> something like sharness [1], that I couldn't easily figure this out.
> 
>  Toby
> 
> [1] http://events.linuxfoundation.org/sites/events/files/slides/Git_and_testing_linuxcon_0.pdf
> 
>> 
>> 
>>> Null vector 2 is likely null vector|| A * v[2] || = 9.22162e-17
>>> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
>>> [0]PETSC ERROR: No support for this operation for this object type
>>> [0]PETSC ERROR: I don't think this makes sense for complex numbers
>>> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
>>> [0]PETSC ERROR: Petsc Development GIT revision: v3.6.1-977-ga9ee8f1  GIT Date: 2015-08-27 23:56:46 -0500
>>> [0]PETSC ERROR: ./ex3 on a arch-opensolaris-cmplx-pkgs-dbg named n-gage by petsc Fri Aug 28 10:53:50 2015
>>> [0]PETSC ERROR: Configure options --with-debugger=/bin/true --with-scalar-type=complex --with-mpi-dir=/export/home/petsc/soft/mpich2-1.2.1p1 --with-c2html=0 --download-cmake=1 --download-metis=1 --download-parmetis=1 --download-triangle=1 --download-fblaslapack=1 --download-scalapack=1 --download-mumps=1 --with-no-output -PETSC_ARCH=arch-opensolaris-cmplx-pkgs-dbg -PETSC_DIR=/export/home/petsc/petsc.clone-3
>>> [0]PETSC ERROR: #1 PetscFVLeastSquaresPseudoInverseSVD_Static() line 2112 in /export/home/petsc/petsc.clone-3/src/dm/dt/interface/dtfv.c
>>> [0]PETSC ERROR: #2 PetscFVComputeGradient_LeastSquares() line 2197 in /export/home/petsc/petsc.clone-3/src/dm/dt/interface/dtfv.c
>>> [0]PETSC ERROR: #3 PetscFVComputeGradient() line 1729 in /export/home/petsc/petsc.clone-3/src/dm/dt/interface/dtfv.c
>>> [0]PETSC ERROR: #4 BuildGradientReconstruction_Internal_Tree() line 1790 in /export/home/petsc/petsc.clone-3/src/dm/impls/plex/plexgeometry.c
>>> [0]PETSC ERROR: #5 DMPlexComputeGradientFVM() line 1845 in /export/home/petsc/petsc.clone-3/src/dm/impls/plex/plexgeometry.c
>>> [0]PETSC ERROR: #6 DMPlexSNESGetGradientDM() line 927 in /export/home/petsc/petsc.clone-3/src/snes/utils/dmplexsnes.c
>>> [0]PETSC ERROR: #7 TestFVGrad() line 523 in /export/home/petsc/petsc.clone-3/src/dm/impls/plex/examples/tests/ex3.c
>>> [0]PETSC ERROR: #8 main() line 889 in /export/home/petsc/petsc.clone-3/src/dm/impls/plex/examples/tests/ex3.c
>>> [0]PETSC ERROR: PETSc Option Table entries:
>>> [0]PETSC ERROR: -dim 2
>>> [0]PETSC ERROR: -display 140.221.10.20:0.0
>>> [0]PETSC ERROR: -dm_plex_max_projection_height 1
>>> [0]PETSC ERROR: -dm_view ascii::ASCII_INFO_DETAIL
>>> [0]PETSC ERROR: -fv_grad
>>> [0]PETSC ERROR: -malloc_dump
>>> [0]PETSC ERROR: -mat_null_space_test_view
>>> [0]PETSC ERROR: -num_comp 2
>>> [0]PETSC ERROR: -petsclimiter_type none
>>> [0]PETSC ERROR: -petscpartitioner_type simple
>>> [0]PETSC ERROR: -petscspace_order 1
>>> [0]PETSC ERROR: -petscspace_poly_tensor
>>> [0]PETSC ERROR: -qorder 1
>>> [0]PETSC ERROR: -simplex 0
>>> [0]PETSC ERROR: -tree
>>> [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint at mcs.anl.gov----------
>>> application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0
>>> [1]PETSC ERROR: --------------------- Error Message -------------------------------------------------------
>>