set ASM subdomains in sub KSP

[Andrew T Barker]

I want to do multiple sweeps of ASM to precondition GMRES.  So I do something like 

-ksp_type gmres -pc_type ksp -ksp_ksp_type richardson -ksp_pc_type asm

which works fine.  Now I want to set my own subdomains with PCASMSetLocalSubdomains(), which isn't accessible from the command line.  

I can use PCKSPGetKSP() to get the KSP (and then the PC) to use in PCASMSetLocalSubdomains(), but I have to call KSPSetUp() on the parent KSP in order to use PCKSPGetKSP().  And then I'm not allowed to call PCASMSetLocalSubdomains() if KSPSetUp() has already been called - object is in wrong state.

I'm at a loss for how to solve this problem - any help would be appreciated.  

Andrew
---
Andrew T. Barker
andrew.barker at colorado.edu
Department of Applied Mathematics
University of Colorado, Boulder
526 UCB, Boulder, CO 80309-0526