Difference between revisions of "Small Scale Heterogeneities Validation"

From SCECpedia
Jump to navigationJump to search
 
(9 intermediate revisions by the same user not shown)
Line 3: Line 3:
 
We would like to validate that the media with the small-scale heterogeneities that are being generated with the new C code are correct. One way to do this would be to take a look at some plots to see if the generated heterogeneities look correct. We present a few different meshes and plots showing the heterogeneities at various depths, as well as difference plots.
 
We would like to validate that the media with the small-scale heterogeneities that are being generated with the new C code are correct. One way to do this would be to take a look at some plots to see if the generated heterogeneities look correct. We present a few different meshes and plots showing the heterogeneities at various depths, as well as difference plots.
  
== Validation by Hurst Value ==
+
Common parameters:
 
 
=== Mesh 0.0 Hurst ===
 
 
 
Specifications:
 
  
Hurst = 0.0<br />
 
L1 = 200<br />
 
 
d = 20<br />
 
d = 20<br />
 
st23 = 5<br />
 
st23 = 5<br />
Line 16: Line 10:
 
n2 = 1000<br />
 
n2 = 1000<br />
 
n3 = 1500<br />
 
n3 = 1500<br />
 +
 +
== Validation by Hurst Value ==
 +
 +
In this test, we vary hurst to be either 0.0, 0.2, or 0.4. L1 is held constant at 200m.
  
 
==== Plots at surface level ====
 
==== Plots at surface level ====
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 0.0 Hurst'''
| '''C SSH Code'''
+
| '''C SSH Code 0.0 Hurst'''
 +
| '''Difference 0.0 Hurst'''
 
|-
 
|-
 
| [[File:matlab_H0_0.png|600px]]
 
| [[File:matlab_H0_0.png|600px]]
 
| [[File:myssh_H0_0.png|600px]]
 
| [[File:myssh_H0_0.png|600px]]
 +
| [[File:compare_H0_0.png|600px]]
 +
|}
 +
 +
{|
 +
| '''Matlab Code 0.2 Hurst'''
 +
| '''C SSH Code 0.2 Hurst'''
 +
| '''Difference 0.2 Hurst'''
 
|-
 
|-
| '''Difference'''
+
| [[File:matlab_H2_0.png|600px]]
|
+
| [[File:myssh_H2_0.png|600px]]
 +
| [[File:compare_H2_0.png|600px]]
 +
|}
 +
 
 +
{|
 +
| '''Matlab Code 0.4 Hurst'''
 +
| '''C SSH Code 0.4 Hurst'''
 +
| '''Difference 0.4 Hurst'''
 
|-
 
|-
| [[File:compare_H0_0.png|600px]]
+
| [[File:matlab_H4_0.png|600px]]
|
+
| [[File:myssh_H4_0.png|600px]]
 +
| [[File:compare_H4_0.png|600px]]
 
|}
 
|}
  
Line 36: Line 50:
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 0.0 Hurst'''
| '''C SSH Code'''
+
| '''C SSH Code 0.0 Hurst'''
 +
| '''Difference 0.0 Hurst'''
 
|-
 
|-
 
| [[File:matlab_H0_250.png|600px]]
 
| [[File:matlab_H0_250.png|600px]]
 
| [[File:myssh_H0_250.png|600px]]
 
| [[File:myssh_H0_250.png|600px]]
 +
| [[File:compare_H0_250.png|600px]]
 +
|}
 +
 +
{|
 +
| '''Matlab Code 0.2 Hurst'''
 +
| '''C SSH Code 0.2 Hurst'''
 +
| '''Difference 0.2 Hurst'''
 
|-
 
|-
| '''Difference'''
+
| [[File:matlab_H2_250.png|600px]]
|
+
| [[File:myssh_H2_250.png|600px]]
 +
| [[File:compare_H2_250.png|600px]]
 +
|}
 +
 
 +
{|
 +
| '''Matlab Code 0.4 Hurst'''
 +
| '''C SSH Code 0.4 Hurst'''
 +
| '''Difference 0.4 Hurst'''
 
|-
 
|-
| [[File:compare_H0_250.png|600px]]
+
| [[File:matlab_H4_250.png|600px]]
|
+
| [[File:myssh_H4_250.png|600px]]
 +
| [[File:compare_H4_250.png|600px]]
 
|}
 
|}
  
Line 52: Line 82:
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 0.0 Hurst'''
| '''C SSH Code'''
+
| '''C SSH Code 0.0 Hurst'''
 +
| '''Difference 0.0 Hurst'''
 
|-
 
|-
 
| [[File:matlab_H0_500.png|600px]]
 
| [[File:matlab_H0_500.png|600px]]
 
| [[File:myssh_H0_500.png|600px]]
 
| [[File:myssh_H0_500.png|600px]]
 +
| [[File:compare_H0_500.png|600px]]
 +
|}
 +
 +
{|
 +
| '''Matlab Code 0.2 Hurst'''
 +
| '''C SSH Code 0.2 Hurst'''
 +
| '''Difference 0.2 Hurst'''
 
|-
 
|-
| '''Difference'''
+
| [[File:matlab_H2_500.png|600px]]
|
+
| [[File:myssh_H2_500.png|600px]]
 +
| [[File:compare_H2_500.png|600px]]
 +
|}
 +
 
 +
{|
 +
| '''Matlab Code 0.4 Hurst'''
 +
| '''C SSH Code 0.4 Hurst'''
 +
| '''Difference 0.4 Hurst'''
 
|-
 
|-
| [[File:compare_H0_499.png|600px]]
+
| [[File:matlab_H4_500.png|600px]]
|
+
| [[File:myssh_H4_500.png|600px]]
 +
| [[File:compare_H4_500.png|600px]]
 
|}
 
|}
  
=== Mesh 0.2 Hurst ===
+
== Validation by Correlation Length ==
  
Specifications:
+
In this condition, we hold hurst to be a constant 0.2, but vary correlation length to be either 100m, 200m, or 300m.
  
Hurst = 0.2<br />
+
==== Plots at surface level ====
L1 = 200<br />
 
d = 20<br />
 
st23 = 5<br />
 
n1 = 500<br />
 
n2 = 1000<br />
 
n3 = 1500<br />
 
  
==== Plots at surface level ====
+
{|
 +
| '''Matlab Code 100m L1'''
 +
| '''C SSH Code 100m L1'''
 +
| '''Difference 100m L1'''
 +
|-
 +
| [[File:matlab_H0_0.png|600px]]
 +
| [[File:myssh_H0_0.png|600px]]
 +
| [[File:compare_H0_0.png|600px]]
 +
|}
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 200m L1'''
| '''C SSH Code'''
+
| '''C SSH Code 200m L1'''
 +
| '''Difference 200m L1'''
 
|-
 
|-
 
| [[File:matlab_H2_0.png|600px]]
 
| [[File:matlab_H2_0.png|600px]]
 
| [[File:myssh_H2_0.png|600px]]
 
| [[File:myssh_H2_0.png|600px]]
 +
| [[File:compare_H2_0.png|600px]]
 +
|}
 +
 +
{|
 +
| '''Matlab Code 300m L1'''
 +
| '''C SSH Code 300m L1'''
 +
| '''Difference 300m L1'''
 
|-
 
|-
| '''Difference'''
+
| [[File:matlab_L300_0.png|600px]]
|
+
| [[File:myssh_L300_0.png|600px]]
|-
+
| [[File:compare_L300_0.png|600px]]
| [[File:compare_H2_0.png|600px]]
 
|
 
 
|}
 
|}
  
Line 96: Line 150:
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 100m L1'''
| '''C SSH Code'''
+
| '''C SSH Code 100m L1'''
 +
| '''Difference 100m L1'''
 +
|-
 +
| [[File:matlab_L100_250.png|600px]]
 +
| [[File:myssh_L100_250.png|600px]]
 +
| [[File:compare_L100_250.png|600px]]
 +
|}
 +
 
 +
{|
 +
| '''Matlab Code 200m L1'''
 +
| '''C SSH Code 200m L1'''
 +
| '''Difference 200m L1'''
 
|-
 
|-
 
| [[File:matlab_H2_250.png|600px]]
 
| [[File:matlab_H2_250.png|600px]]
 
| [[File:myssh_H2_250.png|600px]]
 
| [[File:myssh_H2_250.png|600px]]
|-
 
| '''Difference'''
 
|
 
|-
 
 
| [[File:compare_H2_250.png|600px]]
 
| [[File:compare_H2_250.png|600px]]
|
 
 
|}
 
|}
 
==== Plots for bottom slice ====
 
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 300m L1'''
| '''C SSH Code'''
+
| '''C SSH Code 300m L1'''
|-
+
| '''Difference 300m L1'''
| [[File:matlab_H2_500.png|600px]]
 
| [[File:myssh_H2_500.png|600px]]
 
|-
 
| '''Difference'''
 
|
 
 
|-
 
|-
| [[File:compare_H2_500.png|600px]]
+
| [[File:matlab_L300_250.png|600px]]
|
+
| [[File:myssh_L300_250.png|600px]]
 +
| [[File:compare_L300_250.png|600px]]
 
|}
 
|}
  
=== Mesh 0.4 Hurst ===
+
==== Plots for bottom slice ====
 
 
Specifications:
 
 
 
Hurst = 0.4<br />
 
L1 = 200<br />
 
d = 20<br />
 
st23 = 5<br />
 
n1 = 500<br />
 
n2 = 1000<br />
 
n3 = 1500<br />
 
 
 
==== Plots at surface level ====
 
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 100m L1'''
| '''C SSH Code'''
+
| '''C SSH Code 100m L1'''
|-
+
| '''Difference 100m L1'''
| [[File:matlab_H4_0.png|600px]]
 
| [[File:myssh_H4_0.png|600px]]
 
|-
 
| '''Difference'''
 
|
 
 
|-
 
|-
| [[File:compare_H4_0.png|600px]]
+
| [[File:matlab_L100_500.png|600px]]
|
+
| [[File:myssh_L100_500.png|600px]]
 +
| [[File:compare_L100_500.png|600px]]
 
|}
 
|}
 
==== Plots at 50% depth ====
 
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 200m L1'''
| '''C SSH Code'''
+
| '''C SSH Code 200m L1'''
 +
| '''Difference 200m L1'''
 
|-
 
|-
| [[File:matlab_H4_250.png|600px]]
+
| [[File:matlab_H2_500.png|600px]]
| [[File:myssh_H4_250.png|600px]]
+
| [[File:myssh_H2_500.png|600px]]
|-
+
| [[File:compare_H2_500.png|600px]]
| '''Difference'''
 
|
 
|-
 
| [[File:compare_H4_250.png|600px]]
 
|
 
 
|}
 
|}
 
==== Plots for bottom slice ====
 
  
 
{|
 
{|
| '''Matlab Code'''
+
| '''Matlab Code 300m L1'''
| '''C SSH Code'''
+
| '''C SSH Code 300m L1'''
|-
+
| '''Difference 300m L1'''
| [[File:matlab_H4_500.png|600px]]
 
| [[File:myssh_H4_500.png|600px]]
 
 
|-
 
|-
| '''Difference'''
+
| [[File:matlab_L300_500.png|600px]]
|
+
| [[File:myssh_L300_500.png|600px]]
|-
+
| [[File:compare_L300_500.png|600px]]
| [[File:compare_H4_500.png|600px]]
 
|
 
 
|}
 
|}

Latest revision as of 23:40, 5 March 2014

Introduction

We would like to validate that the media with the small-scale heterogeneities that are being generated with the new C code are correct. One way to do this would be to take a look at some plots to see if the generated heterogeneities look correct. We present a few different meshes and plots showing the heterogeneities at various depths, as well as difference plots.

Common parameters:

d = 20
st23 = 5
n1 = 500
n2 = 1000
n3 = 1500

Validation by Hurst Value

In this test, we vary hurst to be either 0.0, 0.2, or 0.4. L1 is held constant at 200m.

Plots at surface level

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab H0 0.png Myssh H0 0.png Compare H0 0.png
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab H2 0.png Myssh H2 0.png Compare H2 0.png
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst
Matlab H4 0.png Myssh H4 0.png Compare H4 0.png

Plots at 50% depth

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab H0 250.png Myssh H0 250.png Compare H0 250.png
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab H2 250.png Myssh H2 250.png Compare H2 250.png
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst
Matlab H4 250.png Myssh H4 250.png Compare H4 250.png

Plots for bottom slice

Matlab Code 0.0 Hurst C SSH Code 0.0 Hurst Difference 0.0 Hurst
Matlab H0 500.png Myssh H0 500.png Compare H0 500.png
Matlab Code 0.2 Hurst C SSH Code 0.2 Hurst Difference 0.2 Hurst
Matlab H2 500.png Myssh H2 500.png Compare H2 500.png
Matlab Code 0.4 Hurst C SSH Code 0.4 Hurst Difference 0.4 Hurst
Matlab H4 500.png Myssh H4 500.png Compare H4 500.png

Validation by Correlation Length

In this condition, we hold hurst to be a constant 0.2, but vary correlation length to be either 100m, 200m, or 300m.

Plots at surface level

Matlab Code 100m L1 C SSH Code 100m L1 Difference 100m L1
Matlab H0 0.png Myssh H0 0.png Compare H0 0.png
Matlab Code 200m L1 C SSH Code 200m L1 Difference 200m L1
Matlab H2 0.png Myssh H2 0.png Compare H2 0.png
Matlab Code 300m L1 C SSH Code 300m L1 Difference 300m L1
Matlab L300 0.png Myssh L300 0.png Compare L300 0.png

Plots at 50% depth

Matlab Code 100m L1 C SSH Code 100m L1 Difference 100m L1
Matlab L100 250.png Myssh L100 250.png Compare L100 250.png
Matlab Code 200m L1 C SSH Code 200m L1 Difference 200m L1
Matlab H2 250.png Myssh H2 250.png Compare H2 250.png
Matlab Code 300m L1 C SSH Code 300m L1 Difference 300m L1
Matlab L300 250.png Myssh L300 250.png Compare L300 250.png

Plots for bottom slice

Matlab Code 100m L1 C SSH Code 100m L1 Difference 100m L1
Matlab L100 500.png Myssh L100 500.png Compare L100 500.png
Matlab Code 200m L1 C SSH Code 200m L1 Difference 200m L1
Matlab H2 500.png Myssh H2 500.png Compare H2 500.png
Matlab Code 300m L1 C SSH Code 300m L1 Difference 300m L1
Matlab L300 500.png Myssh L300 500.png Compare L300 500.png