From 9f79de85cc327e27668efccdb53832ee509723b8 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Wed, 25 Nov 2020 08:25:02 +0100
Subject: [PATCH] Script to run small tests for parallel implementation
---
wk/test_parallel.m | 44 ++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 44 insertions(+)
create mode 100644 wk/test_parallel.m
diff --git a/wk/test_parallel.m b/wk/test_parallel.m
new file mode 100644
index 00000000..cf70fd46
--- /dev/null
+++ b/wk/test_parallel.m
@@ -0,0 +1,44 @@
+clear all;
+addpath(genpath('../matlab')) % ... add
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%% Set Up parameters
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%% PHYSICAL PARAMETERS
+NU = 1e-2; % Collision frequency
+TAU = 1.0; % e/i temperature ratio
+ETAB = 0.5; % Magnetic gradient
+ETAN = 1.0; % Density gradient
+ETAT = 0.0; % Temperature gradient
+MU = 5e-6; % Hyper diffusivity coefficient
+NOISE0 = 1.0e-5;
+%% GRID PARAMETERS
+N = 256; % Frequency gridpoints (Nkr = N/2)
+L = 66; % Size of the squared frequency domain
+PMAXE = 3; % Highest electron Hermite polynomial degree
+JMAXE = 2; % Highest '' Laguerre ''
+PMAXI = 3; % Highest ion Hermite polynomial degree
+JMAXI = 2; % Highest '' Laguerre ''
+%% TIME PARAMETERS
+TMAX = 400; % Maximal time unit
+DT = 1e-2; % Time step
+SPS0D = 2; % Sampling per time unit for profiler
+SPS2D = 2; % Sampling per time unit for 2D arrays
+SPS5D = 0.5; % Sampling per time unit for 5D arrays
+RESTART = 0; % To restart from last checkpoint
+JOB2LOAD= 0;
+%% OPTIONS
+SIMID = 'test_parallel'; % Name of the simulation
+CO = -2; % Collision operator (0 : L.Bernstein, -1 : Full Coulomb, -2 : Dougherty)
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%% unused
+KR0KH = 0; A0KH = 0; % Background phi mode to drive Ray-Tay inst.
+NO_E = 0; % Remove electrons dynamic
+% DK = 0; % Drift kinetic model (put every kernel_n to 0 except n=0 to 1)
+KREQ0 = 0; % put kr = 0
+KPAR = 0.0; % Parellel wave vector component
+LAMBDAD = 0.0;
+NON_LIN = 1 *(1-KREQ0); % activate non-linearity (is cancelled if KREQ0 = 1)
+CANCEL_ODD_P = 0;% Cancels the odd polynomials degree
+
+setup
--
GitLab