Commit 5e4bd22c authored by Alex Breuer's avatar Alex Breuer
Browse files

Added initialization and error normal evaluation to 2D continuity configs (non-fused only).

parent 14f32f94
......@@ -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,8 +40,59 @@
</size>
</mesh>
<velocity_model>
<domain>
<half_space>
<origin>
<x>0</x>
<y>-9999</y>
</origin>
<normal>
<x>0</x>
<y>1</y>
</normal>
</half_space>
<rho>1</rho>
<lambda>2</lambda>
<mu>1</mu>
</domain>
</velocity_model>
<setups>
<type_default>plane_waves</type_default>
<initial_values>
<![CDATA[
var n[2] := { 1.0/sqrt(2.0), 1.0/sqrt(2.0) };
var k := ( 2.0 * pi ) / 25.0;
var scale := sin( x*k + y*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 r2[5] := {};
var r5[5] := {};
r2[0] := -2.0 * mu * n[0] * n[1];
r2[1] := 2.0 * mu * n[0] * n[1];
r2[2] := mu * ( n[0]*n[0] - n[1]*n[1]);
r2[3] := -n[1] * cs;
r2[4] := n[0] * cs;
r5[0] := lambda + 2.0 * mu * n[0] * n[0];
r5[1] := lambda + 2.0 * mu * n[1] * n[1];
r5[2] := 2.0 * mu * n[0] * n[1];
r5[3] := -n[0] * cp;
r5[4] := -n[1] * cp;
q := (r2 + r5) * scale;
]]>
</initial_values>
<end_time>141.421356237309504880168872420969807856967187537694</end_time>
</setups>
......@@ -52,6 +103,39 @@
<int>10.0</int>
</wave_field>
<error_norms>
<reference_values>
<![CDATA[
var n[2] := { 1.0/sqrt(2.0), 1.0/sqrt(2.0) };
var k := ( 2.0 * pi ) / 25.0;
var scale := sin( x*k + y*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 r2[5] := {};
var r5[5] := {};
r2[0] := -2.0 * mu * n[0] * n[1];
r2[1] := 2.0 * mu * n[0] * n[1];
r2[2] := mu * ( n[0]*n[0] - n[1]*n[1]);
r2[3] := -n[1] * cs;
r2[4] := n[0] * cs;
r5[0] := lambda + 2.0 * mu * n[0] * n[0];
r5[1] := lambda + 2.0 * mu * n[1] * n[1];
r5[2] := 2.0 * mu * n[0] * n[1];
r5[3] := -n[0] * cp;
r5[4] := -n[1] * cp;
q := (r2 + r5) * scale;
]]>
</reference_values>
<type>sout_file</type>
<file>errors.xml</file>
</error_norms>
......
......@@ -40,7 +40,39 @@
</mesh>
<setups>
<type_default>plane_waves</type_default>
<initial_values>
<![CDATA[
var n[2] := { 1.0/sqrt(2.0), 1.0/sqrt(2.0) };
var k := ( 2.0 * pi ) / 25.0;
var scale := sin( x*k + y*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 r2[5] := {};
var r5[5] := {};
r2[0] := -2.0 * mu * n[0] * n[1];
r2[1] := 2.0 * mu * n[0] * n[1];
r2[2] := mu * ( n[0]*n[0] - n[1]*n[1]);
r2[3] := -n[1] * cs;
r2[4] := n[0] * cs;
r5[0] := lambda + 2.0 * mu * n[0] * n[0];
r5[1] := lambda + 2.0 * mu * n[1] * n[1];
r5[2] := 2.0 * mu * n[0] * n[1];
r5[3] := -n[0] * cp;
r5[4] := -n[1] * cp;
q := (r2 + r5) * scale;
]]>
</initial_values>
<end_time>141.421356237309504880168872420969807856967187537694</end_time>
</setups>
......@@ -51,6 +83,39 @@
<int>10.0</int>
</wave_field>
<error_norms>
<reference_values>
<![CDATA[
var n[2] := { 1.0/sqrt(2.0), 1.0/sqrt(2.0) };
var k := ( 2.0 * pi ) / 25.0;
var scale := sin( x*k + y*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 r2[5] := {};
var r5[5] := {};
r2[0] := -2.0 * mu * n[0] * n[1];
r2[1] := 2.0 * mu * n[0] * n[1];
r2[2] := mu * ( n[0]*n[0] - n[1]*n[1]);
r2[3] := -n[1] * cs;
r2[4] := n[0] * cs;
r5[0] := lambda + 2.0 * mu * n[0] * n[0];
r5[1] := lambda + 2.0 * mu * n[1] * n[1];
r5[2] := 2.0 * mu * n[0] * n[1];
r5[3] := -n[0] * cp;
r5[4] := -n[1] * cp;
q := (r2 + r5) * scale;
]]>
</reference_values>
<type>sout_file</type>
<file>errors.xml</file>
</error_norms>
......
......@@ -41,7 +41,39 @@
</mesh>
<setups>
<type_default>plane_waves</type_default>
<initial_values>
<![CDATA[
var n[2] := { 1.0/sqrt(2.0), 1.0/sqrt(2.0) };
var k := ( 2.0 * pi ) / 25.0;
var scale := sin( x*k + y*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 r2[5] := {};
var r5[5] := {};
r2[0] := -2.0 * mu * n[0] * n[1];
r2[1] := 2.0 * mu * n[0] * n[1];
r2[2] := mu * ( n[0]*n[0] - n[1]*n[1]);
r2[3] := -n[1] * cs;
r2[4] := n[0] * cs;
r5[0] := lambda + 2.0 * mu * n[0] * n[0];
r5[1] := lambda + 2.0 * mu * n[1] * n[1];
r5[2] := 2.0 * mu * n[0] * n[1];
r5[3] := -n[0] * cp;
r5[4] := -n[1] * cp;
q := (r2 + r5) * scale;
]]>
</initial_values>
<end_time>141.421356237309504880168872420969807856967187537694</end_time>
</setups>
......@@ -52,6 +84,39 @@
<int>10.0</int>
</wave_field>
<error_norms>
<reference_values>
<![CDATA[
var n[2] := { 1.0/sqrt(2.0), 1.0/sqrt(2.0) };
var k := ( 2.0 * pi ) / 25.0;
var scale := sin( x*k + y*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 r2[5] := {};
var r5[5] := {};
r2[0] := -2.0 * mu * n[0] * n[1];
r2[1] := 2.0 * mu * n[0] * n[1];
r2[2] := mu * ( n[0]*n[0] - n[1]*n[1]);
r2[3] := -n[1] * cs;
r2[4] := n[0] * cs;
r5[0] := lambda + 2.0 * mu * n[0] * n[0];
r5[1] := lambda + 2.0 * mu * n[1] * n[1];
r5[2] := 2.0 * mu * n[0] * n[1];
r5[3] := -n[0] * cp;
r5[4] := -n[1] * cp;
q := (r2 + r5) * 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