Commit c6ab425f authored by Alex Breuer's avatar Alex Breuer
Browse files

Adjusted continuity configs to runtime initialization and error reference values.

parent 5e4bd22c
......@@ -5,7 +5,7 @@
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2016-2017, Regents of the University of California
Copyright (c) 2016-2018, Regents of the University of California
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
......@@ -42,8 +42,72 @@
</size>
</mesh>
<velocity_model>
<domain>
<half_space>
<origin>
<x>0</x>
<y>0</y>
<z>-9999</z>
</origin>
<normal>
<x>0</x>
<y>0</y>
<z>1</z>
</normal>
</half_space>
<rho>1</rho>
<lambda>2</lambda>
<mu>1</mu>
</domain>
</velocity_model>
<setups>
<type_default>plane_waves</type_default>
<!-- Sets up plane waves as a combination of eigenvectors in diagonal direction n=(1,1,1).
The choice of r8 is non-unique, since cs has two eigenvalues.
-->
<initial_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</initial_values>
<end_time>173.205080756887729352744634150587236694280525381038</end_time>
</setups>
......@@ -54,6 +118,48 @@
<int>10.0</int>
</wave_field>
<error_norms>
<!-- Reference values are identical to initial values-->
<reference_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</reference_values>
<type>sout_file</type>
<file>errors.xml</file>
</error_norms>
......
......@@ -5,7 +5,7 @@
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2016-2017, Regents of the University of California
Copyright (c) 2016-2018, Regents of the University of California
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
......@@ -42,8 +42,72 @@
</size>
</mesh>
<velocity_model>
<domain>
<half_space>
<origin>
<x>0</x>
<y>0</y>
<z>-9999</z>
</origin>
<normal>
<x>0</x>
<y>0</y>
<z>1</z>
</normal>
</half_space>
<rho>1</rho>
<lambda>2</lambda>
<mu>1</mu>
</domain>
</velocity_model>
<setups>
<type_default>plane_waves</type_default>
<!-- Sets up plane waves as a combination of eigenvectors in diagonal direction n=(1,1,1).
The choice of r8 is non-unique, since cs has two eigenvalues.
-->
<initial_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</initial_values>
<end_time>173.205080756887729352744634150587236694280525381038</end_time>
</setups>
......@@ -54,6 +118,48 @@
<int>10</int>
</wave_field>
<error_norms>
<!-- Reference values are identical to initial values-->
<reference_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</reference_values>
<type>sout_file</type>
<file>errors.xml</file>
</error_norms>
......
......@@ -5,7 +5,7 @@
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2016-2017, Regents of the University of California
Copyright (c) 2016-2018, Regents of the University of California
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
......@@ -39,7 +39,50 @@
</mesh>
<setups>
<type_default>plane_waves</type_default>
<!-- Sets up plane waves as a combination of eigenvectors in diagonal direction n=(1,1,1).
The choice of r8 is non-unique, since cs has two eigenvalues.
-->
<initial_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</initial_values>
<end_time>173.205080756887729352744634150587236694280525381038</end_time>
</setups>
......@@ -50,6 +93,48 @@
<int>10.0</int>
</wave_field>
<error_norms>
<!-- Reference values are identical to initial values-->
<reference_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</reference_values>
<type>sout_file</type>
<file>errors.xml</file>
</error_norms>
......
......@@ -5,7 +5,7 @@
@author Alexander Breuer (anbreuer AT ucsd.edu)
@section LICENSE
Copyright (c) 2016-2017, Regents of the University of California
Copyright (c) 2016-2018, Regents of the University of California
All rights reserved.
Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
......@@ -40,7 +40,50 @@
</mesh>
<setups>
<type_default>plane_waves</type_default>
<!-- Sets up plane waves as a combination of eigenvectors in diagonal direction n=(1,1,1).
The choice of r8 is non-unique, since cs has two eigenvalues.
-->
<initial_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</initial_values>
<end_time>173.205080756887729352744634150587236694280525381038</end_time>
</setups>
......@@ -51,6 +94,48 @@
<int>10.0</int>
</wave_field>
<error_norms>
<!-- Reference values are identical to initial values-->
<reference_values>
<![CDATA[
var sqrtI := { 1.0/sqrt(3.0) };
var k := ( 2.0 * pi ) / 50.0;
var scale := sin( x*k + y*k + z*k );
var mu := 1.0;
var lambda := 2.0;
var rho := 1.0;
var cs := sqrt( mu / rho );
var cp := sqrt( (lambda + 2*mu) / rho );
var r1[9] := {};
var r8[9] := {};
r1[0] := 4.0 * sqrtI;
r1[1] := 4.0 * sqrtI;
r1[2] := 4.0 * sqrtI;
r1[3] := 1.0 * sqrtI;
r1[4] := 1.0 * sqrtI;
r1[5] := 1.0 * sqrtI;
r1[6] := 1.0;
r1[7] := 1.0;
r1[8] := 1.0;
r8[0] := 2.0 * sqrtI;
r8[1] := 0.0;
r8[2] := -2.0 * sqrtI;
r8[3] := 1.0 * sqrtI;
r8[4] := -1.0 * sqrtI;
r8[5] := 0.0;
r8[6] := -1.0;
r8[7] := 0.0;
r8[8] := 1.0;
q := (r1 + r8) * scale;
]]>
</reference_values>
<type>sout_file</type>
<file>errors.xml</file>
</error_norms>
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment