<< Click to Display Table of Contents >> splinetable 

{ SPLINETABLE.PDE
This example solves the same system as TABLE.PDE, using a Spline interpretation of the
data in the table file 'TABLE.TBL'.
The file format is the same for TABLE or SPLINE TABLE input.
The SPLINE TABLE operator can be used to build spline tables of one or two dimensions.
The resulting interpolation is third order in the coordinates, with continuous values
and derivatives. First or second derivatives of the interpolated function may be computed.
Here the table is used as source and diffusivity in a fictitious heat equation, merely to
show the use of the table variable.
The SAVE function is used to construct a Finite Element interpolation of the data from the
spline table, for comparison of derivatives. Cubic FEM basis is used so that the second
derivative is meaningful.
}
title 'Spline Table Input Test'
select
regrid=off
variables
u
definitions
alpha = spline table('table.tbl') ! construct spline fit of table:
beta = 1/alpha
femalpha = save(alpha) ! save a FEM interpolation of table:
equations
U: div(alpha*grad(u)) + beta = 0
boundaries region 1 start(0,10) value(u) = 0 line to (0,0) to (10,0) to (10,10) to close monitors contour(u)
plots grid(x,y) contour(alpha) as 'table' contour(dx(alpha)) as 'dx(table)' contour(dy(alpha)) as 'dy(table)' vector(grad(alpha)) as 'grad(table)' surface(alpha) as 'table' contour(dxx(alpha)) as 'dxx(table)' contour(dxy(alpha)) as 'dxy(table)' contour(dyy(alpha)) as 'dyy(table)' contour(dxx(alpha)+dyy(alpha)) as "Table Curvature" contour(div(grad(femalpha))) as "FEM Curvature" surface(beta) as "table reciprocal" contour(u) as "temperature solution" surface(u) as "temperature solution"
end
