aboutsummaryrefslogtreecommitdiff
path: root/Examples
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2012-03-21 21:58:39 +0100
committerIan Hinder <ian.hinder@aei.mpg.de>2012-03-21 21:59:52 +0100
commita7128aa0f728af4ce8c9db5d69eeafeec67d4cf4 (patch)
tree64773080c4751559203338b462019da8d27750b5 /Examples
parent74d42c41c16ac853edd4ec817348e9a478a5aeea (diff)
Add version of WaveCaKernel in script language (untested)
Diffstat (limited to 'Examples')
-rw-r--r--Examples/WaveCaKernelScript.kranc24
1 files changed, 24 insertions, 0 deletions
diff --git a/Examples/WaveCaKernelScript.kranc b/Examples/WaveCaKernelScript.kranc
new file mode 100644
index 0000000..12ef712
--- /dev/null
+++ b/Examples/WaveCaKernelScript.kranc
@@ -0,0 +1,24 @@
+begin thorn WaveCaKernelScript
+
+use cakernel
+
+begin variables
+ phi pi
+end variables
+
+begin calculation initial_gaussian_calc scheduled at initial
+ phi = exp(-1*((t+x)/0.1)**2)
+ pi = -2*(t+x)/0.1**2 * exp(-1*((t+x)/0.1)**2)
+end calculation
+
+begin calculation calc_rhs scheduled at MoL_CalcRHS
+ D_t phi = pi
+ D_t pi = Euc^ij*D_ij phi
+end calculation
+
+begin calculation calc_bound_rhs scheduled at MoL_RHSBoundaries
+ D_t phi = 2*(t+x)/0.1**2 * exp(-1*((t+x)/0.1)**2)
+ D_t pi = -2/(exp((t+x)**2/0.1**2) * 0.1**2) + (4*(t+x)**2)/(exp((t+x)**2/0.1**2) * 0.1**4)
+end calculation
+
+end thorn