~fluidity-core/fluidity/adjoint

<?xml version='1.0' encoding='utf-8'?>

<testproblem>

  <name>Shallow Water Equation MMS Optimisation</name>

  <owner userid="sf1409"/>

  <tags>swml adjoint optimisation revolve</tags>

  <problem_definition length="medium" nprocs="1">

    <command_line>time ../../bin/optimality optimiser_A.oml; time ../../bin/optimality optimiser_B.oml; time ../../bin/optimality optimiser_C.oml; time ../../bin/optimality optimiser_D.oml;</command_line>

  </problem_definition>

  <variables>

    <variable name="u_error_initial" language="python">from fluidity_tools import stat_parser as stat

import numpy      

import math

sims = ["_A", "_B", "_C", "_D"]

u_error_initial = numpy.zeros(len(sims))

for i, sim in enumerate(sims):

  u_error_initial[i] = stat("wave%s.stat" % sim)["Fluid"]["VelocityError%magnitude"]["max"][0]

</variable>

    <variable name="u_error_order_final" language="python">from fluidity_tools import stat_parser as stat

import numpy      

import math

sims = ["_A", "_B", "_C", "_D"]

u_error = numpy.zeros(len(sims))

for i, sim in enumerate(sims):

  u_error[i] = stat("wave%s.stat" % sim)["Fluid"]["VelocityError%magnitude"]["max"][-1]

 

u_error_order_final = numpy.zeros(len(sims)-1) 

for i in range(len(sims)-1):

  u_error_order_final[i] = math.log(u_error[i]/u_error[i+1], 2)

</variable>

    <variable name="eta_error_order_initial" language="python">from fluidity_tools import stat_parser as stat

import numpy      

import math

sims = ["_A", "_B", "_C", "_D"]

eta_error = numpy.zeros(len(sims))

for i, sim in enumerate(sims):

  eta_error[i] = stat("wave%s.stat" % sim)["Fluid"]["LayerThicknessError"]["max"][0]

 

eta_error_order_initial = numpy.zeros(len(sims)-1)

for i in range(len(sims)-1):

  eta_error_order_initial[i] = math.log(eta_error[i]/eta_error[i+1], 2)

</variable>

    <variable name="eta_error_order_final" language="python">from fluidity_tools import stat_parser as stat

import numpy      

import math

sims = ["_A", "_B", "_C", "_D"]

eta_error = numpy.zeros(len(sims))

for i, sim in enumerate(sims):

  eta_error[i] = stat("wave%s.stat" % sim)["Fluid"]["LayerThicknessError"]["max"][-1]

 

eta_error_order_final = numpy.zeros(len(sims)-1)

for i in range(len(sims)-1):

  eta_error_order_final[i] = math.log(eta_error[i]/eta_error[i+1], 2)

</variable>

  </variables>

  <pass_tests>

    <test name="u_error_initial_test" language="python">assert min(u_error_initial) == 0.0</test>

    <test name="u_error_order_final_test" language="python">assert min(u_error_order_final) &gt; 1.8</test>

    <test name="eta_error_order_initial_test" language="python">assert min(eta_error_order_initial) &gt; 1.2</test>

    <test name="eta_error_order_final_test" language="python">assert min(eta_error_order_final) &gt; 1.9</test>

  </pass_tests>

  <warn_tests/>

</testproblem>