Hi Folks, I have a basic question regarding matrix size as created by subdomains in PCASM.
In the PCASMSetLocalSubdomains(PC pc,PetscInt n,IS is[],IS is_local[]) routine, I am defining is[] to be the local+overlapping DOFs and is_local [] to be the nonoverlapping local DOFs only. Now if I grab the subksp's from PCASM and print out the matrices, would the matrix size correspond to size of IS in is[] or in is_local[] dofs? Put in another words, do the subdomain matrices are defined for just local DOFs or with local+overlapping DOFs? Thanks, -- --Amneet
