This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Is someone working on 2d interpolation?
- From: Achim Gädke <Achim dot Gaedke at physik dot tu-darmstadt dot de>
- To: gsl-discuss at sources dot redhat dot com
- Date: Tue, 24 Jun 2003 20:54:53 +0200
- Subject: Is someone working on 2d interpolation?
Hi everyone!
I'm using the gsl_matrix type to store values of a 2d magnetic field:
(x,z) -> H
For graphical evaluation, I need to interpolate these values.
Furthermore I plot isolines for H. I've already implemented an
interpolation formula for rectangular sample values (0,0),
(0,1),(1,0),(1,1) . Inside this square I use the interpolation
H=a*x+b*y+c*x*y+d with a=H(1,0)-H(0,0), b=H(0,1)-H(0,0),
c=H(1,1)+H(0,0)-H(0,1)-H(1,0), d=H(0,0)
This results in a continuous surface with hyperbola isolines.
But I'd like to have a smoother solution. Is anyone working on such code?
Thanks, Achim