Within a program the keyword plot(u)
will display u
.
Instruction save('filename',u) will save the data array u
on
disk. If u
is a scalar variable then the (single) value of u
is
appended to the file (this is useful for time dependent problems or any
problem with iteration loop.).
Instruction savemesh('filename') will save the triangulation on disk.
Similarly for reading data with load('filename',u) and loadmesh('filename') .
The file must be in the default directory, else it won't be found. The
file format is best seen by opening them with a text editor. For a data
array f
it is:
ns f[0] .... f[ns-1](
ns
is the number of vertices)
If f
is a constant, its single value is appended to the end of the file;
this is useful for time dependent problems or any problem with iteration
loop.
If precise
is set still the function stored by save
is interpolated on the vertices as the continuous function given by mass lumping (see above).
For triangulations the file format is (nt
= number of triangles):
ns nt q[0].x q[0].y ib[i] ... q[n-1].x q[n-1].y ib[n-1] me[0][0] me[0][1] me[0][2] ngt[0] ... me[n-1][0] me[n-1][1] me[n-1][2] ngt[n-1]