%% sinewaves5.ldf
%% Created by Laurence D. Finston (LDF) Fr 18. Aug 22:12:03 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";

bool do_png;
do_png := false; % true; % 

if do_png:
  verbatim_metapost   "outputformat:=\"png\";outputtemplate := \"%j_%3c.png\";"
                    & "outputformatoptions := \"format=rgb antialias=none\";";
else:
  verbatim_metapost "prologues := 3;outputtemplate := \"%j_%3c.eps\";dotlabeldiam:=2mm;";
fi;




%% * (1)

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

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

%% (* (/ 30.0 16) 9) 16.875


picture v[];

pen huge_pen;

huge_pen := pensquare scaled (2mm, 2mm, 2mm);


%% * (1)

pickup medium_pen;


sinewave s[];

pickup medium_pen;

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

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

set s0 with_amplitude 8 with_frequency 1 with_phase 0 with_increment .1
  with_period 8pi;


r0 := get_rectangle s0 with_scaled (1, 1.025);


% message "r0:";
% show r0;
% pause;

p11 := get_center r0;


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


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

rectangle frame;

n := magnitude(p2 - p1);

frame_ht := n;

frame_wd := n * 16/9;

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

p10 := get_center frame;

shift frame by (p11 - p10);


%% * (1)

beginfig(300);
  draw frame with_pen big_pen;
  s := "\setbox0=\hbox{{\huge Sine Waves Nr.~5}}"
       &  "\vbox{\copy0\vskip2cm"
       & "\hbox to \wd0{\hss{\Large By}\hss}\vskip1cm"
       & "\hbox to \wd0{\hss{\Large Laurence Finston}\hss}\vskip2.5cm"
       & "\hbox to \wd0{\hss{\large Copyright (C) 2023 Laurence Finston}\hss}"
       & "}";
  label(s, p11);
  %dotlabel.top("$p_{10}$", p11);
endfig with_projection parallel_x_y;

beginfig(400);
  draw frame with_pen big_pen;
  s :=   "\hbox{\hfil{\huge The End}\hfil}";
  label(s, p11);
  %dotlabel.top("$p_{11}$", p11);
endfig with_projection parallel_x_y;

%% * (1)

k := 7.5;
for i = 1 upto 15:
  set s[i] with_amplitude k with_frequency 1 with_phase 0 with_increment .1
    with_period 8pi;
  k-=.5;
endfor;

% message "n:";
% show n;
%pause;

for i = 0 upto 6:
  n := length s15;
  n /= 8;
  n *= (i+1);
  n := floor n;
  p[12+i] := get_point (n) s15;
endfor;

%% * (1)

%% This works.  LDF 2023.08.18.

q15 := s15;

n := length s15;
n /= 8;
n := floor n;

q115 := subpath(0, n) of q15;

q14 := s14;

n := length s14;
n /= 8;
n := floor n;

q114 := subpath(0, n) of q14;


q200 := p7 & q115 & p12 & reversed q114 & p7;
q200 += cycle;

% message "q200:";
% show q200;
%pause;

%% * (1)

m := 30;
k := 0;
for i = 15 downto 0:
  beginfig(k);
    if do_png:
      filldraw frame with_pen huge_pen;
      unfill r0;
    fi;
    draw r0 with_pen big_pen;
    draw frame with_pen big_pen;

    if do_labels:
      dotlabel.top("$p_{12}$", p12);
      dotlabel.top("$p_{13}$", p13);
      dotlabel.top("$p_{14}$", p14);
      dotlabel.top("$p_{15}$", p15);
      dotlabel.top("$p_{16}$", p16);
      dotlabel.top("$p_{17}$", p17);
      dotlabel.top("$p_{18}$", p18);
      dotlabel.lft("$p_{7}$", p7);
      dotlabel.rt("$p_{5}$", p5);
    fi;
    for j = 15 downto i:
      draw s[j] with_pen big_pen;
    endfor;

    if false:
      % This works.  LDF 2023.08.18.
      fill q200 with_color dark_green;
      draw q115 with_color red with_pen big_pen;
      draw q114 with_color blue with_pen big_pen;
    fi;
    clear v0;
    v0 := current_picture;

  endfig with_projection parallel_x_y;

  if k <> 15:
    beginfig(m);
      current_picture := v0;
    endfig with_projection parallel_x_y;
    m-= 1;
  fi;
  
  k+=1;

endfor;

j := .5;
for i = 0 upto 15:
  set s[31+i] with_amplitude j with_frequency 1 with_phase pi with_increment .1
    with_period 8pi;
  j+=.5;
endfor;

for i = 0 upto 6:
  n := length s31;
  n /= 8;
  n *= (i+1);
  n := floor n;
  p[19+i] := get_point (n) s31;
endfor;



k := 61;
for i = 0 upto 15:
  %message "i+31: " & decimal (i+31);
  beginfig(31+i);
    if do_png:
      filldraw frame with_pen huge_pen;
      unfill r0;
    fi;
    draw frame with_pen big_pen;
    draw r0 with_pen big_pen;
    draw s[31+i] with_pen big_pen;
    if do_labels:
      dotlabel.rt("$p_5$", p5);
      dotlabel.lft("$p_7$", p7);
      dotlabel.top("$p_{19}$", p19);
      dotlabel.top("$p_{20}$", p20);
      dotlabel.top("$p_{21}$", p21);
      dotlabel.top("$p_{22}$", p22);
      dotlabel.top("$p_{23}$", p23);
      dotlabel.top("$p_{24}$", p24);
      dotlabel.top("$p_{25}$", p25);
    fi;
    n := 31+i;
    for j = n downto 31:
      draw s[j] with_pen big_pen;
    endfor;
    clear v0;
    v0 := current_picture;
  endfig with_projection parallel_x_y;
  if i <> 15:
    %message "k: " & decimal k;
    beginfig(k);
      current_picture := v0;
    endfig with_projection parallel_x_y;
    k -= 1;
  fi;
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:
