RK4 position update weights all stages equally · case 01

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x0, v0, k, c, F, wf, dt, steps):
    def f(t, x, v):
        return v, -k * x - c * v + F * math.cos(wf * t)
    t = 0.0
    x, v = x0, v0
    for _ in range(steps):
        k1x, k1v = f(t, x, v)
        k2x, k2v = f(t + dt / 2, x + dt / 2 * k1x, v + dt / 2 * k1v)
        k3x, k3v = f(t + dt / 2, x + dt / 2 * k2x, v + dt / 2 * k2v)
        k4x, k4v = f(t + dt, x + dt * k3x, v + dt * k3v)
        x = x + dt / 4 * (k1x + k2x + k3x + k4x)
        v = v + dt / 6 * (k1v + 2 * k2v + 2 * k3v + k4v)
        t = t + dt
    return [round(x, 6), round(v, 6), round(t, 6)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[[[0.594, 0.724, 5.574, 1.414, 0.54, 2.544, 0.2, 5], [0.050911, -1.000788, 1.0]], [[0.488, -0.59, 18.046, 0.065, 1.27, 1.284, 0.2, 3], [-0.348659, -0.553321, 0.6]], [[-0.109, -0.682, 16.643, 0.151, 1.208, 3.199, 0.05, 3], [-0.171504, -0.13078, 0.15]], [[0.043, -0.274, 12.514, 0.941, 1.718, 1.943, 0.1, 6], [0.088471, 0.195664, 0.6]], [[0.477, 0.804, 9.463, 0.795, 0.343, 4.866, 0.2, 1], [0.539279, -0.178827, 0.2]], [[0.882, 0.468, 0.596, 1.065, 0.599, 2.931, 0.2, 6], [0.969734, -0.430672, 1.2]], [[-0.156, 0.405, 4.019, 0.13, 0.809, 3.428, 0.1, 5], [0.148937, 0.625786, 0.5]], [[-0.435, 0.18, 4.535, 0.378, 1.241, 0.995, 0.2, 6], [0.675696, 0.34121, 1.2]]], [[[-0.926, 0.952, 5.429, 0.55, 1.365, 4.948, 0.2, 1], [-0.633361, 1.884933, 0.2]], [[-0.032, 0.808, 18.227, 0.749, 0.437, 3.703, 0.1, 1], [0.04836, 0.775101, 0.1]], [[-0.314, -0.856, 15.027, 0.039, 0.971, 3.038, 0.1, 4], [-0.172381, 1.364757, 0.4]], [[-0.876, -0.199, 5.264, 1.176, 0.534, 2.649, 0.1, 3], [-0.724871, 1.082023, 0.3]], [[-0.827, 0.717, 1.853, 0.324, 0.414, 2.766, 0.05, 3], [-0.701343, 0.949832, 0.15]], [[-0.541, -0.377, 16.521, 1.005, 0.304, 1.061, 0.05, 7], [-0.191246, 1.884299, 0.35]], [[0.35, 0.024, 2.459, 0.489, 1.678, 3.184, 0.2, 2], [0.401172, 0.137673, 0.4]], [[-0.492, -0.144, 14.57, 1.424, 1.521, 1.605, 0.05, 3], [-0.421761, 1.012157, 0.15]]], [[[-0.571, 0.019, 10.668, 1.333, 0.813, 4.88, 0.1, 10], [0.238859, -0.217571, 1.0]], [[-0.346, -0.2, 8.452, 1.354, 0.585, 4.532, 0.1, 9], [0.091625, 0.239009, 0.9]], [[0.704, 0.002, 17.554, 0.361, 1.232, 1.599, 0.05, 4], [0.499246, -1.906948, 0.2]], [[0.722, 0.536, 7.744, 0.165, 1.679, 3.18, 0.05, 5], [0.724479, -0.530907, 0.25]], [[0.389, 0.616, 7.403, 1.127, 1.629, 2.121, 0.2, 1], [0.470821, 0.193297, 0.2]], [[0.156, -0.808, 12.089, 0.635, 0.926, 3.013, 0.05, 7], [-0.092465, -0.474045, 0.35]], [[-0.314, 0.869, 7.658, 1.211, 1.937, 0.965, 0.05, 3], [-0.152866, 1.236835, 0.15]], [[-0.394, 0.849, 18.783, 1.27, 1.837, 1.788, 0.1, 10], [0.047492, -1.406422, 1.0]]], [[[0.612, -0.071, 18.951, 0.15, 1.941, 3.737, 0.05, 10], [-0.235385, -2.074863, 0.5]], [[0.28, 0.052, 4.066, 1.361, 0.882, 2.789, 0.2, 1], [0.283811, -0.015909, 0.2]], [[0.561, -0.998, 18.692, 1.273, 1.093, 2.216, 0.2, 6], [0.099116, 0.552827, 1.2]], [[-0.167, -0.626, 6.981, 0.788, 1.493, 3.024, 0.2, 4], [0.018342, 0.36187, 0.8]], [[-0.772, 0.752, 2.823, 0.497, 1.621, 2.92, 0.1, 6], [0.065151, 1.550767, 0.6]], [[-0.198, -0.028, 4.048, 0.079, 1.642, 1.588, 0.2, 8], [0.361524, -1.100738, 1.6]], [[0.962, 0.482, 3.289, 0.436, 1.4, 4.793, 0.05, 12], [0.785086, -1.262878, 0.6]], [[-0.484, -0.908, 18.372, 0.481, 1.964, 2.315, 0.2, 2], [-0.051507, 2.343495, 0.4]]], [[[-0.341, 0.921, 14.484, 0.78, 0.993, 4.954, 0.1, 11], [-0.092121, -0.961694, 1.1]], [[0.0, 0.489, 9.543, 0.742, 1.537, 3.108, 0.1, 13], [-0.277379, -0.761947, 1.3]], [[0.163, 0.02, 4.507, 0.115, 0.982, 3.557, 0.1, 3], [0.175188, 0.032334, 0.3]], [[0.96, 0.368, 6.136, 0.868, 0.713, 4.277, 0.1, 13], [-0.632372, -0.226167, 1.3]], [[0.48, 0.991, 13.504, 0.938, 1.754, 1.99, 0.2, 8], [-0.112907, 0.598276, 1.6]], [[0.949, -0.613, 14.001, 1.335, 0.724, 4.136, 0.1, 10], [-0.470021, 0.99105, 1.0]], [[0.207, 0.017, 8.404, 1.117, 1.656, 2.196, 0.05, 4], [0.207866, -0.012779, 0.2]], [[0.476, 0.606, 16.115, 1.076, 0.293, 2.306, 0.2, 1], [0.443713, -0.858349, 0.2]]]]
for i, (args, expected) in enumerate(fixtures[N-1]):
    check("case %d" % i, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)

SHA-256 / fdf5fd6fccc135fb039aba31af8a0239b905f75a9de5d488a17f60882f2046ce

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x0, v0, k, c, F, wf, dt, steps):
    def f(t, x, v):
        return v, -k * x - c * v + F * math.cos(wf * t)
    t = 0.0
    x, v = x0, v0
    for _ in range(steps):
        k1x, k1v = f(t, x, v)
        k2x, k2v = f(t + dt / 2, x + dt / 2 * k1x, v + dt / 2 * k1v)
        k3x, k3v = f(t + dt / 2, x + dt / 2 * k2x, v + dt / 2 * k2v)
        k4x, k4v = f(t + dt, x + dt * k3x, v + dt * k3v)
        x = x + dt / 6 * (2 * k1x + k2x + k3x + 2 * k4x)
        v = v + dt / 6 * (k1v + 2 * k2v + 2 * k3v + k4v)
        t = t + dt
    return [round(x, 6), round(v, 6), round(t, 6)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[[[0.594, 0.724, 5.574, 1.414, 0.54, 2.544, 0.2, 5], [0.050911, -1.000788, 1.0]], [[0.488, -0.59, 18.046, 0.065, 1.27, 1.284, 0.2, 3], [-0.348659, -0.553321, 0.6]], [[-0.109, -0.682, 16.643, 0.151, 1.208, 3.199, 0.05, 3], [-0.171504, -0.13078, 0.15]], [[0.043, -0.274, 12.514, 0.941, 1.718, 1.943, 0.1, 6], [0.088471, 0.195664, 0.6]], [[0.477, 0.804, 9.463, 0.795, 0.343, 4.866, 0.2, 1], [0.539279, -0.178827, 0.2]], [[0.882, 0.468, 0.596, 1.065, 0.599, 2.931, 0.2, 6], [0.969734, -0.430672, 1.2]], [[-0.156, 0.405, 4.019, 0.13, 0.809, 3.428, 0.1, 5], [0.148937, 0.625786, 0.5]], [[-0.435, 0.18, 4.535, 0.378, 1.241, 0.995, 0.2, 6], [0.675696, 0.34121, 1.2]]], [[[-0.926, 0.952, 5.429, 0.55, 1.365, 4.948, 0.2, 1], [-0.633361, 1.884933, 0.2]], [[-0.032, 0.808, 18.227, 0.749, 0.437, 3.703, 0.1, 1], [0.04836, 0.775101, 0.1]], [[-0.314, -0.856, 15.027, 0.039, 0.971, 3.038, 0.1, 4], [-0.172381, 1.364757, 0.4]], [[-0.876, -0.199, 5.264, 1.176, 0.534, 2.649, 0.1, 3], [-0.724871, 1.082023, 0.3]], [[-0.827, 0.717, 1.853, 0.324, 0.414, 2.766, 0.05, 3], [-0.701343, 0.949832, 0.15]], [[-0.541, -0.377, 16.521, 1.005, 0.304, 1.061, 0.05, 7], [-0.191246, 1.884299, 0.35]], [[0.35, 0.024, 2.459, 0.489, 1.678, 3.184, 0.2, 2], [0.401172, 0.137673, 0.4]], [[-0.492, -0.144, 14.57, 1.424, 1.521, 1.605, 0.05, 3], [-0.421761, 1.012157, 0.15]]], [[[-0.571, 0.019, 10.668, 1.333, 0.813, 4.88, 0.1, 10], [0.238859, -0.217571, 1.0]], [[-0.346, -0.2, 8.452, 1.354, 0.585, 4.532, 0.1, 9], [0.091625, 0.239009, 0.9]], [[0.704, 0.002, 17.554, 0.361, 1.232, 1.599, 0.05, 4], [0.499246, -1.906948, 0.2]], [[0.722, 0.536, 7.744, 0.165, 1.679, 3.18, 0.05, 5], [0.724479, -0.530907, 0.25]], [[0.389, 0.616, 7.403, 1.127, 1.629, 2.121, 0.2, 1], [0.470821, 0.193297, 0.2]], [[0.156, -0.808, 12.089, 0.635, 0.926, 3.013, 0.05, 7], [-0.092465, -0.474045, 0.35]], [[-0.314, 0.869, 7.658, 1.211, 1.937, 0.965, 0.05, 3], [-0.152866, 1.236835, 0.15]], [[-0.394, 0.849, 18.783, 1.27, 1.837, 1.788, 0.1, 10], [0.047492, -1.406422, 1.0]]], [[[0.612, -0.071, 18.951, 0.15, 1.941, 3.737, 0.05, 10], [-0.235385, -2.074863, 0.5]], [[0.28, 0.052, 4.066, 1.361, 0.882, 2.789, 0.2, 1], [0.283811, -0.015909, 0.2]], [[0.561, -0.998, 18.692, 1.273, 1.093, 2.216, 0.2, 6], [0.099116, 0.552827, 1.2]], [[-0.167, -0.626, 6.981, 0.788, 1.493, 3.024, 0.2, 4], [0.018342, 0.36187, 0.8]], [[-0.772, 0.752, 2.823, 0.497, 1.621, 2.92, 0.1, 6], [0.065151, 1.550767, 0.6]], [[-0.198, -0.028, 4.048, 0.079, 1.642, 1.588, 0.2, 8], [0.361524, -1.100738, 1.6]], [[0.962, 0.482, 3.289, 0.436, 1.4, 4.793, 0.05, 12], [0.785086, -1.262878, 0.6]], [[-0.484, -0.908, 18.372, 0.481, 1.964, 2.315, 0.2, 2], [-0.051507, 2.343495, 0.4]]], [[[-0.341, 0.921, 14.484, 0.78, 0.993, 4.954, 0.1, 11], [-0.092121, -0.961694, 1.1]], [[0.0, 0.489, 9.543, 0.742, 1.537, 3.108, 0.1, 13], [-0.277379, -0.761947, 1.3]], [[0.163, 0.02, 4.507, 0.115, 0.982, 3.557, 0.1, 3], [0.175188, 0.032334, 0.3]], [[0.96, 0.368, 6.136, 0.868, 0.713, 4.277, 0.1, 13], [-0.632372, -0.226167, 1.3]], [[0.48, 0.991, 13.504, 0.938, 1.754, 1.99, 0.2, 8], [-0.112907, 0.598276, 1.6]], [[0.949, -0.613, 14.001, 1.335, 0.724, 4.136, 0.1, 10], [-0.470021, 0.99105, 1.0]], [[0.207, 0.017, 8.404, 1.117, 1.656, 2.196, 0.05, 4], [0.207866, -0.012779, 0.2]], [[0.476, 0.606, 16.115, 1.076, 0.293, 2.306, 0.2, 1], [0.443713, -0.858349, 0.2]]]]
for i, (args, expected) in enumerate(fixtures[N-1]):
    check("case %d" % i, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)

SHA-256 / 56f9b4b54b4f06dccb32a2d8cd26bdad2693b58b2b87ecdced797bf5e38ff956

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x0, v0, k, c, F, wf, dt, steps):
    def f(t, x, v):
        return v, -k * x - c * v + F * math.cos(wf * t)
    t = 0.0
    x, v = x0, v0
    for _ in range(steps):
        k1x, k1v = f(t, x, v)
        k2x, k2v = f(t + dt / 2, x + dt / 2 * k1x, v + dt / 2 * k1v)
        k3x, k3v = f(t + dt / 2, x + dt / 2 * k2x, v + dt / 2 * k2v)
        k4x, k4v = f(t + dt, x + dt * k3x, v + dt * k3v)
        x = x + dt / 6 * (k1x + 2 * k2x + 2 * k3x + k4x)
        v = v + dt / 6 * (k1v + 2 * k2v + 2 * k3v + k4v)
        t = t + dt
    return [round(x, 6), round(v, 6), round(t, 6)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[[[0.594, 0.724, 5.574, 1.414, 0.54, 2.544, 0.2, 5], [0.050911, -1.000788, 1.0]], [[0.488, -0.59, 18.046, 0.065, 1.27, 1.284, 0.2, 3], [-0.348659, -0.553321, 0.6]], [[-0.109, -0.682, 16.643, 0.151, 1.208, 3.199, 0.05, 3], [-0.171504, -0.13078, 0.15]], [[0.043, -0.274, 12.514, 0.941, 1.718, 1.943, 0.1, 6], [0.088471, 0.195664, 0.6]], [[0.477, 0.804, 9.463, 0.795, 0.343, 4.866, 0.2, 1], [0.539279, -0.178827, 0.2]], [[0.882, 0.468, 0.596, 1.065, 0.599, 2.931, 0.2, 6], [0.969734, -0.430672, 1.2]], [[-0.156, 0.405, 4.019, 0.13, 0.809, 3.428, 0.1, 5], [0.148937, 0.625786, 0.5]], [[-0.435, 0.18, 4.535, 0.378, 1.241, 0.995, 0.2, 6], [0.675696, 0.34121, 1.2]]], [[[-0.926, 0.952, 5.429, 0.55, 1.365, 4.948, 0.2, 1], [-0.633361, 1.884933, 0.2]], [[-0.032, 0.808, 18.227, 0.749, 0.437, 3.703, 0.1, 1], [0.04836, 0.775101, 0.1]], [[-0.314, -0.856, 15.027, 0.039, 0.971, 3.038, 0.1, 4], [-0.172381, 1.364757, 0.4]], [[-0.876, -0.199, 5.264, 1.176, 0.534, 2.649, 0.1, 3], [-0.724871, 1.082023, 0.3]], [[-0.827, 0.717, 1.853, 0.324, 0.414, 2.766, 0.05, 3], [-0.701343, 0.949832, 0.15]], [[-0.541, -0.377, 16.521, 1.005, 0.304, 1.061, 0.05, 7], [-0.191246, 1.884299, 0.35]], [[0.35, 0.024, 2.459, 0.489, 1.678, 3.184, 0.2, 2], [0.401172, 0.137673, 0.4]], [[-0.492, -0.144, 14.57, 1.424, 1.521, 1.605, 0.05, 3], [-0.421761, 1.012157, 0.15]]], [[[-0.571, 0.019, 10.668, 1.333, 0.813, 4.88, 0.1, 10], [0.238859, -0.217571, 1.0]], [[-0.346, -0.2, 8.452, 1.354, 0.585, 4.532, 0.1, 9], [0.091625, 0.239009, 0.9]], [[0.704, 0.002, 17.554, 0.361, 1.232, 1.599, 0.05, 4], [0.499246, -1.906948, 0.2]], [[0.722, 0.536, 7.744, 0.165, 1.679, 3.18, 0.05, 5], [0.724479, -0.530907, 0.25]], [[0.389, 0.616, 7.403, 1.127, 1.629, 2.121, 0.2, 1], [0.470821, 0.193297, 0.2]], [[0.156, -0.808, 12.089, 0.635, 0.926, 3.013, 0.05, 7], [-0.092465, -0.474045, 0.35]], [[-0.314, 0.869, 7.658, 1.211, 1.937, 0.965, 0.05, 3], [-0.152866, 1.236835, 0.15]], [[-0.394, 0.849, 18.783, 1.27, 1.837, 1.788, 0.1, 10], [0.047492, -1.406422, 1.0]]], [[[0.612, -0.071, 18.951, 0.15, 1.941, 3.737, 0.05, 10], [-0.235385, -2.074863, 0.5]], [[0.28, 0.052, 4.066, 1.361, 0.882, 2.789, 0.2, 1], [0.283811, -0.015909, 0.2]], [[0.561, -0.998, 18.692, 1.273, 1.093, 2.216, 0.2, 6], [0.099116, 0.552827, 1.2]], [[-0.167, -0.626, 6.981, 0.788, 1.493, 3.024, 0.2, 4], [0.018342, 0.36187, 0.8]], [[-0.772, 0.752, 2.823, 0.497, 1.621, 2.92, 0.1, 6], [0.065151, 1.550767, 0.6]], [[-0.198, -0.028, 4.048, 0.079, 1.642, 1.588, 0.2, 8], [0.361524, -1.100738, 1.6]], [[0.962, 0.482, 3.289, 0.436, 1.4, 4.793, 0.05, 12], [0.785086, -1.262878, 0.6]], [[-0.484, -0.908, 18.372, 0.481, 1.964, 2.315, 0.2, 2], [-0.051507, 2.343495, 0.4]]], [[[-0.341, 0.921, 14.484, 0.78, 0.993, 4.954, 0.1, 11], [-0.092121, -0.961694, 1.1]], [[0.0, 0.489, 9.543, 0.742, 1.537, 3.108, 0.1, 13], [-0.277379, -0.761947, 1.3]], [[0.163, 0.02, 4.507, 0.115, 0.982, 3.557, 0.1, 3], [0.175188, 0.032334, 0.3]], [[0.96, 0.368, 6.136, 0.868, 0.713, 4.277, 0.1, 13], [-0.632372, -0.226167, 1.3]], [[0.48, 0.991, 13.504, 0.938, 1.754, 1.99, 0.2, 8], [-0.112907, 0.598276, 1.6]], [[0.949, -0.613, 14.001, 1.335, 0.724, 4.136, 0.1, 10], [-0.470021, 0.99105, 1.0]], [[0.207, 0.017, 8.404, 1.117, 1.656, 2.196, 0.05, 4], [0.207866, -0.012779, 0.2]], [[0.476, 0.606, 16.115, 1.076, 0.293, 2.306, 0.2, 1], [0.443713, -0.858349, 0.2]]]]
for i, (args, expected) in enumerate(fixtures[N-1]):
    check("case %d" % i, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)

SHA-256 / 08e6d5e26821c2314ee75beda264d09368c008e926e549df08481120ac5974cf