蒙特卡诺近似与PBM

介绍蒙特卡诺近似的例子代码

#include<fstream>
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<string>
#include<ctime>
#include<memory>
using std::ifstream;
using std::ofstream;
using std::string;
using std::cout;
using std::endl;
using std::ios;
using std::unique_ptr;
double randx()
{
    return rand() %  / (double);
}

int main(void)
{
    ifstream ifs("D:/tex.pbm");
    string header;
    uint32_t w, h, l;
    ifs >> header;
    ifs >> w >> h >> l;
    cout << "w: " << w << " ,h: " << h << " ,l: " << l << endl;
    ifs.ignore();
    unsigned char* pixles = new unsigned char[w * h * ];
    ifs.read((char*)pixles, w * h * );    

    int nsamples = ;
    srand(time());
    float avgr = , avgg = , avgb = ;
    float sumr = , sumg = , sumb = ;
    for (int i = ; i < nsamples; i++)
    {
        float x = randx() * w;
        float y = randx() * h;
        int n = ((int)(y * w) + int(x)) * ;
        sumr += pixles[n];
        sumg += pixles[n + ];
        sumb += pixles[n + ];
    }
    sumr /= nsamples;
    sumg /= nsamples;
    sumb /= nsamples;
    for (int i = ; i < h; i++)
    {
        for (int j = ; j < w; j++)
        {
            int n = (i * w + j) * ;
            avgr += pixles[n];
            avgg += pixles[n + ];
            avgb += pixles[n + ];
        }
    }
    avgr /= w * h;
    avgg /= w * h;
    avgb /= w * h;

    printf("Average: %0.2f %0.2f %0.2f \n", avgr, avgg, avgb);
    printf("Aproximation: %0.2f %0.2f %0.2f \n", sumr, sumg, sumb);
    delete pixles;

    system("pause");
    return ;
}