example_complex_fractal.cpp

 
#include <Pds/Ca>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Complex Z;
    Pds::Complex C;
    
    unsigned int L=512;
    unsigned int i;
    
    Pds::Vector X=Pds::LinSpace(-2,2,L);
    Pds::Vector Y=Pds::LinSpace(-2,2,L);
    
    Pds::Matrix W(L,L);         
    
    for(unsigned int lin=0;lin<L;lin++)
    for(unsigned int col=0;col<L;col++)
    {
        Z=0;
        C.Rect(X.Get(lin),Y.Get(col));
        for(i=0; (i<255)&&(std::isfinite(Z.Module())) ; i++)
        {
            Z=Z*Z+C;
        }
        if(std::isfinite(Z.Module()))   W.At(lin,col)=255;
        else                            W.At(lin,col)=i;
    }
    
    W.ExportBmpFile (Pds::Colormap::Jet,"mandelbrot.bmp");
    return 0;
}