<div dir="ltr"><div>Hello PETSc Users/Developers,</div><div><br></div><div>A collegue of mine is looking into implementing an adaptive implicit method (AIM) over</div><div>PETSc in our simulator. This has led to some interesting questions about what can</div><div>be done with blocked matrices, which I'm not able to answer myself - does anyone have</div><div>any insight?</div><div><br></div><div>Apparently it would be ideal if we could find a matrix (and vector) type that supports a kind</div><div>of mixed block size:</div><div><br></div><div>
"For AIM [...] we will have matrix elements of various 
shapes: 1x1, 1xN, Nx1 and NxN. [...]. The 
solution and residual will be a mix of 1 and N variable/cell block"</div><div><br></div><div>There are ideas for how to implement what we want using the fixed-block-size objects we</div><div>understand well, but if anything like the above exists it would be very exciting.</div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Daniel<br>

</div></div>