Hi, I meet some trouble about my code. I'd like to use MatZeroRows to perform Dirichlet boundary condition. However, it requires MatAssembly, which packs the sparse matrix. Then the none zero pattern seems to be freezed. The continued MatSetValues may be very low efficient if the item not in the previous none zero pattern..
Is there any way to use MatZeroRows without freeze the none zero pattern? Or I have to redesign my code.... Regards Gong Ding
