Difference between revisions of "File:VorosMap150d.jpg"

Complex map of the Voros function Voros(z)=2z+z^2

is shown with
lines \(u=\Re(\mathrm{Voros}(z)\) and
lines \(v=\Im(\mathrm{Voros}(z)\)
in plane \(z=x\!+\!\mathrm i y\)

C++ generator of curves

Files conto.cin and ado.cin (with lowercase letters) should be loaded in order to compile the code below.

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#define DB double
#define DO(x,y) for(x=0;x<y;x++)
#include <complex>
typedef std::complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
#include "conto.cin"
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
  int M=201,M1=M+1;
  int N=201,N1=N+1;
DB X[M1],Y[N1], g[M1*N1],f[M1*N1], w[M1*N1]; // w is working array.
char v[M1*N1]; // v is working array
FILE *o;o=fopen("VorosMa.eps","w");ado(o,608,608);
fprintf(o,"304 304 translate\n 100 100 scale\n");
DO(m,M1) X[m]=-3.+.03*(m-.5);
DO(n,N1) Y[n]=-3.+.03*(n-.5);

for(m=-3;m<4;m++){if(m==0){M(m,-3.04)L(m,3.04)} else{M(m,-3)L(m,3)}}
for(n=-3;n<4;n++){ M( -3 ,n)L(3,n)}
fprintf(o,".008 W 0 0 0 RGB 2 setlinecap S\n");

// maq(3.);
DO(m,M1)DO(n,N1){g[m*N1+n]=9999; f[m*N1+n]=9999;}
DO(m,M1){x=X[m]; //printf("%5.2f\n",x);
DO(n,N1){y=Y[n]; z=z_type(x,y);
// c=F(.5+E(z));
//        c=-1.8*exp(sqrt(2.)*log(-z));
        c=2.*z+z*z;
        p=Re(c);q=Im(c);
//        if(p>-4.9 && p<4.9) 
        {g[m*N1+n]=p;}
//        if(q>-4.9 && q<4.9 && fabs(q)>1.e-11 ) 
        {f[m*N1+n]=q;}
                        }}

fprintf(o,"1 setlinejoin 2 setlinecap\n");
//p=.8;q=.4;
p=2.;q=.5;
//#include"plof.cin"
for(m=-4;m<4;m++)
for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M,N, (m+.1*n),-q, q);fprintf(o,".005 W 0 .6 0 RGB S\n");
for(m=0;m<4;m++)
for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q, q);fprintf(o,".005 W .9 0 0 RGB S\n");
for(m=0;m<4;m++)
for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q, q);fprintf(o,".005 W 0 0 .9 RGB S\n");

for(m= 1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p); fprintf(o,".011 W .9 0 0 RGB S\n");
for(m= 1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p); fprintf(o,".011 W 0 0 .9 RGB S\n");
for(m=-8;m<9;m++) conto(o,g,w,v,X,Y,M,N, (0.+m),-p,p); fprintf(o,".010 W 0 0 0 RGB S\n");

                  conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".011 W .6 0 .6 RGB S\n");
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
        system("epstopdf VorosMa.eps"); // both linux and mac
        system("convert VorosMa.eps VorosMa.jpg"); // both linux and mac
//        system( "xpdf VorosMa.pdf"); // for linux
//        system( "open VorosMa.pdf"); // for mac
}

Latex generator of labels

File VorosMa.pdf should be generated with code above in order to compile the latex document below:

\documentclass[12pt]{article}
\usepackage{geometry}
\paperwidth 640pt
\paperheight 630pt
\textwidth 700pt
\textheight 700pt
\usepackage{graphics}
\usepackage{rotating}
\pagestyle{empty}
\parindent 0pt
\topmargin -108pt
\oddsidemargin -72pt
\newcommand \ing {\includegraphics}
\newcommand \sx {\scalebox}
\newcommand \rot {\begin{rotate}}
\newcommand \ero {\end{rotate}}
\begin{document}% 
\begin{picture}(638,628)
\put(30,20){\ing{VorosMa.pdf}}
\put(14,615){\sx{2.4}{$y$}}
%\put(14,616){\sx{2.4}{$3$}}
\put(14,516){\sx{2.4}{$2$}}
\put(14,416){\sx{2.4}{$1$}}
\put(14,316){\sx{2.4}{$0$}}
\put(-4,216){\sx{2.4}{$-1$}}
\put(-4,116){\sx{2.4}{$-2$}}
\put(-5,16){\sx{2.4}{$-3$}}
\put(11,0){\sx{2.4}{$-3$}}
\put(111,0){\sx{2.4}{$-2$}}
\put(211,0){\sx{2.4}{$-1$}}
\put(330,0){\sx{2.4}{$0$}}
\put(430,0){\sx{2.4}{$1$}}
\put(530,0){\sx{2.4}{$2$}}
%\put(630,0){\sx{2.4}{$3$}}
\put(624,3){\sx{2.4}{$x$}}

\put(412,562){\rot{37}\sx{2.5}{$u\!=\!-4$}\ero}
\put(460,517){\rot{47}\sx{2.5}{$u\!=\!0$}\ero}
\put(514,482){\rot{54}\sx{2.5}{$u\!=\!4$}\ero}
\put(568,454){\rot{63}\sx{2.5}{$u\!=\!8$}\ero}

\put(566,438){\rot{-18}\sx{2.5}{$v\!=\!8$}\ero}
\put(557,378){\rot{-8}\sx{2.5}{$v\!=\!4$}\ero}
\put(554,317){\sx{2.5}{$v\!=\!0$}}
\put(556,256){\rot{7}\sx{2.5}{$v\!=\!-4$}\ero}
\put(562,195){\rot{16}\sx{2.5}{$v\!=\!-8$}\ero}

\end{picture}
\end{document}

References