Commit e1d7c7b2 authored by Alex Breuer's avatar Alex Breuer

Added initial and reference data to continuous tests for the advection...

Added initial and reference data to continuous tests for the advection equation using line elements.
parent f91e13ad
......@@ -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:
......@@ -31,16 +31,138 @@
<cfr>
<mesh>
<n_elements>
<x>177</x>
<x>250</x>
</n_elements>
<size>
<x>13</x>
<x>5</x>
</size>
</mesh>
<setups>
<type_default>sine</type_default>
<end_time>13</end_time>
<initial_values>
<!-- smooth initial data, 2 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*2*pi;
q[0] := sin(xf*2);
]]>
<!-- smooth initial data, 4 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*3*pi;
q[0] := sin(xf*2);
]]>
<!-- smooth initial data, 8 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*3*pi;
q[0] := sin(xf*4);
]]>
<!-- smooth initial data, 16 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*3*pi;
q[0] := sin(xf*4);
]]>
<!-- discontinuous initial data, 1 sine wave, 1 jump -->
<![CDATA[
var k := 1;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
<!-- discontinuous initial data, 2 sine waves, 2 jumps -->
<![CDATA[
var k := 2;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
<!-- discontinuous initial data, 4 sine waves, 4 jumps -->
<![CDATA[
var k := 4;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
<!-- discontinuous initial data, 8 sine waves, 8 jumps -->
<![CDATA[
var k := 8;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
</initial_values>
<end_time>25</end_time>
</setups>
<output>
......@@ -50,6 +172,129 @@
<int>0.5</int>
</wave_field>
<error_norms>
<reference_values>
<!-- smooth initial data, 2 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*2*pi;
q[0] := sin(xf*2);
]]>
<!-- smooth initial data, 4 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*3*pi;
q[0] := sin(xf*2);
]]>
<!-- smooth initial data, 8 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*3*pi;
q[0] := sin(xf*4);
]]>
<!-- smooth initial data, 16 sine waves -->
<![CDATA[
var xf := ( x / 5.0 )*3*pi;
q[0] := sin(xf*4);
]]>
<!-- discontinuous initial data, 1 sine wave, 1 jump -->
<![CDATA[
var k := 1;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
<!-- discontinuous initial data, 2 sine waves, 2 jumps -->
<![CDATA[
var k := 2;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
<!-- discontinuous initial data, 4 sine waves, 4 jumps -->
<![CDATA[
var k := 4;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
<!-- discontinuous initial data, 8 sine waves, 8 jumps -->
<![CDATA[
var k := 8;
var xf := ( x / 5.0 )*2*pi;
if( xf < pi ) {
q[0] := sin(xf*k*2);
}
for( var ju := 0; ju < k; ju += 1) {
var int := 5.0/( 2*(2*k+1) );
var lo := 5.0/2.0 + (ju*2+1)*int;
var up := 5.0/2.0 + (ju*2+2)*int;
if( x > lo and x < up ) {
if( ju%2 == 0 )
q[0] := 1;
else
q[0] := -1;
}
}
]]>
</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