Function test_velocity_and_steps

tests/unittests/test_diffusion_math.cpp:61–82

}

void test_velocity_and_steps() {
    std::vector<float> velocity{3.0F, 4.0F};
    engine::sampling::clamp_velocity_norm(velocity, {6.0F, 8.0F}, 0.5F);
    require_vector_close(velocity, {3.0F, 4.0F}, 1.0e-5F, "velocity clamp at threshold");

    std::vector<float> scaled{6.0F, 8.0F};
    engine::sampling::clamp_velocity_norm(scaled, {3.0F, 4.0F}, 0.5F);
    require_vector_close(scaled, {1.5F, 2.0F}, 1.0e-5F, "velocity clamp scales");

    std::vector<float> zero{0.0F, 0.0F};
    engine::sampling::clamp_velocity_norm(zero, {0.0F, 0.0F}, 0.5F);
    require_vector_close(zero, {0.0F, 0.0F}, 0.0F, "velocity clamp zero");

    const std::vector<float> x{1.0F, 2.0F, 3.0F};
    const std::vector<float> v{0.5F, -1.0F, 4.0F};
    const std::vector<float> v2{1.0F, 1.0F, 0.0F};
    require_vector_close(engine::sampling::euler_step(x, v, 0.25F), {0.875F, 2.25F, 2.0F}, 0.0F, "euler");
    require_vector_close(engine::sampling::denoise_from_velocity(x, v, 0.75F), {0.625F, 2.75F, 0.0F}, 0.0F, "denoise");
    require_vector_close(engine::sampling::renoise(x, {-1.0F, 0.0F, 2.0F}, 0.2F), {0.6F, 1.6F, 2.8F}, 2.0e-6F, "renoise");
    require_vector_close(engine::sampling::heun_combine_velocity(v, v2), {0.75F, 0.0F, 2.0F}, 0.0F, "heun combine");
    require_vector_close(engine::sampling::heun_step(x, v, v2, 0.25F), {0.8125F, 2.0F, 2.5F}, 0.0F, "heun step");
}

void test_masks() {
    require_vector_close(