1D Euler Equations - Sod Shock Tube (Time Windowed DMD)
See 1D Euler Equations - Sod Shock Tube to familiarize yourself with this case.
Description:
Location: hypar/Examples/1D/Euler1D/SodShockTube_libROM_DMD_Train (This directory contains all the input files needed to run this case.)
Governing equations: 1D Euler equations (euler1d.h)
Reduced Order Modeling: This example trains a DMD object and then predicts the solution using the DMD at the same times that the actual HyPar solution is written at.
References:
- G.A. Sod, "A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws," J. Comput. Phys., 27, 1 (1978).
- C. B. Laney, "Computational Gasdynamics", Cambridge University Press, 1998.
Domain: \(0 \le x \le 1.0\), "extrapolate" (_EXTRAPOLATE_) boundary conditions
Initial Solution:
- \( 0 \le x < 0.5\): \(\rho = 1, u = 0, p = 1\)
- \( 0.5 \le x \le 1\): \(\rho = 0.125, u = 0, p = 0.1\)
Numerical Method:
- Spatial discretization (hyperbolic): Characteristic-based 5th order WENO (Interp1PrimFifthOrderWENOChar())
- Time integration: RK4 (TimeRK(), _RK_44_)
Reduced Order Modeling:
- Mode: train (libROMInterface::m_mode, _ROM_MODE_TRAIN_)
- Component Mode: monolithic (libROMInterface::m_comp_mode, _ROM_COMP_MODE_MONOLITHIC_)
- Type: Dynamic Mode Decomposition (DMD) with time windowing (libROMInterface::m_rom_type)
- Latent subspace dimension: 8 (DMDROMObject::m_rdim)
- Number of samples per time window: 40 (DMDROMObject::m_num_window_samples)
Input files required:
librom.inp
begin
rdim 8
sampling_frequency 1
mode train
dmd_num_win_samples 40
end
solver.inp:
begin
ndims 1
nvars 3
size 201
iproc 10
ghost 3
n_iter 400
time_scheme rk
time_scheme_type 44
hyp_space_scheme weno5
hyp_interp_type characteristic
dt 0.0005
screen_op_iter 1
file_op_iter 50
op_file_format binary
op_overwrite no
model euler1d
end
boundary.inp
2
extrapolate 0 1 0 0
extrapolate 0 -1 0 0
physics.inp
begin
gamma 1.4
upwinding rf-char
end
To generate initial.inp, compile and run the following code in the run directory:
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
int main(){
int NI=101,ndims=1;
FILE *in;
char ip_file_type[50];
strcpy(ip_file_type,"ascii");
printf("Reading file \"solver.inp\"...\n");
in = fopen("solver.inp","r");
if (!in) {
printf("Error: Input file \"solver.inp\" not found. Default values will be used.\n");
} else {
char word[500];
fscanf(in,"%s",word);
if (!strcmp(word, "begin")){
while (strcmp(word, "end")){
fscanf(in,"%s",word);
if (!strcmp(word, "ndims")) fscanf(in,"%d",&ndims);
else if (!strcmp(word, "size")) fscanf(in,"%d",&NI);
else if (!strcmp(word, "ip_file_type")) fscanf(in,"%s",ip_file_type);
}
} else {
printf("Error: Illegal format in solver.inp. Crash and burn!\n");
}
}
fclose(in);
if (ndims != 1) {
printf("ndims is not 1 in solver.inp. this code is to generate 1D initial conditions\n");
return(0);
}
printf("Grid:\t\t\t%d\n",NI);
int i;
double dx = 1.0 / ((double)(NI-1));
double *x, *rho,*rhou,*e;
x = (double*) calloc (NI, sizeof(double));
rho = (double*) calloc (NI, sizeof(double));
rhou = (double*) calloc (NI, sizeof(double));
e = (double*) calloc (NI, sizeof(double));
for (i = 0; i < NI; i++){
x[i] = i*dx;
double RHO,U,P;
if (x[i] < 0.5) {
RHO = 1.0;
U = 0.0;
P = 1.0;
} else {
RHO = 0.125;
U = 0;
P = 0.1;
}
rho[i] = RHO;
rhou[i] = RHO*U;
e[i] = P/0.4 + 0.5*RHO*U*U;
}
if (!strcmp(ip_file_type,"ascii")) {
FILE *out;
out = fopen("initial.inp","w");
for (i = 0; i < NI; i++) fprintf(out,"%lf ",x[i]);
fprintf(out,"\n");
for (i = 0; i < NI; i++) fprintf(out,"%lf ",rho[i]);
fprintf(out,"\n");
for (i = 0; i < NI; i++) fprintf(out,"%lf ",rhou[i]);
fprintf(out,"\n");
for (i = 0; i < NI; i++) fprintf(out,"%lf ",e[i]);
fprintf(out,"\n");
fclose(out);
} else if ((!strcmp(ip_file_type,"binary")) || (!strcmp(ip_file_type,"bin"))) {
printf("Error: Writing binary initial solution file not implemented. ");
printf("Please choose ip_file_type in solver.inp as \"ascii\".\n");
}
free(x);
free(rho);
free(rhou);
free(e);
return(0);
}
Output:
After running the code, there should be the following output files:
- 9 output files op_00000.bin, op_00001.bin, ... op_00008.bin; these are the HyPar solutions.
- 9 output files op_rom_00000.bin, op_rom_00001.bin, ... op_rom_00008.bin; these are the predicted solutions from the DMD object(s).
The first of each of these file sets is the solution at \(t=0\) and the final one is the solution at \(t=2\). Since HyPar::op_overwrite is set to no in solver.inp, a separate file is written for solutions at each output time. All the files are binary (HyPar::op_file_format is set to binary in solver.inp).
The provided Python script (plotSolution.py) can be used to generate plots from the binary files that compare the HyPar and DMD solutions. It will plot the 3 conserved variables - density, momentum, and energy. Alternatively, HyPar::op_file_format can be set to text, and GNUPlot or something similar can be used to plot the resulting text files.
The following plot shows the final solution (density) - FOM (full-order model) refers to the HyPar solution, ROM (reduced-order model) refers to the DMD solution.
The L1, L2, and Linf norms of the diff between the HyPar and ROM solution at the final time are calculated and reported on screen (see below) as well as pde_rom_diff.dat:
201 10 5.0000000000000001E-04 5.6218372311378950E-04 2.2766977309887507E-03 1.6120671488108019E-02 8.0859740000000002E+00 8.1111880000000003E+00
The numbers are: number of grid points (HyPar::dim_global), number of processors (MPIVariables::iproc), time step size (HyPar::dt), L1, L2, and L-infinity norms of the diff (HyPar::rom_diff_norms), solver wall time (seconds) (i.e., not accounting for initialization, and cleaning up), and total wall time.
By default, the code will write the trained DMD object(s) to files in a subdirectory (DMDROMObject::m_dirname - default value is "DMD"). If the subdirectory does not exist, the code may not report an error (or give some HDF5 file-writing error); the DMD objects will not be written! If the subdirectory exists, several files will exist after the simulation is complete - they are in a format that is readable by libROM.
Expected screen output:
HyPar - Parallel (MPI) version with 10 processes
Compiled with PETSc time integration.
Allocated simulation object(s).
Reading solver inputs from file "solver.inp".
No. of dimensions : 1
No. of variables : 3
Domain size : 201
Processes along each dimension : 10
Exact solution domain size : 201
No. of ghosts pts : 3
No. of iter. : 400
Restart iteration : 0
Time integration scheme : rk (44)
Spatial discretization scheme (hyperbolic) : weno5
Split hyperbolic flux term? : no
Interpolation type for hyperbolic term : characteristic
Spatial discretization type (parabolic ) : nonconservative-1stage
Spatial discretization scheme (parabolic ) : 2
Time Step : 5.000000E-04
Check for conservation : no
Screen output iterations : 1
File output iterations : 50
Initial solution file type : ascii
Initial solution read mode : serial
Solution file write mode : serial
Solution file format : binary
Overwrite solution file : no
Physical model : euler1d
Partitioning domain and allocating data arrays.
Reading array from ASCII file initial.inp (Serial mode).
Volume integral of the initial solution:
0: 5.6312499999999999E-01
1: 0.0000000000000000E+00
2: 1.3762500000000000E+00
Reading boundary conditions from boundary.inp.
Boundary extrapolate: Along dimension 0 and face +1
Boundary extrapolate: Along dimension 0 and face -1
2 boundary condition(s) read.
Initializing solvers.
Reading WENO parameters from weno.inp.
Initializing physics. Model = "euler1d"
Reading physical model inputs from file "physics.inp".
Setting up time integration.
Setting up libROM interface.
libROM inputs and parameters:
reduced model dimensionality: 8
sampling frequency: 1
mode: train
type: DMD
save to file: true
local vector size: 60
libROM DMD inputs:
number of samples per window: 40
directory name for DMD onjects: DMD
Solving in time (from 0 to 400 iterations)
Writing solution file op_00000.bin.
DMDROMObject::takeSample() - creating new DMD object, t=0.000000 (total: 1).
iter= 1 t=5.000E-04 CFL=1.183E-01 norm=1.3928E-02 wctime: 6.1E-04 (s)
iter= 2 t=1.000E-03 CFL=1.407E-01 norm=1.3304E-02 wctime: 4.7E-04 (s)
iter= 3 t=1.500E-03 CFL=1.617E-01 norm=1.2540E-02 wctime: 4.6E-04 (s)
iter= 4 t=2.000E-03 CFL=1.756E-01 norm=1.1267E-02 wctime: 4.9E-04 (s)
iter= 5 t=2.500E-03 CFL=1.848E-01 norm=1.0102E-02 wctime: 4.9E-04 (s)
iter= 6 t=3.000E-03 CFL=1.912E-01 norm=9.2728E-03 wctime: 4.8E-04 (s)
iter= 7 t=3.500E-03 CFL=1.956E-01 norm=8.8385E-03 wctime: 4.8E-04 (s)
iter= 8 t=4.000E-03 CFL=1.983E-01 norm=8.8084E-03 wctime: 4.7E-04 (s)
iter= 9 t=4.500E-03 CFL=1.996E-01 norm=8.6687E-03 wctime: 4.6E-04 (s)
iter= 10 t=5.000E-03 CFL=1.999E-01 norm=8.4023E-03 wctime: 4.4E-04 (s)
iter= 11 t=5.500E-03 CFL=1.998E-01 norm=8.1570E-03 wctime: 5.1E-04 (s)
iter= 12 t=6.000E-03 CFL=2.010E-01 norm=8.0063E-03 wctime: 4.9E-04 (s)
iter= 13 t=6.500E-03 CFL=2.055E-01 norm=8.0330E-03 wctime: 4.9E-04 (s)
iter= 14 t=7.000E-03 CFL=2.080E-01 norm=8.0151E-03 wctime: 4.4E-04 (s)
iter= 15 t=7.500E-03 CFL=2.091E-01 norm=7.8389E-03 wctime: 4.9E-04 (s)
iter= 16 t=8.000E-03 CFL=2.095E-01 norm=7.6569E-03 wctime: 4.4E-04 (s)
iter= 17 t=8.500E-03 CFL=2.094E-01 norm=7.5528E-03 wctime: 4.9E-04 (s)
iter= 18 t=9.000E-03 CFL=2.088E-01 norm=7.5441E-03 wctime: 4.5E-04 (s)
iter= 19 t=9.500E-03 CFL=2.079E-01 norm=7.6453E-03 wctime: 4.4E-04 (s)
iter= 20 t=1.000E-02 CFL=2.107E-01 norm=7.6219E-03 wctime: 4.4E-04 (s)
iter= 21 t=1.050E-02 CFL=2.124E-01 norm=7.4669E-03 wctime: 5.4E-04 (s)
iter= 22 t=1.100E-02 CFL=2.132E-01 norm=7.3193E-03 wctime: 4.9E-04 (s)
iter= 23 t=1.150E-02 CFL=2.135E-01 norm=7.2489E-03 wctime: 4.8E-04 (s)
iter= 24 t=1.200E-02 CFL=2.134E-01 norm=7.3010E-03 wctime: 4.8E-04 (s)
iter= 25 t=1.250E-02 CFL=2.130E-01 norm=7.3878E-03 wctime: 4.5E-04 (s)
iter= 26 t=1.300E-02 CFL=2.123E-01 norm=7.3133E-03 wctime: 4.4E-04 (s)
iter= 27 t=1.350E-02 CFL=2.136E-01 norm=7.1491E-03 wctime: 4.6E-04 (s)
iter= 28 t=1.400E-02 CFL=2.146E-01 norm=7.0327E-03 wctime: 4.5E-04 (s)
iter= 29 t=1.450E-02 CFL=2.152E-01 norm=7.0047E-03 wctime: 4.4E-04 (s)
iter= 30 t=1.500E-02 CFL=2.154E-01 norm=7.1079E-03 wctime: 4.8E-04 (s)
iter= 31 t=1.550E-02 CFL=2.154E-01 norm=7.1552E-03 wctime: 4.3E-04 (s)
iter= 32 t=1.600E-02 CFL=2.152E-01 norm=7.0472E-03 wctime: 4.7E-04 (s)
iter= 33 t=1.650E-02 CFL=2.148E-01 norm=6.9057E-03 wctime: 4.5E-04 (s)
iter= 34 t=1.700E-02 CFL=2.156E-01 norm=6.8315E-03 wctime: 4.7E-04 (s)
iter= 35 t=1.750E-02 CFL=2.162E-01 norm=6.8546E-03 wctime: 4.5E-04 (s)
iter= 36 t=1.800E-02 CFL=2.165E-01 norm=6.9765E-03 wctime: 4.6E-04 (s)
iter= 37 t=1.850E-02 CFL=2.167E-01 norm=6.9766E-03 wctime: 4.5E-04 (s)
iter= 38 t=1.900E-02 CFL=2.167E-01 norm=6.8460E-03 wctime: 4.6E-04 (s)
iter= 39 t=1.950E-02 CFL=2.166E-01 norm=6.7252E-03 wctime: 4.5E-04 (s)
iter= 40 t=2.000E-02 CFL=2.164E-01 norm=6.6847E-03 wctime: 4.6E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.020000 (total: 2).
iter= 41 t=2.050E-02 CFL=2.168E-01 norm=6.7577E-03 wctime: 4.5E-04 (s)
iter= 42 t=2.100E-02 CFL=2.172E-01 norm=6.8630E-03 wctime: 4.6E-04 (s)
iter= 43 t=2.150E-02 CFL=2.174E-01 norm=6.8159E-03 wctime: 4.7E-04 (s)
iter= 44 t=2.200E-02 CFL=2.175E-01 norm=6.6782E-03 wctime: 4.4E-04 (s)
iter= 45 t=2.250E-02 CFL=2.175E-01 norm=6.5853E-03 wctime: 4.9E-04 (s)
iter= 46 t=2.300E-02 CFL=2.174E-01 norm=6.5785E-03 wctime: 4.4E-04 (s)
iter= 47 t=2.350E-02 CFL=2.174E-01 norm=6.6914E-03 wctime: 4.8E-04 (s)
iter= 48 t=2.400E-02 CFL=2.177E-01 norm=6.7629E-03 wctime: 4.3E-04 (s)
iter= 49 t=2.450E-02 CFL=2.179E-01 norm=6.6776E-03 wctime: 5.2E-04 (s)
iter= 50 t=2.500E-02 CFL=2.180E-01 norm=6.5471E-03 wctime: 4.3E-04 (s)
Writing solution file op_00001.bin.
iter= 51 t=2.550E-02 CFL=2.181E-01 norm=6.4832E-03 wctime: 5.3E-04 (s)
iter= 52 t=2.600E-02 CFL=2.180E-01 norm=6.5124E-03 wctime: 4.5E-04 (s)
iter= 53 t=2.650E-02 CFL=2.180E-01 norm=6.6413E-03 wctime: 4.3E-04 (s)
iter= 54 t=2.700E-02 CFL=2.180E-01 norm=6.6649E-03 wctime: 4.7E-04 (s)
iter= 55 t=2.750E-02 CFL=2.182E-01 norm=6.5513E-03 wctime: 4.3E-04 (s)
iter= 56 t=2.800E-02 CFL=2.183E-01 norm=6.4377E-03 wctime: 5.2E-04 (s)
iter= 57 t=2.850E-02 CFL=2.184E-01 norm=6.4036E-03 wctime: 4.3E-04 (s)
iter= 58 t=2.900E-02 CFL=2.184E-01 norm=6.4767E-03 wctime: 4.9E-04 (s)
iter= 59 t=2.950E-02 CFL=2.184E-01 norm=6.5974E-03 wctime: 4.3E-04 (s)
iter= 60 t=3.000E-02 CFL=2.183E-01 norm=6.5714E-03 wctime: 4.7E-04 (s)
iter= 61 t=3.050E-02 CFL=2.184E-01 norm=6.4429E-03 wctime: 4.6E-04 (s)
iter= 62 t=3.100E-02 CFL=2.185E-01 norm=6.3513E-03 wctime: 4.5E-04 (s)
iter= 63 t=3.150E-02 CFL=2.186E-01 norm=6.3443E-03 wctime: 4.4E-04 (s)
iter= 64 t=3.200E-02 CFL=2.187E-01 norm=6.4562E-03 wctime: 4.5E-04 (s)
iter= 65 t=3.250E-02 CFL=2.187E-01 norm=6.5452E-03 wctime: 4.5E-04 (s)
iter= 66 t=3.300E-02 CFL=2.186E-01 norm=6.4750E-03 wctime: 4.3E-04 (s)
iter= 67 t=3.350E-02 CFL=2.186E-01 norm=6.3469E-03 wctime: 4.7E-04 (s)
iter= 68 t=3.400E-02 CFL=2.186E-01 norm=6.2820E-03 wctime: 4.3E-04 (s)
iter= 69 t=3.450E-02 CFL=2.187E-01 norm=6.3082E-03 wctime: 4.7E-04 (s)
iter= 70 t=3.500E-02 CFL=2.188E-01 norm=6.4425E-03 wctime: 5.1E-04 (s)
iter= 71 t=3.550E-02 CFL=2.188E-01 norm=6.4853E-03 wctime: 5.1E-04 (s)
iter= 72 t=3.600E-02 CFL=2.188E-01 norm=6.3811E-03 wctime: 4.3E-04 (s)
iter= 73 t=3.650E-02 CFL=2.188E-01 norm=6.2647E-03 wctime: 4.5E-04 (s)
iter= 74 t=3.700E-02 CFL=2.188E-01 norm=6.2260E-03 wctime: 4.7E-04 (s)
iter= 75 t=3.750E-02 CFL=2.188E-01 norm=6.2912E-03 wctime: 4.7E-04 (s)
iter= 76 t=3.800E-02 CFL=2.188E-01 norm=6.4219E-03 wctime: 4.5E-04 (s)
iter= 77 t=3.850E-02 CFL=2.188E-01 norm=6.4134E-03 wctime: 4.5E-04 (s)
iter= 78 t=3.900E-02 CFL=2.189E-01 norm=6.2897E-03 wctime: 4.5E-04 (s)
iter= 79 t=3.950E-02 CFL=2.189E-01 norm=6.1938E-03 wctime: 4.5E-04 (s)
iter= 80 t=4.000E-02 CFL=2.189E-01 norm=6.1825E-03 wctime: 4.5E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.040000 (total: 3).
iter= 81 t=4.050E-02 CFL=2.189E-01 norm=6.2883E-03 wctime: 4.3E-04 (s)
iter= 82 t=4.100E-02 CFL=2.188E-01 norm=6.3926E-03 wctime: 4.9E-04 (s)
iter= 83 t=4.150E-02 CFL=2.189E-01 norm=6.3366E-03 wctime: 4.3E-04 (s)
iter= 84 t=4.200E-02 CFL=2.189E-01 norm=6.2077E-03 wctime: 4.7E-04 (s)
iter= 85 t=4.250E-02 CFL=2.189E-01 norm=6.1362E-03 wctime: 4.3E-04 (s)
iter= 86 t=4.300E-02 CFL=2.189E-01 norm=6.1543E-03 wctime: 4.6E-04 (s)
iter= 87 t=4.350E-02 CFL=2.189E-01 norm=6.2873E-03 wctime: 4.6E-04 (s)
iter= 88 t=4.400E-02 CFL=2.189E-01 norm=6.3486E-03 wctime: 4.5E-04 (s)
iter= 89 t=4.450E-02 CFL=2.189E-01 norm=6.2553E-03 wctime: 4.4E-04 (s)
iter= 90 t=4.500E-02 CFL=2.189E-01 norm=6.1349E-03 wctime: 4.7E-04 (s)
iter= 91 t=4.550E-02 CFL=2.189E-01 norm=6.0902E-03 wctime: 4.7E-04 (s)
iter= 92 t=4.600E-02 CFL=2.190E-01 norm=6.1452E-03 wctime: 4.5E-04 (s)
iter= 93 t=4.650E-02 CFL=2.190E-01 norm=6.2833E-03 wctime: 4.7E-04 (s)
iter= 94 t=4.700E-02 CFL=2.190E-01 norm=6.2949E-03 wctime: 4.3E-04 (s)
iter= 95 t=4.750E-02 CFL=2.190E-01 norm=6.1774E-03 wctime: 4.7E-04 (s)
iter= 96 t=4.800E-02 CFL=2.190E-01 norm=6.0745E-03 wctime: 4.3E-04 (s)
iter= 97 t=4.850E-02 CFL=2.190E-01 norm=6.0558E-03 wctime: 4.3E-04 (s)
iter= 98 t=4.900E-02 CFL=2.190E-01 norm=6.1502E-03 wctime: 4.3E-04 (s)
iter= 99 t=4.950E-02 CFL=2.190E-01 norm=6.2693E-03 wctime: 4.6E-04 (s)
iter= 100 t=5.000E-02 CFL=2.190E-01 norm=6.2327E-03 wctime: 4.5E-04 (s)
Writing solution file op_00002.bin.
iter= 101 t=5.050E-02 CFL=2.190E-01 norm=6.1057E-03 wctime: 5.1E-04 (s)
iter= 102 t=5.100E-02 CFL=2.190E-01 norm=6.0269E-03 wctime: 4.3E-04 (s)
iter= 103 t=5.150E-02 CFL=2.190E-01 norm=6.0361E-03 wctime: 4.3E-04 (s)
iter= 104 t=5.200E-02 CFL=2.190E-01 norm=6.1634E-03 wctime: 4.4E-04 (s)
iter= 105 t=5.250E-02 CFL=2.191E-01 norm=6.2457E-03 wctime: 4.5E-04 (s)
iter= 106 t=5.300E-02 CFL=2.191E-01 norm=6.1680E-03 wctime: 4.3E-04 (s)
iter= 107 t=5.350E-02 CFL=2.191E-01 norm=6.0440E-03 wctime: 4.3E-04 (s)
iter= 108 t=5.400E-02 CFL=2.191E-01 norm=5.9906E-03 wctime: 4.3E-04 (s)
iter= 109 t=5.450E-02 CFL=2.191E-01 norm=6.0322E-03 wctime: 4.7E-04 (s)
iter= 110 t=5.500E-02 CFL=2.191E-01 norm=6.1736E-03 wctime: 4.3E-04 (s)
iter= 111 t=5.550E-02 CFL=2.191E-01 norm=6.2092E-03 wctime: 4.3E-04 (s)
iter= 112 t=5.600E-02 CFL=2.191E-01 norm=6.1020E-03 wctime: 4.3E-04 (s)
iter= 113 t=5.650E-02 CFL=2.191E-01 norm=5.9920E-03 wctime: 4.5E-04 (s)
iter= 114 t=5.700E-02 CFL=2.192E-01 norm=5.9645E-03 wctime: 4.4E-04 (s)
iter= 115 t=5.750E-02 CFL=2.192E-01 norm=6.0444E-03 wctime: 4.5E-04 (s)
iter= 116 t=5.800E-02 CFL=2.192E-01 norm=6.1771E-03 wctime: 4.4E-04 (s)
iter= 117 t=5.850E-02 CFL=2.192E-01 norm=6.1633E-03 wctime: 4.3E-04 (s)
iter= 118 t=5.900E-02 CFL=2.192E-01 norm=6.0397E-03 wctime: 4.6E-04 (s)
iter= 119 t=5.950E-02 CFL=2.192E-01 norm=5.9503E-03 wctime: 5.0E-04 (s)
iter= 120 t=6.000E-02 CFL=2.192E-01 norm=5.9483E-03 wctime: 4.6E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.060000 (total: 4).
iter= 121 t=6.050E-02 CFL=2.192E-01 norm=6.0644E-03 wctime: 4.3E-04 (s)
iter= 122 t=6.100E-02 CFL=2.192E-01 norm=6.1675E-03 wctime: 4.4E-04 (s)
iter= 123 t=6.150E-02 CFL=2.193E-01 norm=6.1089E-03 wctime: 4.5E-04 (s)
iter= 124 t=6.200E-02 CFL=2.193E-01 norm=5.9826E-03 wctime: 4.3E-04 (s)
iter= 125 t=6.250E-02 CFL=2.193E-01 norm=5.9186E-03 wctime: 4.3E-04 (s)
iter= 126 t=6.300E-02 CFL=2.193E-01 norm=5.9463E-03 wctime: 4.4E-04 (s)
iter= 127 t=6.350E-02 CFL=2.193E-01 norm=6.0858E-03 wctime: 4.4E-04 (s)
iter= 128 t=6.400E-02 CFL=2.193E-01 norm=6.1458E-03 wctime: 4.3E-04 (s)
iter= 129 t=6.450E-02 CFL=2.193E-01 norm=6.0515E-03 wctime: 4.3E-04 (s)
iter= 130 t=6.500E-02 CFL=2.193E-01 norm=5.9337E-03 wctime: 4.3E-04 (s)
iter= 131 t=6.550E-02 CFL=2.193E-01 norm=5.8951E-03 wctime: 4.5E-04 (s)
iter= 132 t=6.600E-02 CFL=2.193E-01 norm=5.9579E-03 wctime: 4.4E-04 (s)
iter= 133 t=6.650E-02 CFL=2.194E-01 norm=6.0993E-03 wctime: 4.3E-04 (s)
iter= 134 t=6.700E-02 CFL=2.194E-01 norm=6.1110E-03 wctime: 4.3E-04 (s)
iter= 135 t=6.750E-02 CFL=2.193E-01 norm=5.9943E-03 wctime: 4.3E-04 (s)
iter= 136 t=6.800E-02 CFL=2.193E-01 norm=5.8947E-03 wctime: 4.4E-04 (s)
iter= 137 t=6.850E-02 CFL=2.193E-01 norm=5.8818E-03 wctime: 4.3E-04 (s)
iter= 138 t=6.900E-02 CFL=2.194E-01 norm=5.9821E-03 wctime: 4.3E-04 (s)
iter= 139 t=6.950E-02 CFL=2.194E-01 norm=6.1035E-03 wctime: 4.3E-04 (s)
iter= 140 t=7.000E-02 CFL=2.194E-01 norm=6.0676E-03 wctime: 4.4E-04 (s)
iter= 141 t=7.050E-02 CFL=2.194E-01 norm=5.9413E-03 wctime: 4.4E-04 (s)
iter= 142 t=7.100E-02 CFL=2.194E-01 norm=5.8649E-03 wctime: 4.3E-04 (s)
iter= 143 t=7.150E-02 CFL=2.193E-01 norm=5.8786E-03 wctime: 4.3E-04 (s)
iter= 144 t=7.200E-02 CFL=2.193E-01 norm=6.0094E-03 wctime: 4.3E-04 (s)
iter= 145 t=7.250E-02 CFL=2.193E-01 norm=6.0946E-03 wctime: 4.6E-04 (s)
iter= 146 t=7.300E-02 CFL=2.193E-01 norm=6.0186E-03 wctime: 4.3E-04 (s)
iter= 147 t=7.350E-02 CFL=2.193E-01 norm=5.8957E-03 wctime: 4.3E-04 (s)
iter= 148 t=7.400E-02 CFL=2.193E-01 norm=5.8452E-03 wctime: 4.3E-04 (s)
iter= 149 t=7.450E-02 CFL=2.193E-01 norm=5.8903E-03 wctime: 4.4E-04 (s)
iter= 150 t=7.500E-02 CFL=2.193E-01 norm=6.0345E-03 wctime: 4.4E-04 (s)
Writing solution file op_00003.bin.
iter= 151 t=7.550E-02 CFL=2.193E-01 norm=6.0738E-03 wctime: 4.9E-04 (s)
iter= 152 t=7.600E-02 CFL=2.193E-01 norm=5.9680E-03 wctime: 4.3E-04 (s)
iter= 153 t=7.650E-02 CFL=2.193E-01 norm=5.8580E-03 wctime: 4.3E-04 (s)
iter= 154 t=7.700E-02 CFL=2.193E-01 norm=5.8327E-03 wctime: 4.4E-04 (s)
iter= 155 t=7.750E-02 CFL=2.193E-01 norm=5.9142E-03 wctime: 4.4E-04 (s)
iter= 156 t=7.800E-02 CFL=2.193E-01 norm=6.0506E-03 wctime: 4.3E-04 (s)
iter= 157 t=7.850E-02 CFL=2.193E-01 norm=6.0420E-03 wctime: 4.3E-04 (s)
iter= 158 t=7.900E-02 CFL=2.193E-01 norm=5.9199E-03 wctime: 4.3E-04 (s)
iter= 159 t=7.950E-02 CFL=2.193E-01 norm=5.8307E-03 wctime: 4.7E-04 (s)
iter= 160 t=8.000E-02 CFL=2.193E-01 norm=5.8307E-03 wctime: 4.4E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.080000 (total: 5).
iter= 161 t=8.050E-02 CFL=2.193E-01 norm=5.9473E-03 wctime: 4.7E-04 (s)
iter= 162 t=8.100E-02 CFL=2.193E-01 norm=6.0559E-03 wctime: 4.3E-04 (s)
iter= 163 t=8.150E-02 CFL=2.193E-01 norm=6.0020E-03 wctime: 4.4E-04 (s)
iter= 164 t=8.200E-02 CFL=2.193E-01 norm=5.8757E-03 wctime: 4.3E-04 (s)
iter= 165 t=8.250E-02 CFL=2.193E-01 norm=5.8111E-03 wctime: 4.3E-04 (s)
iter= 166 t=8.300E-02 CFL=2.193E-01 norm=5.8394E-03 wctime: 4.3E-04 (s)
iter= 167 t=8.350E-02 CFL=2.193E-01 norm=5.9806E-03 wctime: 4.3E-04 (s)
iter= 168 t=8.400E-02 CFL=2.192E-01 norm=6.0481E-03 wctime: 5.4E-04 (s)
iter= 169 t=8.450E-02 CFL=2.192E-01 norm=5.9579E-03 wctime: 4.7E-04 (s)
iter= 170 t=8.500E-02 CFL=2.192E-01 norm=5.8393E-03 wctime: 4.3E-04 (s)
iter= 171 t=8.550E-02 CFL=2.192E-01 norm=5.8004E-03 wctime: 4.3E-04 (s)
iter= 172 t=8.600E-02 CFL=2.192E-01 norm=5.8624E-03 wctime: 4.3E-04 (s)
iter= 173 t=8.650E-02 CFL=2.192E-01 norm=6.0078E-03 wctime: 4.3E-04 (s)
iter= 174 t=8.700E-02 CFL=2.192E-01 norm=6.0268E-03 wctime: 4.3E-04 (s)
iter= 175 t=8.750E-02 CFL=2.192E-01 norm=5.9117E-03 wctime: 4.3E-04 (s)
iter= 176 t=8.800E-02 CFL=2.192E-01 norm=5.8099E-03 wctime: 4.4E-04 (s)
iter= 177 t=8.850E-02 CFL=2.192E-01 norm=5.7963E-03 wctime: 4.4E-04 (s)
iter= 178 t=8.900E-02 CFL=2.192E-01 norm=5.8961E-03 wctime: 4.3E-04 (s)
iter= 179 t=8.950E-02 CFL=2.192E-01 norm=6.0244E-03 wctime: 4.7E-04 (s)
iter= 180 t=9.000E-02 CFL=2.192E-01 norm=5.9951E-03 wctime: 4.3E-04 (s)
iter= 181 t=9.050E-02 CFL=2.193E-01 norm=5.8688E-03 wctime: 4.3E-04 (s)
iter= 182 t=9.100E-02 CFL=2.193E-01 norm=5.7901E-03 wctime: 4.3E-04 (s)
iter= 183 t=9.150E-02 CFL=2.193E-01 norm=5.8026E-03 wctime: 4.4E-04 (s)
iter= 184 t=9.200E-02 CFL=2.193E-01 norm=5.9346E-03 wctime: 4.4E-04 (s)
iter= 185 t=9.250E-02 CFL=2.193E-01 norm=6.0270E-03 wctime: 4.4E-04 (s)
iter= 186 t=9.300E-02 CFL=2.193E-01 norm=5.9545E-03 wctime: 4.4E-04 (s)
iter= 187 t=9.350E-02 CFL=2.193E-01 norm=5.8295E-03 wctime: 4.3E-04 (s)
iter= 188 t=9.400E-02 CFL=2.192E-01 norm=5.7767E-03 wctime: 4.6E-04 (s)
iter= 189 t=9.450E-02 CFL=2.193E-01 norm=5.8211E-03 wctime: 4.3E-04 (s)
iter= 190 t=9.500E-02 CFL=2.193E-01 norm=5.9698E-03 wctime: 4.3E-04 (s)
iter= 191 t=9.550E-02 CFL=2.193E-01 norm=6.0156E-03 wctime: 4.3E-04 (s)
iter= 192 t=9.600E-02 CFL=2.193E-01 norm=5.9107E-03 wctime: 4.3E-04 (s)
iter= 193 t=9.650E-02 CFL=2.193E-01 norm=5.7980E-03 wctime: 4.3E-04 (s)
iter= 194 t=9.700E-02 CFL=2.193E-01 norm=5.7708E-03 wctime: 4.5E-04 (s)
iter= 195 t=9.750E-02 CFL=2.193E-01 norm=5.8525E-03 wctime: 4.4E-04 (s)
iter= 196 t=9.800E-02 CFL=2.192E-01 norm=5.9944E-03 wctime: 4.3E-04 (s)
iter= 197 t=9.850E-02 CFL=2.192E-01 norm=5.9890E-03 wctime: 4.6E-04 (s)
iter= 198 t=9.900E-02 CFL=2.192E-01 norm=5.8653E-03 wctime: 4.2E-04 (s)
iter= 199 t=9.950E-02 CFL=2.192E-01 norm=5.7735E-03 wctime: 4.3E-04 (s)
iter= 200 t=1.000E-01 CFL=2.192E-01 norm=5.7725E-03 wctime: 4.3E-04 (s)
Writing solution file op_00004.bin.
DMDROMObject::takeSample() - creating new DMD object, t=0.100000 (total: 6).
iter= 201 t=1.005E-01 CFL=2.192E-01 norm=5.8924E-03 wctime: 4.9E-04 (s)
iter= 202 t=1.010E-01 CFL=2.192E-01 norm=6.0060E-03 wctime: 4.8E-04 (s)
iter= 203 t=1.015E-01 CFL=2.192E-01 norm=5.9524E-03 wctime: 4.3E-04 (s)
iter= 204 t=1.020E-01 CFL=2.192E-01 norm=5.8236E-03 wctime: 4.6E-04 (s)
iter= 205 t=1.025E-01 CFL=2.192E-01 norm=5.7570E-03 wctime: 4.4E-04 (s)
iter= 206 t=1.030E-01 CFL=2.192E-01 norm=5.7854E-03 wctime: 4.6E-04 (s)
iter= 207 t=1.035E-01 CFL=2.192E-01 norm=5.9313E-03 wctime: 4.5E-04 (s)
iter= 208 t=1.040E-01 CFL=2.192E-01 norm=6.0005E-03 wctime: 4.5E-04 (s)
iter= 209 t=1.045E-01 CFL=2.192E-01 norm=5.9077E-03 wctime: 4.5E-04 (s)
iter= 210 t=1.050E-01 CFL=2.192E-01 norm=5.7868E-03 wctime: 4.3E-04 (s)
iter= 211 t=1.055E-01 CFL=2.192E-01 norm=5.7472E-03 wctime: 4.7E-04 (s)
iter= 212 t=1.060E-01 CFL=2.192E-01 norm=5.8124E-03 wctime: 4.3E-04 (s)
iter= 213 t=1.065E-01 CFL=2.192E-01 norm=5.9628E-03 wctime: 4.7E-04 (s)
iter= 214 t=1.070E-01 CFL=2.192E-01 norm=5.9801E-03 wctime: 4.3E-04 (s)
iter= 215 t=1.075E-01 CFL=2.192E-01 norm=5.8614E-03 wctime: 4.6E-04 (s)
iter= 216 t=1.080E-01 CFL=2.192E-01 norm=5.7581E-03 wctime: 4.5E-04 (s)
iter= 217 t=1.085E-01 CFL=2.192E-01 norm=5.7446E-03 wctime: 5.2E-04 (s)
iter= 218 t=1.090E-01 CFL=2.192E-01 norm=5.8500E-03 wctime: 4.7E-04 (s)
iter= 219 t=1.095E-01 CFL=2.192E-01 norm=5.9799E-03 wctime: 4.5E-04 (s)
iter= 220 t=1.100E-01 CFL=2.192E-01 norm=5.9455E-03 wctime: 4.5E-04 (s)
iter= 221 t=1.105E-01 CFL=2.192E-01 norm=5.8160E-03 wctime: 4.5E-04 (s)
iter= 222 t=1.110E-01 CFL=2.192E-01 norm=5.7377E-03 wctime: 4.5E-04 (s)
iter= 223 t=1.115E-01 CFL=2.192E-01 norm=5.7530E-03 wctime: 4.3E-04 (s)
iter= 224 t=1.120E-01 CFL=2.192E-01 norm=5.8921E-03 wctime: 4.6E-04 (s)
iter= 225 t=1.125E-01 CFL=2.192E-01 norm=5.9819E-03 wctime: 4.3E-04 (s)
iter= 226 t=1.130E-01 CFL=2.192E-01 norm=5.9028E-03 wctime: 4.6E-04 (s)
iter= 227 t=1.135E-01 CFL=2.192E-01 norm=5.7759E-03 wctime: 4.3E-04 (s)
iter= 228 t=1.140E-01 CFL=2.192E-01 norm=5.7244E-03 wctime: 4.6E-04 (s)
iter= 229 t=1.145E-01 CFL=2.192E-01 norm=5.7742E-03 wctime: 4.4E-04 (s)
iter= 230 t=1.150E-01 CFL=2.192E-01 norm=5.9275E-03 wctime: 4.3E-04 (s)
iter= 231 t=1.155E-01 CFL=2.192E-01 norm=5.9664E-03 wctime: 4.3E-04 (s)
iter= 232 t=1.160E-01 CFL=2.192E-01 norm=5.8555E-03 wctime: 4.4E-04 (s)
iter= 233 t=1.165E-01 CFL=2.192E-01 norm=5.7434E-03 wctime: 4.4E-04 (s)
iter= 234 t=1.170E-01 CFL=2.192E-01 norm=5.7192E-03 wctime: 4.3E-04 (s)
iter= 235 t=1.175E-01 CFL=2.192E-01 norm=5.8097E-03 wctime: 4.8E-04 (s)
iter= 236 t=1.180E-01 CFL=2.192E-01 norm=5.9517E-03 wctime: 4.3E-04 (s)
iter= 237 t=1.185E-01 CFL=2.192E-01 norm=5.9368E-03 wctime: 4.3E-04 (s)
iter= 238 t=1.190E-01 CFL=2.192E-01 norm=5.8086E-03 wctime: 4.3E-04 (s)
iter= 239 t=1.195E-01 CFL=2.192E-01 norm=5.7192E-03 wctime: 4.4E-04 (s)
iter= 240 t=1.200E-01 CFL=2.192E-01 norm=5.7226E-03 wctime: 4.4E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.120000 (total: 7).
iter= 241 t=1.205E-01 CFL=2.192E-01 norm=5.8515E-03 wctime: 4.4E-04 (s)
iter= 242 t=1.210E-01 CFL=2.192E-01 norm=5.9595E-03 wctime: 4.4E-04 (s)
iter= 243 t=1.215E-01 CFL=2.192E-01 norm=5.8960E-03 wctime: 4.4E-04 (s)
iter= 244 t=1.220E-01 CFL=2.192E-01 norm=5.7660E-03 wctime: 5.0E-04 (s)
iter= 245 t=1.225E-01 CFL=2.193E-01 norm=5.7038E-03 wctime: 4.3E-04 (s)
iter= 246 t=1.230E-01 CFL=2.193E-01 norm=5.7397E-03 wctime: 4.7E-04 (s)
iter= 247 t=1.235E-01 CFL=2.193E-01 norm=5.8917E-03 wctime: 4.3E-04 (s)
iter= 248 t=1.240E-01 CFL=2.193E-01 norm=5.9509E-03 wctime: 4.7E-04 (s)
iter= 249 t=1.245E-01 CFL=2.193E-01 norm=5.8495E-03 wctime: 4.5E-04 (s)
iter= 250 t=1.250E-01 CFL=2.193E-01 norm=5.7299E-03 wctime: 4.7E-04 (s)
Writing solution file op_00005.bin.
iter= 251 t=1.255E-01 CFL=2.192E-01 norm=5.6951E-03 wctime: 4.9E-04 (s)
iter= 252 t=1.260E-01 CFL=2.193E-01 norm=5.7704E-03 wctime: 4.6E-04 (s)
iter= 253 t=1.265E-01 CFL=2.193E-01 norm=5.9208E-03 wctime: 4.5E-04 (s)
iter= 254 t=1.270E-01 CFL=2.193E-01 norm=5.9260E-03 wctime: 4.6E-04 (s)
iter= 255 t=1.275E-01 CFL=2.193E-01 norm=5.8021E-03 wctime: 4.6E-04 (s)
iter= 256 t=1.280E-01 CFL=2.193E-01 norm=5.7032E-03 wctime: 4.3E-04 (s)
iter= 257 t=1.285E-01 CFL=2.193E-01 norm=5.6958E-03 wctime: 5.0E-04 (s)
iter= 258 t=1.290E-01 CFL=2.193E-01 norm=5.8120E-03 wctime: 4.3E-04 (s)
iter= 259 t=1.295E-01 CFL=2.193E-01 norm=5.9354E-03 wctime: 4.8E-04 (s)
iter= 260 t=1.300E-01 CFL=2.193E-01 norm=5.8887E-03 wctime: 4.3E-04 (s)
iter= 261 t=1.305E-01 CFL=2.192E-01 norm=5.7570E-03 wctime: 4.6E-04 (s)
iter= 262 t=1.310E-01 CFL=2.193E-01 norm=5.6842E-03 wctime: 4.5E-04 (s)
iter= 263 t=1.315E-01 CFL=2.193E-01 norm=5.7073E-03 wctime: 4.6E-04 (s)
iter= 264 t=1.320E-01 CFL=2.193E-01 norm=5.8542E-03 wctime: 4.7E-04 (s)
iter= 265 t=1.325E-01 CFL=2.193E-01 norm=5.9332E-03 wctime: 4.6E-04 (s)
iter= 266 t=1.330E-01 CFL=2.193E-01 norm=5.8441E-03 wctime: 5.3E-04 (s)
iter= 267 t=1.335E-01 CFL=2.193E-01 norm=5.7190E-03 wctime: 5.0E-04 (s)
iter= 268 t=1.340E-01 CFL=2.193E-01 norm=5.6739E-03 wctime: 4.5E-04 (s)
iter= 269 t=1.345E-01 CFL=2.193E-01 norm=5.7341E-03 wctime: 4.5E-04 (s)
iter= 270 t=1.350E-01 CFL=2.193E-01 norm=5.8887E-03 wctime: 4.5E-04 (s)
iter= 271 t=1.355E-01 CFL=2.192E-01 norm=5.9141E-03 wctime: 4.3E-04 (s)
iter= 272 t=1.360E-01 CFL=2.192E-01 norm=5.7963E-03 wctime: 4.7E-04 (s)
iter= 273 t=1.365E-01 CFL=2.192E-01 norm=5.6885E-03 wctime: 4.3E-04 (s)
iter= 274 t=1.370E-01 CFL=2.193E-01 norm=5.6709E-03 wctime: 4.8E-04 (s)
iter= 275 t=1.375E-01 CFL=2.193E-01 norm=5.7729E-03 wctime: 4.3E-04 (s)
iter= 276 t=1.380E-01 CFL=2.193E-01 norm=5.9096E-03 wctime: 4.4E-04 (s)
iter= 277 t=1.385E-01 CFL=2.193E-01 norm=5.8816E-03 wctime: 4.3E-04 (s)
iter= 278 t=1.390E-01 CFL=2.193E-01 norm=5.7509E-03 wctime: 4.6E-04 (s)
iter= 279 t=1.395E-01 CFL=2.193E-01 norm=5.6678E-03 wctime: 4.5E-04 (s)
iter= 280 t=1.400E-01 CFL=2.192E-01 norm=5.6789E-03 wctime: 4.8E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.140000 (total: 8).
iter= 281 t=1.405E-01 CFL=2.192E-01 norm=5.8167E-03 wctime: 4.8E-04 (s)
iter= 282 t=1.410E-01 CFL=2.192E-01 norm=5.9140E-03 wctime: 4.7E-04 (s)
iter= 283 t=1.415E-01 CFL=2.192E-01 norm=5.8391E-03 wctime: 4.8E-04 (s)
iter= 284 t=1.420E-01 CFL=2.192E-01 norm=5.7099E-03 wctime: 4.5E-04 (s)
iter= 285 t=1.425E-01 CFL=2.192E-01 norm=5.6546E-03 wctime: 4.6E-04 (s)
iter= 286 t=1.430E-01 CFL=2.192E-01 norm=5.7003E-03 wctime: 4.3E-04 (s)
iter= 287 t=1.435E-01 CFL=2.192E-01 norm=5.8556E-03 wctime: 4.4E-04 (s)
iter= 288 t=1.440E-01 CFL=2.192E-01 norm=5.9019E-03 wctime: 4.3E-04 (s)
iter= 289 t=1.445E-01 CFL=2.192E-01 norm=5.7932E-03 wctime: 4.9E-04 (s)
iter= 290 t=1.450E-01 CFL=2.192E-01 norm=5.6776E-03 wctime: 4.5E-04 (s)
iter= 291 t=1.455E-01 CFL=2.192E-01 norm=5.6496E-03 wctime: 4.6E-04 (s)
iter= 292 t=1.460E-01 CFL=2.192E-01 norm=5.7358E-03 wctime: 4.5E-04 (s)
iter= 293 t=1.465E-01 CFL=2.192E-01 norm=5.8824E-03 wctime: 4.4E-04 (s)
iter= 294 t=1.470E-01 CFL=2.192E-01 norm=5.8743E-03 wctime: 4.3E-04 (s)
iter= 295 t=1.475E-01 CFL=2.192E-01 norm=5.7467E-03 wctime: 4.5E-04 (s)
iter= 296 t=1.480E-01 CFL=2.192E-01 norm=5.6537E-03 wctime: 5.0E-04 (s)
iter= 297 t=1.485E-01 CFL=2.192E-01 norm=5.6536E-03 wctime: 4.4E-04 (s)
iter= 298 t=1.490E-01 CFL=2.192E-01 norm=5.7797E-03 wctime: 4.6E-04 (s)
iter= 299 t=1.495E-01 CFL=2.192E-01 norm=5.8944E-03 wctime: 4.3E-04 (s)
iter= 300 t=1.500E-01 CFL=2.192E-01 norm=5.8364E-03 wctime: 4.4E-04 (s)
Writing solution file op_00006.bin.
iter= 301 t=1.505E-01 CFL=2.192E-01 norm=5.7050E-03 wctime: 5.3E-04 (s)
iter= 302 t=1.510E-01 CFL=2.192E-01 norm=5.6387E-03 wctime: 4.4E-04 (s)
iter= 303 t=1.515E-01 CFL=2.192E-01 norm=5.6703E-03 wctime: 4.4E-04 (s)
iter= 304 t=1.520E-01 CFL=2.192E-01 norm=5.8219E-03 wctime: 4.4E-04 (s)
iter= 305 t=1.525E-01 CFL=2.192E-01 norm=5.8890E-03 wctime: 4.3E-04 (s)
iter= 306 t=1.530E-01 CFL=2.192E-01 norm=5.7918E-03 wctime: 4.9E-04 (s)
iter= 307 t=1.535E-01 CFL=2.192E-01 norm=5.6697E-03 wctime: 4.3E-04 (s)
iter= 308 t=1.540E-01 CFL=2.192E-01 norm=5.6314E-03 wctime: 4.6E-04 (s)
iter= 309 t=1.545E-01 CFL=2.192E-01 norm=5.7017E-03 wctime: 4.3E-04 (s)
iter= 310 t=1.550E-01 CFL=2.192E-01 norm=5.8550E-03 wctime: 4.3E-04 (s)
iter= 311 t=1.555E-01 CFL=2.192E-01 norm=5.8685E-03 wctime: 4.3E-04 (s)
iter= 312 t=1.560E-01 CFL=2.192E-01 norm=5.7464E-03 wctime: 4.4E-04 (s)
iter= 313 t=1.565E-01 CFL=2.192E-01 norm=5.6436E-03 wctime: 4.6E-04 (s)
iter= 314 t=1.570E-01 CFL=2.192E-01 norm=5.6324E-03 wctime: 4.3E-04 (s)
iter= 315 t=1.575E-01 CFL=2.192E-01 norm=5.7437E-03 wctime: 4.9E-04 (s)
iter= 316 t=1.580E-01 CFL=2.192E-01 norm=5.8736E-03 wctime: 4.6E-04 (s)
iter= 317 t=1.585E-01 CFL=2.192E-01 norm=5.8348E-03 wctime: 4.3E-04 (s)
iter= 318 t=1.590E-01 CFL=2.192E-01 norm=5.7032E-03 wctime: 4.3E-04 (s)
iter= 319 t=1.595E-01 CFL=2.192E-01 norm=5.6263E-03 wctime: 4.5E-04 (s)
iter= 320 t=1.600E-01 CFL=2.192E-01 norm=5.6448E-03 wctime: 4.3E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.160000 (total: 9).
iter= 321 t=1.605E-01 CFL=2.192E-01 norm=5.7889E-03 wctime: 4.4E-04 (s)
iter= 322 t=1.610E-01 CFL=2.192E-01 norm=5.8766E-03 wctime: 4.4E-04 (s)
iter= 323 t=1.615E-01 CFL=2.192E-01 norm=5.7937E-03 wctime: 4.3E-04 (s)
iter= 324 t=1.620E-01 CFL=2.192E-01 norm=5.6663E-03 wctime: 4.3E-04 (s)
iter= 325 t=1.625E-01 CFL=2.192E-01 norm=5.6168E-03 wctime: 4.3E-04 (s)
iter= 326 t=1.630E-01 CFL=2.192E-01 norm=5.6709E-03 wctime: 4.6E-04 (s)
iter= 327 t=1.635E-01 CFL=2.192E-01 norm=5.8269E-03 wctime: 4.3E-04 (s)
iter= 328 t=1.640E-01 CFL=2.192E-01 norm=5.8627E-03 wctime: 4.5E-04 (s)
iter= 329 t=1.645E-01 CFL=2.192E-01 norm=5.7490E-03 wctime: 4.6E-04 (s)
iter= 330 t=1.650E-01 CFL=2.193E-01 norm=5.6374E-03 wctime: 4.9E-04 (s)
iter= 331 t=1.655E-01 CFL=2.193E-01 norm=5.6155E-03 wctime: 4.5E-04 (s)
iter= 332 t=1.660E-01 CFL=2.193E-01 norm=5.7106E-03 wctime: 4.3E-04 (s)
iter= 333 t=1.665E-01 CFL=2.193E-01 norm=5.8530E-03 wctime: 4.3E-04 (s)
iter= 334 t=1.670E-01 CFL=2.192E-01 norm=5.8351E-03 wctime: 4.4E-04 (s)
iter= 335 t=1.675E-01 CFL=2.192E-01 norm=5.7055E-03 wctime: 4.4E-04 (s)
iter= 336 t=1.680E-01 CFL=2.192E-01 norm=5.6174E-03 wctime: 4.3E-04 (s)
iter= 337 t=1.685E-01 CFL=2.192E-01 norm=5.6233E-03 wctime: 4.9E-04 (s)
iter= 338 t=1.690E-01 CFL=2.192E-01 norm=5.7561E-03 wctime: 4.3E-04 (s)
iter= 339 t=1.695E-01 CFL=2.193E-01 norm=5.8632E-03 wctime: 4.5E-04 (s)
iter= 340 t=1.700E-01 CFL=2.193E-01 norm=5.7972E-03 wctime: 4.3E-04 (s)
iter= 341 t=1.705E-01 CFL=2.193E-01 norm=5.6665E-03 wctime: 4.3E-04 (s)
iter= 342 t=1.710E-01 CFL=2.193E-01 norm=5.6059E-03 wctime: 4.3E-04 (s)
iter= 343 t=1.715E-01 CFL=2.193E-01 norm=5.6447E-03 wctime: 4.4E-04 (s)
iter= 344 t=1.720E-01 CFL=2.193E-01 norm=5.7991E-03 wctime: 4.4E-04 (s)
iter= 345 t=1.725E-01 CFL=2.192E-01 norm=5.8571E-03 wctime: 4.3E-04 (s)
iter= 346 t=1.730E-01 CFL=2.192E-01 norm=5.7542E-03 wctime: 5.0E-04 (s)
iter= 347 t=1.735E-01 CFL=2.192E-01 norm=5.6347E-03 wctime: 4.3E-04 (s)
iter= 348 t=1.740E-01 CFL=2.192E-01 norm=5.6015E-03 wctime: 4.5E-04 (s)
iter= 349 t=1.745E-01 CFL=2.193E-01 norm=5.6796E-03 wctime: 4.3E-04 (s)
iter= 350 t=1.750E-01 CFL=2.193E-01 norm=5.8311E-03 wctime: 4.3E-04 (s)
Writing solution file op_00007.bin.
iter= 351 t=1.755E-01 CFL=2.193E-01 norm=5.8354E-03 wctime: 4.9E-04 (s)
iter= 352 t=1.760E-01 CFL=2.193E-01 norm=5.7104E-03 wctime: 4.5E-04 (s)
iter= 353 t=1.765E-01 CFL=2.193E-01 norm=5.6119E-03 wctime: 4.3E-04 (s)
iter= 354 t=1.770E-01 CFL=2.193E-01 norm=5.6061E-03 wctime: 4.7E-04 (s)
iter= 355 t=1.775E-01 CFL=2.193E-01 norm=5.7247E-03 wctime: 4.6E-04 (s)
iter= 356 t=1.780E-01 CFL=2.192E-01 norm=5.8490E-03 wctime: 4.4E-04 (s)
iter= 357 t=1.785E-01 CFL=2.192E-01 norm=5.8017E-03 wctime: 4.6E-04 (s)
iter= 358 t=1.790E-01 CFL=2.192E-01 norm=5.6695E-03 wctime: 4.3E-04 (s)
iter= 359 t=1.795E-01 CFL=2.192E-01 norm=5.5972E-03 wctime: 4.3E-04 (s)
iter= 360 t=1.800E-01 CFL=2.193E-01 norm=5.6217E-03 wctime: 4.3E-04 (s)
DMDROMObject::takeSample() - creating new DMD object, t=0.180000 (total: 10).
iter= 361 t=1.805E-01 CFL=2.193E-01 norm=5.7703E-03 wctime: 5.1E-04 (s)
iter= 362 t=1.810E-01 CFL=2.193E-01 norm=5.8501E-03 wctime: 4.3E-04 (s)
iter= 363 t=1.815E-01 CFL=2.193E-01 norm=5.7607E-03 wctime: 5.0E-04 (s)
iter= 364 t=1.820E-01 CFL=2.193E-01 norm=5.6350E-03 wctime: 4.6E-04 (s)
iter= 365 t=1.825E-01 CFL=2.193E-01 norm=5.5905E-03 wctime: 4.5E-04 (s)
iter= 366 t=1.830E-01 CFL=2.192E-01 norm=5.6519E-03 wctime: 4.3E-04 (s)
iter= 367 t=1.835E-01 CFL=2.192E-01 norm=5.8082E-03 wctime: 4.4E-04 (s)
iter= 368 t=1.840E-01 CFL=2.192E-01 norm=5.8348E-03 wctime: 4.6E-04 (s)
iter= 369 t=1.845E-01 CFL=2.192E-01 norm=5.7168E-03 wctime: 4.5E-04 (s)
iter= 370 t=1.850E-01 CFL=2.192E-01 norm=5.6084E-03 wctime: 4.4E-04 (s)
iter= 371 t=1.855E-01 CFL=2.193E-01 norm=5.5912E-03 wctime: 5.3E-04 (s)
iter= 372 t=1.860E-01 CFL=2.193E-01 norm=5.6940E-03 wctime: 5.9E-04 (s)
iter= 373 t=1.865E-01 CFL=2.193E-01 norm=5.8326E-03 wctime: 4.8E-04 (s)
iter= 374 t=1.870E-01 CFL=2.193E-01 norm=5.8058E-03 wctime: 4.5E-04 (s)
iter= 375 t=1.875E-01 CFL=2.193E-01 norm=5.6746E-03 wctime: 4.5E-04 (s)
iter= 376 t=1.880E-01 CFL=2.193E-01 norm=5.5909E-03 wctime: 4.5E-04 (s)
iter= 377 t=1.885E-01 CFL=2.192E-01 norm=5.6023E-03 wctime: 4.3E-04 (s)
iter= 378 t=1.890E-01 CFL=2.192E-01 norm=5.7411E-03 wctime: 4.8E-04 (s)
iter= 379 t=1.895E-01 CFL=2.192E-01 norm=5.8409E-03 wctime: 4.3E-04 (s)
iter= 380 t=1.900E-01 CFL=2.192E-01 norm=5.7671E-03 wctime: 4.8E-04 (s)
iter= 381 t=1.905E-01 CFL=2.192E-01 norm=5.6370E-03 wctime: 4.4E-04 (s)
iter= 382 t=1.910E-01 CFL=2.193E-01 norm=5.5809E-03 wctime: 4.6E-04 (s)
iter= 383 t=1.915E-01 CFL=2.193E-01 norm=5.6266E-03 wctime: 4.5E-04 (s)
iter= 384 t=1.920E-01 CFL=2.193E-01 norm=5.7833E-03 wctime: 4.6E-04 (s)
iter= 385 t=1.925E-01 CFL=2.193E-01 norm=5.8324E-03 wctime: 4.5E-04 (s)
iter= 386 t=1.930E-01 CFL=2.193E-01 norm=5.7241E-03 wctime: 4.5E-04 (s)
iter= 387 t=1.935E-01 CFL=2.192E-01 norm=5.6073E-03 wctime: 4.5E-04 (s)
iter= 388 t=1.940E-01 CFL=2.192E-01 norm=5.5788E-03 wctime: 4.3E-04 (s)
iter= 389 t=1.945E-01 CFL=2.192E-01 norm=5.6648E-03 wctime: 4.6E-04 (s)
iter= 390 t=1.950E-01 CFL=2.192E-01 norm=5.8138E-03 wctime: 4.3E-04 (s)
iter= 391 t=1.955E-01 CFL=2.192E-01 norm=5.8084E-03 wctime: 4.9E-04 (s)
iter= 392 t=1.960E-01 CFL=2.193E-01 norm=5.6805E-03 wctime: 4.3E-04 (s)
iter= 393 t=1.965E-01 CFL=2.193E-01 norm=5.5859E-03 wctime: 4.4E-04 (s)
iter= 394 t=1.970E-01 CFL=2.193E-01 norm=5.5852E-03 wctime: 4.3E-04 (s)
iter= 395 t=1.975E-01 CFL=2.193E-01 norm=5.7113E-03 wctime: 4.4E-04 (s)
iter= 396 t=1.980E-01 CFL=2.193E-01 norm=5.8292E-03 wctime: 4.3E-04 (s)
iter= 397 t=1.985E-01 CFL=2.193E-01 norm=5.7732E-03 wctime: 4.5E-04 (s)
iter= 398 t=1.990E-01 CFL=2.192E-01 norm=5.6408E-03 wctime: 4.7E-04 (s)
iter= 399 t=1.995E-01 CFL=2.192E-01 norm=5.5732E-03 wctime: 4.3E-04 (s)
iter= 400 t=2.000E-01 CFL=2.192E-01 norm=5.6040E-03 wctime: 4.7E-04 (s)
Completed time integration (Final time: 0.200000), total wctime: 0.860301 (seconds).
Writing solution file op_00008.bin.
libROM: Training ROM.
DMDRomObject::train() - training DMD object 0 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 1 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 2 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 3 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 4 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 5 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 6 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 7 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 8 with 41 samples.
Using 8 basis vectors out of 40.
DMDRomObject::train() - training DMD object 9 with 40 samples.
Using 8 basis vectors out of 39.
libROM: wallclock time: 0.130851 (seconds).
libROM: Predicting solution at time 0.0000e+00 using ROM.
libROM: wallclock time: 0.000330 (seconds).
Writing solution file op_rom_00000.bin.
libROM: Predicting solution at time 2.5000e-02 using ROM.
libROM: wallclock time: 0.000290 (seconds).
Writing solution file op_rom_00001.bin.
libROM: Predicting solution at time 5.0000e-02 using ROM.
libROM: wallclock time: 0.000291 (seconds).
Writing solution file op_rom_00002.bin.
libROM: Predicting solution at time 7.5000e-02 using ROM.
libROM: wallclock time: 0.000278 (seconds).
Writing solution file op_rom_00003.bin.
libROM: Predicting solution at time 1.0000e-01 using ROM.
libROM: wallclock time: 0.000260 (seconds).
Writing solution file op_rom_00004.bin.
libROM: Predicting solution at time 1.2500e-01 using ROM.
libROM: wallclock time: 0.000313 (seconds).
Writing solution file op_rom_00005.bin.
libROM: Predicting solution at time 1.5000e-01 using ROM.
libROM: wallclock time: 0.000433 (seconds).
Writing solution file op_rom_00006.bin.
libROM: Predicting solution at time 1.7500e-01 using ROM.
libROM: wallclock time: 0.000276 (seconds).
Writing solution file op_rom_00007.bin.
libROM: Predicting solution at time 2.0000e-01 using ROM.
libROM: wallclock time: 0.000296 (seconds).
libROM: Calculating diff between PDE and ROM solutions.
Writing solution file op_rom_00008.bin.
libROM: total prediction/query wallclock time: 0.002767 (seconds).
libROMInterface::saveROM() - saving ROM objects.
Saving DMD object with filename root DMD/dmdobj_0000.
Saving DMD object with filename root DMD/dmdobj_0001.
Saving DMD object with filename root DMD/dmdobj_0002.
Saving DMD object with filename root DMD/dmdobj_0003.
Saving DMD object with filename root DMD/dmdobj_0004.
Saving DMD object with filename root DMD/dmdobj_0005.
Saving DMD object with filename root DMD/dmdobj_0006.
Saving DMD object with filename root DMD/dmdobj_0007.
Saving DMD object with filename root DMD/dmdobj_0008.
Saving DMD object with filename root DMD/dmdobj_0009.
Norms of the diff between ROM and PDE solutions for domain 0:
L1 Norm : 5.6218372311378950E-04
L2 Norm : 2.2766977309887507E-03
Linfinity Norm : 1.6120671488108019E-02
Solver runtime (in seconds): 8.0859740000000002E+00
Total runtime (in seconds): 8.1111880000000003E+00
Deallocating arrays.
Finished.