diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-03-21 21:58:39 +0100 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-03-21 21:59:52 +0100 |
commit | a7128aa0f728af4ce8c9db5d69eeafeec67d4cf4 (patch) | |
tree | 64773080c4751559203338b462019da8d27750b5 /Examples | |
parent | 74d42c41c16ac853edd4ec817348e9a478a5aeea (diff) |
Add version of WaveCaKernel in script language (untested)
Diffstat (limited to 'Examples')
-rw-r--r-- | Examples/WaveCaKernelScript.kranc | 24 |
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 |