|
|
Blitz Support : |
From: Julian Cummings (cummings_at_[hidden])
Date: 2005-06-16 14:15:54
Hello Vladimir,
I have a question about this. How were you planning on representing the
scalar function g(u)? Is u simply a Cartesian coordinate or a more
complicated function itself?
Regards, Julian C.
Dr. Julian C. Cummings Office: PB-111
Caltech/CACR, MC 158-79 Phone: 626-395-2543
1200 E. California Blvd. Fax: 626-584-5917
Pasadena, CA 91125
> -----Original Message-----
> From: blitz-support-bounces_at_[hidden]
> [mailto:blitz-support-bounces_at_[hidden]] On Behalf Of
> Vladimir Chalupecky
> Sent: Tuesday, June 14, 2005 5:18 AM
> To: blitz-support_at_[hidden]
> Subject: [Blitz-support] divgrad stencil operator
>
> Hi all,
>
> I would like to implement a stencil operator for div(g(u)
> grad u) operator in 2D (and 3D) using the classical
> approximation (for simplicity only in 1D)
> div(g(u) grad u)|_i = \frac{g_{i+1/2}(u_{i
> +1}-u_i)-g_{i-1/2}(u_i-u_{i-1})}{h^2}
> where
> g_{i+1/2} = 1/2*(g(u_{i+1})+g(u_i))
> and
> g_{i-1/2} = 1/2*(g(u_{i-1})+g(u_i))
> What would be the best approach? Is it doable in an elegant
> and efficient way?
>
> Thanks, Vladimir
>
> --
> Ing. Vladimir Chalupecky
> Department of Mathematics
> CTU FNSPE Prague
>
> _______________________________________________
> Blitz-support mailing list
> Blitz-support_at_[hidden]
> http://www.oonumerics.org/mailman/listinfo.cgi/blitz-support
>
>