StratifiedSampler.cpp

#include "StratifiedSampler.h"

namespace util {
StratifiedSampler::StratifiedSampler(const std::shared_ptr<Sampler>& sampler, size_t n)
    : sampler_(sampler)
    , n_(std::ceil(std::sqrt(n))) {};

Vec3 StratifiedSampler::color(float x, float y) const
{
    Vec3 color(0, 0, 0);
    for (int xi = 0; xi < n_; xi++) {
        for (int yi = 0; yi < n_; yi++) {
            double rx = dis(gen);
            double ry = dis(gen);
            double xs = x + (xi + rx) / n_;
            double ys = y + (yi + ry) / n_;
            color = color + sampler_->color(xs, ys);
        }
    }
    //if (color != Vec3 (16, 16, 16) && color != Vec3 (0, 0, 8))
    //std::cout << color << " and " << color / (n_ * n_) << n_ << std::endl;
    return color / (n_ * n_);
}
}