summaryrefslogtreecommitdiff
path: root/graphics/epix/samples/vfield.xp
diff options
context:
space:
mode:
Diffstat (limited to 'graphics/epix/samples/vfield.xp')
-rw-r--r--graphics/epix/samples/vfield.xp71
1 files changed, 0 insertions, 71 deletions
diff --git a/graphics/epix/samples/vfield.xp b/graphics/epix/samples/vfield.xp
deleted file mode 100644
index ba6e1fd8ba..0000000000
--- a/graphics/epix/samples/vfield.xp
+++ /dev/null
@@ -1,71 +0,0 @@
-/* -*-ePiX-*- */
-#include "epix.h"
-using namespace ePiX;
-
-// A spatial vector field depending on three variables
-P F(double x, double y, double z)
-{
- return -0.25*P(0.5*x+y, 0.5*y-x, 1+z);
-}
-
-domain R(P(-1,-1,-1), P(1,1,1), mesh(6,6,4), mesh());
-
-int main()
-{
- picture(P(-2,-2),P(2,2), "4x4in");
-
- begin();
- border();
-
- label(P(0, ymax()), P(0,-4), "$F(x,y,z)=-(x+2y, y-2x, 1+z)$",b);
-
- camera.at(P(11,8,5));
-
- // coordinate grids
- const int GRIDS(6);
- plain(Black(0.3));
- grid(P(-1, -1,-1), P(-1, 1, 1), GRIDS, GRIDS);
- grid(P(-1, -1,-1), P( 1,-1, 1), GRIDS, GRIDS);
- grid(P(-1, -1,-1), P(1, 1, -1), GRIDS, GRIDS);
-
- // flow lines
- bold();
- for (int i=0; i<12; ++i)
- {
- double t(M_PI_2+i*M_PI/6);
- pen(RGB(0.25*(3+Sin(t)), 0.25, 0.25*(3+Cos(t))));
- ode_plot(F, cyl(1.4, t, 1), 0, 10, 60);
- }
-
- // draw arrows in plain black borderd by a white edge;
- // vectors will appear to hide vectors behind them
- plain(Black());
- base(White(), "2pt");
-
- // the field is drawn in successive slices x_3 = const
-#ifdef VECT
- vector_field(F, R, 0.75); // draw arrowheads at 75% of default size
-#else
- dart_field(F, R, 0.5); // draw darts at 50% of default length
-#endif
-
- const double DX(0.125); // label padding
- font_size("scriptsize");
- degrees();
-
- label_angle(30); // adjusted visually; could calculate from camera position
- label(P(1+DX, -1, -1), P(-4,-2), "$y=-1$", bl);
- label(P(1+DX, 0, -1), P(-4,-2), "$y=0$", bl);
- label(P(1+DX, 1, -1), P(-4,-2), "$y=1$", bl);
-
- label_angle(-15);
- label(P(-1, 1+DX, -1), P(4,0), "$x=-1$", br);
- label(P( 0, 1+DX, -1), P(4,0), "$x=0$", br);
- label(P( 1, 1+DX, -1), P(4,0), "$x=1$", br);
-
- label(P( 1, -1-DX, -1), P(-4,0), "$z=-1$", tl);
- label(P( 1, -1-DX, 0), P(-4,0), "$z=0$", tl);
- label(P( 1, -1-DX, 1), P(-4,0), "$z=1$", tl);
-
- end();
-}