The restriction is the list of indices in the Outer subdomain:
that you get with -ksp_view
. You can define in your external solver a matrix of dimension (size of local matrix, that you get in the -ksp_view
as well) x (size of global matrix, in your case 5654). Then, you just fill with ones for all lines at each column index given in the -ksp_view
.