How do we define the matrix for GMRES, if we do not want to solve the boundary elements but only the interior ones.I am using pentagonal elements so in a row there are 6 elements (cell itself + 5 fluxes) which means that in a row boundary elements are also accounted. But I need to make the matrix square according to the code here and I do not want to solve for boundary elements? What should I do?
GMRES: Making the matrix square without solving for boundaries
sparse;c++;boundary conditions;matrix;gmres