Padua2D: Lagrange interpolation at Padua points on bivariate domains

There are three drivers: drvrec (for rectangles), drvell (for ellipses) and
drvtri (for triangles). To compile,

$(FC) $(FFLAGS) drvrec.f pdpts.f cheb.f padua2.f pd2val.f tstfn2.f f.f -o drvrec -l$(BLAS)
$(FC) $(FFLAGS) drvell.f pdpts.f cheb.f padua2.f pd2val.f tstfn2.f f.f -o drvell -l$(BLAS)
$(FC) $(FFLAGS) drvtri.f pdpts.f cheb.f padua2.f pd2val.f tstfn2.f f.f -o drvtri -l$(BLAS)

where $(FC) is your Fortran77 compiler (usually f77 or g77), $(FFLAGS) the
compiler options (usually -O3 or -O2) and $(BLAS) your BLAS library (usually blas). 
If the BLAS library is not available on your system, replace -l$(BLAS) with blas/*.f 
in the command lines above. 

tstfn2.f is a subroutine by Robert J. Renka. It is included in

R. J. Renka, Algorithm 790: CSHEP2D: Cubic Shepard Method for Bivariate Interpolation 
of Scattered Data, ACM Trans. Math. Soft., 25(1), 1999, 70--73.

f.f is the Padua2D interface to tstfn2.f.

---------------------------------------------
Copyright (C) 2006 2007 

Marco Caliari and Marco Vianello <{mcaliari,marcov}@math.unipd.it>
Department of Pure and Applied Mathematics
University of Padua (Italy)

Stefano De Marchi <stefano.demarchi@univr.it>
Department of Computer Science
University of Verona (Italy)

