[petsc-dev] (no subject)

[Stefano Zampini]

3.10 works fine

Il giorno lun 22 lug 2019 alle ore 15:28 Stefano Zampini <
stefano.zampini at gmail.com> ha scritto:

> I just checked, and this is a bug also in maint.
>
> Il giorno lun 22 lug 2019 alle ore 13:23 Stefano Zampini via petsc-dev <
> petsc-dev at mcs.anl.gov> ha scritto:
>
>> Junchao,
>>
>> I found an issue with PetscSFReduceBegin/End. It seems that we can no
>> longer do (always worked before)
>>
>>     ierr =
>> PetscSFReduceBegin(sf,MPIU_INT,leaf,root,MPI_MAX);CHKERRQ(ierr);
>>     ierr =
>> PetscSFReduceBegin(sf,MPIU_INT,leaf2,root2,MPI_MAX);CHKERRQ(ierr);
>>     ierr = PetscSFReduceEnd(sf,MPIU_INT,leaf,root,MPI_MAX);CHKERRQ(ierr);
>>     ierr =
>> PetscSFReduceEnd(sf,MPIU_INT,leaf2,root2,MPI_MAX);CHKERRQ(ierr);
>>
>> You can reproduce with the currrent master
>> $ cd $PETSC_DIR/src/ksp/ksp/examples/tutorials
>> $ make ex71
>> $ mpiexec -n 4  ./ex71 -pde_type Elasticity -cells 7,9 -dim 2
>> -pc_bddc_levels 1 -pc_bddc_coarsening_ratio 2 -ksp_error_if_not_converged
>> -pc_bddc_monolithic -pc_bddc_use_faces
>> -pc_bddc_coarse_pc_bddc_corner_selection
>> -pc_bddc_coarse_l1_pc_bddc_corner_selection -mat_partitioning_type average
>>  -pc_bddc_coarse_pc_bddc_use_deluxe_scaling
>> -pc_bddc_coarse_sub_schurs_mat_solver_type petsc
>>
>> The attached patch shows that if we instead do
>>
>>     ierr =
>> PetscSFReduceBegin(sf,MPIU_INT,leaf,root,MPI_MAX);CHKERRQ(ierr);
>>     ierr = PetscSFReduceEnd(sf,MPIU_INT,leaf,root,MPI_MAX);CHKERRQ(ierr);
>>     ierr =
>> PetscSFReduceBegin(sf,MPIU_INT,leaf2,root2,MPI_MAX);CHKERRQ(ierr);
>>     ierr =
>> PetscSFReduceEnd(sf,MPIU_INT,leaf2,root2,MPI_MAX);CHKERRQ(ierr);
>>
>> everything works properly.
>> I think the issue is with the SF GetPack mechanism, but I couldn't figure
>> out where.
>>
>> --
>> Stefano
>>
>
>
> --
> Stefano
>


-- 
Stefano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20190722/0c043c2e/attachment.html>