Commit 9a04ce33 authored by Dave Moxey's avatar Dave Moxey
Browse files

Add MultiRegions pyramid regression tests

parent 24f44507
......@@ -125,6 +125,8 @@ ADD_NEKTAR_TEST(Helmholtz3D_CG_Tet)
ADD_NEKTAR_TEST(Helmholtz3D_CG_Prism)
ADD_NEKTAR_TEST(Helmholtz3D_CG_Prism_iter_ml)
ADD_NEKTAR_TEST(Helmholtz3D_CG_Prism_Deformed)
ADD_NEKTAR_TEST(Helmholtz3D_CG_Pyr)
ADD_NEKTAR_TEST(Helmholtz3D_CG_Pyr_Deformed)
ADD_NEKTAR_TEST(Helmholtz3D_CG_Homo1D)
ADD_NEKTAR_TEST_LENGTHY(Helmholtz3D_HDG_Homo1D)
ADD_NEKTAR_TEST(Helmholtz3D_HDG_Prism)
......
<?xml version="1.0" encoding="utf-8" ?>
<test>
<description>Helmholtz 3D CG for Pyramid</description>
<executable>Helmholtz3D</executable>
<parameters>Helmholtz3D_Pyr.xml</parameters>
<files>
<file description="Session File">Helmholtz3D_Pyr.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value tolerance="1e-12">5.36868e-07</value>
</metric>
<metric type="Linf" id="2">
<value tolerance="1e-12">3.78881e-06</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8" ?>
<test>
<description>Helmholtz 3D CG for Pyramid (deformed)</description>
<executable>Helmholtz3D</executable>
<parameters>Helmholtz3D_Pyr_Deformed.xml</parameters>
<files>
<file description="Session File">Helmholtz3D_Pyr_Deformed.xml</file>
</files>
<metrics>
<metric type="L2" id="1">
<value tolerance="1e-12">2.76145e-05</value>
</metric>
<metric type="Linf" id="2">
<value tolerance="1e-12">0.000161114</value>
</metric>
</metrics>
</test>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="7" TYPE="MODIFIED" FIELDS="u" />
</EXPANSIONS>
<CONDITIONS>
<PARAMETERS>
<P> Lambda = 1 </P>
</PARAMETERS>
<VARIABLES>
<V ID="0"> u </V>
</VARIABLES>
<BOUNDARYREGIONS>
<B ID="0"> C[1-6] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="u" VALUE="sin(PI/2*x)*sin(PI/2*y)*sin(PI/2*z)" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="Forcing">
<E VAR="u" VALUE="-(Lambda+3*PI*PI/4)*sin(PI/2*x)*sin(PI/2*y)*sin(PI/2*z)" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="sin(PI/2*x)*sin(PI/2*y)*sin(PI/2*z)" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="3" SPACE="3">
<VERTEX>
<V ID="0">0.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="1">1.00000000e+00 0.00000000e+00 0.00000000e+00</V>
<V ID="2">1.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="3">0.00000000e+00 1.00000000e+00 0.00000000e+00</V>
<V ID="4">5.00000000e-01 5.00000000e-01 5.00000000e-01</V>
<V ID="5">0.00000000e+00 0.00000000e+00 1.00000000e+00</V>
<V ID="6">1.00000000e+00 0.00000000e+00 1.00000000e+00</V>
<V ID="7">1.00000000e+00 1.00000000e+00 1.00000000e+00</V>
<V ID="8">0.00000000e+00 1.00000000e+00 1.00000000e+00</V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 3 2 </E>
<E ID="3"> 0 3 </E>
<E ID="4"> 0 4 </E>
<E ID="5"> 1 4 </E>
<E ID="6"> 2 4 </E>
<E ID="7"> 3 4 </E>
<E ID="8"> 0 5 </E>
<E ID="9"> 5 6 </E>
<E ID="10"> 1 6 </E>
<E ID="11"> 5 4 </E>
<E ID="12"> 6 4 </E>
<E ID="13"> 6 7 </E>
<E ID="14"> 2 7 </E>
<E ID="15"> 7 4 </E>
<E ID="16"> 7 8 </E>
<E ID="17"> 3 8 </E>
<E ID="18"> 8 4 </E>
<E ID="19"> 8 5 </E>
</EDGE>
<FACE>
<Q ID="0"> 0 1 2 3</Q>
<T ID="1"> 0 5 4</T>
<T ID="2"> 1 6 5</T>
<T ID="3"> 2 6 7</T>
<T ID="4"> 3 7 4</T>
<Q ID="5"> 8 9 10 0</Q>
<T ID="6"> 8 11 4</T>
<T ID="7"> 9 12 11</T>
<T ID="8"> 10 12 5</T>
<Q ID="9"> 10 13 14 1</Q>
<T ID="10"> 13 15 12</T>
<T ID="11"> 14 15 6</T>
<Q ID="12"> 14 16 17 2</Q>
<T ID="13"> 16 18 15</T>
<T ID="14"> 17 18 7</T>
<Q ID="15"> 17 19 8 3</Q>
<T ID="16"> 19 11 18</T>
<Q ID="17"> 19 16 13 9</Q>
</FACE>
<ELEMENT>
<P ID="0"> 0 1 2 3 4 </P>
<P ID="1"> 5 6 7 8 1 </P>
<P ID="2"> 9 8 10 11 2 </P>
<P ID="3"> 12 11 13 14 3 </P>
<P ID="4"> 15 14 16 6 4 </P>
<P ID="5"> 17 16 13 10 7 </P>
</ELEMENT>
<COMPOSITE>
<C ID="0"> P[0-5] </C>
<C ID="1"> F[0] </C>
<C ID="2"> F[5] </C>
<C ID="3"> F[9] </C>
<C ID="4"> F[12] </C>
<C ID="5"> F[15] </C>
<C ID="6"> F[17] </C>
</COMPOSITE>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
<?xml version="1.0" encoding="utf-8" ?>
<NEKTAR>
<EXPANSIONS>
<E COMPOSITE="C[0]" NUMMODES="7" TYPE="MODIFIED" FIELDS="u" />
</EXPANSIONS>
<CONDITIONS>
<PARAMETERS>
<P> Lambda = 1 </P>
</PARAMETERS>
<VARIABLES>
<V ID="0"> u </V>
</VARIABLES>
<BOUNDARYREGIONS>
<B ID="0"> C[1] </B>
</BOUNDARYREGIONS>
<BOUNDARYCONDITIONS>
<REGION REF="0">
<D VAR="u" VALUE="sin(PI/2*x)*sin(PI/2*y)*sin(PI/2*z)" />
</REGION>
</BOUNDARYCONDITIONS>
<FUNCTION NAME="Forcing">
<E VAR="u" VALUE="-(Lambda+3*PI*PI/4)*sin(PI/2*x)*sin(PI/2*y)*sin(PI/2*z)" />
</FUNCTION>
<FUNCTION NAME="ExactSolution">
<E VAR="u" VALUE="sin(PI/2*x)*sin(PI/2*y)*sin(PI/2*z)" />
</FUNCTION>
</CONDITIONS>
<GEOMETRY DIM="3" SPACE="3">
<VERTEX>
<V ID="0"> 0.0 0.0 0.0 </V>
<V ID="1"> 1.0 0.0 0.0 </V>
<V ID="2"> 1.0 1.0 0.0 </V>
<V ID="3"> 0.0 1.0 0.0 </V>
<V ID="4"> 0.5 0.5 1.0 </V>
</VERTEX>
<EDGE>
<E ID="0"> 0 1 </E>
<E ID="1"> 1 2 </E>
<E ID="2"> 3 2 </E>
<E ID="3"> 0 3 </E>
<E ID="4"> 0 4 </E>
<E ID="5"> 1 4 </E>
<E ID="6"> 2 4 </E>
<E ID="7"> 3 4 </E>
</EDGE>
<FACE>
<Q ID="0"> 0 1 2 3 </Q>
<T ID="1"> 0 5 4 </T>
<T ID="2"> 1 6 5 </T>
<T ID="3"> 2 6 7 </T>
<T ID="4"> 3 7 4 </T>
</FACE>
<ELEMENT>
<P ID="0"> 0 1 2 3 4 </P>
</ELEMENT>
<COMPOSITE>
<C ID="0"> P[0] </C>
<C ID="1"> F[0-4] </C>
</COMPOSITE>
<CURVED>
<E ID="0" EDGEID="0" TYPE="PolyEvenlySpaced" NUMPOINTS="3"> 0 0 0 0.5 -0.1 -0.3 1 0 0 </E>
</CURVED>
<DOMAIN> C[0] </DOMAIN>
</GEOMETRY>
</NEKTAR>
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