Add files via upload

Parameters.md

@@ -9,6 +9,117 @@ function animation. Some simulations use variants of the published code. The lis
 updated gradually. For older simulations, see files `Parameters_July21.md` etc. 
 
 
+### 3 August 21 - Wave protection comparison 3: regular vs random square grid ###
+
+**Program:** `wave_comparison.c` (variant of `wave_billiard.c`) 
+
+**Initial condition in function `animation()`:** `init_planar_wave(XMIN + 0.01, 0.0, phi, psi, xy_in);`
+
+```
+#define MOVIE 1         /* set to 1 to generate movie */
+
+/* General geometrical parameters */
+
+#define WINWIDTH 	1280  /* window width */
+#define WINHEIGHT 	720   /* window height */
+
+#define NX 1280          /* number of grid points on x axis */
+#define NY 720          /* number of grid points on y axis */
+
+/* setting NX to WINWIDTH and NY to WINHEIGHT increases resolution */
+/* but will multiply run time by 4                                 */
+
+#define XMIN -2.0
+#define XMAX 2.0	/* x interval */
+#define YMIN -1.125
+#define YMAX 1.125	/* y interval for 9/16 aspect ratio */
+
+#define JULIA_SCALE 1.0 /* scaling for Julia sets */
+
+/* Choice of the billiard table */
+
+#define B_DOMAIN_A 20      /* choice of domain shape, see list in global_pdes.c */
+#define B_DOMAIN_B 20      /* choice of domain shape, see list in global_pdes.c */
+
+#define CIRCLE_PATTERN_A 0    /* pattern of circles, see list in global_pdes.c */
+#define CIRCLE_PATTERN_B 2    /* pattern of circles, see list in global_pdes.c */
+
+#define P_PERCOL 0.25       /* probability of having a circle in C_RAND_PERCOL arrangement */
+#define NPOISSON 300        /* number of points for Poisson C_RAND_POISSON arrangement */
+
+#define LAMBDA 0.75	    /* parameter controlling the dimensions of domain */
+#define MU 0.03	    /* parameter controlling the dimensions of domain */
+#define NPOLY 3             /* number of sides of polygon */
+#define APOLY 1.0           /* angle by which to turn polygon, in units of Pi/2 */ 
+#define MDEPTH 4            /* depth of computation of Menger gasket */
+#define MRATIO 3            /* ratio defining Menger gasket */
+#define MANDELLEVEL 1000      /* iteration level for Mandelbrot set */
+#define MANDELLIMIT 10.0     /* limit value for approximation of Mandelbrot set */
+#define FOCI 1              /* set to 1 to draw focal points of ellipse */
+#define NGRIDX 15            /* number of grid point for grid of disks */
+#define NGRIDY 20           /* number of grid point for grid of disks */
+
+/* You can add more billiard tables by adapting the functions */
+/* xy_in_billiard and draw_billiard below */
+
+/* Physical parameters of wave equation */
+
+#define TWOSPEEDS 0         /* set to 1 to replace hardcore boundary by medium with different speed */
+
+#define OMEGA 0.9           /* frequency of periodic excitation */
+#define COURANT 0.01       /* Courant number */
+#define COURANTB 0.0075      /* Courant number in medium B */
+#define GAMMA 0.0      /* damping factor in wave equation */
+#define GAMMA_SIDES 1.0e-4      /* damping factor on boundary */
+#define GAMMA_TOPBOT 1.0e-6      /* damping factor on boundary */
+#define KAPPA 0.0       /* "elasticity" term enforcing oscillations */
+#define KAPPA_SIDES 5.0e-4       /* "elasticity" term on absorbing boundary */
+#define KAPPA_TOPBOT 0.0       /* "elasticity" term on absorbing boundary */
+/* The Courant number is given by c*DT/DX, where DT is the time step and DX the lattice spacing */
+/* The physical damping coefficient is given by GAMMA/(DT)^2 */
+/* Increasing COURANT speeds up the simulation, but decreases accuracy */
+/* For similar wave forms, COURANT^2*GAMMA should be kept constant */
+
+/* Boundary conditions, see list in global_pdes.c  */
+
+#define B_COND 3
+
+/* Parameters for length and speed of simulation */
+
+#define NSTEPS 3000      /* number of frames of movie */
+#define NVID 40          /* number of iterations between images displayed on screen */
+#define NSEG 100         /* number of segments of boundary */
+#define INITIAL_TIME 250    /* time after which to start saving frames */
+
+#define PAUSE 1000         /* number of frames after which to pause */
+#define PSLEEP 1         /* sleep time during pause */
+#define SLEEP1  1        /* initial sleeping time */
+#define SLEEP2  1   /* final sleeping time */
+
+/* Plot type, see list in global_pdes.c  */
+
+#define PLOT 1
+
+/* Color schemes */
+
+#define BLACK 1          /* background */
+
+#define COLOR_SCHEME 1   /* choice of color scheme, see list in global_pdes.c  */
+
+#define SCALE 0          /* set to 1 to adjust color scheme to variance of field */
+#define SLOPE 1.0        /* sensitivity of color on wave amplitude */
+#define ATTENUATION 0.0  /* exponential attenuation coefficient of contrast with time */
+#define E_SCALE 750.0     /* scaling factor for energy representation */
+
+#define COLORHUE 260     /* initial hue of water color for scheme C_LUM */
+#define COLORDRIFT 0.0   /* how much the color hue drifts during the whole simulation */
+#define LUMMEAN 0.5      /* amplitude of luminosity variation for scheme C_LUM */
+#define LUMAMP 0.3       /* amplitude of luminosity variation for scheme C_LUM */
+#define HUEMEAN 220.0    /* mean value of hue for color scheme C_HUE */
+#define HUEAMP -220.0      /* amplitude of variation of hue for color scheme C_HUE */
+
+```
+
 ### 2 August 21 - Avoiding a laser in a room of mirrors, long version with close-up and particle counter ###
 
 **Program:** Variant of `particle_billiard.c` (to be published later on)
@@ -33,7 +144,7 @@ updated gradually. For older simulations, see files `Parameters_July21.md` etc.
 
 #define B_DOMAIN 21      /* choice of domain shape */
 
-#define CIRCLE_PATTERN 3    /* pattern of circles */
+#define CIRCLE_PATTERN 11    /* pattern of circles */
 
 #define ABSORBING_CIRCLES 1 /* set to 1 for circular scatterers to be absorbing */