diff --git a/regression/README.md b/regression/README.md new file mode 100644 index 00000000..703fda22 --- /dev/null +++ b/regression/README.md @@ -0,0 +1,309 @@ +# MrHyDE Regression Testing + +This document describes the regression testing framework for MrHyDE. + +## Quick Start + +```bash +# From MrHyDE build directory +cd MrHyDE/regression +ln -sf /path/to/mrhyde.exe mrhyde +python3 scripts/runtests.py +``` + +## Directory Structure + +``` +regression/ + scripts/ + runtests.py # Main test runner + mrhyde_test_support.py # Test utilities for .tst scripts + burgers/ + 1D_bump/ + input.yaml # Problem configuration + mrhyde.tst # Test script (executable) + mrhyde.gold # Expected output + maxwell/ + PlaneWave/ + input.yaml, input_mesh.yaml, input_functions.yaml + mrhyde.tst + mrhyde.gold + maxwell_hybrid/ + input.yaml, ... + run.sh # Manual run script (no automated testing) + ref/mrhyde.ocs # Reference output for comparison + hpc/ # HPC tests (run with -k hpc) + scripts/ + parse_mrhyde_log.py # Log file parser + scaling_analysis.py # Scaling analysis utilities + strong_scaling/ + maxwell_cavity/ # Strong scaling test for Maxwell + ... +``` + +## Test Types + +### Type A: Automated Tests (.tst + .gold) + +These tests are discovered and run automatically by `runtests.py`. + +**Required files:** +- `input.yaml` - Problem configuration (mesh, physics, solver parameters) +- `mrhyde.tst` - Python test script with `#TESTING` directives (uses `mrhyde_test_support.py` utilities) +- `mrhyde.gold` - Expected output for diff comparison + +**How it works:** +1. `runtests.py` finds all `.tst` files recursively +2. Parses `#TESTING` directives to determine if test is active, keywords, etc. +3. Executes the `.tst` script +4. The script runs MrHyDE, cleans the log, and diffs against `.gold` +5. Exit status 0 = pass, non-zero = fail + +**Example .tst file** (`burgers/1D_bump/mrhyde.tst`): +```python +#!/usr/bin/env python3 +import sys, os +sys.path.append("../../scripts") +from mrhyde_test_support import * + +desc = '''Burgers 1D bump test''' +its = mrhyde_test_support(desc) +its.opts.verbose = True + +root = 'mrhyde' + +# Test metadata parsed by runtests.py +#TESTING active +#TESTING -n 1 +#TESTING -k Burgers,transient,1D + +status = 0 +status += its.call('mpiexec -n 1 ../../mrhyde >& mrhyde.log') +status += its.clean_log() +status += its.call('diff -y %s.log %s.gold' % (root, root)) + +if status == 0: print('Success.') +else: print('Failure.') +sys.exit(status) +``` + +### Type B: Manual Tests (run.sh only) + +These exist for development or manual verification. They are **not** run by `runtests.py`. + +**Files:** +- `input*.yaml` - Problem configuration +- `run.sh` - Shell script to run MrHyDE +- `ref/` - Reference output for manual comparison + +**Example** (`maxwell/maxwell_hybrid/run.sh`): +```bash +#!/bin/bash +mpiexec -n 1 ../../mrhyde >& mrhyde.log +``` + +## Test Script Directives (#TESTING) + +Place these comments in `.tst` files to control test behavior: + +| Directive | Description | +|-----------|-------------| +| `#TESTING active` | **Required** to run the test | +| `#TESTING -n N` | Number of MPI processes | +| `#TESTING -k key1,key2` | Keywords for filtering (comma = AND) | +| `#TESTING -K key1,key2` | Exclude keywords | +| `#TESTING -m machine1` | Run only on specific machines | +| `#TESTING -M machine1` | Exclude specific machines | +| `#TESTING -t hours` | Expected runtime in hours | + +## Running Tests + +### Basic Usage + +```bash +cd MrHyDE/regression +ln -sf /path/to/mrhyde.exe mrhyde + +# Run all active tests +python3 scripts/runtests.py + +# List tests without running +python3 scripts/runtests.py -i + +# Verbose output +python3 scripts/runtests.py -v +``` + +### Filtering Tests + +```bash +# Run only Maxwell tests +python3 scripts/runtests.py -k Maxwells + +# Run tests that are BOTH transient AND 3D +python3 scripts/runtests.py -k transient,3D + +# Exclude long-running tests +python3 scripts/runtests.py -K long + +# Run specific test directory +python3 scripts/runtests.py -d maxwell/PlaneWave +``` + +### Parallel Execution (-S flag) + +The `-S CORES` flag runs multiple *tests* concurrently, not MPI processes. Each test still launches its own MPI job (as specified by `#TESTING -n`). + +The scheduler tracks available cores and only launches a test when `test.nprocs <= available_cores`. Tests are sorted largest-first to optimize packing. + +```bash +# If you have 8 cores and tests use 1-4 MPI procs each +python3 scripts/runtests.py -S 8 +``` +If all tests use 4 MPI procs, only 2 tests run concurrently (8/4=2). + +Without `-S`, tests run sequentially (one at a time). + +### Other Options + +```bash +# Override processor count for all tests +python3 scripts/runtests.py -n 2 + +# Non-recursive (current directory only) +python3 scripts/runtests.py -N + +# List all keywords used across tests +python3 scripts/runtests.py --list-keywords +``` + +## runtests.py Reference + +### Working Options + +| Flag | Description | +|------|-------------| +| `-d DIR` | Starting directory (default: current) | +| `-e EXT` | Test extension (default: `.tst`) | +| `-f` | Force run inactive tests | +| `-i` | Info mode: list tests, don't run | +| `-k KEYWORDS` | Include matching keywords | +| `-K KEYWORDS` | Exclude matching keywords | +| `-n NPROCS` | Override processor count | +| `--nrange MIN,MAX` | Filter by processor range | +| `-N` | Non-recursive search | +| `-S [CORES]` | Parallel execution | +| `-t TESTS` | Comma-separated test list | +| `-T FILE` | File with test list | +| `-v` | Verbose output | +| `-s` | Simple reporting | +| `--list-keywords` | List all keywords | +| `--print-keywords` | Show keywords in results | +| `--include-all` | Include inactive in report | + +### Legacy Options (potentially broken) + +These options were inherited from the DGM project and have hardcoded values from ~2015: + +| Flag | Status | +|------|--------| +| `-b` | **Broken**: References non-existent scripts | +| `-q` | **Legacy**: Only works on old SNL clusters (redsky, glory, skybridge, curie) | +| `-c` | Passed to tests but doesn't affect runtests.py | +| `--32/--64` | Legacy bit mode | + +## Running on SLURM Clusters + +The built-in cluster support (`-b`, `-q`) is outdated. Use these interactive or batch scripts to run these tests on your local cluster. + +## HPC Regression Testing + +HPC regression tests verify parallel scaling behavior and are located in the `hpc/` directory. +These tests are not run by default. The test driver automatically excludes tests whose +keywords include `hpc` unless you explicitly request them with `-k hpc`. + +### Running HPC Tests + +```bash +cd MrHyDE/regression + +# Run all HPC tests +python3 scripts/runtests.py -k hpc +``` + +### HPC Test Structure + +``` +regression/ + hpc/ + scripts/ + parse_mrhyde_log.py # Log file parser + scaling_analysis.py # Scaling analysis utilities + strong_scaling/ + maxwell_cavity/ + input.yaml # Problem configuration + mrhyde.tst # Test script +``` + +### HPC Helper Scripts (standalone use) + +The two helper scripts in `hpc/scripts/` can be used directly on any MrHyDE log files +that contain compatible TimeMonitor output (`verbose: 10`). + +**1. Log parser: `parse_mrhyde_log.py`** + +From the `regression/` directory: + +```bash +# Parse a single log and print workload and timing tables +python3 hpc/scripts/parse_mrhyde_log.py \ + hpc/strong_scaling/maxwell_cavity/mrhyde_np1.log +``` +This script extracts: +- Per-processor workload information (elements, memory usage) +- Timing data from the TimeMonitor section + +**2. Scaling analysis: `scaling_analysis.py`** + +From the `regression/` directory: + +```bash +# Analyze scaling across multiple logs (for example 1, 2, 4 MPI processes) +python3 hpc/scripts/scaling_analysis.py \ + hpc/strong_scaling/maxwell_cavity/mrhyde_np1.log \ + hpc/strong_scaling/maxwell_cavity/mrhyde_np2.log \ + hpc/strong_scaling/maxwell_cavity/mrhyde_np4.log +This script: +- Loads timing data from each log (via `parse_mrhyde_log.py`) +- Computes speedup and parallel efficiency versus processor count + +## Creating a New Test + +1. Create directory: `regression/mymodule/mytest/` +2. Copy an existing test (e.g., `maxwell/PlaneWave/`) as template +3. Modify `input.yaml` and update `#TESTING` directives in `mrhyde.tst` +4. Run once, then `cp mrhyde.log mrhyde.gold` to create baseline +5. Verify: `./mrhyde.tst` should exit 0 + +As a guideline, design tests that run on the order of seconds to allow for frequent and fast regression testing. + +## Output + +Tests produce XML output (`TEST-package.xml`) compatible with Jenkins CI written to working directory from where `runtests.py` is run. +The summary shows: + +``` +Test Results from Directory: /path/MrHyDE/regression +Total number of test(s): 139 +----------------------------------------------------------------------------------------------- + 1/139 pass 14.53s np=4 maxwell/Planewave-LeapFrog1 + 2/139 pass 1.17s np=4 maxwell/PlaneWave + 3/139 pass 1.05s np=4 maxwell/Planewave-LeapFrog2 + 4/139 pass 1.79s np=4 phasefield/2d-3phi + ... + 139/139 pass 1.38s np=1 levelSet/wildfire +----------------------------------------------------------------------------------------------- + Pass: 139 Fail: 0 Skipped: 0 Total: 139 + +Total Runtime: 129.42s +``` diff --git a/sandbox/testCases/hcurl_verification/2d_sinusoidal/input.yaml b/regression/boundary_conditions/hcurl_2d_sinusoidal/input.yaml similarity index 96% rename from sandbox/testCases/hcurl_verification/2d_sinusoidal/input.yaml rename to regression/boundary_conditions/hcurl_2d_sinusoidal/input.yaml index 161f133f..e08034b8 100644 --- a/sandbox/testCases/hcurl_verification/2d_sinusoidal/input.yaml +++ b/regression/boundary_conditions/hcurl_2d_sinusoidal/input.yaml @@ -20,8 +20,8 @@ ANONYMOUS: xmax: 1.0 ymin: 0.0 ymax: 1.0 - NX: 32 - NY: 32 + NX: 16 + NY: 16 blocknames: eblock-0_0 Physics: @@ -79,8 +79,7 @@ ANONYMOUS: analysis type: forward Postprocess: - write solution: true - write frequency: 1 + write solution: false compute errors: true eblock-0_0: True solutions: diff --git a/sandbox/testCases/hcurl_verification/2d_sinusoidal/mrhyde.gold b/regression/boundary_conditions/hcurl_2d_sinusoidal/mrhyde.gold similarity index 63% rename from sandbox/testCases/hcurl_verification/2d_sinusoidal/mrhyde.gold rename to regression/boundary_conditions/hcurl_2d_sinusoidal/mrhyde.gold index f890dcbf..1f60e0d6 100644 --- a/sandbox/testCases/hcurl_verification/2d_sinusoidal/mrhyde.gold +++ b/regression/boundary_conditions/hcurl_2d_sinusoidal/mrhyde.gold @@ -2,16 +2,12 @@ STK Meta data: Element blocks = "eblock-0_0" Sidesets = "bottom" "left" "right" "top" Nodesets = "lower_left" "origin" - Fields = "coordinates" "edges" "PROC_ID" "LOAD_BAL" "Ex" "Ey" "B" + Fields = "coordinates" "edges" "PROC_ID" "LOAD_BAL" DOFManager Field Information: Element Block = eblock-0_0 "E" is field ID 0 "B" is field ID 1 -********************************************************* -***** Writing the solution to output.exo -********************************************************* - ******************************************************* @@ -24,21 +20,21 @@ DOFManager Field Information: ********************************************************* ***** Iteration: 0 -***** Norm of nonlinear residual: 0.00612665 +***** Norm of nonlinear residual: 0.0243957 ***** Scaled Norm of nonlinear residual: 1 ********************************************************* ********************************************************* ***** Iteration: 1 -***** Norm of nonlinear residual: 1.12645e-06 -***** Scaled Norm of nonlinear residual: 0.000183861 +***** Norm of nonlinear residual: 2.17702e-06 +***** Scaled Norm of nonlinear residual: 8.9238e-05 ********************************************************* Norm of solution: 1.01551 (overlapped vector so results may differ on multiple procs) ********************************************************* ***** Computing errors ****** -***** L2 norm of the error for E = 4.33856e-06 (time = 0) -***** L2 norm of the error for E = 4.33856e-06 (time = 1e-12) +***** L2 norm of the error for E = 3.28164e-05 (time = 0) +***** L2 norm of the error for E = 3.28164e-05 (time = 1e-12) ***** L2 norm of the error for B = 0 (time = 0) -***** L2 norm of the error for B = 3.14038e-12 (time = 1e-12) +***** L2 norm of the error for B = 3.13693e-12 (time = 1e-12) diff --git a/regression/boundary_conditions/hcurl_2d_sinusoidal/mrhyde.tst b/regression/boundary_conditions/hcurl_2d_sinusoidal/mrhyde.tst new file mode 100755 index 00000000..e79efed1 --- /dev/null +++ b/regression/boundary_conditions/hcurl_2d_sinusoidal/mrhyde.tst @@ -0,0 +1,42 @@ +#!/usr/bin/env python3 +#------------------------------------------------------------------------------- + +import sys, os +sys.path.append("../../scripts") +from mrhyde_test_support import * + +# ============================================================================== +# Parsing input + +desc = '''HCURL 2D Inhomogeneous Dirichlet BC verification + E = (sin(pi*y), sin(pi*x)) with component-based BCs + ''' + +its = mrhyde_test_support(desc) + +print('Because of the diff test on the log file, this test needs ') +print('to run with "-v". There is a buffering issue.') +print('Setting the verbosity to True.') +its.opts.verbose = True + +#------------------------------------------------------------------------------- +# Problem Parameters + +root = 'mrhyde' + +# These comments are for testing with the runtest.py utility. +#TESTING active +#TESTING -n 1 +#TESTING -k maxwell,HCURL,Dirichlet,2D,inhomogeneous + +# ============================================================================== +status = 0 + +status += its.call('mpiexec -n 2 ../../mrhyde >& mrhyde.log') +status += its.clean_log() +status += its.call('diff -y %s.log %s.gold' % (root, root)) + +# ============================================================================== +if status == 0: print('Success.') +else: print('Failure.') +sys.exit(status) diff --git a/sandbox/testCases/hcurl_verification/rotating_linear/input.yaml b/regression/boundary_conditions/hcurl_rotating_linear/input.yaml similarity index 97% rename from sandbox/testCases/hcurl_verification/rotating_linear/input.yaml rename to regression/boundary_conditions/hcurl_rotating_linear/input.yaml index 607d7980..5abbb319 100644 --- a/sandbox/testCases/hcurl_verification/rotating_linear/input.yaml +++ b/regression/boundary_conditions/hcurl_rotating_linear/input.yaml @@ -79,8 +79,7 @@ ANONYMOUS: analysis type: forward Postprocess: - write solution: true - write frequency: 1 + write solution: false compute errors: true eblock-0_0_0: True solutions: diff --git a/sandbox/testCases/hcurl_verification/rotating_linear/mrhyde.gold b/regression/boundary_conditions/hcurl_rotating_linear/mrhyde.gold similarity index 60% rename from sandbox/testCases/hcurl_verification/rotating_linear/mrhyde.gold rename to regression/boundary_conditions/hcurl_rotating_linear/mrhyde.gold index 97d674d6..93a7b7bd 100644 --- a/sandbox/testCases/hcurl_verification/rotating_linear/mrhyde.gold +++ b/regression/boundary_conditions/hcurl_rotating_linear/mrhyde.gold @@ -3,16 +3,12 @@ STK Meta data: Element blocks = "eblock-0_0_0" Sidesets = "back" "bottom" "front" "left" "right" "top" Nodesets = "origin" - Fields = "coordinates" "edges" "faces" "PROC_ID" "LOAD_BAL" "Ex" "Ey" "Ez" "Bx" "By" "Bz" + Fields = "coordinates" "edges" "faces" "PROC_ID" "LOAD_BAL" DOFManager Field Information: Element Block = eblock-0_0_0 "E" is field ID 0 "B" is field ID 1 -********************************************************* -***** Writing the solution to output.exo -********************************************************* - ******************************************************* @@ -31,23 +27,23 @@ DOFManager Field Information: ********************************************************* ***** Iteration: 1 -***** Norm of nonlinear residual: 3.57288e-08 -***** Scaled Norm of nonlinear residual: 1.78644e-08 +***** Norm of nonlinear residual: 3.44976e-08 +***** Scaled Norm of nonlinear residual: 1.72488e-08 ********************************************************* ********************************************************* ***** Iteration: 2 -***** Norm of nonlinear residual: 3.57288e-08 -***** Scaled Norm of nonlinear residual: 1.78644e-08 +***** Norm of nonlinear residual: 3.44976e-08 +***** Scaled Norm of nonlinear residual: 1.72488e-08 ********************************************************* -Norm of solution: 1.20221 (overlapped vector so results may differ on multiple procs) +Norm of solution: 1.08253 (overlapped vector so results may differ on multiple procs) ********************************************************* ***** Computing errors ****** -***** L2 norm of the error for E = 2.43416e-16 (time = 0) -***** L2 norm of the error for E = 2.48556e-16 (time = 1e-10) -***** L2 norm of the error for curl(E) = 2.1725e-15 (time = 0) -***** L2 norm of the error for curl(E) = 2.1978e-15 (time = 1e-10) +***** L2 norm of the error for E = 2.54667e-16 (time = 0) +***** L2 norm of the error for E = 2.53146e-16 (time = 1e-10) +***** L2 norm of the error for curl(E) = 2.30949e-15 (time = 0) +***** L2 norm of the error for curl(E) = 2.25341e-15 (time = 1e-10) ***** L2 norm of the error for B = 0 (time = 0) ***** L2 norm of the error for B = 2e-10 (time = 1e-10) diff --git a/regression/boundary_conditions/hcurl_rotating_linear/mrhyde.tst b/regression/boundary_conditions/hcurl_rotating_linear/mrhyde.tst new file mode 100755 index 00000000..12cfa5ec --- /dev/null +++ b/regression/boundary_conditions/hcurl_rotating_linear/mrhyde.tst @@ -0,0 +1,43 @@ +#!/usr/bin/env python3 +#------------------------------------------------------------------------------- + +import sys, os +sys.path.append("../../scripts") +from mrhyde_test_support import * + +# ============================================================================== +# Parsing input + +desc = '''HCURL 3D Component-based Dirichlet BC verification + E = (-y, x, 0) with curl(E) = (0, 0, 2) + L2 error expected at machine precision + ''' + +its = mrhyde_test_support(desc) + +print('Because of the diff test on the log file, this test needs ') +print('to run with "-v". There is a buffering issue.') +print('Setting the verbosity to True.') +its.opts.verbose = True + +#------------------------------------------------------------------------------- +# Problem Parameters + +root = 'mrhyde' + +# These comments are for testing with the runtest.py utility. +#TESTING active +#TESTING -n 1 +#TESTING -k maxwell,HCURL,Dirichlet,3D,component-BC + +# ============================================================================== +status = 0 + +status += its.call('mpiexec -n 1 ../../mrhyde >& mrhyde.log') +status += its.clean_log() +status += its.call('diff -y %s.log %s.gold' % (root, root)) + +# ============================================================================== +if status == 0: print('Success.') +else: print('Failure.') +sys.exit(status) diff --git a/regression/boundary_conditions/hgrad_dirichlet/input.yaml b/regression/boundary_conditions/hgrad_dirichlet/input.yaml new file mode 100644 index 00000000..180673b7 --- /dev/null +++ b/regression/boundary_conditions/hgrad_dirichlet/input.yaml @@ -0,0 +1,45 @@ +%YAML 1.1 +--- +# HGRAD Dirichlet BC verification test +# Inhomogeneous Dirichlet BCs: T = x*y on all boundaries +# Manufactured solution: T = x*y, source = 0 (\Delta T = 0) +ANONYMOUS: + verbosity: 5 + debug level: 0 + Mesh: + dimension: 2 + element type: quad + xmin: 0.0 + xmax: 1.0 + ymin: 0.0 + ymax: 1.0 + NX: 10 + NY: 10 + Functions: + thermal source: '0.0' + Physics: + modules: thermal + assemble boundary terms: false + build face terms: true + Dirichlet conditions: + T: + all boundaries: x*y + Initial conditions: + T: '0.0' + Discretization: + order: + T: 1 + quadrature: 2 + Solver: + solver: steady-state + nonlinear TOL: 1.0e-08 + max nonlinear iters: 2 + use strong DBCs: true + Analysis: + analysis type: forward + Postprocess: + compute errors: true + write solution: false + True solutions: + T: x*y +... diff --git a/regression/boundary_conditions/hgrad_dirichlet/mrhyde.gold b/regression/boundary_conditions/hgrad_dirichlet/mrhyde.gold new file mode 100644 index 00000000..6f62812d --- /dev/null +++ b/regression/boundary_conditions/hgrad_dirichlet/mrhyde.gold @@ -0,0 +1,26 @@ +STK Meta data: + Element blocks = "eblock-0_0" + Sidesets = "bottom" "left" "right" "top" + Nodesets = "lower_left" "origin" + Fields = "coordinates" "edges" "PROC_ID" "LOAD_BAL" +DOFManager Field Information: + Element Block = eblock-0_0 + "T" is field ID 0 + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.46667 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.88178e-16 +***** Scaled Norm of nonlinear residual: 6.05576e-16 +********************************************************* +Norm of solution: 3.85 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Computing errors ****** + +***** L2 norm of the error for T = 5.11314e-08 (time = 0) diff --git a/regression/boundary_conditions/hgrad_dirichlet/mrhyde.tst b/regression/boundary_conditions/hgrad_dirichlet/mrhyde.tst new file mode 100755 index 00000000..aea009d3 --- /dev/null +++ b/regression/boundary_conditions/hgrad_dirichlet/mrhyde.tst @@ -0,0 +1,44 @@ +#!/usr/bin/env python3 +#------------------------------------------------------------------------------- + +import sys, os +sys.path.append("../../scripts") +from mrhyde_test_support import * + +# ============================================================================== +# Parsing input + +desc = '''HGRAD Dirichlet BC verification - polynomial solution + T = x*y on all boundaries with zero source + Solution is exactly representable with p=1 elements + Verifies inhomogeneous Dirichlet BC enforcement for HGRAD basis + ''' + +its = mrhyde_test_support(desc) + +print('Because of the diff test on the log file, this test needs ') +print('to run with "-v". There is a buffering issue.') +print('Setting the verbosity to True.') +its.opts.verbose = True + +#------------------------------------------------------------------------------- +# Problem Parameters + +root = 'mrhyde' + +# These comments are for testing with the runtest.py utility. +#TESTING active +#TESTING -n 1 +#TESTING -k thermal,HGRAD,Dirichlet,steady-state + +# ============================================================================== +status = 0 + +status += its.call('mpiexec -n 1 ../../mrhyde >& mrhyde.log') +status += its.clean_log() +status += its.call('diff -y %s.log %s.gold' % (root, root)) + +# ============================================================================== +if status == 0: print('Success.') +else: print('Failure.') +sys.exit(status) diff --git a/regression/boundary_conditions/pec_3d_cavity/input.yaml b/regression/boundary_conditions/pec_3d_cavity/input.yaml new file mode 100644 index 00000000..bb91f708 --- /dev/null +++ b/regression/boundary_conditions/pec_3d_cavity/input.yaml @@ -0,0 +1,109 @@ +%YAML 1.1 +--- +# Test: 3D PEC Cavity TM_110 Mode - Energy Conservation +# +# Tests Perfect Electric Conductor (PEC) boundary conditions in 3D: n x E = 0. +# Uses TM_110 mode which naturally satisfies PEC on all walls. +# +# Domain: [0,1]x[0,1]x[0,0.5] +# Mode: TM_110 with Ez = sin(pi*x)*sin(pi*y), Ex = Ey = 0 +# +# PEC verification via energy conservation: +# - In a lossless PEC cavity, total EM energy U = U_E + U_B is constant (within numerical tolerance) +# - Paraview (qualitative test - should be near machine precision) +# - On x=0,1: max|Ey|, max|Ez| +# - On y=0,1: max|Ex|, max|Ez| +# - On z=0,0.5: max|Ex|, max|Ey| + + +ANONYMOUS: + debug level: 0 + verbosity: 5 + + Mesh: + dimension: 3 + shape: hex + xmin: 0.0 + xmax: 1.0 + ymin: 0.0 + ymax: 1.0 + zmin: 0.0 + zmax: 0.5 + NX: 10 + NY: 10 + NZ: 4 + blocknames: eblock-0_0_0 + + Physics: + Dirichlet conditions: + scalar data: false + eblock-0_0_0: + modules: maxwell + assemble boundary terms: false + Dirichlet conditions: + # PEC on all boundaries: tangential E = 0 + # For 3D, we specify all components as 0 on all boundaries + # HCURL basis will enforce only tangential components + Ex: + all boundaries: '0.0' + Ey: + all boundaries: '0.0' + Ez: + all boundaries: '0.0' + Initial conditions: + # TM_110 mode at t=0: Ez = sin(pi*x)*sin(pi*y), B = 0 + 'E[x]': '0.0' + 'E[y]': '0.0' + 'E[z]': 'sin(pi*x)*sin(pi*y)' + 'B[x]': '0.0' + 'B[y]': '0.0' + 'B[z]': '0.0' + + Discretization: + eblock-0_0_0: + order: + E: 1 + B: 1 + quadrature: 2 + + Functions: + current x: '0.0' + current y: '0.0' + current z: '0.0' + permittivity: '1.0' + permeability: '1.0' + + Solver: + solver: transient + transient BDF order: 1 + transient Butcher tableau: CN + workset size: 100 + nonlinear TOL: 1.0e-12 + linear TOL: 1.0e-14 + max nonlinear iters: 3 + max linear iters: 500 + final time: 0.50 + number of steps: 10 + initial type: L2-projection + use preconditioner: true + use direct solver: true + + Analysis: + analysis type: forward + + Postprocess: + write solution: false + compute errors: false + compute integrated quantities: true + eblock-0_0_0: + Integrated quantities: + electric_energy: + integrand: '0.5*(E[x]^2 + E[y]^2 + E[z]^2)' + location: 'volume' + magnetic_energy: + integrand: '0.5*(B[x]^2 + B[y]^2 + B[z]^2)' + location: 'volume' + total_energy: + integrand: '0.5*(E[x]^2 + E[y]^2 + E[z]^2) + 0.5*(B[x]^2 + B[y]^2 + B[z]^2)' + location: 'volume' +... diff --git a/regression/boundary_conditions/pec_3d_cavity/mrhyde.gold b/regression/boundary_conditions/pec_3d_cavity/mrhyde.gold new file mode 100644 index 00000000..84367215 --- /dev/null +++ b/regression/boundary_conditions/pec_3d_cavity/mrhyde.gold @@ -0,0 +1,392 @@ +p=0 | CubeHexMesh: Building sub cells +STK Meta data: + Element blocks = "eblock-0_0_0" + Sidesets = "back" "bottom" "front" "left" "right" "top" + Nodesets = "origin" + Fields = "coordinates" "edges" "faces" "PROC_ID" "LOAD_BAL" +DOFManager Field Information: + Element Block = eblock-0_0_0 + "E" is field ID 0 + "B" is field ID 1 +Time = 0, Integrated quantities: electric_energy = 0.0624974 magnetic_energy = 0 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 1 +**** Current time is 0 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.23553 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 3.86524e-16 +***** Scaled Norm of nonlinear residual: 3.1284e-16 +********************************************************* +Norm of solution: 0.635362 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.23553 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 1.15957e-15 +***** Scaled Norm of nonlinear residual: 9.3852e-16 +********************************************************* +Norm of solution: 0.61975 (overlapped vector so results may differ on multiple procs) +Time = 0.05, Integrated quantities: electric_energy = 0.0594637 magnetic_energy = 0.00303363 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 2 +**** Current time is 0.05 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.20538 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 4.833e-16 +***** Scaled Norm of nonlinear residual: 4.00953e-16 +********************************************************* +Norm of solution: 0.604217 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.17522 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 9.77243e-16 +***** Scaled Norm of nonlinear residual: 8.31537e-16 +********************************************************* +Norm of solution: 0.589371 (overlapped vector so results may differ on multiple procs) +Time = 0.1, Integrated quantities: electric_energy = 0.0509518 magnetic_energy = 0.0115455 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 3 +**** Current time is 0.1 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.11618 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 5.97926e-16 +***** Scaled Norm of nonlinear residual: 5.35689e-16 +********************************************************* +Norm of solution: 0.543473 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 1.05714 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 1.05936e-15 +***** Scaled Norm of nonlinear residual: 1.0021e-15 +********************************************************* +Norm of solution: 0.530118 (overlapped vector so results may differ on multiple procs) +Time = 0.15, Integrated quantities: electric_energy = 0.0386144 magnetic_energy = 0.023883 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 4 +**** Current time is 0.15 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.971932 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 6.07311e-16 +***** Scaled Norm of nonlinear residual: 6.24849e-16 +********************************************************* +Norm of solution: 0.456149 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.886727 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.57801e-16 +***** Scaled Norm of nonlinear residual: 9.6738e-16 +********************************************************* +Norm of solution: 0.444941 (overlapped vector so results may differ on multiple procs) +Time = 0.2, Integrated quantities: electric_energy = 0.0248468 magnetic_energy = 0.0376506 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 5 +**** Current time is 0.2 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.7796 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 6.97846e-16 +***** Scaled Norm of nonlinear residual: 8.95134e-16 +********************************************************* +Norm of solution: 0.346635 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.672473 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 7.96678e-16 +***** Scaled Norm of nonlinear residual: 1.1847e-15 +********************************************************* +Norm of solution: 0.338118 (overlapped vector so results may differ on multiple procs) +Time = 0.25, Integrated quantities: electric_energy = 0.0123222 magnetic_energy = 0.0501752 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 6 +**** Current time is 0.25 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.548854 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.04636e-16 +***** Scaled Norm of nonlinear residual: 1.46603e-15 +********************************************************* +Norm of solution: 0.220645 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.425234 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 6.64141e-16 +***** Scaled Norm of nonlinear residual: 1.56182e-15 +********************************************************* +Norm of solution: 0.215223 (overlapped vector so results may differ on multiple procs) +Time = 0.3, Integrated quantities: electric_energy = 0.00347238 magnetic_energy = 0.059025 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 7 +**** Current time is 0.3 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.29126 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.57195e-16 +***** Scaled Norm of nonlinear residual: 2.94306e-15 +********************************************************* +Norm of solution: 0.0867944 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.157286 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 7.48874e-16 +***** Scaled Norm of nonlinear residual: 4.76122e-15 +********************************************************* +Norm of solution: 0.0846616 (overlapped vector so results may differ on multiple procs) +Time = 0.35, Integrated quantities: electric_energy = 1.56227e-05 magnetic_energy = 0.0624817 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 8 +**** Current time is 0.35 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.0877938 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 9.15082e-16 +***** Scaled Norm of nonlinear residual: 1.04231e-14 +********************************************************* +Norm of solution: 0.0685084 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.118336 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.37796e-16 +***** Scaled Norm of nonlinear residual: 7.07978e-15 +********************************************************* +Norm of solution: 0.066825 (overlapped vector so results may differ on multiple procs) +Time = 0.4, Integrated quantities: electric_energy = 0.00262309 magnetic_energy = 0.0598743 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 9 +**** Current time is 0.4 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.25327 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.98315e-16 +***** Scaled Norm of nonlinear residual: 3.54687e-15 +********************************************************* +Norm of solution: 0.201782 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.388203 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 9.12117e-16 +***** Scaled Norm of nonlinear residual: 2.34959e-15 +********************************************************* +Norm of solution: 0.196824 (overlapped vector so results may differ on multiple procs) +Time = 0.45, Integrated quantities: electric_energy = 0.0107885 magnetic_energy = 0.0517088 total_energy = 0.0624974 + + +******************************************************* + +**** Beginning Time Step 10 +**** Current time is 0.45 + +******************************************************* + + + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.513598 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 6.16634e-16 +***** Scaled Norm of nonlinear residual: 1.20062e-15 +********************************************************* +Norm of solution: 0.329555 (overlapped vector so results may differ on multiple procs) + +********************************************************* +***** Iteration: 0 +***** Norm of nonlinear residual: 0.638993 +***** Scaled Norm of nonlinear residual: 1 +********************************************************* + +********************************************************* +***** Iteration: 1 +***** Norm of nonlinear residual: 8.90771e-16 +***** Scaled Norm of nonlinear residual: 1.39402e-15 +********************************************************* +Norm of solution: 0.321457 (overlapped vector so results may differ on multiple procs) +Time = 0.5, Integrated quantities: electric_energy = 0.0229265 magnetic_energy = 0.0395708 total_energy = 0.0624974 + +********************************************************* +****** Storing Integrated Quantities ****** +********************************************************* + +********************************************************* +****** Integrated Quantities on block : eblock-0_0_0 ****** +********************************************************* +electric_energy : 0.0229265 +magnetic_energy : 0.0395708 +total_energy : 0.0624974 diff --git a/regression/boundary_conditions/pec_3d_cavity/mrhyde.tst b/regression/boundary_conditions/pec_3d_cavity/mrhyde.tst new file mode 100755 index 00000000..96affc0d --- /dev/null +++ b/regression/boundary_conditions/pec_3d_cavity/mrhyde.tst @@ -0,0 +1,43 @@ +#!/usr/bin/env python3 +#------------------------------------------------------------------------------- + +import sys, os +sys.path.append("../../scripts") +from mrhyde_test_support import * + +# ============================================================================== +# Parsing input + +desc = '''3D PEC Cavity TM_110 Mode Test + Verifies 3D PEC (n x E = 0) on all 6 faces using TM_110 mode. + Ez = sin(pi*x)*sin(pi*y), Ex = Ey = 0, all B components = 0 + ''' + +its = mrhyde_test_support(desc) + +print('Because of the diff test on the log file, this test needs ') +print('to run with "-v". There is a buffering issue.') +print('Setting the verbosity to True.') +its.opts.verbose = True + +#------------------------------------------------------------------------------- +# Problem Parameters + +root = 'mrhyde' + +# These comments are for testing with the runtest.py utility. +#TESTING active +#TESTING -n 1 +#TESTING -k maxwell,HCURL,PEC,3D,cavity,TM110 + +# ============================================================================== +status = 0 + +status += its.call('mpiexec -n 1 ../../mrhyde >& mrhyde.log') +status += its.clean_log() +status += its.call('diff -y %s.log %s.gold' % (root, root)) + +# ============================================================================== +if status == 0: print('Success.') +else: print('Failure.') +sys.exit(status) diff --git a/regression/hpc/scripts/parse_mrhyde_log.py b/regression/hpc/scripts/parse_mrhyde_log.py new file mode 100644 index 00000000..907ca307 --- /dev/null +++ b/regression/hpc/scripts/parse_mrhyde_log.py @@ -0,0 +1,159 @@ +#!/usr/bin/env python3 +""" +Parse MrHyDE log files to extract processor workload and timing tables. +""" + +import re +import pandas as pd + + +def parse_workload(text): + """Extract processor workload info into a DataFrame.""" + + # patterns for each data type + elem_pat = r"Processor (\d+) has (\d+) elements" + belem_pat = r"Processor (\d+) has (\d+) boundary elements" + vol_pat = r"Processor (\d+) is using ([\d.]+) MB to store volumetric data" + face_pat = r"Processor (\d+) is using ([\d.]+) MB to store face data" + bnd_pat = r"Processor (\d+) is using ([\d.]+) MB to store boundary data" + + # data per processor + data = {} + for match in re.finditer(elem_pat, text): + proc = int(match.group(1)) + data.setdefault(proc, {})["elements"] = int(match.group(2)) + + for match in re.finditer(belem_pat, text): + proc = int(match.group(1)) + data.setdefault(proc, {})["boundary_elements"] = int(match.group(2)) + + for match in re.finditer(vol_pat, text): + proc = int(match.group(1)) + data.setdefault(proc, {})["volumetric_MB"] = float(match.group(2)) + + for match in re.finditer(face_pat, text): + proc = int(match.group(1)) + data.setdefault(proc, {})["face_MB"] = float(match.group(2)) + + for match in re.finditer(bnd_pat, text): + proc = int(match.group(1)) + data.setdefault(proc, {})["boundary_MB"] = float(match.group(2)) + + if not data: + return pd.DataFrame() + + rows = [] + for proc in sorted(data.keys()): + row = {"processor": proc} + row.update(data[proc]) + rows.append(row) + + return pd.DataFrame(rows) + + +def parse_timing(text): + """Extract TimeMonitor table into a DataFrame.""" + + # Find the TimeMonitor section + section_pat = r"={50,}\n(.*?TimeMonitor.*?)\n={50,}" + match = re.search(section_pat, text, re.DOTALL) + + if not match: + return pd.DataFrame() + + section = match.group(0) + lines = section.strip().split("\n") + + # format check: single-processor vs multi-processor + is_single_proc = "Global time" in section + + data_lines = [] + for line in lines: + if line.startswith("MrHyDE::") or line.startswith("STK_") or \ + line.startswith("panzer") or line.startswith("Solver") or \ + line.startswith("Time ") or line.startswith("Total Time"): + data_lines.append(line) + + if not data_lines: + return pd.DataFrame() + + def parse_value_calls(s): + s = s.strip() + m = re.match(r"([\d.e+-]+)\s*\((\d+\.?\d*)\)", s) + if m: + return float(m.group(1)), float(m.group(2)) + return 0.0, 0 + + rows = [] + for line in data_lines: + parts = re.split(r"\s{2,}", line.strip()) + + if is_single_proc: + if len(parts) < 2: + continue + timer_name = parts[0] + val, calls = parse_value_calls(parts[1]) + rows.append({ + "timer_name": timer_name, + "GlobalTime": val, + "GlobalTime_calls": int(calls), + }) + else: + if len(parts) < 5: + continue + timer_name = parts[0] + min_val, min_calls = parse_value_calls(parts[1]) + mean_val, mean_calls = parse_value_calls(parts[2]) + max_val, max_calls = parse_value_calls(parts[3]) + mean_per_call, mean_per_call_n = parse_value_calls(parts[4]) + rows.append({ + "timer_name": timer_name, + "MinOverProcs": min_val, + "MinOverProcs_calls": int(min_calls), + "MeanOverProcs": mean_val, + "MeanOverProcs_calls": int(mean_calls), + "MaxOverProcs": max_val, + "MaxOverProcs_calls": int(max_calls), + "MeanOverCallCounts": mean_per_call, + "MeanOverCallCounts_calls": int(mean_per_call_n), + }) + + return pd.DataFrame(rows) + + +def parse_log(filepath): + """ + Parse a MrHyDE log file. + + Returns: + tuple: (workload_df, timing_df) + """ + with open(filepath, "r") as f: + text = f.read() + + workload_df = parse_workload(text) + timing_df = parse_timing(text) + + return workload_df, timing_df + + +if __name__ == "__main__": + import sys + + if len(sys.argv) < 2: + print("Usage: python parse_mrhyde_log.py ") + sys.exit(1) + + logfile = sys.argv[1] + workload, timing = parse_log(logfile) + + print("=" * 60) + print("PROCESSOR WORKLOAD") + print("=" * 60) + print(workload.to_string(index=False)) + + print("\n") + print("=" * 60) + print("TIMING DATA") + print("=" * 60) + print(timing.to_string(index=False)) diff --git a/regression/hpc/scripts/scaling_analysis.py b/regression/hpc/scripts/scaling_analysis.py new file mode 100644 index 00000000..1fe99a04 --- /dev/null +++ b/regression/hpc/scripts/scaling_analysis.py @@ -0,0 +1,360 @@ +#!/usr/bin/env python3 +""" +Strong scaling analysis for MrHyDE HPC regression testing. + +Analyzes how solver routines scale with increasing MPI tasks. +Provides pass/fail assessment based on qualitative efficiency thresholds. +""" + +import sys +from pathlib import Path + +import numpy as np +import pandas as pd + +from parse_mrhyde_log import parse_log + + +def get_num_procs(workload_df): + """Infer number of processors from workload DataFrame.""" + if workload_df.empty: + return 1 + return workload_df["processor"].nunique() + + +def get_time_column(timing_df): + """Determine which time column to use based on available columns.""" + if "GlobalTime" in timing_df.columns: + return "GlobalTime" + elif "MeanOverProcs" in timing_df.columns: + return "MeanOverProcs" + else: + raise ValueError("No recognized time column in timing data") + + +def load_all_logs(log_files): + """ + Load timing data from all log files. + + Returns: + dict: {num_procs: timing_df} mapping processor count to timing data + """ + results = {} + for filepath in log_files: + workload_df, timing_df = parse_log(filepath) + num_procs = get_num_procs(workload_df) + + if timing_df.empty: + print(f"Warning: No timing data in {filepath}, skipping") + continue + + time_col = get_time_column(timing_df) + timing_df = timing_df[["timer_name", time_col]].copy() + timing_df.columns = ["timer_name", "time"] + + timing_df = timing_df[timing_df["time"] > 0] + + if num_procs in results: + print(f"Warning: Duplicate processor count {num_procs}, " + f"overwriting with {filepath}") + results[num_procs] = timing_df + + return results + + +def get_total_time(timing_df): + """Get total simulation time from timing DataFrame.""" + total_row = timing_df[timing_df["timer_name"].str.contains("Total Time|MrHyDE::driver", case=False)] + if not total_row.empty: + return total_row["time"].values[0] + # Fallback: sum of top-level timers + return timing_df["time"].max() + + +def get_all_timer_names(all_data): + """Union of timer names across all logs.""" + names = set() + for timing_df in all_data.values(): + names.update(timing_df["timer_name"].tolist()) + return sorted(names) + + +def get_top_routines(timing_table, fraction=0.25, baseline_col=None): + """ + Return routine names that account for the top fraction of time at baseline. + + Args: + timing_table: DataFrame index=timer_name, columns=proc counts + fraction: take top 25% by count (0.25) + baseline_col: column for ranking (default: first column) + + Returns: + list of timer_name (top fraction by time at baseline) + """ + if baseline_col is None: + baseline_col = timing_table.columns[0] + series = timing_table[baseline_col].dropna().sort_values(ascending=False) + n = max(1, int(np.ceil(fraction * len(series)))) + return series.head(n).index.tolist() + + +def build_timing_table(all_data, timer_names=None): + """ + Build a table with routines as rows and processor counts as columns. + + Args: + all_data: dict mapping num_procs to timing DataFrame + timer_names: list of timer names to include (None = use total time only) + + Returns: + pd.DataFrame with timer_name index and proc counts as columns + """ + proc_counts = sorted(all_data.keys()) + + if timer_names is None: + # Just use total time + timer_names = ["Total"] + rows = [{"timer_name": "Total"}] + for nprocs in proc_counts: + rows[0][nprocs] = get_total_time(all_data[nprocs]) + else: + rows = [] + for timer in timer_names: + row = {"timer_name": timer} + for nprocs in proc_counts: + df = all_data[nprocs] + match = df[df["timer_name"] == timer] + if not match.empty: + row[nprocs] = match["time"].values[0] + else: + row[nprocs] = np.nan + rows.append(row) + + table = pd.DataFrame(rows).set_index("timer_name") + return table + + +def write_scaling_table(timing_table, routine_names, savepath): + """ + Write time vs nprocs scaling data to a text file. + + Args: + timing_table: DataFrame index=timer_name, columns=proc counts + routine_names: list of timer_name to include + savepath: path to save text file + """ + routine_names = [r for r in routine_names if r in timing_table.index] + if not routine_names: + return + + proc_counts = timing_table.columns.tolist() + baseline_col = proc_counts[0] + baseline_procs = baseline_col + + with open(savepath, 'w') as f: + f.write("=" * 80 + "\n") + f.write("Top Routines Scaling Analysis\n") + f.write("=" * 80 + "\n\n") + + for name in routine_names: + row = timing_table.loc[name] + times = [row[c] if not np.isnan(row[c]) else None for c in proc_counts] + baseline_time = times[0] + + if baseline_time is None or baseline_time <= 0: + continue + + f.write(f"{name}\n") + f.write("-" * 80 + "\n") + f.write(f"{'MPI procs':<12} {'Time (s)':<15} {'Speedup':<15} {'Efficiency (%)':<15}\n") + f.write("-" * 80 + "\n") + + for i, nprocs in enumerate(proc_counts): + time_val = times[i] + if time_val is None or time_val <= 0: + f.write(f"{nprocs:<12} {'N/A':<15} {'N/A':<15} {'N/A':<15}\n") + else: + speedup = baseline_time / time_val + scale = nprocs / baseline_procs + efficiency = (speedup / scale) * 100 + f.write(f"{nprocs:<12} {time_val:<15.4f} {speedup:<15.4f} {efficiency:<15.2f}\n") + + f.write("\n") + + f.write("=" * 80 + "\n") + + +def compute_speedup(timing_table): + """Compute speedup relative to smallest processor count.""" + baseline_col = timing_table.columns[0] + speedup = timing_table.copy() + for col in timing_table.columns: + speedup[col] = timing_table[baseline_col] / timing_table[col] + return speedup + + +def compute_efficiency(timing_table): + """Compute parallel efficiency: T1 / (N * TN) * 100.""" + baseline_col = timing_table.columns[0] + baseline_procs = baseline_col + efficiency = timing_table.copy() + for col in timing_table.columns: + nprocs = col + scale = nprocs / baseline_procs + efficiency[col] = (timing_table[baseline_col] / + (scale * timing_table[col])) * 100 + return efficiency + + +def check_scaling_quality(timing_table, min_efficiency=50, verbose=True): + """ + Check if scaling is acceptable. + + Args: + timing_table: DataFrame with timing data + min_efficiency: minimum acceptable parallel efficiency (%) + verbose: print detailed output + + Returns: + (passed, message) tuple + """ + speedup_table = compute_speedup(timing_table) + efficiency_table = compute_efficiency(timing_table) + + proc_counts = list(timing_table.columns) + baseline_procs = proc_counts[0] + + if verbose: + print("\nStrong Scaling Analysis") + print("=" * 60) + print(f"Processor counts: {proc_counts}") + print("\nTiming Summary:") + + all_pass = True + messages = [] + + for timer in timing_table.index: + if verbose: + print(f"\n {timer}:") + + for i, nprocs in enumerate(proc_counts): + time_val = timing_table.loc[timer, nprocs] + speedup = speedup_table.loc[timer, nprocs] + eff = efficiency_table.loc[timer, nprocs] + + if np.isnan(time_val): + if verbose: + print(f" {nprocs} procs: N/A") + continue + + if i == 0: + if verbose: + print(f" {nprocs} proc(s): {time_val:.2f}s (baseline)") + else: + if verbose: + print(f" {nprocs} procs: {time_val:.2f}s " + f"(speedup {speedup:.2f}x, efficiency {eff:.0f}%)") + + # Check efficiency threshold (skip baseline) + if eff < min_efficiency: + all_pass = False + messages.append( + f"{timer} at {nprocs} procs: efficiency {eff:.0f}% < {min_efficiency}%" + ) + + if verbose: + print("\n" + "=" * 60) + if all_pass: + print("Assessment: PASS") + print(" Scaling is within acceptable range.") + else: + print("Assessment: FAIL") + for msg in messages: + print(f" - {msg}") + + return all_pass, messages + + +def analyze_scaling(log_files, min_efficiency=50, verbose=True, + top_fraction=None, output_path=None): + """ + Main entry point for scaling analysis. + + Args: + log_files: list of log file paths + min_efficiency: minimum acceptable parallel efficiency (%) + verbose: print detailed output + top_fraction: if set (e.g. 0.25), write side-by-side scaling table for + top fraction of most consuming routines + output_path: path to save the top-routines text file (used if + top_fraction is set) + + Returns: + 0 if passed, 1 if failed + """ + if verbose: + print(f"Loading {len(log_files)} log files...") + + all_data = load_all_logs(log_files) + + if len(all_data) < 2: + print("Error: Need at least 2 valid log files for scaling analysis") + return 1 + + if verbose: + print(f"Found processor counts: {sorted(all_data.keys())}") + + # Pass/fail uses total time only + timing_table_total = build_timing_table(all_data) + passed, messages = check_scaling_quality(timing_table_total, min_efficiency, verbose) + + # Full timing DataFrame and optional top-fraction output + if top_fraction is not None and output_path and 0 < top_fraction <= 1: + all_timers = get_all_timer_names(all_data) + full_timing_table = build_timing_table(all_data, timer_names=all_timers) + top_routines = get_top_routines(full_timing_table, fraction=top_fraction) + if top_routines and verbose: + print(f"\nTop {int(100 * top_fraction)}% routines by time: {len(top_routines)}") + write_scaling_table(full_timing_table, top_routines, output_path) + if verbose and output_path: + print(f"Scaling table saved to {output_path}") + + return 0 if passed else 1 + + +if __name__ == "__main__": + import argparse + + parser = argparse.ArgumentParser( + description="Strong scaling analysis for MrHyDE HPC regression testing" + ) + parser.add_argument( + "logfiles", + nargs="+", + help="Log files to analyze (one per processor count)" + ) + parser.add_argument( + "--min-efficiency", + type=float, + default=50, + help="Minimum acceptable parallel efficiency in %% (default: 50)" + ) + parser.add_argument( + "-q", "--quiet", + action="store_true", + help="Quiet mode (only print pass/fail)" + ) + + args = parser.parse_args() + + for f in args.logfiles: + if not Path(f).exists(): + print(f"Error: File not found: {f}") + sys.exit(1) + + status = analyze_scaling( + args.logfiles, + min_efficiency=args.min_efficiency, + verbose=not args.quiet + ) + sys.exit(status) diff --git a/regression/hpc/strong_scaling/maxwell_cavity/expected_scaling.yaml b/regression/hpc/strong_scaling/maxwell_cavity/expected_scaling.yaml new file mode 100644 index 00000000..7b0451b8 --- /dev/null +++ b/regression/hpc/strong_scaling/maxwell_cavity/expected_scaling.yaml @@ -0,0 +1,50 @@ +# Expected Scaling Behavior for Maxwell Cavity Strong Scaling Test +# +# This file documents qualitative expectations for the strong scaling test. +# These are NOT hard thresholds - they provide guidance on what "reasonable" +# scaling looks like for this problem. + +# Minimum parallel efficiency threshold (percent) +# The test will FAIL if efficiency drops below this value. +# This is intentionally set low to accommodate various hardware configurations. +minimum_parallel_efficiency: 50 + +# Problem characteristics +problem: + type: "Maxwell 3D cavity (TM_110 mode)" + mesh_size: "20x20x10 hexahedral elements (4000 total)" + solver: "Iterative with domain decomposition preconditioner" + time_steps: 20 + +# Expected behavior notes +notes: | + Strong scaling efficiency depends heavily on: + + 1. Problem Size: + - Small problems have less parallel work, leading to lower efficiency + - This test uses 4000 elements which is modest; efficiency may be lower + than production-scale runs + + 2. Hardware Configuration: + - Shared-memory systems typically show 70-90% efficiency at low core counts + - Distributed memory systems may show lower efficiency due to communication + - NUMA effects can reduce efficiency on multi-socket nodes + + 3. Solver Characteristics: + - Iterative solvers with domain decomposition depend on subdomain coupling + - Communication overhead increases with processor count + - Preconditioner setup cost may not scale perfectly + + 4. What to expect: + - 1 to 2 procs: 80-95% efficiency (minimal communication) + - 2 to 4 procs: 60-85% efficiency (increasing communication) + - Above 4 procs: efficiency depends strongly on problem size + + If this test fails: + - First check if the run completed successfully + - Review timing breakdown to identify bottlenecks + - Consider if hardware issues (oversubscription, memory pressure) affected results + - For unusual failures, run manually and inspect log files + +# Processor counts tested +processor_counts: [1, 2, 4] diff --git a/regression/hpc/strong_scaling/maxwell_cavity/input.yaml b/regression/hpc/strong_scaling/maxwell_cavity/input.yaml new file mode 100644 index 00000000..fe8d6459 --- /dev/null +++ b/regression/hpc/strong_scaling/maxwell_cavity/input.yaml @@ -0,0 +1,120 @@ +%YAML 1.1 +--- +# HPC Strong Scaling Test: 3D Maxwell Cavity with Iterative Solver +# +# This test is designed for strong scaling analysis. It uses: +# - A larger mesh (40x40x20) for meaningful parallel work +# - DIRK-1,2 with Ifpack2 Jacobi preconditioner and BiCGStab (aligned with cavity-pot-mono implicit) +# - More time steps to increase solve time +# +# The problem is a TM_110 cavity mode which naturally satisfies PEC +# boundary conditions on all walls. + +ANONYMOUS: + debug level: 0 + verbosity: 10 + + Mesh: + dimension: 3 + shape: hex + xmin: 0.0 + xmax: 1.0 + ymin: 0.0 + ymax: 1.0 + zmin: 0.0 + zmax: 0.5 + NX: 80 + NY: 80 + NZ: 80 + blocknames: eblock-0_0_0 + + Physics: + Dirichlet conditions: + scalar data: false + eblock-0_0_0: + modules: maxwell + assemble boundary terms: false + Dirichlet conditions: + # PEC on all boundaries: tangential E = 0 + Ex: + all boundaries: '0.0' + Ey: + all boundaries: '0.0' + Ez: + all boundaries: '0.0' + Initial conditions: + # TM_110 mode at t=0: Ez = sin(pi*x)*sin(pi*y), B = 0 + 'E[x]': '0.0' + 'E[y]': '0.0' + 'E[z]': 'sin(pi*x)*sin(pi*y)' + 'B[x]': '0.0' + 'B[y]': '0.0' + 'B[z]': '0.0' + + Discretization: + eblock-0_0_0: + order: + E: 1 + B: 1 + quadrature: 2 + + Functions: + current x: '0.0' + current y: '0.0' + current z: '0.0' + permittivity: '1.0' + permeability: '1.0' + + Solver: + solver: transient + transient BDF order: 1 + transient Butcher tableau: DIRK-1,2 + workset size: 100 + final time: 0.50 + number of steps: 20 + initial type: L2-projection + assembly partitioning: sequential + use basis database: false + use mass database: false + database TOL: 1.0e-13 + linear TOL: 1.0e-8 + max linear iters: 500 + nonlinear TOL: 1.0e-07 + max nonlinear iters: 1 + use direct solver: false + use preconditioner: true + preconditioner type: Ifpack2 + preconditioner variant: RELAXATION + right preconditioner: true + reuse preconditioner: true + reuse Jacobian: true + preconditioner reuse type: full + Belos implicit residual scaling: Norm of Initial Residual + storage proportion: 1.0 + Belos solver: BiCGStab + enable autotune: false + Preconditioner Settings: + "relaxation: type": "Jacobi" + "relaxation: damping factor": 0.5 + "relaxation: backward mode": false + "relaxation: sweeps": 2 + + Analysis: + analysis type: forward + + Postprocess: + write solution: false + compute errors: false + compute integrated quantities: true + eblock-0_0_0: + Integrated quantities: + electric_energy: + integrand: '0.5*(E[x]^2 + E[y]^2 + E[z]^2)' + location: 'volume' + magnetic_energy: + integrand: '0.5*(B[x]^2 + B[y]^2 + B[z]^2)' + location: 'volume' + total_energy: + integrand: '0.5*(E[x]^2 + E[y]^2 + E[z]^2) + 0.5*(B[x]^2 + B[y]^2 + B[z]^2)' + location: 'volume' +... diff --git a/regression/hpc/strong_scaling/maxwell_cavity/mrhyde.tst b/regression/hpc/strong_scaling/maxwell_cavity/mrhyde.tst new file mode 100755 index 00000000..f4f2d71f --- /dev/null +++ b/regression/hpc/strong_scaling/maxwell_cavity/mrhyde.tst @@ -0,0 +1,140 @@ +#!/usr/bin/env python3 +""" +HPC Strong Scaling Test: Maxwell 3D Cavity + +Runs the same problem with 1, 2, and 4 MPI processes and verifies +that parallel efficiency is within acceptable bounds. + +This test only runs when explicitly requested with -k hpc. +""" + +import sys +import os + +# Get the directory containing this test file +test_dir = os.path.dirname(os.path.abspath(__file__)) +# hpc/scripts is two levels up from maxwell_cavity +hpc_scripts_dir = os.path.join(test_dir, "..", "..", "scripts") +hpc_scripts_abs = os.path.abspath(hpc_scripts_dir) +sys.path.insert(0, hpc_scripts_abs) # Use insert(0) to prioritize this path + +# Also add regression/scripts for mrhyde_test_support +regression_scripts_dir = os.path.join(test_dir, "..", "..", "..", "scripts") +regression_scripts_abs = os.path.abspath(regression_scripts_dir) +sys.path.insert(0, regression_scripts_abs) + +from mrhyde_test_support import * + +# Test metadata parsed by runtests.py +#TESTING active +#TESTING -n 4 +#TESTING -k hpc,strong_scaling,maxwell + +desc = '''HPC Strong Scaling Test: Maxwell 3D Cavity''' +its = mrhyde_test_support(desc) +its.opts.verbose = True + +# Configuration +PROC_COUNTS = [1, 2, 4] +MIN_EFFICIENCY = 50 # percent + +status = 0 + +print("=" * 70) +print("HPC Strong Scaling Test: Maxwell 3D Cavity") +print("=" * 70) + +# Run MrHyDE with different processor counts +log_files = [] +for nprocs in PROC_COUNTS: + logfile = f"mrhyde_np{nprocs}.log" + log_files.append(logfile) + print(f"\nRunning with {nprocs} MPI process(es)...") + run_status = its.call(f'mpiexec -n {nprocs} ../../../mrhyde >& {logfile}') + if run_status != 0: + print(f" ERROR: MrHyDE failed with {nprocs} processes") + status += 1 + +if status != 0: + print("\nFailed to run MrHyDE, cannot perform scaling analysis") + sys.exit(status) + +# Perform scaling analysis +print("\n" + "=" * 70) +print("Scaling Analysis") +print("=" * 70) + +try: + print(f"Attempting to import scaling_analysis from: {hpc_scripts_abs}") + from scaling_analysis import analyze_scaling + print(f"Successfully imported scaling_analysis") + print(f"Calling analyze_scaling with {len(log_files)} log files...") + analysis_status = analyze_scaling( + log_files, + min_efficiency=MIN_EFFICIENCY, + top_fraction=0.50, + output_path="scaling_top50.txt", + ) + print(f"Scaling analysis completed with status: {analysis_status}") + status += analysis_status +except (ImportError, ModuleNotFoundError) as e: + print(f"\nWarning: Could not import scaling_analysis: {e}") + print(f" Searched in: {hpc_scripts_abs}") + print("Falling back to simplified analysis...") + # Fallback: simple timing extraction without pandas + import re + + def extract_total_time(logfile): + """Extract total time from log file.""" + with open(logfile, 'r') as f: + text = f.read() + # Look for Total Time in TimeMonitor output + match = re.search(r"Total Time[:\s]+([\d.]+)", text) + if match: + return float(match.group(1)) + # Fallback: look for MrHyDE::driver time + match = re.search(r"MrHyDE::driver\s+([\d.]+)", text) + if match: + return float(match.group(1)) + return None + + times = {} + for i, nprocs in enumerate(PROC_COUNTS): + t = extract_total_time(log_files[i]) + if t: + times[nprocs] = t + + if len(times) < 2: + print("ERROR: Could not extract timing data from log files") + status = 1 + else: + baseline_procs = PROC_COUNTS[0] + baseline_time = times[baseline_procs] + print(f"\nTiming Summary:") + print(f" {baseline_procs} proc(s): {baseline_time:.2f}s (baseline)") + + for nprocs in PROC_COUNTS[1:]: + if nprocs in times: + t = times[nprocs] + speedup = baseline_time / t + scale = nprocs / baseline_procs + efficiency = (speedup / scale) * 100 + print(f" {nprocs} procs: {t:.2f}s " + f"(speedup {speedup:.2f}x, efficiency {efficiency:.0f}%)") + if efficiency < MIN_EFFICIENCY: + print(f" WARNING: Efficiency below threshold ({MIN_EFFICIENCY}%)") + status = 1 +except Exception as e: + print(f"\nError during scaling analysis: {e}") + import traceback + traceback.print_exc() + status = 1 + +print("\n" + "=" * 70) +if status == 0: + print("Success: Strong scaling test PASSED") +else: + print("Failure: Strong scaling test FAILED") +print("=" * 70) + +sys.exit(status) diff --git a/regression/scripts/runtests.py b/regression/scripts/runtests.py index 678c7f32..0ea92b97 100755 --- a/regression/scripts/runtests.py +++ b/regression/scripts/runtests.py @@ -115,12 +115,12 @@ def __init__(self, fullpath, startingDir): hostname = hostname.rstrip('1234567890') if not os.path.exists(self.fullpath): - print('%s does not exist!') %(self.fullpath) + print('%s does not exist!' % self.fullpath) self.status = 1 self.statusStr = '!exist' self.skipped = 1 elif not os.access(self.fullpath, os.X_OK): - print('%s not executable!') %(self.fullpath) + print('%s not executable!' % self.fullpath) self.status = 1 self.statusStr = '!exec' self.skipped = 1 @@ -340,7 +340,7 @@ def testOutput(self,test,runtime): % (test.index+1, self.list_length, \ test.statusStr, runtime, test.nprocs, fname) print(stmt2) - + #=============================================================================== @@ -369,8 +369,8 @@ def serial_testing(opts,listOfTests,doc): # default is 32 bit so only add if 64 bit if opts.mode_64: cmd += ['--64'] if test.test_args: cmd += shlex.split(test.test_args) - if opts.verbose==1: print('Executing the command %s' % (cmd)) - # launch command and wait for completion + if opts.verbose==1: print('Executing the command %s' % (cmd)) + # launch command and wait for completion p = Popen(cmd, stdout=PIPE, stderr=PIPE, cwd=head) test.stdout, test.stderr = p.communicate() test.status = p.wait() @@ -459,16 +459,16 @@ def smp_testing(opts,listOfTests,doc): if fd not in subprocs: continue test = subprocs[fd] if fd == test.p.stdout.fileno(): - chunk = os.read(fd,n) + chunk = os.read(fd,n).decode('utf-8', errors='replace') test.stdout += chunk while len(chunk) == n: - chunk = os.read(fd,n) + chunk = os.read(fd,n).decode('utf-8', errors='replace') test.stdout += chunk elif fd == test.p.stderr.fileno(): - chunk = os.read(fd,n) + chunk = os.read(fd,n).decode('utf-8', errors='replace') test.stderr += chunk while len(chunk) == n: - chunk = os.read(fd,n) + chunk = os.read(fd,n).decode('utf-8', errors='replace') test.stderr += chunk done_tests = [] for key in subprocs: @@ -833,7 +833,7 @@ def main(): p.add_option("-c", dest="computer", default=None, \ action="store", \ - help='''Computer to run on (e.g. SandiaCray, + help='''Computer to run on (e.g. SandiaCray, SandiaSkybridge, MacPro, Linux)''') p.add_option("-d", dest="startingDir", default='.', \ @@ -881,9 +881,9 @@ def main(): Composite example showing both comma separated and multiple -k options: - e.g. -k acoustic,short -k Trilinos,elastic + e.g. -k acoustic,short -k Trilinos,elastic selects tests that have both the acoustic and - short keywords, and selects additional tests that + short keywords, and selects additional tests that have both the elastic and Trilinos keywords.''') p.add_option("-K", dest="exclude_keywords", default=[], \ @@ -897,14 +897,14 @@ def main(): either Trilinos tests or short tests. Example with comma separated -K option: - e.g. -K Trilinos,short only excludes a test if - that test has both the Trilinos and short + e.g. -K Trilinos,short only excludes a test if + that test has both the Trilinos and short keywords. Composite example with both -k and -K options: e.g. -k acoustic,short -K Trilinos,Rol will first - select tests that are both acoustic and short, - and will then remove any of those tests that + select tests that are both acoustic and short, + and will then remove any of those tests that have both the Trilinos and Rol keywords.''') p.add_option("--list-keywords", dest="listKeywords", default=False, \ @@ -1006,11 +1006,11 @@ def main(): p.add_option("-v", dest="verbose", default=False, action="store_true", \ help='''Echo the commands to output.''') - + p.add_option("--print-keywords", dest="printKeywords", default=False, \ action="store_true", \ help='''Print the keywords with results of each test.''') - + # Parse the command line options opts, args = p.parse_args() @@ -1027,6 +1027,19 @@ def main(): if opts.startingDir[-1] == "/": opts.startingDir = opts.startingDir[:-1] + def _keyword_requested(keyword_list, keyword): + for item in keyword_list: + for tok in item.split(','): + if tok.strip() == keyword: + return True + return False + + # By default, do not run HPC tests unless explicitly requested. + # (HPC tests are typically environment-dependent and expensive.) + if not _keyword_requested(opts.include_keywords, "hpc"): + if not _keyword_requested(opts.exclude_keywords, "hpc"): + opts.exclude_keywords.append("hpc") + # If test file names are specified, change to full relative path if opts.listOfTestFileNames != []: tmpList = [ os.path.join(opts.startingDir,fname) for diff --git a/sandbox/testCases/hcurl_dirichlet_test/debug.out b/sandbox/testCases/hcurl_dirichlet_test/debug.out deleted file mode 100644 index ad621140..00000000 --- a/sandbox/testCases/hcurl_dirichlet_test/debug.out +++ /dev/null @@ -1,6 +0,0 @@ -[setDirichletData DEBUG] Block=eblock-0_0_0, Var=E, Sideset=bottom, sideEntities.size()=4 -[setDirichletData DEBUG] side_output.size()=4 (elements touching sideset) -[setDirichletData DEBUG] Element 0: local_side_Id=0, getGIDFieldOffsets_closure returned 4 offsets: [0,4,8,9] -[setDirichletData DEBUG] Element 1: local_side_Id=0, getGIDFieldOffsets_closure returned 4 offsets: [0,4,8,9] -[setDirichletData DEBUG] Element 2: local_side_Id=0, getGIDFieldOffsets_closure returned 4 offsets: [0,4,8,9] -[setDirichletData DEBUG] FINAL: Block=eblock-0_0_0, Var=E, unique Dirichlet DOFs marked=12 diff --git a/sandbox/testCases/hcurl_dirichlet_test/input.yaml b/sandbox/testCases/hcurl_dirichlet_test/input.yaml deleted file mode 100644 index 52c4607c..00000000 --- a/sandbox/testCases/hcurl_dirichlet_test/input.yaml +++ /dev/null @@ -1,49 +0,0 @@ -%YAML 1.1 ---- -# Simple test case for HCURL Dirichlet boundary conditions (PEC) -# E field is HCURL, we apply E x n = 0 on bottom face by setting E DOFs to zero -ANONYMOUS: - debug level: 0 - verbosity: 10 # High verbosity to enable debug output in setDirichletData - Mesh input file: input_mesh.yaml - Physics: - eblock-0_0_0: - modules: maxwell - assemble boundary terms: false - Dirichlet conditions: - E: - bottom: '0.0' - Initial conditions: - scalar data: true - eblock-0_0_0: - E: 0.0 - B: 0.0 - Discretization: - eblock-0_0_0: - order: - E: 1 - B: 1 - quadrature: 2 - Functions input file: input_functions.yaml - Solver: - solver: 'transient' - transient BDF order: 1 - transient Butcher tableau: BWE - workset size: 100 - nonlinear TOL: 1.0e-07 - linear TOL: 1.0e-09 - max nonlinear iters: 1 - max linear iters: 200 - final time: 1.0e-2 - number of steps: 1 - Initial type: none - use preconditioner: true - use direct solver: true - allow backtracking: false - Analysis: - analysis type: forward - Postprocess: - write solution: false - compute errors: false -... - diff --git a/sandbox/testCases/hcurl_dirichlet_test/input_functions.yaml b/sandbox/testCases/hcurl_dirichlet_test/input_functions.yaml deleted file mode 100644 index be3e196a..00000000 --- a/sandbox/testCases/hcurl_dirichlet_test/input_functions.yaml +++ /dev/null @@ -1,11 +0,0 @@ -%YAML 1.1 ---- -ANONYMOUS: - Functions: - current x: '0.0' - current y: '0.0' - current z: '0.0' - permittivity: '1.0' - permeability: '1.0' -... - diff --git a/sandbox/testCases/hcurl_dirichlet_test/input_mesh.yaml b/sandbox/testCases/hcurl_dirichlet_test/input_mesh.yaml deleted file mode 100644 index a0fafc37..00000000 --- a/sandbox/testCases/hcurl_dirichlet_test/input_mesh.yaml +++ /dev/null @@ -1,21 +0,0 @@ -%YAML 1.1 ---- -ANONYMOUS: - Mesh: - dimension: 3 - shape: hex - xmin: 0.0 - xmax: 1.0 - ymin: 0.0 - ymax: 1.0 - zmin: 0.0 - zmax: 1.0 - NX: 2 - NY: 2 - NZ: 2 - Xprocs: 1 - Yprocs: 1 - Zprocs: 1 - blocknames: eblock-0_0_0 -... - diff --git a/sandbox/testCases/hcurl_verification/pec_driven/input.yaml b/sandbox/testCases/hcurl_verification/pec_driven/input.yaml deleted file mode 100644 index 4e92f9c4..00000000 --- a/sandbox/testCases/hcurl_verification/pec_driven/input.yaml +++ /dev/null @@ -1,78 +0,0 @@ -%YAML 1.1 ---- -# Test: PEC Cavity with Interior Source -# -# What it tests: -# PEC walls with interior current source Jz = sin(pi*x)*sin(pi*y) driving -# a TM_110-like cavity mode. Longer time evolution than other PEC tests. -# -# Analytical solution: None (visualize in Paraview) -# -# Expected qualitative result: Ez peaks at center, zero at x/y edges, Ex/Ey small -ANONYMOUS: - debug level: 0 - verbosity: 5 - Mesh: - dimension: 3 - shape: hex - xmin: 0.0 - xmax: 1.0 - ymin: 0.0 - ymax: 1.0 - zmin: 0.0 - zmax: 1.0 - NX: 16 - NY: 16 - NZ: 8 - blocknames: eblock-0_0_0 - Physics: - Dirichlet conditions: - scalar data: false - eblock-0_0_0: - modules: maxwell - assemble boundary terms: false - Dirichlet conditions: - Ex: - all boundaries: '0.0' - Ey: - all boundaries: '0.0' - Ez: - all boundaries: '0.0' - Initial conditions: - scalar data: true - E: 0.0 - B: 0.0 - Discretization: - eblock-0_0_0: - order: - E: 1 - B: 1 - quadrature: 3 - Functions: - current x: '0.0' - current y: '0.0' - current z: 'sin(pi*x)*sin(pi*y)' - permittivity: '1.0' - permeability: '1.0' - Solver: - solver: transient - transient BDF order: 1 - transient Butcher tableau: BWE - workset size: 100 - nonlinear TOL: 1.0e-10 - linear TOL: 1.0e-12 - max nonlinear iters: 10 - max linear iters: 500 - final time: 3.0 - number of steps: 100 - use preconditioner: true - use direct solver: true - Analysis: - analysis type: forward - Postprocess: - write solution: true - write frequency: 1 - exodus write frequency: 1 - compute errors: false -... - diff --git a/sandbox/testCases/hcurl_verification/pec_driven/mrhyde.gold b/sandbox/testCases/hcurl_verification/pec_driven/mrhyde.gold deleted file mode 100644 index cbcc15e4..00000000 --- a/sandbox/testCases/hcurl_verification/pec_driven/mrhyde.gold +++ /dev/null @@ -1,2314 +0,0 @@ -p=0 | CubeHexMesh: Building sub cells -STK Meta data: - Element blocks = "eblock-0_0_0" - Sidesets = "back" "bottom" "front" "left" "right" "top" - Nodesets = "origin" - Fields = "coordinates" "edges" "faces" "PROC_ID" "LOAD_BAL" "Ex" "Ey" "Ez" "Bx" "By" "Bz" -DOFManager Field Information: - Element Block = eblock-0_0_0 - "E" is field ID 0 - "B" is field ID 1 - -********************************************************* -***** Writing the solution to output.exo -********************************************************* - - -******************************************************* - -**** Beginning Time Step 1 -**** Current time is 0 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00776246 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.54092e-16 -***** Scaled Norm of nonlinear residual: 3.27335e-14 -********************************************************* -Norm of solution: 0.0469044 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 2 -**** Current time is 0.03 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0230009 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.48936e-16 -***** Scaled Norm of nonlinear residual: 1.51705e-14 -********************************************************* -Norm of solution: 0.092167 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 3 -**** Current time is 0.06 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0451963 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.26418e-16 -***** Scaled Norm of nonlinear residual: 7.22222e-15 -********************************************************* -Norm of solution: 0.135025 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 4 -**** Current time is 0.09 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0662116 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.3494e-16 -***** Scaled Norm of nonlinear residual: 3.54832e-15 -********************************************************* -Norm of solution: 0.17477 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 5 -**** Current time is 0.12 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0856996 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.53203e-16 -***** Scaled Norm of nonlinear residual: 2.95454e-15 -********************************************************* -Norm of solution: 0.210761 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 6 -**** Current time is 0.15 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.103346 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.80762e-16 -***** Scaled Norm of nonlinear residual: 2.71673e-15 -********************************************************* -Norm of solution: 0.242435 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 7 -**** Current time is 0.18 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.118873 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.92726e-16 -***** Scaled Norm of nonlinear residual: 3.30374e-15 -********************************************************* -Norm of solution: 0.269314 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 8 -**** Current time is 0.21 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.132047 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.19266e-16 -***** Scaled Norm of nonlinear residual: 2.41782e-15 -********************************************************* -Norm of solution: 0.291011 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 9 -**** Current time is 0.24 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.142678 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 4.20338e-16 -***** Scaled Norm of nonlinear residual: 2.94605e-15 -********************************************************* -Norm of solution: 0.307238 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 10 -**** Current time is 0.27 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.150625 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 4.17182e-16 -***** Scaled Norm of nonlinear residual: 2.76967e-15 -********************************************************* -Norm of solution: 0.317808 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 11 -**** Current time is 0.3 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.155795 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 4.85249e-16 -***** Scaled Norm of nonlinear residual: 3.11466e-15 -********************************************************* -Norm of solution: 0.322637 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 12 -**** Current time is 0.33 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.158147 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 5.71396e-16 -***** Scaled Norm of nonlinear residual: 3.61307e-15 -********************************************************* -Norm of solution: 0.321741 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 13 -**** Current time is 0.36 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.157688 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 4.83246e-16 -***** Scaled Norm of nonlinear residual: 3.06457e-15 -********************************************************* -Norm of solution: 0.31524 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 14 -**** Current time is 0.39 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.154476 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.37694e-16 -***** Scaled Norm of nonlinear residual: 2.18606e-15 -********************************************************* -Norm of solution: 0.303347 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 15 -**** Current time is 0.42 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.148615 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.32502e-16 -***** Scaled Norm of nonlinear residual: 2.23734e-15 -********************************************************* -Norm of solution: 0.28637 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 16 -**** Current time is 0.45 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.140255 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.36561e-16 -***** Scaled Norm of nonlinear residual: 2.39964e-15 -********************************************************* -Norm of solution: 0.2647 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 17 -**** Current time is 0.48 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.129585 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.03948e-16 -***** Scaled Norm of nonlinear residual: 2.34556e-15 -********************************************************* -Norm of solution: 0.238806 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 18 -**** Current time is 0.51 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.116833 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.18507e-16 -***** Scaled Norm of nonlinear residual: 2.72618e-15 -********************************************************* -Norm of solution: 0.209231 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 19 -**** Current time is 0.54 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.102258 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.44293e-16 -***** Scaled Norm of nonlinear residual: 3.3669e-15 -********************************************************* -Norm of solution: 0.176578 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 20 -**** Current time is 0.57 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.086148 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.80775e-16 -***** Scaled Norm of nonlinear residual: 3.25922e-15 -********************************************************* -Norm of solution: 0.141523 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 21 -**** Current time is 0.6 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0688116 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.65599e-16 -***** Scaled Norm of nonlinear residual: 3.85979e-15 -********************************************************* -Norm of solution: 0.104831 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 22 -**** Current time is 0.63 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0505738 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.42943e-16 -***** Scaled Norm of nonlinear residual: 4.80373e-15 -********************************************************* -Norm of solution: 0.0675302 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 23 -**** Current time is 0.66 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0317699 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.02178e-16 -***** Scaled Norm of nonlinear residual: 6.36384e-15 -********************************************************* -Norm of solution: 0.0322778 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 24 -**** Current time is 0.69 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0127389 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.1863e-16 -***** Scaled Norm of nonlinear residual: 1.71624e-14 -********************************************************* -Norm of solution: 0.0228921 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 25 -**** Current time is 0.72 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00627368 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.22396e-16 -***** Scaled Norm of nonlinear residual: 3.5449e-14 -********************************************************* -Norm of solution: 0.0537313 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 26 -**** Current time is 0.75 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0246632 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.29991e-16 -***** Scaled Norm of nonlinear residual: 9.32527e-15 -********************************************************* -Norm of solution: 0.0884126 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 27 -**** Current time is 0.78 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.042389 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.53072e-16 -***** Scaled Norm of nonlinear residual: 5.97023e-15 -********************************************************* -Norm of solution: 0.121792 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 28 -**** Current time is 0.81 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0590621 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.86945e-16 -***** Scaled Norm of nonlinear residual: 4.85836e-15 -********************************************************* -Norm of solution: 0.152742 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 29 -**** Current time is 0.84 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0744092 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.29604e-16 -***** Scaled Norm of nonlinear residual: 3.0857e-15 -********************************************************* -Norm of solution: 0.18061 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 30 -**** Current time is 0.87 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0881845 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.08126e-16 -***** Scaled Norm of nonlinear residual: 2.36012e-15 -********************************************************* -Norm of solution: 0.204902 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 31 -**** Current time is 0.9 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.100174 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.78399e-16 -***** Scaled Norm of nonlinear residual: 2.77914e-15 -********************************************************* -Norm of solution: 0.22523 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 32 -**** Current time is 0.93 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.1102 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.88873e-16 -***** Scaled Norm of nonlinear residual: 2.62134e-15 -********************************************************* -Norm of solution: 0.24129 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 33 -**** Current time is 0.96 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.118121 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.61089e-16 -***** Scaled Norm of nonlinear residual: 3.05694e-15 -********************************************************* -Norm of solution: 0.252869 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 34 -**** Current time is 0.99 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.123835 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.16255e-16 -***** Scaled Norm of nonlinear residual: 2.55384e-15 -********************************************************* -Norm of solution: 0.259836 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 35 -**** Current time is 1.02 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.12728 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.92118e-16 -***** Scaled Norm of nonlinear residual: 3.08074e-15 -********************************************************* -Norm of solution: 0.262145 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 36 -**** Current time is 1.05 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.128436 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.24248e-16 -***** Scaled Norm of nonlinear residual: 2.52458e-15 -********************************************************* -Norm of solution: 0.259836 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 37 -**** Current time is 1.08 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.127323 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 4.23755e-16 -***** Scaled Norm of nonlinear residual: 3.32818e-15 -********************************************************* -Norm of solution: 0.253027 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 38 -**** Current time is 1.11 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.124 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.63781e-16 -***** Scaled Norm of nonlinear residual: 2.12726e-15 -********************************************************* -Norm of solution: 0.241915 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 39 -**** Current time is 1.14 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.118564 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.58184e-16 -***** Scaled Norm of nonlinear residual: 2.17759e-15 -********************************************************* -Norm of solution: 0.226769 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 40 -**** Current time is 1.17 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.111147 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.75388e-16 -***** Scaled Norm of nonlinear residual: 2.47768e-15 -********************************************************* -Norm of solution: 0.207924 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 41 -**** Current time is 1.2 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.101914 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.86235e-16 -***** Scaled Norm of nonlinear residual: 2.80859e-15 -********************************************************* -Norm of solution: 0.185773 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 42 -**** Current time is 1.23 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0910577 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.0766e-16 -***** Scaled Norm of nonlinear residual: 3.37874e-15 -********************************************************* -Norm of solution: 0.160764 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 43 -**** Current time is 1.26 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0787971 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.74432e-16 -***** Scaled Norm of nonlinear residual: 2.21369e-15 -********************************************************* -Norm of solution: 0.133383 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 44 -**** Current time is 1.29 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0653714 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.63397e-16 -***** Scaled Norm of nonlinear residual: 2.49952e-15 -********************************************************* -Norm of solution: 0.104155 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 45 -**** Current time is 1.32 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0510361 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.8685e-16 -***** Scaled Norm of nonlinear residual: 3.66113e-15 -********************************************************* -Norm of solution: 0.0736293 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 46 -**** Current time is 1.35 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0360582 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.39567e-16 -***** Scaled Norm of nonlinear residual: 3.87061e-15 -********************************************************* -Norm of solution: 0.0423914 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 47 -**** Current time is 1.38 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0207112 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.04413e-16 -***** Scaled Norm of nonlinear residual: 9.86967e-15 -********************************************************* -Norm of solution: 0.0113334 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 48 -**** Current time is 1.41 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00527027 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.51922e-16 -***** Scaled Norm of nonlinear residual: 2.88262e-14 -********************************************************* -Norm of solution: 0.0207064 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 49 -**** Current time is 1.44 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00999259 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.48177e-16 -***** Scaled Norm of nonlinear residual: 1.48286e-14 -********************************************************* -Norm of solution: 0.0507481 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 50 -**** Current time is 1.47 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0248132 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.81657e-16 -***** Scaled Norm of nonlinear residual: 7.32099e-15 -********************************************************* -Norm of solution: 0.0795176 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 51 -**** Current time is 1.5 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0389398 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.12941e-16 -***** Scaled Norm of nonlinear residual: 5.46847e-15 -********************************************************* -Norm of solution: 0.106419 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 52 -**** Current time is 1.53 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0521373 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.26252e-16 -***** Scaled Norm of nonlinear residual: 4.33955e-15 -********************************************************* -Norm of solution: 0.130999 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 53 -**** Current time is 1.56 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0641906 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.08716e-16 -***** Scaled Norm of nonlinear residual: 3.2515e-15 -********************************************************* -Norm of solution: 0.152864 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 54 -**** Current time is 1.59 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.074909 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.85065e-16 -***** Scaled Norm of nonlinear residual: 2.47053e-15 -********************************************************* -Norm of solution: 0.171675 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 55 -**** Current time is 1.62 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0841279 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.96501e-16 -***** Scaled Norm of nonlinear residual: 2.33574e-15 -********************************************************* -Norm of solution: 0.187158 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 56 -**** Current time is 1.65 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0917123 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.37813e-16 -***** Scaled Norm of nonlinear residual: 2.59303e-15 -********************************************************* -Norm of solution: 0.199101 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 57 -**** Current time is 1.68 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.097558 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.06686e-16 -***** Scaled Norm of nonlinear residual: 2.1186e-15 -********************************************************* -Norm of solution: 0.207356 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 58 -**** Current time is 1.71 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.101593 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.90902e-16 -***** Scaled Norm of nonlinear residual: 2.8634e-15 -********************************************************* -Norm of solution: 0.211844 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 59 -**** Current time is 1.74 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.103779 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.89927e-16 -***** Scaled Norm of nonlinear residual: 2.79371e-15 -********************************************************* -Norm of solution: 0.212556 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 60 -**** Current time is 1.77 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.104109 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.4595e-16 -***** Scaled Norm of nonlinear residual: 2.36244e-15 -********************************************************* -Norm of solution: 0.209546 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 61 -**** Current time is 1.8 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.10261 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.42167e-16 -***** Scaled Norm of nonlinear residual: 2.36007e-15 -********************************************************* -Norm of solution: 0.202935 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 62 -**** Current time is 1.83 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0993406 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.40083e-16 -***** Scaled Norm of nonlinear residual: 3.4234e-15 -********************************************************* -Norm of solution: 0.192906 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 63 -**** Current time is 1.86 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0943891 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.14833e-16 -***** Scaled Norm of nonlinear residual: 2.27604e-15 -********************************************************* -Norm of solution: 0.179698 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 64 -**** Current time is 1.89 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0878715 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.01142e-16 -***** Scaled Norm of nonlinear residual: 2.28905e-15 -********************************************************* -Norm of solution: 0.163607 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 65 -**** Current time is 1.92 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0799294 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.61798e-16 -***** Scaled Norm of nonlinear residual: 3.27537e-15 -********************************************************* -Norm of solution: 0.144977 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 66 -**** Current time is 1.95 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0707268 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.97139e-16 -***** Scaled Norm of nonlinear residual: 2.78734e-15 -********************************************************* -Norm of solution: 0.124199 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 67 -**** Current time is 1.98 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0604468 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.08827e-16 -***** Scaled Norm of nonlinear residual: 3.45472e-15 -********************************************************* -Norm of solution: 0.101714 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 68 -**** Current time is 2.01 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0492885 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.89647e-16 -***** Scaled Norm of nonlinear residual: 3.8477e-15 -********************************************************* -Norm of solution: 0.0780325 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 69 -**** Current time is 2.04 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0374624 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.65286e-16 -***** Scaled Norm of nonlinear residual: 4.41204e-15 -********************************************************* -Norm of solution: 0.0538296 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 70 -**** Current time is 2.07 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0251874 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.84639e-16 -***** Scaled Norm of nonlinear residual: 7.33061e-15 -********************************************************* -Norm of solution: 0.0305296 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 71 -**** Current time is 2.1 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0126864 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.75992e-16 -***** Scaled Norm of nonlinear residual: 1.38725e-14 -********************************************************* -Norm of solution: 0.0162162 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 72 -**** Current time is 2.13 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00414612 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.48172e-16 -***** Scaled Norm of nonlinear residual: 3.57374e-14 -********************************************************* -Norm of solution: 0.0294816 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 73 -**** Current time is 2.16 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0121064 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.59872e-16 -***** Scaled Norm of nonlinear residual: 1.32055e-14 -********************************************************* -Norm of solution: 0.0514037 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 74 -**** Current time is 2.19 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0239678 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.41664e-16 -***** Scaled Norm of nonlinear residual: 5.9106e-15 -********************************************************* -Norm of solution: 0.0734681 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 75 -**** Current time is 2.22 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0352018 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.61394e-16 -***** Scaled Norm of nonlinear residual: 4.58483e-15 -********************************************************* -Norm of solution: 0.0942799 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 76 -**** Current time is 2.25 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0456227 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.81379e-16 -***** Scaled Norm of nonlinear residual: 3.97564e-15 -********************************************************* -Norm of solution: 0.113259 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 77 -**** Current time is 2.28 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0550622 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.74278e-16 -***** Scaled Norm of nonlinear residual: 4.98124e-15 -********************************************************* -Norm of solution: 0.130021 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 78 -**** Current time is 2.31 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0633723 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.10968e-16 -***** Scaled Norm of nonlinear residual: 3.32902e-15 -********************************************************* -Norm of solution: 0.144274 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 79 -**** Current time is 2.34 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0704272 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.40504e-16 -***** Scaled Norm of nonlinear residual: 3.41493e-15 -********************************************************* -Norm of solution: 0.155792 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 80 -**** Current time is 2.37 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0761254 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.7453e-16 -***** Scaled Norm of nonlinear residual: 3.60629e-15 -********************************************************* -Norm of solution: 0.16441 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 81 -**** Current time is 2.4 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0803908 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.9963e-16 -***** Scaled Norm of nonlinear residual: 2.48324e-15 -********************************************************* -Norm of solution: 0.170022 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 82 -**** Current time is 2.43 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0831738 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.63887e-16 -***** Scaled Norm of nonlinear residual: 3.17271e-15 -********************************************************* -Norm of solution: 0.172576 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 83 -**** Current time is 2.46 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0844517 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.46532e-16 -***** Scaled Norm of nonlinear residual: 2.91921e-15 -********************************************************* -Norm of solution: 0.172079 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 84 -**** Current time is 2.49 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0842285 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.84938e-16 -***** Scaled Norm of nonlinear residual: 3.38292e-15 -********************************************************* -Norm of solution: 0.168591 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 85 -**** Current time is 2.52 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0825343 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 3.09879e-16 -***** Scaled Norm of nonlinear residual: 3.75454e-15 -********************************************************* -Norm of solution: 0.162223 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 86 -**** Current time is 2.55 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0794245 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 2.34926e-16 -***** Scaled Norm of nonlinear residual: 2.95785e-15 -********************************************************* -Norm of solution: 0.153137 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 87 -**** Current time is 2.58 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0749785 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.96854e-16 -***** Scaled Norm of nonlinear residual: 2.62548e-15 -********************************************************* -Norm of solution: 0.141539 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 88 -**** Current time is 2.61 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0692974 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.48187e-16 -***** Scaled Norm of nonlinear residual: 2.13842e-15 -********************************************************* -Norm of solution: 0.127677 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 89 -**** Current time is 2.64 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0625023 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.35141e-16 -***** Scaled Norm of nonlinear residual: 2.16217e-15 -********************************************************* -Norm of solution: 0.111834 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 90 -**** Current time is 2.67 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0547318 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.25855e-16 -***** Scaled Norm of nonlinear residual: 2.29949e-15 -********************************************************* -Norm of solution: 0.0943253 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 91 -**** Current time is 2.7 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.046139 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.20247e-16 -***** Scaled Norm of nonlinear residual: 2.60618e-15 -********************************************************* -Norm of solution: 0.0754956 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 92 -**** Current time is 2.73 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0368888 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.36133e-16 -***** Scaled Norm of nonlinear residual: 3.69037e-15 -********************************************************* -Norm of solution: 0.0557173 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 93 -**** Current time is 2.76 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0271546 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.32667e-16 -***** Scaled Norm of nonlinear residual: 4.88561e-15 -********************************************************* -Norm of solution: 0.03542 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 94 -**** Current time is 2.79 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0171153 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.17368e-16 -***** Scaled Norm of nonlinear residual: 6.85748e-15 -********************************************************* -Norm of solution: 0.0153866 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 95 -**** Current time is 2.82 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00695217 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.34846e-16 -***** Scaled Norm of nonlinear residual: 1.93963e-14 -********************************************************* -Norm of solution: 0.00880062 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 96 -**** Current time is 2.85 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.00335121 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.56556e-16 -***** Scaled Norm of nonlinear residual: 4.67163e-14 -********************************************************* -Norm of solution: 0.0272635 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 97 -**** Current time is 2.88 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0130295 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.14362e-16 -***** Scaled Norm of nonlinear residual: 8.77717e-15 -********************************************************* -Norm of solution: 0.0463182 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 98 -**** Current time is 2.91 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0225031 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.01975e-16 -***** Scaled Norm of nonlinear residual: 4.53159e-15 -********************************************************* -Norm of solution: 0.0643993 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 99 -**** Current time is 2.94 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0314169 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.2317e-16 -***** Scaled Norm of nonlinear residual: 3.9205e-15 -********************************************************* -Norm of solution: 0.0810939 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 100 -**** Current time is 2.97 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 0.0396244 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.03651e-16 -***** Scaled Norm of nonlinear residual: 2.61584e-15 -********************************************************* -Norm of solution: 0.0961051 (overlapped vector so results may differ on multiple procs) diff --git a/sandbox/testCases/hcurl_verification/pec_free/input.yaml b/sandbox/testCases/hcurl_verification/pec_free/input.yaml deleted file mode 100644 index 7df0ffa1..00000000 --- a/sandbox/testCases/hcurl_verification/pec_free/input.yaml +++ /dev/null @@ -1,84 +0,0 @@ -%YAML 1.1 ---- -# Test: PEC Cavity -# -# What it tests: -# Homogeneous PEC boundary conditions (n x E = 0) with non-zero interior IC. -# Verifies that boundary DOFs are clamped to zero while interior evolves. -# -# Analytical solution: None (complex cavity resonance modes) -# -# Expected qualitative result: Tangential E = 0 at walls, non-zero interior field - -ANONYMOUS: - debug level: 0 - verbosity: 5 - - Mesh: - dimension: 3 - shape: hex - xmin: 0.0 - xmax: 1.0 - ymin: 0.0 - ymax: 1.0 - zmin: 0.0 - zmax: 1.0 - NX: 8 - NY: 8 - NZ: 8 - blocknames: eblock-0_0_0 - - Physics: - Dirichlet conditions: - scalar data: false - eblock-0_0_0: - modules: maxwell - assemble boundary terms: false - Dirichlet conditions: - Ex: - all boundaries: '0.0' - Ey: - all boundaries: '0.0' - Ez: - all boundaries: '0.0' - Initial conditions: - scalar data: true - E: 1.0 - B: 0.0 - - Discretization: - eblock-0_0_0: - order: - E: 1 - B: 1 - quadrature: 2 - - Functions: - current x: '0.0' - current y: '0.0' - current z: '0.0' - permittivity: '1.0' - permeability: '1.0' - - Solver: - solver: transient - transient BDF order: 1 - transient Butcher tableau: BWE - workset size: 100 - nonlinear TOL: 1.0e-06 - linear TOL: 1.0e-10 - max nonlinear iters: 3 - max linear iters: 500 - final time: 1.0 - number of steps: 100 - use preconditioner: true - use direct solver: true - - Analysis: - analysis type: forward - - Postprocess: - write solution: true - write frequency: 1 - compute errors: false -... diff --git a/sandbox/testCases/hcurl_verification/pec_free/mrhyde.gold b/sandbox/testCases/hcurl_verification/pec_free/mrhyde.gold deleted file mode 100644 index 6fd2f310..00000000 --- a/sandbox/testCases/hcurl_verification/pec_free/mrhyde.gold +++ /dev/null @@ -1,2314 +0,0 @@ -p=0 | CubeHexMesh: Building sub cells -STK Meta data: - Element blocks = "eblock-0_0_0" - Sidesets = "back" "bottom" "front" "left" "right" "top" - Nodesets = "origin" - Fields = "coordinates" "edges" "faces" "PROC_ID" "LOAD_BAL" "Ex" "Ey" "Ez" "Bx" "By" "Bz" -DOFManager Field Information: - Element Block = eblock-0_0_0 - "E" is field ID 0 - "B" is field ID 1 - -********************************************************* -***** Writing the solution to output.exo -********************************************************* - - -******************************************************* - -**** Beginning Time Step 1 -**** Current time is 0 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 128 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.25913e-13 -***** Scaled Norm of nonlinear residual: 9.83695e-16 -********************************************************* -Norm of solution: 38.7432 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 2 -**** Current time is 0.01 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 126.722 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.06775e-13 -***** Scaled Norm of nonlinear residual: 8.42598e-16 -********************************************************* -Norm of solution: 38.4677 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 3 -**** Current time is 0.02 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 124.239 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.58062e-13 -***** Scaled Norm of nonlinear residual: 1.27224e-15 -********************************************************* -Norm of solution: 38.0695 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 4 -**** Current time is 0.03 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 120.657 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.47156e-13 -***** Scaled Norm of nonlinear residual: 1.21962e-15 -********************************************************* -Norm of solution: 37.5632 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 5 -**** Current time is 0.04 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 116.103 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.01473e-13 -***** Scaled Norm of nonlinear residual: 8.7399e-16 -********************************************************* -Norm of solution: 36.9646 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 6 -**** Current time is 0.05 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 110.724 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.07197e-13 -***** Scaled Norm of nonlinear residual: 9.68154e-16 -********************************************************* -Norm of solution: 36.2891 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 7 -**** Current time is 0.06 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 104.672 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.09939e-13 -***** Scaled Norm of nonlinear residual: 1.05032e-15 -********************************************************* -Norm of solution: 35.5503 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 8 -**** Current time is 0.07 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 98.1076 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.12174e-14 -***** Scaled Norm of nonlinear residual: 9.29769e-16 -********************************************************* -Norm of solution: 34.7596 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 9 -**** Current time is 0.08 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 91.1848 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.26851e-13 -***** Scaled Norm of nonlinear residual: 1.39115e-15 -********************************************************* -Norm of solution: 33.9255 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 10 -**** Current time is 0.09 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 84.0501 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.78936e-14 -***** Scaled Norm of nonlinear residual: 1.16471e-15 -********************************************************* -Norm of solution: 33.0541 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 11 -**** Current time is 0.1 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 76.8372 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.23239e-13 -***** Scaled Norm of nonlinear residual: 1.60389e-15 -********************************************************* -Norm of solution: 32.1499 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 12 -**** Current time is 0.11 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 69.6635 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.16126e-14 -***** Scaled Norm of nonlinear residual: 1.31507e-15 -********************************************************* -Norm of solution: 31.2163 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 13 -**** Current time is 0.12 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 65.3206 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.44826e-13 -***** Scaled Norm of nonlinear residual: 2.21716e-15 -********************************************************* -Norm of solution: 30.2566 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 14 -**** Current time is 0.13 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 69.1513 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.03528e-13 -***** Scaled Norm of nonlinear residual: 1.49712e-15 -********************************************************* -Norm of solution: 29.275 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 15 -**** Current time is 0.14 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 71.9251 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.03672e-13 -***** Scaled Norm of nonlinear residual: 1.44138e-15 -********************************************************* -Norm of solution: 28.2764 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 16 -**** Current time is 0.15 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 73.5963 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.16183e-13 -***** Scaled Norm of nonlinear residual: 1.57866e-15 -********************************************************* -Norm of solution: 27.2675 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 17 -**** Current time is 0.16 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 74.1558 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.19471e-13 -***** Scaled Norm of nonlinear residual: 1.61108e-15 -********************************************************* -Norm of solution: 26.2561 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 18 -**** Current time is 0.17 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 73.6527 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.32091e-14 -***** Scaled Norm of nonlinear residual: 9.93978e-16 -********************************************************* -Norm of solution: 25.2515 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 19 -**** Current time is 0.18 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 72.4737 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.23848e-13 -***** Scaled Norm of nonlinear residual: 1.70887e-15 -********************************************************* -Norm of solution: 24.2633 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 20 -**** Current time is 0.19 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 70.4829 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.686e-14 -***** Scaled Norm of nonlinear residual: 1.23236e-15 -********************************************************* -Norm of solution: 23.3013 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 21 -**** Current time is 0.2 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 67.7889 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.70544e-14 -***** Scaled Norm of nonlinear residual: 1.43172e-15 -********************************************************* -Norm of solution: 22.3744 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 22 -**** Current time is 0.21 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 64.5078 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.45377e-14 -***** Scaled Norm of nonlinear residual: 1.3105e-15 -********************************************************* -Norm of solution: 21.4901 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 23 -**** Current time is 0.22 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 60.7575 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.36907e-14 -***** Scaled Norm of nonlinear residual: 1.54204e-15 -********************************************************* -Norm of solution: 20.653 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 24 -**** Current time is 0.23 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 58.3734 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.49542e-14 -***** Scaled Norm of nonlinear residual: 1.62667e-15 -********************************************************* -Norm of solution: 19.8645 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 25 -**** Current time is 0.24 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 61.8575 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.04888e-14 -***** Scaled Norm of nonlinear residual: 1.3012e-15 -********************************************************* -Norm of solution: 19.1222 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 26 -**** Current time is 0.25 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 64.6989 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.19653e-13 -***** Scaled Norm of nonlinear residual: 1.84938e-15 -********************************************************* -Norm of solution: 18.4195 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 27 -**** Current time is 0.26 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 66.8693 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.10378e-13 -***** Scaled Norm of nonlinear residual: 1.65066e-15 -********************************************************* -Norm of solution: 17.7463 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 28 -**** Current time is 0.27 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 68.3571 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.09123e-13 -***** Scaled Norm of nonlinear residual: 1.59637e-15 -********************************************************* -Norm of solution: 17.0898 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 29 -**** Current time is 0.28 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 69.1646 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.96468e-14 -***** Scaled Norm of nonlinear residual: 1.44072e-15 -********************************************************* -Norm of solution: 16.4357 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 30 -**** Current time is 0.29 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 69.3046 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.14001e-13 -***** Scaled Norm of nonlinear residual: 1.64492e-15 -********************************************************* -Norm of solution: 15.7706 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 31 -**** Current time is 0.3 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 68.7965 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.78579e-14 -***** Scaled Norm of nonlinear residual: 1.13171e-15 -********************************************************* -Norm of solution: 15.0835 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 32 -**** Current time is 0.31 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 67.6633 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.06014e-13 -***** Scaled Norm of nonlinear residual: 1.56679e-15 -********************************************************* -Norm of solution: 14.3685 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 33 -**** Current time is 0.32 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 65.9284 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.46916e-14 -***** Scaled Norm of nonlinear residual: 1.43628e-15 -********************************************************* -Norm of solution: 13.6273 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 34 -**** Current time is 0.33 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 63.6135 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.14905e-13 -***** Scaled Norm of nonlinear residual: 1.80629e-15 -********************************************************* -Norm of solution: 12.8711 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 35 -**** Current time is 0.34 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 60.7371 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.14845e-13 -***** Scaled Norm of nonlinear residual: 1.89086e-15 -********************************************************* -Norm of solution: 12.1238 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 36 -**** Current time is 0.35 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 57.3142 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.22977e-13 -***** Scaled Norm of nonlinear residual: 2.14567e-15 -********************************************************* -Norm of solution: 11.424 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 37 -**** Current time is 0.36 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 53.3558 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.86966e-14 -***** Scaled Norm of nonlinear residual: 1.84978e-15 -********************************************************* -Norm of solution: 10.8258 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 38 -**** Current time is 0.37 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 48.8708 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.98453e-14 -***** Scaled Norm of nonlinear residual: 1.63381e-15 -********************************************************* -Norm of solution: 10.3962 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 39 -**** Current time is 0.38 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 43.8667 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.00109e-14 -***** Scaled Norm of nonlinear residual: 2.05192e-15 -********************************************************* -Norm of solution: 10.2054 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 40 -**** Current time is 0.39 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 38.3525 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.56888e-14 -***** Scaled Norm of nonlinear residual: 2.23424e-15 -********************************************************* -Norm of solution: 10.3094 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 41 -**** Current time is 0.4 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 32.3403 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.29435e-14 -***** Scaled Norm of nonlinear residual: 2.56471e-15 -********************************************************* -Norm of solution: 10.7319 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 42 -**** Current time is 0.41 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 32.8161 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.72341e-14 -***** Scaled Norm of nonlinear residual: 2.65827e-15 -********************************************************* -Norm of solution: 11.4573 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 43 -**** Current time is 0.42 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 35.9322 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.22427e-14 -***** Scaled Norm of nonlinear residual: 2.28883e-15 -********************************************************* -Norm of solution: 12.44 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 44 -**** Current time is 0.43 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 38.7843 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.73821e-14 -***** Scaled Norm of nonlinear residual: 1.99519e-15 -********************************************************* -Norm of solution: 13.6205 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 45 -**** Current time is 0.44 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 41.3224 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.08458e-14 -***** Scaled Norm of nonlinear residual: 1.95646e-15 -********************************************************* -Norm of solution: 14.9395 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 46 -**** Current time is 0.45 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 43.5042 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.55952e-14 -***** Scaled Norm of nonlinear residual: 1.96751e-15 -********************************************************* -Norm of solution: 16.345 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 47 -**** Current time is 0.46 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 45.2962 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.36697e-14 -***** Scaled Norm of nonlinear residual: 1.84717e-15 -********************************************************* -Norm of solution: 17.7939 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 48 -**** Current time is 0.47 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 46.6745 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.6342e-14 -***** Scaled Norm of nonlinear residual: 2.06412e-15 -********************************************************* -Norm of solution: 19.2518 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 49 -**** Current time is 0.48 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 47.6262 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.62816e-14 -***** Scaled Norm of nonlinear residual: 2.02161e-15 -********************************************************* -Norm of solution: 20.6922 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 50 -**** Current time is 0.49 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 48.1493 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.63049e-14 -***** Scaled Norm of nonlinear residual: 2.00013e-15 -********************************************************* -Norm of solution: 22.0944 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 51 -**** Current time is 0.5 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 48.2537 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.87854e-14 -***** Scaled Norm of nonlinear residual: 2.04721e-15 -********************************************************* -Norm of solution: 23.4428 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 52 -**** Current time is 0.51 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 54.1449 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.58056e-14 -***** Scaled Norm of nonlinear residual: 1.76943e-15 -********************************************************* -Norm of solution: 24.726 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 53 -**** Current time is 0.52 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 61.7658 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.97238e-14 -***** Scaled Norm of nonlinear residual: 1.61455e-15 -********************************************************* -Norm of solution: 25.9355 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 54 -**** Current time is 0.53 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 68.8905 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.23248e-13 -***** Scaled Norm of nonlinear residual: 1.78905e-15 -********************************************************* -Norm of solution: 27.0657 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 55 -**** Current time is 0.54 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 75.4156 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.70612e-14 -***** Scaled Norm of nonlinear residual: 1.15442e-15 -********************************************************* -Norm of solution: 28.1127 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 56 -**** Current time is 0.55 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 81.2471 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.16756e-13 -***** Scaled Norm of nonlinear residual: 1.43705e-15 -********************************************************* -Norm of solution: 29.0743 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 57 -**** Current time is 0.56 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 86.3041 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.23788e-13 -***** Scaled Norm of nonlinear residual: 1.43433e-15 -********************************************************* -Norm of solution: 29.9491 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 58 -**** Current time is 0.57 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 90.5204 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.20885e-14 -***** Scaled Norm of nonlinear residual: 1.01732e-15 -********************************************************* -Norm of solution: 30.7368 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 59 -**** Current time is 0.58 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 93.8469 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.70692e-14 -***** Scaled Norm of nonlinear residual: 1.03434e-15 -********************************************************* -Norm of solution: 31.4376 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 60 -**** Current time is 0.59 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 96.2523 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.94957e-14 -***** Scaled Norm of nonlinear residual: 1.0337e-15 -********************************************************* -Norm of solution: 32.0521 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 61 -**** Current time is 0.6 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 97.7238 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.00575e-13 -***** Scaled Norm of nonlinear residual: 1.02918e-15 -********************************************************* -Norm of solution: 32.5811 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 62 -**** Current time is 0.61 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 98.2669 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.47353e-14 -***** Scaled Norm of nonlinear residual: 8.62297e-16 -********************************************************* -Norm of solution: 33.0261 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 63 -**** Current time is 0.62 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 97.9045 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.64847e-14 -***** Scaled Norm of nonlinear residual: 8.83358e-16 -********************************************************* -Norm of solution: 33.3885 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 64 -**** Current time is 0.63 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 96.6752 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.11089e-14 -***** Scaled Norm of nonlinear residual: 8.38984e-16 -********************************************************* -Norm of solution: 33.6705 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 65 -**** Current time is 0.64 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 94.6319 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.53039e-14 -***** Scaled Norm of nonlinear residual: 7.95756e-16 -********************************************************* -Norm of solution: 33.8744 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 66 -**** Current time is 0.65 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 91.839 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.37041e-14 -***** Scaled Norm of nonlinear residual: 1.02031e-15 -********************************************************* -Norm of solution: 34.0032 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 67 -**** Current time is 0.66 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 88.3702 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 6.75482e-14 -***** Scaled Norm of nonlinear residual: 7.64378e-16 -********************************************************* -Norm of solution: 34.0604 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 68 -**** Current time is 0.67 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 85.9181 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.83291e-14 -***** Scaled Norm of nonlinear residual: 1.14445e-15 -********************************************************* -Norm of solution: 34.05 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 69 -**** Current time is 0.68 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 85.9137 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 6.8748e-14 -***** Scaled Norm of nonlinear residual: 8.00199e-16 -********************************************************* -Norm of solution: 33.9765 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 70 -**** Current time is 0.69 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 85.2081 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.4375e-13 -***** Scaled Norm of nonlinear residual: 1.68705e-15 -********************************************************* -Norm of solution: 33.8447 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 71 -**** Current time is 0.7 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 83.7966 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.91437e-14 -***** Scaled Norm of nonlinear residual: 1.18315e-15 -********************************************************* -Norm of solution: 33.66 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 72 -**** Current time is 0.71 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 81.6939 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.07687e-13 -***** Scaled Norm of nonlinear residual: 1.31818e-15 -********************************************************* -Norm of solution: 33.4274 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 73 -**** Current time is 0.72 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 78.9333 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.45341e-14 -***** Scaled Norm of nonlinear residual: 1.07096e-15 -********************************************************* -Norm of solution: 33.1521 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 74 -**** Current time is 0.73 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 75.5667 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.80286e-14 -***** Scaled Norm of nonlinear residual: 1.03258e-15 -********************************************************* -Norm of solution: 32.8385 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 75 -**** Current time is 0.74 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 71.6625 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.29822e-14 -***** Scaled Norm of nonlinear residual: 1.15796e-15 -********************************************************* -Norm of solution: 32.4906 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 76 -**** Current time is 0.75 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 67.3032 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 6.19578e-14 -***** Scaled Norm of nonlinear residual: 9.20578e-16 -********************************************************* -Norm of solution: 32.1114 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 77 -**** Current time is 0.76 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 63.9366 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.31331e-14 -***** Scaled Norm of nonlinear residual: 1.30024e-15 -********************************************************* -Norm of solution: 31.7028 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 78 -**** Current time is 0.77 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 67.8378 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.54033e-14 -***** Scaled Norm of nonlinear residual: 1.40634e-15 -********************************************************* -Norm of solution: 31.2654 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 79 -**** Current time is 0.78 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 71.5615 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.65754e-14 -***** Scaled Norm of nonlinear residual: 1.2098e-15 -********************************************************* -Norm of solution: 30.7987 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 80 -**** Current time is 0.79 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 75.0148 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.85548e-14 -***** Scaled Norm of nonlinear residual: 1.04719e-15 -********************************************************* -Norm of solution: 30.301 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 81 -**** Current time is 0.8 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 78.1096 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 6.97871e-14 -***** Scaled Norm of nonlinear residual: 8.93451e-16 -********************************************************* -Norm of solution: 29.7696 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 82 -**** Current time is 0.81 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 80.7645 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.08927e-13 -***** Scaled Norm of nonlinear residual: 1.3487e-15 -********************************************************* -Norm of solution: 29.2008 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 83 -**** Current time is 0.82 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 82.9074 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.88222e-14 -***** Scaled Norm of nonlinear residual: 1.19196e-15 -********************************************************* -Norm of solution: 28.5907 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 84 -**** Current time is 0.83 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 84.4768 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.10744e-13 -***** Scaled Norm of nonlinear residual: 1.31094e-15 -********************************************************* -Norm of solution: 27.9346 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 85 -**** Current time is 0.84 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 85.4233 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.23295e-14 -***** Scaled Norm of nonlinear residual: 8.46718e-16 -********************************************************* -Norm of solution: 27.2281 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 86 -**** Current time is 0.85 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 85.7108 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.63429e-14 -***** Scaled Norm of nonlinear residual: 8.90703e-16 -********************************************************* -Norm of solution: 26.4669 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 87 -**** Current time is 0.86 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 85.3169 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 8.50603e-14 -***** Scaled Norm of nonlinear residual: 9.96993e-16 -********************************************************* -Norm of solution: 25.6474 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 88 -**** Current time is 0.87 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 84.2332 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.75182e-14 -***** Scaled Norm of nonlinear residual: 1.15772e-15 -********************************************************* -Norm of solution: 24.7662 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 89 -**** Current time is 0.88 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 82.465 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.06276e-13 -***** Scaled Norm of nonlinear residual: 1.28874e-15 -********************************************************* -Norm of solution: 23.8211 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 90 -**** Current time is 0.89 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 80.0307 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.26184e-14 -***** Scaled Norm of nonlinear residual: 1.15729e-15 -********************************************************* -Norm of solution: 22.8106 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 91 -**** Current time is 0.9 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 76.9611 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.35901e-14 -***** Scaled Norm of nonlinear residual: 1.21607e-15 -********************************************************* -Norm of solution: 21.734 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 92 -**** Current time is 0.91 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 73.2977 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.04835e-13 -***** Scaled Norm of nonlinear residual: 1.43026e-15 -********************************************************* -Norm of solution: 20.5918 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 93 -**** Current time is 0.92 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 69.0916 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.68545e-14 -***** Scaled Norm of nonlinear residual: 1.40183e-15 -********************************************************* -Norm of solution: 19.3851 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 94 -**** Current time is 0.93 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 64.4014 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.39985e-14 -***** Scaled Norm of nonlinear residual: 1.14902e-15 -********************************************************* -Norm of solution: 18.1163 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 95 -**** Current time is 0.94 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 59.2919 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 9.45555e-14 -***** Scaled Norm of nonlinear residual: 1.59475e-15 -********************************************************* -Norm of solution: 16.7883 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 96 -**** Current time is 0.95 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 53.8318 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.01409e-13 -***** Scaled Norm of nonlinear residual: 1.88382e-15 -********************************************************* -Norm of solution: 15.405 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 97 -**** Current time is 0.96 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 48.0921 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 7.8183e-14 -***** Scaled Norm of nonlinear residual: 1.62569e-15 -********************************************************* -Norm of solution: 13.9713 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 98 -**** Current time is 0.97 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 42.6908 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.03987e-13 -***** Scaled Norm of nonlinear residual: 2.43581e-15 -********************************************************* -Norm of solution: 12.493 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 99 -**** Current time is 0.98 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 37.1893 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.19822e-13 -***** Scaled Norm of nonlinear residual: 3.22194e-15 -********************************************************* -Norm of solution: 10.9774 (overlapped vector so results may differ on multiple procs) - - -******************************************************* - -**** Beginning Time Step 100 -**** Current time is 0.99 - -******************************************************* - - - -********************************************************* -***** Iteration: 0 -***** Norm of nonlinear residual: 31.5601 -***** Scaled Norm of nonlinear residual: 1 -********************************************************* - -********************************************************* -***** Iteration: 1 -***** Norm of nonlinear residual: 1.3165e-13 -***** Scaled Norm of nonlinear residual: 4.1714e-15 -********************************************************* -Norm of solution: 9.43367 (overlapped vector so results may differ on multiple procs) diff --git a/sandbox/testCases/hgrad_dirichlet_test/debug.out b/sandbox/testCases/hgrad_dirichlet_test/debug.out deleted file mode 100644 index db1534a4..00000000 --- a/sandbox/testCases/hgrad_dirichlet_test/debug.out +++ /dev/null @@ -1,6 +0,0 @@ -[setDirichletData DEBUG] Block=eblock-0_0_0, Var=T, Sideset=bottom, sideEntities.size()=4 -[setDirichletData DEBUG] side_output.size()=4 (elements touching sideset) -[setDirichletData DEBUG] Element 0: local_side_Id=0, getGIDFieldOffsets_closure returned 4 offsets: [0,1,4,5] -[setDirichletData DEBUG] Element 1: local_side_Id=0, getGIDFieldOffsets_closure returned 4 offsets: [0,1,4,5] -[setDirichletData DEBUG] Element 2: local_side_Id=0, getGIDFieldOffsets_closure returned 4 offsets: [0,1,4,5] -[setDirichletData DEBUG] FINAL: Block=eblock-0_0_0, Var=T, unique Dirichlet DOFs marked=9 diff --git a/sandbox/testCases/hgrad_dirichlet_test/input.yaml b/sandbox/testCases/hgrad_dirichlet_test/input.yaml deleted file mode 100644 index a32cd5db..00000000 --- a/sandbox/testCases/hgrad_dirichlet_test/input.yaml +++ /dev/null @@ -1,43 +0,0 @@ -%YAML 1.1 ---- -# Simple test case for HGRAD Dirichlet boundary conditions -# T field is HGRAD, we apply T = 0 on bottom face -ANONYMOUS: - debug level: 0 - verbosity: 10 # High verbosity to enable debug output in setDirichletData - Mesh: - dimension: 3 - element type: hex - xmin: 0.0 - xmax: 1.0 - ymin: 0.0 - ymax: 1.0 - zmin: 0.0 - zmax: 1.0 - NX: 2 - NY: 2 - NZ: 2 - Physics: - eblock-0_0_0: - modules: thermal - Dirichlet conditions: - T: - bottom: '0.0' - Initial conditions: - T: '0.0' - Discretization: - eblock-0_0_0: - order: - T: 1 - quadrature: 2 - Functions: - thermal source: '0.0' - Solver: - solver: steady-state - Analysis: - analysis type: forward - Postprocess: - compute errors: false - write solution: false -... - diff --git a/sandbox/testCases/hgrad_dirichlet_test/input_mesh.yaml b/sandbox/testCases/hgrad_dirichlet_test/input_mesh.yaml deleted file mode 100644 index a0fafc37..00000000 --- a/sandbox/testCases/hgrad_dirichlet_test/input_mesh.yaml +++ /dev/null @@ -1,21 +0,0 @@ -%YAML 1.1 ---- -ANONYMOUS: - Mesh: - dimension: 3 - shape: hex - xmin: 0.0 - xmax: 1.0 - ymin: 0.0 - ymax: 1.0 - zmin: 0.0 - zmax: 1.0 - NX: 2 - NY: 2 - NZ: 2 - Xprocs: 1 - Yprocs: 1 - Zprocs: 1 - blocknames: eblock-0_0_0 -... - diff --git a/src/managers/postprocess/postprocessManager_integrated_quantities.hpp b/src/managers/postprocess/postprocessManager_integrated_quantities.hpp index ba5a5ef3..c4c699b9 100644 --- a/src/managers/postprocess/postprocessManager_integrated_quantities.hpp +++ b/src/managers/postprocess/postprocessManager_integrated_quantities.hpp @@ -249,6 +249,18 @@ void PostprocessManager::computeIntegratedQuantities(vector &c integratedQuantities[iLocal][iIQ].val(0) = allsums[iIQ]; } // end loop over integrated quantities + // Print integrated quantities at each time step for debugging + if (Comm->getRank() == 0 && verbosity >= 5) + { + std::cout << "Time = " << currenttime << ", Integrated quantities:"; + for (size_t k = 0; k < integratedQuantities[iLocal].size(); ++k) + { + std::cout << " " << integratedQuantities[iLocal][k].name << " = " + << integratedQuantities[iLocal][k].val(0); + } + std::cout << std::endl; + } + // need to put in the right place now (accessible to the residual) and // update any parameters which depend on the IQs // TODO :: BWR this ultimately is an "explicit" idea but doing things implicitly