First, added da object for aux. Now each of q and aux has da object (q_da, aux_da), local vector (lqVec, lauxVec) and global vector (gqVec, gauxVec). Second, removed creating da from grid.__init__() and added methods init_q_petsc_structures, init_aux_petsc_structures, fill_q_petsc_structures, fill_aux_petsc_structures for setting up petsc structures. Third, Modified grid method deepcopy to copy the petsc structures as well.
First, added da object for aux. Now each of q and aux has da object (q_da, aux_da), local vector (lqVec, lauxVec) and global vector (gqVec, gauxVec). Second, removed creating da from grid.__init__() and added methods init_q_petsc_structures, init_aux_petsc_structures, fill_q_petsc_structures, fill_aux_petsc_structures for setting up petsc structures. Third, Modified grid method deepcopy to copy the petsc structures as well.

unterweg/peanoclaw,unterweg/peanoclaw,unterweg/peanoclaw,unterweg/peanoclaw,unterweg/peanoclaw,unterweg/peanoclaw