On Wed, Nov 2, 2011 at 8:38 PM, xiaohl <span dir="ltr"><<a href="mailto:xiaohl@ices.utexas.edu">xiaohl@ices.utexas.edu</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi<br>
<br>
I am going to implement cell center difference method for<br>
u = - K grad p<br>
div u = f<br>
where p is the pressure , u is the velocity, f is the source term.<br>
<br>
my goal is to assemble the matrix and test the performance of different linear solvers in parallel.<br>
<br>
my question is how can I read the input file for K where K is n*n tensor.<br></blockquote><div><br></div><div>MatLoad()</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
second one is that do you have any similar examples?<br></blockquote><div><br></div><div>Nothing with the mixed-discretization of the Laplacian.</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hailong<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>