From b64e957fff8b5cc22de43b41f0569bb5d02a3fea Mon Sep 17 00:00:00 2001 From: Karl Berry Date: Mon, 17 Aug 2009 23:53:54 +0000 Subject: new metapost package garrigues (11aug09) git-svn-id: svn://tug.org/texlive/trunk@14742 c570f23f-e606-0410-a88d-b1316a301751 --- Master/texmf-dist/metapost/garrigues/garrigues.mp | 1392 +++++++++++++++++++++ 1 file changed, 1392 insertions(+) create mode 100644 Master/texmf-dist/metapost/garrigues/garrigues.mp (limited to 'Master/texmf-dist/metapost') diff --git a/Master/texmf-dist/metapost/garrigues/garrigues.mp b/Master/texmf-dist/metapost/garrigues/garrigues.mp new file mode 100644 index 00000000000..babeb87841e --- /dev/null +++ b/Master/texmf-dist/metapost/garrigues/garrigues.mp @@ -0,0 +1,1392 @@ +%-*- coding: utf-8 -*- + +% Damien Garrigues' ``nomogramme,'' reproduced from: +% +% @Article{garrigues:chronometrie1939, +% author = "Damien Garrigues", +% title = "{Généralisation de la formule pascale de Gauss: +% Nomogramme du Calendrier perpétuel}", +% journal = "Annales françaises de chronométrie", +% year = "1939", +% date = "1er trimestre 1939", +% volume = "", +% number = "", +% pages = "47-60", +% month = "", +% language = "french", +% } +% +% +% Denis Roegel +% 10-11 June 2001: initial version (complete, but not meant to be +% described in an article) +% +% 16 October 2008: improved version using latexmp +% 19 October 2008: some changes to make it more logical +% (in particular, col4 has been renumbered from the bottom) +% 20-24 October 2008: more changes and restructuration of the code, +% improvements in labels and braces +% 3-7 November 2008: addition of the special case of epact 25 which +% had previously been forgotten (footnote (a) in +% the Easter table area) +% +% +% +% Only part of the ``nomogramme'' is reproduced, +% and this version is not completely identical to the original. +% Some of the idiosyncrasies of the original drawing are difficult +% to reproduce algorithmically and the purpose is not to create an +% exact copy. +% +% It should also be noted that the code is much more complex than +% it needs to be for the pure reproduction of Garrigues' drawing, +% because we wanted in addition to show how the drawing is constructed, +% and provide a number of partial views. As a consequence, the code +% has been broken into many different pieces, and some related pieces +% had to be put in different macros. +% +% For detailed explanations, see the accompanying article. + +input latexmp % for flexible latex labels +input rboxes % for round boxes in column 10 + +numeric u; +u=1mm; + +labeloffset:=7bp; + +numeric th_ex; +th_ex=.5pt; % thickness of example dashed lines + +% label shifts in column 8: +numeric col_shift_eight_a,col_shift_eight_b; +col_shift_eight_a=6u; +col_shift_eight_b=10u; + +def whateverpair=(whatever,whatever) enddef; + +% this macro formats the epact so that it fits on two characters +% (using a phantom 0 on the left) +def ep_st(expr i)= + if i=0: textext("\phantom{0}$\star$") + elseif i<10: + textext("\phantom{0}"&decimal(i)) + else: + textext(decimal(i)) + fi +enddef; + +% this macro formats a number on two digits, adding a phantom 0 when +% necessary +def gstring(expr i)= + if i<10: + textext("\vphantom{(}\phantom{0}"&decimal(i)) + else: + textext("\vphantom{(}"&decimal(i)) + fi +enddef; + +def gpstring(expr i)= + if i<10: + textext("(\phantom{0}"&decimal(i)&")") + else: + textext("("&decimal(i)&")") + fi +enddef; + +% this macro adds a label within the Easter area at position (x,y) +def label_easter_slot(expr x,y,lab)= + label(lab,.5[s[x]r[y],s[x+1]l[y+1]]); +enddef; + +% transform 1->A, 2->B, 3->C, 4->D, 5->E, 6->F, 7->G +% for the correct display of the dominical letter +def DL(expr i)= + char(64+i) +enddef; + +% i is the point number in column 2 (1 at the bottom), +% and this macro returns a pair: +% a) the Golden number associated to this point +% b) the Julian epact associated to this point +vardef gn_epact(expr i)= + save G,JE; + G=1+((9-11i) mod 19); + JE=(11G-3) mod 30; + (G,JE) % value returned +enddef; + +% i is the point number in column 2 (1 at the bottom), +% and this macro returns a pair: +% a) the Golden number associated to this point +% b) the corresponding point in column 4 (1 at the bottom) +% (col. 4 = Julian epact) +vardef gn_epactl(expr i)= + save G,JE,JEL; + G=1+((9-11i) mod 19); + JE=(11G-3) mod 30; + JEL=30-((JE+12) mod 30); + (G,JEL) % value returned +enddef; + +vardef sval(text sec)(text f)= + save s;string s; + % we build a string with space-separated values: + for $=sec: + if unknown s: + s=f $; + else: + s:=s & " " & f $; + fi; + endfor; + s +enddef; + +% draw the values of the Gregorian secular parts in column 8; +% each call corresponds to one line +% i=circle from the bottom (i=1 to 4) +% j=line within the circle (j=1 at the top) +vardef secular_year(expr i,j)(text sec)= + save vd; + vd=4u; % vertical shift of the first line + label(textext(sval(sec)(decimal)) scaled .5, + col8[2i-1]+(10u,vd-(j-1)*.5vd)); +enddef; + +def tddec expr i= + if i<10: "0" & decimal(i) + else: + decimal(i) + fi +enddef; + +% i=box number, from 1 (at the bottom) to 7 (at the top) +% j=1 or 3 (1 for the line at the top of a box, 3 for the line at the bottom) +vardef yn(expr i,j)(text y)= + save vd; + vd=2u; % vertical shift of the first line + label(textext(sval(y)(tddec)) scaled .5, + col10[i]+(0,vd-(j-1)*vd)); +enddef; + +% these numbers always appear at a fixed position: +def yn_left(expr i)(text y)= + label.rt(textext(sval(y)(tddec)) scaled .5,col10[i]+(-10u,0)); +enddef; + +% these numbers always appear at a fixed position: +def yn_right(expr i)(text y)= + label.lft(textext(sval(y)(tddec)) scaled .5,col10[i]+(10u,0)); +enddef; + +% list is a string with comma-separated values +% l is the line number (1 at the bottom) +vardef col_six_f(expr list,l)= + save n,i;n=0; + for $=scantokens(list): + n:=n+1; + endfor; + % n is now the number of values in the list + i=0; + for $=scantokens(list): + i:=i+1; + if n>1: + % if there is more than 1 value, the extreme values + % are put at 2u below and above the center, + % and the other values (if any) are spread evenly + % in-between + label(textext(decimal $) scaled .7, + col6a[l]+(0,-2u+(i-1)*(4u/(n-1)))); + else: + % if there is only one value, it is centered; + % there is only one such case + label(textext(decimal $) scaled .7,col6a[l]); + fi; + endfor; +enddef; + +vardef col_one_three_f(expr list,l,c)= + save n,i;n=0; + for $=scantokens(list): + n:=n+1; + endfor; + i=0; + for $=scantokens(list): + i:=i+1; + label(textext(if c=3: (tddec $) else: decimal $ fi) + scaled .6, + col[c][l]+((if c=3: 2.5u else: 2u fi,0) + rotated (180-(i-1)*360/n))); + endfor; +enddef; + +% variables for the regular columns and circle radiuses +pair col[][],col[]a[]; +numeric height;height=150u; % vertical dimension of the drawing +numeric diam[]; + +% variables for labels +string col[][]st; + +% variables for the Easter table: +pair B,C; +pair s[]r[],s[]l[]; + +% define the points in columns 1-7 +vardef define_first_cols= + save E; + diam1=(height-18*.5)/18; + for i:=1 upto 19: + col1[i]=(0,(i-1)*height/18); + endfor; + for i:=1 upto 19: + col3[i]=(40u,(i-1)*height/18); + endfor; + for i:=1 upto 37: + col2[i]=((xpart(col1[1])+xpart(col3[1]))/2, + .5*(i-1)*height/18); + endfor; + for i:=1 upto 37: + E:=ypart(gn_epactl(i)); + col4[E]=(60u,(E-1)*height/29); + endfor; + for i:=1 upto 30: + col6[i]=(110u,(i-1)*height/29); + endfor; + for i:=1 upto 58: + col5[i]=((xpart(col4[1])+xpart(col6[1]))/2,.5*(i-1)*height/29); + endfor; + for i:=1 upto 30: + % the circles for column 6 appear alternately + % on the right and on the left of the center corresponding to col6 + if i mod 2=0: + col6a[i]=col6[i]-(8u,0); + else: + col6a[i]=col6[i]+(8u,0); + fi; + endfor; + for i:=1 upto 30: + % the values in column 7 are shifted vertically by 20u: + col7[i]=(130u,20u+(i-1)*(height-20u)/29); + endfor; +enddef; + +% define strings containing the numbers in the circles of column 1 +vardef define_col_one_labels= + save p; + for i=0 upto 84: % only S=0 to 84 are considered + p:=1+(9+3i) mod 19; + if unknown col1[p]st: + col1[p]st=decimal(i); + else: + col1[p]st + :=col1[p]st & "," & decimal(i); + fi; + endfor; +enddef; + +% define strings containing the numbers in the circles of column 3 +vardef define_col_three_labels= + save p; + for i=0 upto 99: + p:=1+(15+12i) mod 19; + if unknown col3[p]st: + col3[p]st=decimal(i); + else: + col3[p]st + :=col3[p]st & "," & decimal(i); + fi; + endfor; +enddef; + +% define strings containing the numbers in the circles of column 6 +vardef define_col_six_labels= + save p,M; + for i=15 upto 84: % all secular parts of column 6 + % find the line number where it goes + M:=(15+i-floor(i/4)-floor((8i+13)/25)) mod 30; + p:=1+(M+21) mod 30; + if unknown col6[p]st: + col6[p]st=decimal(i); + else: + col6[p]st:=col6[p]st & "," & decimal(i); + fi; + endfor; +enddef; + +def draw_col_one_points= + for i:=1 upto 19: + drawdot col1[i]; + endfor; +enddef; + +def draw_col_two_points= + for i:=1 upto 37: + drawdot col2[i]; + endfor; +enddef; + +def draw_col_three_points= + for i:=1 upto 19: + drawdot col3[i]; + endfor; +enddef; + +vardef draw_col_four_points= + save GNE;pair GNE; + for i:=1 upto 37: + GNE:=whateverpair;GNE=gn_epactl(i); + drawdot col4[ypart(GNE)]; + endfor; +enddef; + +def draw_col_five_points= + for i:=1 upto 58: + drawdot col5[i]; + endfor; +enddef; + +def draw_col_six_points= + for i:=1 upto 30: + drawdot col6[i]; + endfor; +enddef; + +def draw_col_seven_points= + for i:=1 upto 30: + drawdot col7[i]; + endfor; +enddef; + +def draw_col_eight_points= + for i=1 upto 7: + drawdot col8[i]; + endfor; +enddef; + +def draw_col_nine_points= + for i=1 upto 13: + drawdot col9[i]; + endfor; +enddef; + +def draw_col_ten_points= + for i=1 upto 7: + drawdot col10[i]; + endfor; +enddef; + +def draw_col_eleven_points= + for i=1 upto 7: + drawdot col11[i]; + endfor; +enddef; + +vardef draw_block_one_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_col_one_points; + draw_col_two_points; + draw_col_three_points; + pickup oldpen; +enddef; + +vardef draw_block_two_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_col_two_points; + %draw_col_three_points; + draw_col_four_points; + pickup oldpen; +enddef; + +vardef draw_block_three_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_col_four_points; + draw_col_five_points; + draw_col_six_points; + pickup oldpen; +enddef; + +vardef draw_block_four_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_col_five_points; + %draw_col_six_points; + draw_col_seven_points; + pickup oldpen; +enddef; + +vardef draw_block_five_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_col_eight_points; + draw_col_nine_points; + draw_col_ten_points; + pickup oldpen; +enddef; + +vardef draw_block_six_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + %draw_col_eight_points; + draw_col_nine_points; + draw_col_eleven_points; + pickup oldpen; +enddef; + +vardef draw_block_seven_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_col_seven_points; + draw_col_eleven_points; + pickup oldpen; +enddef; + +% draw the various points in the columns: +vardef draw_all_points= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw_block_one_points; + draw_block_three_points; + draw_block_five_points; + draw_block_seven_points; + pickup oldpen; +enddef; + +vardef draw_block_two_connections= + for i:=1 upto 37: + draw col2[i]--col4[ypart(gn_epactl(i))]; + endfor; +enddef; + +def draw_block_four_connections= + for i:=1 upto 58: + draw col5[i]--col7[1+((i-1) mod 30)]; + endfor; +enddef; + +def draw_block_six_connections= + for i=1 upto 13: + draw col9[i]--col11[1+(13-i) mod 7]; + endfor; +enddef; + +% draw the lines connecting the points +def draw_all_connections= + draw_block_two_connections; + draw_block_four_connections; + draw_block_six_connections; +enddef; + +% connecting the secular part to the last two digits of the year: +vardef draw_block_one_example= + save oldpen; + oldpen=savepen; + draw col1[10]--col3[9] withcolor white; % erasing the old connection + pickup pencircle scaled th_ex; + draw col1[10]--col3[9] dashed evenly; + pickup oldpen; +enddef; + +% connecting the Golden number to the Julian epact: +vardef draw_block_two_example= + save oldpen; + oldpen=savepen; + draw col2[18]--col4[29] withcolor white; % erasing the old connection + pickup pencircle scaled th_ex; + draw col2[18]--col4[29] dashed evenly; + pickup oldpen; +enddef; + +% connecting the Julian epact to M: +vardef draw_block_three_example= + save oldpen; + oldpen=savepen; + draw col4[29]--col6[16] withcolor white; % erasing the old connection + pickup pencircle scaled th_ex; + draw col4[29]--col6[16] dashed evenly; + pickup oldpen; +enddef; + +% connecting the Gregorian epacts: +vardef draw_block_four_example= + save oldpen; + oldpen=savepen; + draw col5[44]--col7[14] withcolor white; % erasing the old connection + pickup pencircle scaled th_ex; + draw col5[44]--col7[14] dashed evenly; + pickup oldpen; +enddef; + +% connecting the secular part of the year to the +% last two digits, in order to obtain the dominical letter: +vardef draw_block_five_example= + save oldpen; + oldpen=savepen; + draw col8[7]--col10[5] withcolor white; % erasing the old connection + pickup pencircle scaled th_ex; + draw col8[7]--col10[5] dashed evenly; + pickup oldpen; +enddef; + +% connecting the dominical letters together: +vardef draw_block_six_example= + save oldpen; + oldpen=savepen; + draw col9[11]--col11[3] withcolor white; % erasing the old connection + pickup pencircle scaled th_ex; + draw col9[11]--col11[3] dashed evenly; + pickup oldpen; +enddef; + +vardef draw_block_seven_example= + save oldpen; + oldpen=savepen; + pickup pencircle scaled th_ex; + % connecting all Easter dates with the same dominical letter A: + draw col11[3]--C dashed evenly; + % connecting all Easter dates with the same Gregorian epact value: + draw col7[14]--B dashed evenly; + % showing the Easter date: + pickup pencircle scaled 3pt; + draw (col11[3]--C) intersectionpoint (col7[14]--B); + pickup oldpen; +enddef; + +% example for 1939: +def draw_all_example= + draw_block_one_example; + draw_block_two_example; + draw_block_three_example; + draw_block_four_example; + draw_block_five_example; + draw_block_six_example; + draw_block_seven_example; +enddef; + +def draw_dbl_circle(expr da,db,c) text options= + unfill fullcircle scaled db shifted c; + draw fullcircle scaled db shifted c; + draw fullcircle scaled da shifted c options; +enddef; + +% circles in column 1: +def draw_col_one_circles= + for i:=1 upto 19: + draw_dbl_circle(.9diam1,diam1,col1[i]); + endfor; +enddef; + +% extra labels for block one: +def draw_col_one_circle_labels= + for i:=1 upto 19: + label.lft(textext(decimal((4+8*(i-1)) mod 19)),col1[i]+7u*left); + endfor; +enddef; + +% circles in column 3: +def draw_col_three_circles= + for i:=1 upto 19: + % this produces a simple circle: + draw_dbl_circle(diam1,diam1,col3[i]); + endfor; +enddef; + +% extra labels for block one: +def draw_col_three_circle_labels= + for i:=1 upto 19: + label.lft(textext(decimal((13+8*(i-1)) mod 19)),col3[i]+13u*right); + endfor; +enddef; + +def draw_col_six_circles= + % circles around column 6: + for i:=1 upto 30: + draw_dbl_circle(.9diam1,diam1,col6a[i]); + endfor; + % special case of the Julian calendar: + draw_dbl_circle(.9diam1,diam1,col6[7]-(17u,0)) dashed evenly; + drawarrow (col6[7]-(17u,0))--col6[7] + cutbefore (fullcircle scaled diam1 shifted (col6[7]-(17u,0))); + label(textext("\vbox{\halign{\hfil#\hfil\cr Jul.\cr Cal.\cr}}") scaled .7, + col6[7]-(17u,0)); +enddef; + +% circles around column 8: +def draw_col_eight_circles= + for i=1 upto 7: + draw_dbl_circle(.9diam2,diam2,col8[i]-(col_shift_eight_a,0)); + if i mod 2=1: + draw_dbl_circle(.9diam3,diam3,col8[i]+(col_shift_eight_b,0)); + fi; + endfor; +enddef; + +def draw_block_one_circles= + draw_col_one_circles; + draw_col_one_circle_labels; + draw_col_three_circles; + draw_col_three_circle_labels; +enddef; + +def draw_block_two_circles= + draw_col_three_circles; +enddef; + +def draw_block_three_circles= + draw_col_six_circles; +enddef; + +def draw_block_four_circles= + draw_col_six_circles; +enddef; + +def draw_block_five_circles= + draw_col_eight_circles; +enddef; + +def draw_block_six_circles= + draw_col_eight_circles; +enddef; + +def draw_all_circles= + draw_col_one_circles; + draw_col_three_circles; + draw_col_six_circles; + draw_col_eight_circles; +enddef; + +% display values in circles of column 1: +def draw_labels_col_one= + for i=1 upto 19: + col_one_three_f(col1[i]st,i,1); + endfor; + label(textext("I"),col1[1]-(0,10u)); + label(textext("$S$") scaled 1.5,col1[19]+(col1[19]-col1[18])); +enddef; + +def draw_labels_col_two_a= + label(textext("II"),col2[1]-(0,10u)); +enddef; + +% this is an extended display +vardef draw_labels_col_two_b= + save GNE;pair GNE; + for i:=1 upto 37: + GNE:=whateverpair;GNE=gn_epact(i); + % some of the values here are written twice, but this causes no harm: + label.rt(gstring(xpart(GNE)),col2[i]); + label.lft(gstring(ypart(GNE)),col2[i]); + endfor; + label.lft(textext("$E_J$"),col2[37]+(col2[37]-col2[36])); + label.rt(textext("$G$"),col2[37]+(col2[37]-col2[36])); +enddef; + +% this is an extended display +vardef draw_labels_col_two_c= + save GNE;pair GNE; + for i:=1 upto 37: + GNE:=whateverpair;GNE=gn_epact(i); + % some of the values here are written twice, but this causes no harm: + label.lft(gstring(xpart(GNE)),col2[i]+15u*left); + label.lft(gpstring(ypart(GNE)),col2[i]+5u*left); + endfor; + label.lft(textext("$G$"),col2[37]+(col2[37]-col2[36])+15u*left); + label.lft(textext("$E_J$"),col2[37]+(col2[37]-col2[36])+5u*left); +enddef; + +% display values in circles of column 3: +def draw_labels_col_three= + for i=1 upto 19: + col_one_three_f(col3[i]st,i,3); + endfor; + label(textext("III"),col3[1]-(0,10u)); + label(textext("$A$") scaled 1.5,col3[19]+(col3[19]-col3[18])); +enddef; + +vardef draw_labels_col_four_base= + save GNE;pair GNE; + for i:=1 upto 37: + GNE:=whateverpair;GNE=gn_epactl(i); + % some of the values here are written twice, but this causes no harm: + label.rt(gstring(xpart(GNE)),col4[ypart(GNE)]); + endfor; +enddef; + +% display labels in column 4 (values of the Golden number) +def draw_labels_col_four_a= + draw_labels_col_four_base; + label.lft(textext("IV"),col4[1]-(0,10u)); + label.rt(textext("Gold. N.") rotated 90,col4[1]-(0,10u)); +enddef; + +def draw_labels_col_four_b= + draw_labels_col_four_base; + label(textext("IV"),col4[1]-(0,10u)); +enddef; + +vardef draw_labels_col_four_c= + save GNE,GNEL;pair GNE,GNEL; + for i:=1 upto 37: + GNEL:=whateverpair;GNEL=gn_epactl(i); + GNE:=whateverpair;GNE=gn_epact(i); + label.rt(gpstring(ypart(GNE)),col4[ypart(GNEL)]+5u*right); + endfor; + label.rt(textext("$G$"),col4[29]+(col4[29]-col4[28])); + label.rt(textext("$E_J$"),col4[29]+(col4[29]-col4[28])+5u*right); +enddef; + +vardef draw_labels_col_four_d= + save GNE,GNEL;pair GNE,GNEL; + for i:=1 upto 37: + GNEL:=whateverpair;GNEL=gn_epactl(i); + GNE:=whateverpair;GNE=gn_epact(i); + label.lft(gstring(xpart(GNEL)),col4[ypart(GNEL)]+8u*left); + label.lft(gpstring(ypart(GNE)),col4[ypart(GNEL)]); + endfor; + label.lft(textext("$G$"),col4[29]+(col4[29]-col4[28])+8u*left); + label.lft(textext("$E_J$"),col4[29]+(col4[29]-col4[28])); + label.lft(textext("IV"),col4[1]-(0,10u)); +enddef; + +def draw_labels_col_five= + for i=1 upto 58: + label.lft(ep_st((24-i) mod 30) scaled .7,col5[i]); + endfor; + label.lft(textext("$E_G$"),col5[58]+2(col5[58]-col5[57])); + label(textext("V"),col5[1]-(0,10u)); +enddef; + +def draw_labels_col_five_a= + label(textext("V"),col5[1]-(0,10u)); +enddef; + +% display values in circles of column 6: +def draw_labels_col_six_a= + for i=1 upto 30: + col_six_f(col6[i]st,i); + endfor; + label(textext("VI"),col6[1]-(0,10u)); + label(textext("$S$") scaled 1.5,col6[30]+2(col6[30]-col6[29])); +enddef; + +vardef draw_labels_col_six_b= + save lo; + lo=labeloffset; + labeloffset:=3bp; + for i=1 upto 30: + label.rt(gstring((7-i) mod 30) scaled .5,col6[i]); + endfor; + labeloffset:=lo; + label(textext("$15-M$") scaled .7,col6[30]+.5(col6[30]-col6[29])+2u*right); +enddef; + +% display values of the Gregorian epact +def draw_labels_col_seven_a= + for i:=1 upto 30: + label.rt(ep_st((24-i) mod 30),col7[i]); + endfor; + label(textext("VII"),col7[1]-(0,20u)); + label.rt(textext("Epact:") rotated 90,col7[1]-(0,10u)); +enddef; + +def draw_labels_col_seven_b= + draw_labels_col_seven_a; + label.rt(textext("$E_G$"),col7[30]+(col7[30]-col7[29])); +enddef; + +def draw_labels_col_seven_c= + for i:=1 upto 30: + label.lft(ep_st((24-i) mod 30),col7[i]); + endfor; + label.lft(textext("$E_G$"),col7[30]+2(col7[30]-col7[29])); + label(textext("VII"),col7[1]-(0,20u)); + label.lft(textext("Epact:") rotated 90,col7[1]-(0,10u)); +enddef; + +def draw_block_one_labels= + draw_labels_col_one; + draw_labels_col_two_a; + draw_labels_col_three; +enddef; + +def draw_block_two_labels= + draw_labels_col_two_a; + draw_labels_col_four_b; +enddef; + +def draw_block_three_labels= + draw_labels_col_four_d; + draw_labels_col_five; + draw_labels_col_six_a; + draw_labels_col_six_b; +enddef; + +def draw_block_four_labels= + draw_labels_col_five; + draw_labels_col_seven_b; +enddef; + +% draw the labels in columns 1-7 +def draw_first_labels= + draw_labels_col_one; + draw_labels_col_three; + draw_labels_col_four_a; + draw_labels_col_five_a; + draw_labels_col_six_a; + draw_labels_col_seven_a; +enddef; + +def define_corners= + corner1=whatever[col7[0],B] + =C+whatever*up; + corner3=.3[B,col7[31]]; % can be positioned freely + corner2=(C--corner3) intersectionpoint (B--corner1); + corner4=whatever[B,col7[31]] + =C+whatever*up; +enddef; + +vardef define_easter_table= + save corner,p; + pair corner[]; + % B and C can be positioned freely: + C=(xpart(col7[1])+10u,-10u); + B=(xpart(C)+150u,ypart(col7[5])); + % we define two additional points in column 7, + % one above the 30th, and one below the first: + col7[1]-col7[0]=col7[31]-col7[30]=col7[2]-col7[1]; + % these four corners define the shape of area XI: + define_corners; + % area XI is divided into 8 slices (9 boundaries, + % only 8 of which are drawn) + for i=1 upto 9: + s[i]l0=s[i]r0=(((i-1)/8))[corner4,corner3]; + s[i]l5=s[i]r5 + =whatever[corner1,corner2]=whatever[s[i]l0,C]; + endfor; + % each of the eight slices is itself divided four times: + for i=1 upto 8: % i varies on the ``vertical lines,'' + % from left to right; 8 vertical lines enclose the + % 35 easter slots; + for j=1 upto 4: % j varies on the horizontal inner divisions + p:=30-i-(j-1)*7; + % from top to bottom + if i<8: + % division leaving to the right of vertical line i + s[i]r[j]=(s[i]l0--s[i]l5) + intersectionpoint (B--.5[col7[p],col7[p-1]]); + fi; + if i>1: + % division leaving to the left of vertical line i + s[i]l[j]=(s[i]l0--s[i]l5) + intersectionpoint + (B--.5[col7[p+1],col7[p]]); + fi; + endfor; + endfor; +enddef; + +% draw points B and C +vardef draw_B_and_C= + save oldpen; + oldpen=savepen; + draw_dbl_circle(.8*3u,3u,B); + draw_dbl_circle(.8*3u,3u,C); + pickup pencircle scaled 2pt; + drawdot B;drawdot C; + pickup oldpen; + label.ulft(textext("$B$"),B); + label.ulft(textext("$C$"),C); +enddef; + +% draw the slices in the Easter table: +vardef draw_easter_table_slices= + save oldpen; + oldpen=savepen; + for i=1 upto 8: + draw s[i]l0--s[i]l5; + endfor; + draw s8l5--s8l0--s1l0--s1l5--cycle; + for i=1 upto 7: + for j=1 upto 4: + draw s[i]r[j]--s[i+1]l[j]; + endfor; + endfor; + pickup pencircle scaled 2pt; + % March/April divisions: + draw s8l3--s7r3--s7l3--s6r3--s6l3--s5r3-- + s5l4--s4r4--s4l4--s3r4--s3l4-- + s2r4--s2l4--s1r4; + pickup oldpen; +enddef; + +vardef draw_arrow_label@#(expr lab,c,s)= + save oldpen; + oldpen=savepen; + pickup pencircle scaled 2pt; + draw (c); + pickup oldpen; + drawarrow (c+s)--(c+u*unitvector(s)); + label@#(lab,c+s); +enddef; + + +vardef draw_easter_table_details= + save corner; + pair corner[]; + define_corners; + draw_arrow_label.top(textext("\tt corner4"),corner4,(0,10u)); + draw_arrow_label.urt(textext("\tt corner3"),corner3,10u*dir(20)); + draw_arrow_label.lrt(textext("\tt corner2"),corner2,10u*dir(-20)); + draw_arrow_label.llft(textext("\tt corner1"),corner1,25u*dir(-132)); + + draw_arrow_label.top(textext("\tt s2l0"),s[2]l[0],30u*dir(80)); + draw_arrow_label.rt(textext("\tt s2r0"),s[2]r[0],40u*dir(10)); + draw_arrow_label.top(textext("\tt s4l0"),s[4]l[0],12u*dir(90)); + draw_arrow_label.rt(textext("\tt s4r0"),s[4]r[0],20u*dir(5)); + draw_arrow_label.lrt(textext("\tt s7l3"),s[7]l[3],40u*dir(-10)); + draw_arrow_label.lrt(textext("\tt s7r3"),s[7]r[3],40u*dir(-40)); + draw_arrow_label.rt(textext("\tt s7r2"),s[7]r[2],60u*dir(5)); + + draw corner1--corner2 dashed evenly; + draw corner4--corner3 dashed evenly; + draw corner2--corner3 dashed evenly; +enddef; + +% draw the labels in the Easter table: +vardef draw_easter_table_labels= + save laban,march,april,sl,j; + string march,april; + for i=1 upto 35: % 35 dates from March 22 till April 25 + sl:=1+(7-(i mod 7)) mod 7; + j:=4-floor((i-1)/7); + label_easter_slot(sl,j,textext(if i=35:"25 (a)" else: + decimal(if i>10:i-10 else: i+21 fi) fi) + if i<8: scaled .7 fi); + endfor; + laban=angle(s8l0-s8l5); + march="$\underbrace{\kern" + & decimal(arclength(s8l3--s8l5)-5) + & "bp}_{\hbox{MARCH}}$"; + april="$\underbrace{\kern" + & decimal(arclength(s8l0--s8l3)-5) + & "bp}_{\hbox{APRIL}}$"; + label(textext(march) rotated laban, + .5[s8l3,s8l5]+3u*unitvector((s8l0-s8l5) rotated -90)); + label(textext(april) rotated laban, + .5[s8l0,s8l3]+3u*unitvector((s8l0-s8l5) rotated -90)); +enddef; + +def draw_easter_table= + draw_easter_table_slices; + draw_B_and_C; + draw_easter_table_labels; +enddef; + +vardef define_dominical_letters= + save shift; + pair shift; + shift=(3u,3u); + for i=1 upto 8: + col11[i]=whatever[C,.5[s[i]r0,s[i+1]l0]] + =whatever[s1r0+shift,s8l0+shift]; + endfor; +enddef; + +def draw_labels_col_eleven= + for i=1 upto 7: + label.ulft(textext(DL(1+(10-i) mod 7)),col11[i]); + endfor; + label(textext("XI") rotated (angle(s1r0-s8l0)-90), + .6[col11[8],col11[7]]); +enddef; + +% define columns 8-11 +def define_last_cols= + diam2=7u; + diam3=13u; + for i=1 upto 7: + col8[i]=s8l0+(15u,10u+(i-1)*ypart(s1l0-s8l0)/7); + col10[i]=col8[i]+(50u,0); + endfor; + for i=1 upto 13: + col9[i]=(xpart(col8[1]+col10[1])/2, + ypart(col8[1])+(i-1)*(ypart(col8[7]-col8[1]))/12); + endfor; + define_dominical_letters; +enddef; + +vardef draw_labels_col_eight_julian= + save v; + for i=1 upto 7: + for j=1 upto 3: + v:=((4+i) mod 7)+(j-1)*7; + label(textext(decimal(v)) scaled .5, + col8[i]-(col_shift_eight_a,0) + +((0,1.4u) rotated ((j-1)*120))); + endfor; + endfor; +enddef; + +def draw_labels_col_eight_gregorian= + % the following lines could be parameterized, but it's not worth it + secular_year(1,1)(16,20,24); + secular_year(1,2)(28,32,36,40); + secular_year(1,3)(44,48,52,56); + secular_year(1,4)(60,64,68,72); + secular_year(1,5)(76,80,84); + + secular_year(2,1)(17,21); + secular_year(2,2)(25,29,33,37); + secular_year(2,3)(41,45,49,53); + secular_year(2,4)(57,61,65,69); + secular_year(2,5)(73,77,81); + + secular_year(3,1)(18,22); + secular_year(3,2)(26,30,34,38); + secular_year(3,3)(42,46,50,54); + secular_year(3,4)(58,62,66,70); + secular_year(3,5)(74,78,82); + + secular_year(4,1)(15,19,23); + secular_year(4,2)(27,31,35,39); + secular_year(4,3)(43,47,51,55); + secular_year(4,4)(59,63,67,71); + secular_year(4,5)(75,79,83); +enddef; + +vardef draw_labels_col_ten= + save rbr; + rbr=rbox_radius; + rbox_radius:=15pt; + for i=1 upto 7: + rboxit.rb[i](""); + rb[i].c=col10[i]; + rb[i].dx=9u;rb[i].dy=3.3u; + unfill bpath(rb[i]); + drawboxes(rb[i]); + endfor; + rbox_radius:=rbr; + + yn(1,1)(3,8,14,25,31); + yn_left(1)(36,42);yn_right(1)(53,59); + yn(1,3)(64,70,81,87,92,98); + + yn(2,1)(2,13,19,24,30); + yn_left(2)(41,47);yn_right(2)(52,58); + yn(2,3)(69,75,80,86,97); + + yn(3,1)(1,7,12,18,29); + yn_left(3)(35,40);yn_right(3)(46,57); + yn(3,3)(63,68,74,85,91,96); + + yn(4,1)(0,6,17,23,28); + yn_left(4)(34,45);yn_right(4)(51,56); + yn(4,3)(62,73,79,84,90); + + yn(5,1)(5,11,16,22,33); + yn_left(5)(39,44);yn_right(5)(50,61); + yn(5,3)(67,72,78,89,95); + + yn(6,1)(4,10,21,27,32); + yn_left(6)(38,49);yn_right(6)(55,60); + yn(6,3)(66,77,83,88,94); + + yn(7,1)(9,15,20,26,37); + yn_left(7)(43,48);yn_right(7)(54,65); + yn(7,3)(71,76,82,93,99); + + label(textext("X"),col10[1]-(0,10u)); + label(textext("$A$") scaled 1.5,col10[7]+2(col8[7]-col8[6])); +enddef; + +def draw_labels_col_ten_a= + draw_labels_col_ten; + for i=1 upto 7: + label.rt(textext(decimal((i+3) mod 7)),col10[i]+10u*right); + endfor; +enddef; + +def draw_labels_col_eight_a= + draw_labels_col_eight_julian; + draw_labels_col_eight_gregorian; + label(textext("VIII"),col8[1]-(0,10u)); + label(textext("J.C.") scaled 1.5, + col8[7]+(-col_shift_eight_a,0)+(col8[7]-col8[6])); + label(textext("G.C.") scaled 1.5, + col8[7]+(col_shift_eight_b,0)+(col8[7]-col8[6])); +enddef; + +def draw_labels_col_eight_b= + for i=1 upto 7: + label.top(decimal(i-1),col8[i]); + endfor; +enddef; + +def draw_labels_col_nine= + for i=1 upto 13: + label.rt(textext(DL(1+(3+i) mod 7)),col9[i]); + endfor; + label(textext("IX"),col9[1]-(0,10u)); + label.rt(textext("Dom. L.") rotated 90,col9[1]-(0,10u)); +enddef; + +def draw_labels_col_nine_a= + for i=1 upto 13: + label.rt(textext(DL(1+(3+i) mod 7) + & " ("&decimal((i+3) mod 7)&")"),col9[i]); + endfor; + label(textext("IX"),col9[1]-(0,10u)); + label.rt(textext("Dom. L.") rotated 90,col9[1]-(0,10u)); +enddef; + +def draw_block_five_labels= + draw_labels_col_eight_a; + draw_labels_col_eight_b; + draw_labels_col_nine_a; + draw_labels_col_ten_a; +enddef; + +def draw_block_six_labels= + draw_labels_col_nine; + draw_labels_col_eleven; +enddef; + +def draw_block_seven_labels= + draw_labels_col_seven_c; + draw_labels_col_eleven; +enddef; + +def draw_all_labels= + draw_first_labels; + draw_labels_col_eight_julian; + draw_labels_col_eight_gregorian; + label(textext("VIII"),col8[1]-(0,10u)); + draw_labels_col_nine; + draw_labels_col_ten; + draw_labels_col_eleven; +enddef; + +def draw_block_one_brace= + label(textext("$\overbrace{\kern85bp}^{\hbox{Year}}$") scaled 1.5, + col2[37]+4(col2[37]-col2[36])); +enddef; + +def draw_col_eight_brace= + label(textext("$\overbrace{\kern50bp}^{\hbox{$S$}}$") scaled 1.5, + col8[7]+2(col8[7]-col8[6])+(.5(col_shift_eight_b-col_shift_eight_a),0)); +enddef; + +def draw_block_five_braces= + draw_col_eight_brace; + label(textext("$\overbrace{\kern100bp}^{\hbox{Year}}$") scaled 1.5, + col9[13]+6(col9[13]-col9[12]) + +(.25(col_shift_eight_b-col_shift_eight_a),0)); +enddef; + +def draw_block_six_brace= + draw_col_eight_brace; +enddef; + +def draw_all_braces= + draw_block_one_brace; + draw_block_five_braces; +enddef; + +def draw_block_one= + draw_block_one_example; + draw_block_one_circles; + draw_block_one_points; % after the circles + draw_block_one_labels; + draw_block_one_brace; + draw_labels_col_two_b; +enddef; + +def draw_block_two= + draw_block_two_connections; + draw_block_two_example; + draw_block_two_points; % after the circles + draw_block_two_labels; + draw_labels_col_four_c; + draw_labels_col_two_c; +enddef; + +def draw_block_three= + draw_block_three_example; + draw_block_three_circles; + draw_block_three_points; % after the circles + draw_block_three_labels; +enddef; + +def draw_block_four= + draw_block_four_connections; + draw_block_four_example; + draw_block_four_points; % after the circles + draw_block_four_labels; +enddef; + +def draw_block_five= + draw_block_five_example; + draw_block_five_circles; + draw_block_five_labels; + draw_block_five_points; % must be after the labels + draw_block_five_braces; +enddef; + +def draw_block_six= + draw_block_six_connections; + draw_block_six_example; + draw_block_six_points; % after the example + draw_block_six_labels; +enddef; + +def draw_block_seven= + draw_block_seven_example; + draw_easter_table; % after the example + draw_block_seven_labels; + draw_block_seven_points; % after the example + draw_easter_table_details; +enddef; + +define_first_cols; +define_easter_table; +define_last_cols; +define_col_six_labels; +define_col_one_labels; +define_col_three_labels; + +beginfig(1); + draw_all_connections; + draw_all_example; + draw_easter_table; + label(textext("\footnotesize "& + "\parbox{4cm}{\raggedright "& + "Note (a). The Gregorian years which have " & + "epact $=$ 25, dominical letter $=$ C, and " & + "Golden Number $>$ 11 have their Easter on " & + "April 18 instead of April 25. The only such " & + "years are 1954, 2049, 2106, 3165, 3260, 3317, "& + "3852, 3909, 4004, 6399, etc." & + "}"),C+35u*right); + + draw_all_circles; + draw_all_labels; + draw_all_points; % after the circles and the labels + draw_all_braces; +endfig; + +beginfig(11); + draw_block_one; +endfig; + +beginfig(12); + draw_block_two; +endfig; + +beginfig(13); + draw_block_three; +endfig; + +beginfig(14); + draw_block_four; +endfig; + +beginfig(15); + draw_block_five; +endfig; + +beginfig(16); + draw_block_six; +endfig; + +beginfig(17); + draw_block_seven; +endfig; + +beginfig(21); + numeric v,w,last; + v=10u; + w=20u; + last=5; + pair P[],Q[],R[]; + pickup pencircle scaled 2pt; + for i=0 upto last: + P[i]=(0,i*v);draw P[i]; + Q[i]=(2w,i*v);draw Q[i]; + endfor; + for i=0 upto 2last: + R[i]=(w,i*v/2);draw R[i]; + endfor; + pickup pencircle scaled .5; + + label.bot(textext("$p_i$"),P0); + label.bot(textext("$q_i$"),Q0); + label.bot(textext("$r_i$"),R0); + draw P2--Q3; + for i=0 upto 5: + label.lft(textext(decimal i) scaled .7,P[i]); + endfor; + for i=0 upto 5: + label.rt(textext(decimal i) scaled .7,Q[i]); + endfor; + for i=0 upto 4: + label.rt(textext(decimal i) scaled .7,R[i]); + endfor; + unfill bbox(thelabel.rt(textext("5") scaled .7,R5)); + label.rt(textext("5") scaled .7,R5); + for i=6 upto 10: + label.rt(textext(decimal i) scaled .7,R[i]); + endfor; + +endfig; + +beginfig(25); + numeric v,w,last; + v=10u; + w=20u; + last=5; + pair P[],Q[],R[]; + pickup pencircle scaled 2pt; + for i=0 upto last: + P[i]=(0,i*v);draw P[i]; + Q[i]=(2w,i*v);draw Q[i]; + endfor; + for i=0 upto 2last: + R[i]=(w,i*v/2);draw R[i]; + endfor; + pickup pencircle scaled .5; + + label.bot(textext("$p_i$"),P0); + label.bot(textext("$q_i$"),Q0); + label.bot(textext("$r_i$"),R0); + draw P2--Q3; + for i=0 upto 5: + label.lft(textext("$x_" & decimal i & "$") scaled .7,P[i]); + endfor; + for i=0 upto 5: + label.rt(textext("$y_" & decimal i & "$") scaled .7,Q[i]); + endfor; + for i=0 upto 4: + label.rt(textext("$z_" & decimal i & "$") scaled .7,R[i]); + endfor; + unfill bbox(thelabel.rt(textext("$z_5$") scaled .7,R5)); + label.rt(textext("$z_5$") scaled .7,R5); + for i=6 upto 10: + label.rt(textext("$z_{" & decimal i & "}$") scaled .7,R[i]); + endfor; + +endfig; + + + +end -- cgit v1.2.3