<div dir="auto">mat/examples/tests/ex1.c<div dir="auto"><br></div><div dir="auto">Yes, it raises an error for the matrix types aij and mffd, but the -h output is still relevant. I'll open a PR tomorrow<br><div dir="auto"><br></div><div dir="auto"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">Il Dom 9 Set 2018, 21:00 Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
> On Sep 9, 2018, at 4:54 AM, Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com" target="_blank" rel="noreferrer">stefano.zampini@gmail.com</a>> wrote:<br>
> <br>
> I just noticed a strange behaviour in master. Take mat/examples/ex1.c <br>
<br>
   There is no such example and the examples in the tutorials and tests subdirectories don't work with the given arguments so where is this example?<br>
<br>
<br>
> <br>
> ./ex1 -h -mat_type mpiaij <br>
> <br>
> does not print the help for the specific MPIAIJ options, only the MPIAIJ options header<br>
> <br>
> ....<br>
> Matrix (Mat) options -------------------------------------------------<br>
>   -mat_block_size <-1>: Set the blocksize used to store the matrix (MatSetBlockSize)<br>
>   -mat_type <aij>: Matrix type (one of) mffd mpimaij seqmaij maij is shell composite mpiaij seqaij mpiaijperm seqaijperm seqaijcrl mpiaijcrl mpibaij seqbaij mpisbaij seqsbaij mpidense seqdense mpiadj scatter blockmat nest mpisell seqsell preallocator dummy hypre (MatSetType)<br>
>   -mat_is_symmetric: Checks if mat is symmetric on MatAssemblyEnd() (MatIsSymmetric)<br>
>   -mat_is_symmetric <0.>: Checks if mat is symmetric on MatAssemblyEnd() (MatIsSymmetric)<br>
>   -mat_null_space_test: <FALSE> Checks if provided null space is correct in MatAssemblyEnd() (MatSetNullSpaceTest)<br>
>   MPIAIJ options<br>
>   -mat_new_nonzero_location_err: <FALSE> Generate an error if new nonzeros are created in the matrix structure (useful to test preallocation) (MatSetOption)<br>
> <br>
> Instead, if you run e.g. MFFD, it works <br>
> <br>
> Matrix (Mat) options -------------------------------------------------<br>
>   -mat_block_size <-1>: Set the blocksize used to store the matrix (MatSetBlockSize)<br>
>   -mat_type <aij>: Matrix type (one of) mffd mpimaij seqmaij maij is shell composite mpiaij seqaij mpiaijperm seqaijperm seqaijcrl mpiaijcrl mpibaij seqbaij mpisbaij seqsbaij mpidense seqdense mpiadj scatter blockmat nest mpisell seqsell preallocator dummy hypre (MatSetType)<br>
>   -mat_is_symmetric: Checks if mat is symmetric on MatAssemblyEnd() (MatIsSymmetric)<br>
>   -mat_is_symmetric <0.>: Checks if mat is symmetric on MatAssemblyEnd() (MatIsSymmetric)<br>
>   -mat_null_space_test: <FALSE> Checks if provided null space is correct in MatAssemblyEnd() (MatSetNullSpaceTest)<br>
> Matrix-free Finite Differencing (MatMFFD) options -------------------------------------------------<br>
>   -mat_mffd_type <(null)>: Matrix free type (one of) ds wp (MatMFFDSetType)<br>
>   -mat_mffd_err <1.49012e-08>: set sqrt relative error in function (MatMFFDSetFunctionError)<br>
>   -mat_mffd_period <1>: how often h is recomputed (MatMFFDSetPeriod)<br>
>   -mat_mffd_check_positivity: <FALSE> Insure that U + h*a is nonnegative (MatMFFDSetCheckh)<br>
>   -mat_new_nonzero_location_err: <FALSE> Generate an error if new nonzeros are created in the matrix structure (useful to test preallocation) (MatSetOption)<br>
>   -mat_new_nonzero_allocation_err: <FALSE> Generate an error if new nonzeros are allocated in the matrix structure (useful to test preallocation) (MatSetOption)<br>
> <br>
> What is the proper fix? remove the PetscOptionsBegin call and replace PetscOptionsEnd() with PetscOptionsTail()?<br>
> <br>
> -- <br>
> Stefano<br>
<br>
</blockquote></div>