fdaPDE (version 1.1-21)

image.FEM.time: Image plot of a 2D FEM.time object at a given time

Description

Image plot of a FEM.time object, generated by the function FEM.time or returned by smooth.FEM.time. Only FEM objects defined over a 2D mesh can be plotted with this method.

Usage

# S3 method for FEM.time
image(x,t,lambdaS=1,lambdaT=1,num_refinements=NULL,...)

Value

No return value

Arguments

x

A 2D-mesh FEM.time object.

t

time at which do the plot

lambdaS

index of the space penalization parameter to use for the plot, useful when FEM.time returned by smooth.FEM.time using GCV

lambdaT

index of the time penalization parameter to use for the plot, useful when FEM.time returned by smooth.FEM.time using GCV

num_refinements

A natural number specifying how many bisections should by applied to each triangular element for plotting purposes. This functionality is useful where a discretization with 2nd order Finite Element is applied.

...

Arguments representing graphical options to be passed to plot3d.

See Also

FEM.time image.FEM.time

Examples

library(fdaPDE)
## Upload the horseshoe2D data
data(horseshoe2D)
boundary_nodes = horseshoe2D$boundary_nodes
boundary_segments = horseshoe2D$boundary_segments
locations = horseshoe2D$locations

## Create the 2D mesh
mesh = create.mesh.2D(nodes = rbind(boundary_nodes, locations), segments = boundary_segments)
## Create the FEM basis
FEMbasis = create.FEM.basis(mesh)
## Compute the coeff vector evaluating the desired function at the mesh nodes
## In this case we consider the fs.test() function introduced by Wood et al. 2008
time = 1:5
coeff = rep(fs.test(mesh$nodes[,1], mesh$nodes[,2]),5)*time
## Create the FEM.time object
FEM_time_function = FEM.time(coeff=coeff, time_mesh=1:5,FEMbasis=FEMbasis,FLAG_PARABOLIC=TRUE)

## Plot the FEM function
t = c(1.2,1.5,3.6,2.4,4.5)
image(FEM_time_function,t)