src/rcrt/materials/WoodMaterial.cpp

Go to the documentation of this file.

#include "WoodMaterial.h"

namespace rcrt
{

WoodMaterial::WoodMaterial(const Point3D& org, const Vec3D& dir, const RGBColor& light, const RGBColor& dark):wood(0.03,1,0.4,org,dir,light,dark)
{
        
}

WoodMaterial::~WoodMaterial()
{
}


RGBColor WoodMaterial::sample(const Vec3D& wOut, const Vec3D& wInc,
                                Intersection& is)
{
        //std::cout << "sample wood material" << std::endl;
        const Point3D& p = is.getPosition();
        RGBColor color = wood.getColor(p).getRGB();
        //std::cout << "color: " << color << std::endl;
        return color/M_PI;
}

RGBColor WoodMaterial::sampleDiffuse(const Vec3D& wOut, const Vec3D& wInc, Intersection& is)
{
        return sample(wOut,wInc,is);
}

RGBColor WoodMaterial::sampleSpecular(const Vec3D& wOut, const Vec3D& wInc, Intersection& is)
{
        return 0;
}

RGBColor WoodMaterial::getEmitted(const Vec3D& wOut, Intersection& is)
{
        return 0;
}

bool WoodMaterial::hasDiffuse() const
{
        return true;
}

bool WoodMaterial::hasSpecular() const
{
        return false;
}

bool WoodMaterial::hasTransmissive() const
{
        return false;
}

bool WoodMaterial::refracts() const
{
        return false;
}

ScatterEvent WoodMaterial::scatterPhoton(Intersection& is, Photon* photon)
{
        Vec3D normal(getShadingNormal(is));
        const RGBColor Kd(wood.getColor(is.getPosition()).getRGB());
        //const RGBColor& P = photon->getPower();
        const float Pd = Kd.max();
        const float xi = float(rand()) / float(RAND_MAX); 
        if(xi <= Pd){
                photon->setPos(is.getPosition());
                //diffuse reflection
                photon->scalePower(Kd/Pd);
                //cout << " scaling " << Kd/Pd << endl;
                //get a random direction in the upper hemisphere
                float xi1 = float(rand()) / float(RAND_MAX);
                float xi2 = float(rand()) / float(RAND_MAX);
                float z2 = xi2 * 2.0f * M_PI;
                Vec3D u;
                Vec3D v;
                normal.getCS(u, v);
                Vec3D outDir(u * cos(z2) + v * sin(z2) * sqrt(1-xi1) + normal * sqrt(xi1));
                outDir.normalize();
                photon->setDir(outDir);
                return DIFFUSE;
        } else
                return STORE;
}

Vec3D WoodMaterial::getShadingNormal(Intersection& is) const
{
        return wood.getBump(is.getPosition(),is.getSNormalW()) * (is.backSide() ? -1 : 1);
}

}

---