Visualization

 1real[int] xx(10), yy(10);
 2
 3mesh Th = square(5,5);
 4
 5fespace Vh(Th, P1);
 6
 7//plot scalar and vectorial FE function
 8Vh uh=x*x+y*y, vh=-y^2+x^2;
 9plot(Th, uh, [uh, vh], value=true, ps="three.eps", wait=true);
10
11//zoom on box defined by the two corner points [0.1,0.2] and [0.5,0.6]
12plot(uh, [uh, vh], bb=[[0.1, 0.2], [0.5, 0.6]],
13   wait=true, grey=true, fill=true, value=true, ps="threeg.eps");
14
15//compute a cut
16for (int i = 0; i < 10; i++){
17   x = i/10.;
18   y = i/10.;
19   xx[i] = i;
20   yy[i] = uh; //value of uh at point (i/10., i/10.)
21}
22plot([xx, yy], ps="likegnu.eps", wait=true);

 1// from: \url{http://en.wikipedia.org/wiki/HSV_color_space}
 2// The HSV (Hue, Saturation, Value) model defines a color space
 3// in terms of three constituent components:
 4// HSV color space as a color wheel
 5// Hue, the color type (such as red, blue, or yellow):
 6// Ranges from 0-360 (but normalized to 0-100% in some applications, like here)
 7// Saturation, the "vibrancy" of the color: Ranges from 0-100%
 8// The lower the saturation of a color, the more "grayness" is present
 9// and the more faded the color will appear.
10// Value, the brightness of the color: Ranges from 0-100%
11
12mesh Th = square(10, 10, [2*x-1, 2*y-1]);
13
14fespace Vh(Th, P1);
15Vh uh=2-x*x-y*y;
16
17real[int] colorhsv=[ // color hsv model
18   4./6., 1 , 0.5, // dark blue
19   4./6., 1 , 1, // blue
20   5./6., 1 , 1, // magenta
21   1, 1. , 1, // red
22   1, 0.5 , 1 // light red
23   ];
24 real[int] viso(31);
25
26 for (int i = 0; i < viso.n; i++)
27   viso[i] = i*0.1;
28
29 plot(uh, viso=viso(0:viso.n-1), value=true, fill=true, wait=true, hsv=colorhsv);

1{// file for gnuplot
2   ofstream gnu("plot.gp");
3   for (int i = 0; i < n; i++)
4      gnu << xx[i] << " " << yy[i] << endl;
5}
6
7// to call gnuplot command and wait 5 second (due to the Unix command)
8// and make postscript plot
9exec("echo 'plot \"plot.gp\" w l \n pause 5 \n set term postscript \n set output \"gnuplot.eps\" \n replot \n quit' | gnuplot");

1load "medit"
2
3mesh Th = square(10, 10, [2*x-1, 2*y-1]);
4
5fespace Vh(Th, P1);
6Vh u=2-x*x-y*y;
7
8medit("u", Th, u);

 1mesh Th = square(10, 10, [2*x-1, 2*y-1]);
 2
 3fespace Vh(Th, P1);
 4Vh u=2-x*x-y*y;
 5
 6savemesh(Th, "u", [x, y, u*.5]); //save u.points and u.faces file
 7// build a u.bb file for medit
 8{
 9   ofstream file("u.bb");
10   file << "2 1 1 " << u[].n << " 2 \n";
11   for (int j = 0; j < u[].n; j++)
12      file << u[][j] << endl;
13}
14//call medit command
15exec("ffmedit u");
16//clean files on unix-like OS
17exec("rm u.bb u.faces u.points");

 1load "iovtk"
 2
 3mesh Th = square(10, 10, [2*x-1, 2*y-1]);
 4
 5fespace Vh(Th, P1);
 6Vh u=2-x*x-y*y;
 7
 8int[int] Order = [1];
 9string DataName = "u";
10savevtk("u.vtu", Th, u, dataname=DataName, order=Order);

1include "ffmatlib.idp"
2
3mesh Th = square(10, 10, [2*x-1, 2*y-1]);
4fespace Vh(Th, P1);
5Vh u=2-x*x-y*y;
6
7savemesh(Th,"export_mesh.msh");
8ffSaveVh(Th,Vh,"export_vh.txt");
9ffSaveData(u,"export_data.txt");