[petsc-dev] NULL arg for AIJ field with MATIS

[Mark Adams]

I added a case to the SeqAIJ options switch:

 case MAT_COO_BLOCKED_ASSEMBLY:
    A->coo_blocked_assembly = flg;
    break;

And added the flag for this to mat/tests/ex123* and when we have '-mat_type
is', I get matis calling *MatSetOption* with a NULL argument (a->A)

What is wrong here? Not being able to use blocked assembly for matis is
fine, but this is not a useful error message.

x* thread #1, queue = 'com.apple.main-thread', stop reason = breakpoint 1.1
    frame #0: 0x0000000128c58223 libpetsc.3.023.dylib`PetscError(comm=0,
line=-1, func=0x0000000000000000, file="", n=PETSC_SUCCESS,
p=PETSC_ERROR_INITIAL, mess="*Null Pointer: Parameter # %d*") at err.c:387
    frame #1: 0x000000011fa3df99
libpetsc.3.023.dylib`*MatSetOption(mat=0x0000000000000000,
op=MAT_COO_BLOCKED_ASSEMBLY, flg=PETSC_TRUE) a*t matrix.c:6048:3
  * frame #2: 0x000000011ebdcfef
libpetsc.3.023.dylib`MatSetOption_IS(A=0x0000621000004350,
op=MAT_COO_BLOCKED_ASSEMBLY, flg=PETSC_TRUE) at matis.c:3286:3
    frame #3: 0x000000011fa4a7e0
libpetsc.3.023.dylib`MatSetOption(mat=0x0000621000004350,
op=MAT_COO_BLOCKED_ASSEMBLY, flg=PETSC_TRUE) at matrix.c:6127:3
    frame #4: 0x000000011fd87f66
libpetsc.3.023.dylib`MatSetFromOptions(B=0x0000621000004350) at
gcreate.c:320:12
    frame #5: 0x00000001000039d3 ex123`main(argc=11,
args=0x000000030412d1f0) at ex123.c:44:3
    frame #6: 0x0000000200018530 dyld`start + 3056
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20250703/76fc394c/attachment.html>