%% colorflow2.ldf
%% Created by Laurence D. Finston (LDF) Do 17. Aug 13:20:08 CEST 2023

%% * (1) Copyright and License.

%%%% Copyright (C) 2023 Laurence Finston

%%%% This file is free software; you can redistribute it and/or modify 
%%%% it under the terms of the GNU General Public License as published by 
%%%% the Free Software Foundation; either version 3 of the License, or 
%%%% (at your option) any later version.  

%%%% This file is distributed in the hope that it will be useful, 
%%%% but WITHOUT ANY WARRANTY; without even the implied warranty of 
%%%% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the 
%%%% GNU General Public License for more details.  

%%%% You should have received a copy of the GNU General Public License 
%%%% along with this file; if not, write to the Free Software 
%%%% Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA 

%%%% Please send bug reports to Laurence.Finston@gmx.de

%%%% The author can be contacted at: 

%%%% Laurence Finston 
%%%% Laurence.Finston@gmx.de

input "plainldf.lmc";

medium_pen := pencircle scaled (.333mm, .333mm, .333mm);

verbatim_tex "\magnification=\magstep2\font\huge=cmbx12 scaled 4000"
  & "\font\Large=cmbx12 scaled 3000\font\large=cmbx12 scaled 1600";

verbatim_metapost "prologues := 3;outputtemplate := \"%j_%3c.eps\";dotlabeldiam:=2mm;";

% verbatim_metapost   "outputformat:=\"png\";outputtemplate := \"%j_%3c.png\";"
%                   & "outputformatoptions := \"format=rgb antialias=none\";";

%% * (1)

point p[];
path q[];
rectangle r[];
transform t[];
picture v[];

triangle T[];

circle c[];

string s;

focus f;

set f with_position (0, 10, -40) with_direction (0, 10, 10) 
  with_distance 45;

numeric frame_wd;
numeric frame_ht;

% \special{papersize=400mm, 225mm} %% Frame size, proportion 16x9
% \hsize=400mm
% \vsize=225mm

%% (* (/ 300 16.0) 9) 168.75

frame_wd := 300mm;
frame_ht := 168.75mm;

rectangle frame;

frame := ((-.5frame_wd, -.5frame_ht), (.5frame_wd, -.5frame_ht),
          (.5frame_wd, .5frame_ht), (-.5frame_wd, .5frame_ht));

p10 := get_center frame;

%% * (1)

pickup medium_pen;

pickup medium_pen;

dot_pen := pencircle scaled (2mm, 2mm, 2mm);

bool do_labels;
do_labels := false; % true; % 


color C[];

%% * (1)

for i = 0 upto 3:
  p[i] := get_point (i) frame;
endfor;

for i = 0 upto 3:
  p[i+4] := mediate(get_point (i) frame, get_point((i+1) mod 4) frame);
endfor;

%% * (1)

beginfig(300);
  unfill frame;
  draw frame with_pen big_pen;
  
  s :=   "\vbox to 0pt{\vss\hbox to 0pt{"
       & "\hss\vbox{\hbox to 300mm{\hfil{\huge Color Flow Test Nr.~2}\hfil}\vskip2cm"
       & "\hbox to 300mm{\hfil{\Large By}\hfil}\vskip1cm"
       & "\hbox to 300mm{\hfil{\Large Laurence Finston}\hfil}\vskip2.5cm"
       & "\hbox to 300mm{\hfil{\large Copyright (C) 2023 Laurence Finston}\hfil}"
       & "}\hss}\vss}";

  label(s, p10);
endfig with_projection parallel_x_y;

beginfig(301);
  unfill frame;
  s :=   "\vbox to 0pt{\vss\hbox to 0pt{\hss\hbox to 300mm{\hfil{\huge Screen Wipes}\hfil}\hss}\vss}";
  
  draw frame with_pen big_pen;
  label(s, p10);
endfig with_projection parallel_x_y;



beginfig(400);
  unfill frame;
  s :=   "\vbox to 0pt{\vss\hbox to 0pt{\hss\hbox to 300mm{\hfil{\huge The End}\hfil}\hss}\vss}";
  
  draw frame with_pen big_pen;
  label(s, p10);
endfig with_projection parallel_x_y;

%% * (1)

point a[];

beginfig(0);
  unfill frame;
  draw frame with_pen big_pen;

  if do_labels:
    dotlabel.bot("$p_0$", p0) shifted (0, -2);
    dotlabel.bot("$p_1$", p1) shifted (0, -2);
    dotlabel.top("$p_2$", p2) shifted (0, 2);
    dotlabel.top("$p_3$", p3) shifted (0, 2);
    dotlabel.bot("$p_4$", p4) shifted (0, -2);
    dotlabel.rt("$p_5$", p5) shifted (3, 0);;
    dotlabel.top("$p_6$", p6) shifted (0, 2);
    dotlabel.lft("$p_7$", p7) shifted (-2, 0);
  fi;
  
  for i = 0 upto 20:
    a[i] := mediate(p0, p1, i*1/20);
    s := "$\scriptstyle a_{" & decimal (i) & "}$";
    if do_labels:
      dotlabel.bot(s, a[i]);
    fi;

    a[i+31] := mediate(p3, p2, i*1/20);
    s := "$\scriptstyle a_{" & decimal (i+31) & "}$";
    if do_labels:
      label.top(s, a[i+31]) shifted (0, .25);
      drawdot a[i+31] with_pen pencircle scaled (2mm, 2mm, 2mm);
    fi;
  endfor;

  for i = 0 upto 9:
    a[i+21] := mediate(p0, p3, i*1/10);
    s := "$\scriptstyle a_{" & decimal (i+21) & "}$";

    if do_labels:
      dotlabel.lft(s, a[i+21]);
    fi;


    a[i+52] := mediate(p1, p2, i*1/10);
    s := "$\scriptstyle a_{" & decimal (i+52) & "}$";

    if do_labels:
      label.rt(s, a[i+52]) shifted (.25, 0);
      drawdot a[i+52] with_pen pencircle scaled (2mm, 2mm, 2mm);
    fi;
  endfor;

  q0 := a0 -- a1 -- a22 -- cycle;
  draw q0 on_picture v100;
  
  
  for i = 1 upto 19:
    q[i] := a[i] -- a[i+1] -- a[22+i] -- a[21+i] -- cycle;
    draw q[i] on_picture v100;
  endfor;

  for i = 0 upto 8:
    q[20+i] := a[52+i] -- a[53+i] -- a[42+i] -- a[41+i] -- cycle;
    draw q[20+i] on_picture v100;
  endfor;

  q29 := a61 -- a50 -- a51 -- cycle;

  draw q29 on_picture v100;

  if do_labels:
    label.bot("$o$", origin) shifted (-4pt, -2pt);
    drawdot origin with_pen dot_pen;
  fi;

  for i = 0 upto 29:
    q[30+i] := q[i] rotated (0, 180);
    draw q[30+i] with_color blue on_picture v100;
  endfor;

  v0 := current_picture;

  current_picture += v100;
  
endfig with_projection parallel_x_y;

%% * (1) Red

n := 1/27;

for i = 1 upto 30:
  beginfig(i);
    current_picture := v0;
    for j = (i-1) downto 0:
      C0 := (.3+n*j, 0, 0);
      filldraw q[j] with_color C0 with_pen big_pen;
    endfor;
    if i == 30:
      v1 := current_picture;
    fi;
  endfig with_projection parallel_x_y;
endfor;


%% * (1) Blue

for i = 1 upto 30:
  m := 30 - i;
  beginfig(i+30);
    current_picture := v1;
    j := 0;
    for k = 29 downto m:
      C0 := (0, 0, .3+n*j);
      filldraw q[k] with_color C0 with_pen big_pen;
      j += 1;
    endfor;
    if i == 30:
      v2 := current_picture;
    fi;
  endfig with_projection parallel_x_y;
endfor;


%% * (1) Green

for i = 1 upto 30:
  m := 60 - i;
  beginfig(i+60);
    current_picture := v2;
    j := 0;
    for k = 59 downto m:
      C0 := (0, .3+n*j, 0);
      filldraw q[k] with_color C0 with_pen big_pen;
      j += 1;
    endfor;
    if i == 30:
      v3 := current_picture;
    fi;
  endfig with_projection parallel_x_y;
endfor;

%% * (1)

color_vector cv;

for i = 0 upto 30:
  C0 := (.3+i*n, .3+i*n, 0);
  cv += C0;
endfor;
    
%% * (1) Yellow

for i = 1 upto 30:
  beginfig(i+90);
    current_picture := v3;

    for j = 0 upto (i-1):
      filldraw q[30+j] with_color cv[j] with_pen big_pen;
    endfor;
    if i == 30:
      v4 := current_picture;
    fi;
    
  endfig with_projection parallel_x_y;
endfor;

%% * (1) Magenta

for i = 1 upto 30:
  beginfig(i+120);
    current_picture := v4;
    for j = (i-1) downto 0:
      C0 := (.3+n*j, 0, .3+n*j);
      filldraw q[j] with_color C0 with_pen big_pen;
    endfor;
    if i == 30:
      v5 := current_picture;
    fi;
  endfig with_projection parallel_x_y;
endfor;


%% * (1) Cyan

for i = 1 upto 30:
  m := 30 - i;
  beginfig(i+150);
    current_picture := v5;
    j := 0;
    for k = 29 downto m:
      C0 := (0, .3+n*j, .3+n*j);
      filldraw q[k] with_color C0 with_pen big_pen;
      j += 1;
    endfor;
    if i == 30:
      v6 := current_picture;
    fi;
  endfig with_projection parallel_x_y;
endfor;

bye;

%% ** (2) 

%% * (1) Emacs-Lisp code for use in indirect buffers when using the          
%%       GNU Emacs editor.  The local variable list is not evaluated when an 
%%       indirect buffer is visited, so it's necessary to evaluate the       
%%       following s-expression in order to use the facilities normally      
%%       accessed via the local variables list.                              
%%       LDF 2004.02.12.

%% (progn (metafont-mode) (outline-minor-mode t) (setq fill-column 80) (ignore '(  
%% )) (setq outline-regexp "%% [*\f]+"))

%% * (1) Local variables for Emacs.

%% Local Variables:
%% mode:Metapost
%% eval:(outline-minor-mode t)
%% eval:(read-abbrev-file abbrev-file-name)
%% outline-regexp:"%% [*\f]+"
%% End: