I'm interested in trying to use multi-grid techniques to solve my problem, but it's unclear to me how to define the right-hand side and system matrix for my case where the grids are not uniform (ie, the grid lengths x,y,z are variable, and the grid properties are variable).
Right now I'm solving my equations using a DA. So, on the finest scale I know how to compute the RHS and the Jacobian. It's unclear to me how PETSc does the multi-gridding, and therefore how I would write my code. Any pointers would be appreciated. Thanks, Randy -- Randall Mackie GSY-USA, Inc. PMB# 643 2261 Market St., San Francisco, CA 94114-1600 Tel (415) 469-8649 Fax (415) 469-5044 California Registered Geophysicist License No. GP 1034
