diff options
| author | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-01-27 15:27:52 -0600 |
|---|---|---|
| committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2012-01-27 15:27:52 -0600 |
| commit | 5a381ef764da7d29fcaa88ab9b787d3cc2b01377 (patch) | |
| tree | d1a1ec027efc0cb7990dbc5402f21ebdc4aa9b91 /Examples | |
| parent | aaa5eb199c4cc9bc31dfde4bec69e29ab1b4fad7 (diff) | |
SimpleWaveCaKernel.m: Switch to Gaussian initial data
If this is close enough to zero at the boundaries, we don't need to worry about a boundary condition there.
Diffstat (limited to 'Examples')
| -rw-r--r-- | Examples/SimpleWaveCaKernel.m | 15 |
1 files changed, 14 insertions, 1 deletions
diff --git a/Examples/SimpleWaveCaKernel.m b/Examples/SimpleWaveCaKernel.m index f6908c7..3560710 100644 --- a/Examples/SimpleWaveCaKernel.m +++ b/Examples/SimpleWaveCaKernel.m @@ -21,6 +21,19 @@ initialSineCalc = } }; +f[x_] := Exp[-(x/0.1)^2]; + +initialGaussianCalc = +{ + Name -> "initial_gaussian", + Schedule -> {"AT INITIAL"}, + Equations -> + { + phi -> f[t+x], + pi -> D[f[t+x],t] + } +}; + evolveCalc = { Name -> "calc_rhs", @@ -48,7 +61,7 @@ integrateCalc = CreateKrancThornTT[groups, ".", "SimpleWaveCaKernel", - Calculations -> {initialSineCalc, evolveCalc, integrateCalc}, + Calculations -> {initialGaussianCalc, evolveCalc, integrateCalc}, PartialDerivatives -> derivatives, UseCaKernel -> True, EvolutionTimelevels -> 1, |