<?xml version="1.0" encoding="UTF-8"?> <simulation xmds-version="2"> <name>lorenz_single</name> <author>Graham Dennis</author> <description> The Lorenz Attractor, an example of chaos. </description> <features> <auto_vectorise /> <benchmark /> <error_check /> <bing /> <precision>single</precision> <diagnostics /> <globals> <![CDATA[ real sigma = 10.0; real b = 8.0/3.0; real r = 28.0; ]]> </globals> </features> <geometry> <propagation_dimension> t </propagation_dimension> </geometry> <vector name="main" type="real"> <components> x y z </components> <initialisation> <![CDATA[ x = y = z = 1.0; ]]> </initialisation> </vector> <sequence> <integrate algorithm="ARK45" interval="100.0" steps="25000" tolerance="3e-5"> <samples>5000</samples> <operators> <integration_vectors>main</integration_vectors> <![CDATA[ dx_dt = sigma*(y-x); dy_dt = r*x - y - x*z; dz_dt = x*y - b*z; ]]> </operators> </integrate> </sequence> <output> <sampling_group initial_sample="yes"> <moments>xR yR zR</moments> <dependencies>main</dependencies> <![CDATA[ xR = x; yR = y; zR = z; ]]> </sampling_group> </output> </simulation>