Lisandro : I fixed bug in branch hzhang/fix-ao2petscis: https://bitbucket.org/petsc/petsc/commits/260dbea8cadd3eeca7e59945a235ff92ceae1333
I tested "-iga_dof 3" as well. Let me know if you still encounter problems, or have any comments on this fix. Hong On 20 June 2017 at 01:46, Hong <[email protected]> wrote: > > Thanks for detailed instruction. I built PetIGA and reproduced error. > > > > Did you permute IS used as input for MatCreateSubMatrix? ISGetMinMax() > gives > > incorrect output min and max, thus delegates to incorrect private > > implementation routine. > > > > I don't think so, but look at the implementations: > > https://github.com/dalcinl/PetIGA/blob/master/src/petigamat.c#L34 > > The test fails for the case of bs=1 (default unless you pass let say > "-iga_dof 3" for bs=3 in the command line), so... > > Maybe OApplicationToPetscIS() is not handling appropriately stride index > sets? > > > > > > -- > Lisandro Dalcin > ============ > Research Scientist > Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) > Extreme Computing Research Center (ECRC) > King Abdullah University of Science and Technology (KAUST) > http://ecrc.kaust.edu.sa/ > > 4700 King Abdullah University of Science and Technology > al-Khawarizmi Bldg (Bldg 1), Office # 0109 > Thuwal 23955-6900, Kingdom of Saudi Arabia > http://www.kaust.edu.sa > > Office Phone: +966 12 808-0459 >
