Re: [Libmesh-users] matrix free scheme

2008-02-01 Thread Benjamin Kirk
> This could be useful to solve Euler equation with a DG method that > employs a Riemann solver to compute the numerical fluxes > (discontinuous at element faces-edges) without dealing with Godunov > fluxes to compute the jacobian... Am I wrong? Absolutely. And a common trick from the finite volu

Re: [Libmesh-users] matrix free scheme

2008-02-01 Thread Lorenzo Botti
This could be useful to solve Euler equation with a DG method that employs a Riemann solver to compute the numerical fluxes (discontinuous at element faces-edges) without dealing with Godunov fluxes to compute the jacobian... Am I wrong? Lorenzo 2008/2/1, Benjamin Kirk <[EMAIL PROTECTED]>: > > >

Re: [Libmesh-users] matrix free scheme

2008-02-01 Thread Benjamin Kirk
> Let me rewrite the expression you wrote as > > J*v = (F(u+epsilon*v)-F(u)) / epsilon > > Where epsilon is a small perturbation and F(u) is the nonlinear residual > function and J is the jacobian matrix of the nonlinear system. The above > formula compute the action of a Jacobian on a given vec

Re: [Libmesh-users] matrix free scheme

2008-01-31 Thread Vijay M
ahadevan Cc: [email protected] Subject: Re: [Libmesh-users] matrix free scheme thanx a lot! Can you write me the exact expression of the problem you're solving? Recently, I read the paper of Knoll about JFNK method (Jacobian-free Newton–Krylov methods: a survey of approaches

Re: [Libmesh-users] matrix free scheme

2008-01-31 Thread li pan
thanx a lot! Can you write me the exact expression of the problem you're solving? Recently, I read the paper of Knoll about JFNK method (Jacobian-free Newton–Krylov methods: a survey of approaches and applications). I'm not sure whether I understood well. Please correct me. If r_0 is the first resi

Re: [Libmesh-users] matrix free scheme

2008-01-30 Thread Vijay S. Mahadevan
I've attached the message from Ben Kirk. Also attached is the example program that he sent out. See below for his original mail. On a side note, I have written a working code with Matrix-free technique to solve a nonlinear, transient diffusion-reaction problem. It still needs few refinements and I

[Libmesh-users] matrix free scheme

2008-01-30 Thread li pan
Dear all, I remember that there was a discussion about matrix free scheme with libmesh before X'mas. I'd like to ask if somebody has got a example code for this. thanx pan Never miss a thing. Make Yahoo