[petsc-users] questions

[Barry Smith]

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