[petsc-users] questions
Barry Smith
bsmith at mcs.anl.gov
Wed Nov 2 15:53:27 CDT 2011
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
More information about the petsc-users
mailing list