On Nov 2, 2011, at 3:42 PM, Matthew Knepley wrote: > On Wed, Nov 2, 2011 at 8:38 PM, xiaohl <xiaohl at ices.utexas.edu> wrote: > Hi > > I am going to implement cell center difference method for > u = - K grad p > div u = f > where p is the pressure , u is the velocity, f is the source term. > > my goal is to assemble the matrix and test the performance of different > linear solvers in parallel. > > my question is how can I read the input file for K where K is n*n tensor. > > MatLoad()
Hm, I think you should use a DMDA with n*n size dof and then use VecLoad() to load the entries of K. Barry > > second one is that do you have any similar examples? > > Nothing with the mixed-discretization of the Laplacian. > > Matt > > Hailong > > > > -- > What most experimenters take for granted before they begin their experiments > is infinitely more interesting than any results to which their experiments > lead. > -- Norbert Wiener
