Update: Displacement fields the 3rd

To make the plotting of displacement fields more comprehensive, I deprecated plotDisplacementField in favor of a plot-method for objects of class DisplacementField, returning (and optionally plotting) an object of the (newly introduced) class DisplacementPlot, that in turn has its own plot-method.

For better visualization, and as rgl has no rendering options for arrows, I decided to render the starting positions as points, leading to lollipop graphs (Fig. 1 shows the resulting plot, using the example from yesterday). The code for the plot below is simply:

plot(dispfield)
wire3d(dummytrans)

example 1
Fig. 1: Pimped displacement field rendering - with reference mesh added as wireframe

And finally, there is a new complementing function invertDisplacementField to invert an existing displacement field.