I am trying to do Monte Carlo simulation for a large problem which requires eigensolution of a matrix for each sample. The matrix itself is quite large so much so that I want the eigensolution itself to be parallel as well.
This will probably require dividing the processors in subgroups and then within each group this eigensolution will be parallelized. I am trying to use PETSc and SLEPC for this.
Simple parallelization was possible for me but I am unable make use of such subgrouping. Does anybody know, how this can be done?