diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2023-03-07 15:32:39 +0100 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2023-03-07 15:32:39 +0100 |
commit | b8fc822e94382883c8a772344f3d8f77686a6282 (patch) | |
tree | d97509f036b36ea1e2496de69a9e18bc3f9fcd6a /IsingScalingFunction.wl | |
parent | 0c9f4d522b557a0c9f1436e80bf666d9d2630592 (diff) | |
download | paper-b8fc822e94382883c8a772344f3d8f77686a6282.tar.gz paper-b8fc822e94382883c8a772344f3d8f77686a6282.tar.bz2 paper-b8fc822e94382883c8a772344f3d8f77686a6282.zip |
Fixes and cleaning of example file.
Diffstat (limited to 'IsingScalingFunction.wl')
-rw-r--r-- | IsingScalingFunction.wl | 11 |
1 files changed, 8 insertions, 3 deletions
diff --git a/IsingScalingFunction.wl b/IsingScalingFunction.wl index 8545567..4069a5d 100644 --- a/IsingScalingFunction.wl +++ b/IsingScalingFunction.wl @@ -25,6 +25,9 @@ DScriptFPlusMinusDξList::usage = DScriptF0DηList::usage = "DScriptF0DηList computes the first m derivatives of the scaling function F_0." +DScriptF0Dη::usage = + "DScriptF0Dη computes the mth derivative of the scaling function F_0." + DScriptMCasDξList::usage = "Computes the first m derivatives of the scaling function M given by Caselle et al." uf::usage = "uf computes the singular free energy u_f." @@ -45,6 +48,8 @@ PrepareArgument::usage = "Converts scaling function data into appropriate argume θ0Cas::usage = "" +h0Cas::usage = "" + gsCas::usage = "" Φs::usage = "" @@ -252,7 +257,7 @@ DScriptF0Dη[params__][m_, θ_] := Last[DScriptF0DηList[params][m, θ]] DufDut[θ0_, θYL_, B_, C0_, CYL_, Gs_, gs_][m_][R_, θ_] := m! RealAbs[uh[θ0, gs][R, θ]]^(2 / Δ - m / Δ) DScriptF0Dη[θ0, θYL, B, C0, CYL, Gs, gs][m, θ] + Log[uh[θ0, gs][R, θ]^2] / (8 π Δ) Derivative[m][Function[utp, utp^2]][ut[R, θ]] -DufDuh[θ0_, θYL_, B_, C0_, CYL_, Gs_, gs_][m_][R_, θ_] := m! RealAbs[ut[R, θ]]^(2-m Δ) DScriptFPlusMinusDξ[θ0, θYL, B, C0, CYL, Gs, gs][m, θ] + Derivative[m][Identity][θ] ut[R, θ]^2 / (8 π) Log[ut[R, θ]^2] +DufDuh[θ0_, θYL_, B_, C0_, CYL_, Gs_, gs_][m_][R_, θ_] := m! RealAbs[ut[R, θ]]^(2-m Δ) DScriptFPlusMinusDξ[θ0, θYL, B, C0, CYL, Gs, gs][m, θ] + Derivative[m][1&][θ] ut[R, θ]^2 / (8 π) Log[ut[R, θ]^2] ruleB[θ0_, gs_] := (2 * OverlineS / π) * (- g[θ0, gs]'[θ0] / t[θ0]^Δ) @@ -265,8 +270,8 @@ h0Cas := a b ρ /. { b -> (-c4/c2^3)^(1/2), ρ -> 2.00881 } /. { - c2->Ghs[[3]]2!, - c4->Ghs[[5]]["Value"]4! + c2->Ghs[[3]] 2!, + c4->Ghs[[5]]["Value"] 4! } gsCas := h0Cas * { |