Differences

This shows you the differences between two versions of the page.

--- playground:playground [2022/11/07 16:54] – johnsoevans
+++ playground:playground [2024/01/23 07:45] – johnsoevans
@@ Line 1: / Line 1: @@
-dummy page to create a page about some code to generate .inp code for a TLS refinement.
+====== PlayGround ======
-[[:TLS_refinement|TLS_refinement]]
+Play with any wiki page tests you like here.
-[[:Mathematica_scripts|Mathematica Scripts]]
-can i write the word print
-<code> print </code>
 <code topas>
-test
+/* this is a comment */
-max 0
+' this is a quotation comment
-hello
+   TOF_XYE....
-Sin
+    local !two_theta 145
-for
+    local !Lf 25
-print
+    local mu 0 min 0
-</code>
+    Abs_Lobanov
-<code>
+   TOF_XYE....
-#include <iostream>
+    local !two_theta 90
-using namespace std;
+    local !Lf 26
+    local mu 0 min 0
-int main() {
+    Abs_Lobanov
-  cout << "Hello World!";
-  return 0;
-  hello
-}
+    macro Abs_Lobanov
+    {
+    prm !Lf 60.183 /* POWGEN 90 degree flightpath */
+    prm !h_bar 6.626176e-34 /* Plank's constant */
+    prm !m_n 1.67495e-27 /* neutron mass */
+    prm !td_const = 10^4 h_bar/(m_n Lf);
+    local !sintheta = Sin(Deg_on_2 two_theta);
+    local !lambda = Xo 10^4 h_bar/(m_n Lf);
+    local !k0 = 1.697653;
+    local !k1 = (25.99978 - (0.01911 sintheta^0.5)Exp(-0.024514 sintheta^2)) + 0.109561 sintheta - 26.0456;
+    local !k2 = -0.02489 - (0.39499 sintheta^2) + (1.219077 sintheta^3) - (1.31268 sintheta^4) + (0.871081 sintheta^5) - (0.2327 sintheta^6);
+    local !k3 = 0.003045 + (0.018167 sintheta^2) - (0.03305 sintheta^4);
+    local !k4 = 1.433902 + (11.07504 sintheta^2) - (8.77629 sintheta^4) + (10.02088 sintheta^6) - (3.36778 sintheta^8);
+    local !k5 = ((0.013869 - 0.01249 sintheta^2)Exp(3.27094 sintheta^2)) + ((0.337894 + (13.77317 sintheta^2))/((1 + 11.53544 sintheta^2)^1.555039));
+    local !k6 = (1.933433 / (1 + 23.12967 sintheta^2)^1.686715) - 0.13576 sintheta + 1.163198;
+    local !k7 = 0.044365 - (0.04259/((1 + 0.41051 sintheta^2)^148.4202));
+    scale_pks = If( (mu*lambda) < 3, Exp(-k0 (mu*lambda) - k1 (mu*lambda)^2 - k2 (mu*lambda)^3 - k3 (mu*lambda)^4),
+				Exp((k4 - k7/(1 + k5 ((mu*lambda)-3))^k6)+k7)
+				);
+    }
 </code>
-  hello
-  this is indented
-  print
-  Print