source: GTP/trunk/App/Demos/Illum/PathMap/PathMapEffect.cpp @ 896

#include "dxstdafx.h"
#include ".\pathmapeffect.h"
#include "Material.hpp"
#include "TexturedMaterial.h"
#include "WoodMaterial.hpp"
#include "TriangleMesh.h"
#include "Transformed.h"
#include "shlwapi.h"
#include <queue>
#include <fstream>

#define TEST_CLUST 0
//#define GENERATE_PATH_MAPS

//! D3D vector-matrix multiplication
D3DVECTOR operator *(const D3DVECTOR& v, const D3DMATRIX& m)
{
        D3DVECTOR r;
        r.x = v.x * m._11 + v.y * m._21 + v.z * m._31 + m._41;
        r.y = v.x * m._12 + v.y * m._22 + v.z * m._32 + m._42;
        r.z = v.x * m._13 + v.y * m._23 + v.z * m._33 + m._43;
        float h = v.x * m._14 + v.y * m._24 + v.z * m._34 + m._44;
        if(h > 0.00001f || h < 0.00001f)
        {
                r.x /= h;
                r.y /= h;
                r.z /= h;
        }
        return r;
}

//! method name strings to be displayed on screen
const wchar_t* PathMapEffect::Method::methodNames[10] =
{
        L"PRM",
        L"PRM texture",
        L"_",
        L"_"
        , NULL, NULL, NULL, NULL, NULL, NULL
};

//! method constants
const PathMapEffect::Method PathMapEffect::Method::PRM(0);
const PathMapEffect::Method PathMapEffect::Method::SHOWTEX(1);
const PathMapEffect::Method PathMapEffect::Method::LAST(2);

//! constructor
PathMapEffect::PathMapEffect(LPDIRECT3DDEVICE9 device)
:method(Method::PRM)
{
        rayId = 1;
        this->device = device;

        //first person camera allows more freedom of movement
        camera = new CFirstPersonCamera();
        lightCamera = new CFirstPersonCamera();

        HRESULT hr;
        DWORD effectCompileFlag=0;
        
        LPD3DXBUFFER compilationErrors;
        if(FAILED(
                hr = 
                        D3DXCreateEffectFromFile(
                                device,
                                L"pathMap.fx",
                                NULL,
                                NULL,
                                0,
                                NULL,
                                &effect,
                                &compilationErrors) )){
                MessageBoxA( NULL, (LPSTR)compilationErrors->GetBufferPointer(), "Failed to load effect file!", MB_OK);
                exit(-1);
        }

        //store buffers so we can reset them after rendering to texture render targets
        device->GetRenderTarget(0, &frameColorBuffer);
        device->GetDepthStencilSurface(&frameDepthStencilBuffer);
        
        //load empty texture
        D3DXCreateTextureFromFile(device, L"media\\empty.bmp", &emptyTexture);

        //set up a scene

        //load hangar mesh
        loadMesh( L"media\\fulltexchamber.x" );

        //create hangar entity
        Entity e1;
        e1.owner = this;
        e1.renderMesh = renderMeshes.at(0);
        D3DXMatrixIdentity(&e1.modelWorldTransform);
//      D3DXMatrixScaling(&e1.modelWorldTransform, 0.5f, 1.0f, 0.5f);
        D3DXMatrixIdentity(&e1.inverseTransposedModelWorldTransform);
//      D3DXMatrixScaling(&e1.inverseTransposedModelWorldTransform, 1.0f/0.5f, 1.0f, 1.0f/0.5f);

        e1.createRayTraceEntity();
        
        entities.push_back(e1);


        //load 'steps' mesh
/*      loadMesh( L"media\\steps.x" );

        //create step entities
        //for(int eitor=0; eitor < 6; eitor++)
        {
                Entity e2;
                e2.owner = this;
                e2.renderMesh = renderMeshes.at(1);
                D3DXMatrixTranslation(&e2.modelWorldTransform, 0.0f, -14.0f, 0.0f);
                D3DXMatrixInverse(&e2.inverseTransposedModelWorldTransform, NULL, &e2.modelWorldTransform);
                e2.createRayTraceEntity();
                entities.push_back(e2);

                D3DXMATRIX trasi, roti;
                D3DXMatrixTranslation(&trasi, 7.0f, 14.0f, -7.0f);//7.25f, 13.5, -7.25); //rot 90 -90 90 trans 7.25 13.5 -7.25
                D3DXMatrixRotationYawPitchRoll(&roti, -Vector::PI * 0.5, Vector::PI, 0.0);
//              D3DXMatrixIdentity(&roti);
                D3DXMatrixMultiply(&e2.modelWorldTransform, &roti, &trasi);
                D3DXMatrixInverse(&e2.inverseTransposedModelWorldTransform, NULL, &e2.modelWorldTransform);
                e2.createRayTraceEntity();
                entities.push_back(e2);
        }*/

/*      for(int eitor=0; eitor < 6; eitor++)
        {
                Entity e2;
                e2.owner = this;
                e2.renderMesh = renderMeshes.at(1);
                D3DXMatrixIdentity(&e2.modelWorldTransform);
                D3DXMatrixScaling(&e2.modelWorldTransform, -0.5, 0.5, -0.5);
                D3DXMATRIX trmx;
                D3DXMatrixTranslation(&trmx, -8.0 + eitor * 2.8f, -3.2f, -8.6f + eitor * 0.6f);
                e2.modelWorldTransform *= trmx;
                D3DXMatrixIdentity(&e2.inverseTransposedModelWorldTransform);
                D3DXMatrixScaling(&e2.inverseTransposedModelWorldTransform, -2.0, 2.0, -2.0);
                D3DXMatrixTranslation(&trmx, +8.0 - eitor * 2.8f, 3.2f, 8.6f - eitor * 0.6f);
                e2.inverseTransposedModelWorldTransform *= trmx;
                e2.createRayTraceEntity();
                entities.push_back(e2);
        }*/

        //compute the surface area of the complete geometry. useful for random sampling.
        sumSurfaceArea = 0.0;
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                sumSurfaceArea += entityIterator->rayTraceEntity->getSurfaceArea();
                entityIterator++;
        }

        //initialize camera
        camera->SetViewParams( &D3DXVECTOR3(0.0f, 0.0f, 4.0f), &D3DXVECTOR3(0.0f, 0.0f, 0.0f));
        D3DSURFACE_DESC fbdesc;
        frameColorBuffer->GetDesc(&fbdesc);
        camera->SetProjParams( D3DX_PI/2, (float)fbdesc.Width / fbdesc.Height, 0.1f, 300.0f );

        //create resources for entities
        createPRMTextures();

        //set up spotlight
        lightCamera->SetViewParams( &D3DXVECTOR3(5.0f, 0.0f, 0.0f), &D3DXVECTOR3(0.0f, 0.0f, 0.0f));
        lightCamera->SetProjParams( D3DX_PI/1.9, 1.0f, 0.1f, 300.0f );

        //create ray tracing kd-tree containg ray-traceable versions of entities 
        Intersectable** objs = new Intersectable*[entities.size()];
        for(int u=0; u<entities.size(); u++)
                objs[u] = entities.at(u).rayTraceEntity;
        kdtree = new KDTree(objs, entities.size());
        delete [] objs;

        //create dummy render target texture for depth map rendering
        device->CreateTexture(DEPTHMAPRES, DEPTHMAPRES, 1, D3DUSAGE_RENDERTARGET, D3DFMT_R32F, D3DPOOL_DEFAULT, &fakeTexture, NULL);
        fakeTexture->GetSurfaceLevel(0, &fakeSurface);

        //create a depthstencil texture for depth map rendering
        device->CreateTexture(DEPTHMAPRES, DEPTHMAPRES, 1,
                D3DUSAGE_DEPTHSTENCIL, D3DFMT_D24X8, D3DPOOL_DEFAULT, &depthMapTexture, NULL);
        depthMapTexture->GetSurfaceLevel(0, &depthMapDepthStencilBuffer);

        //create a stencil surface for PRM rendering (to avoid blend-adding to the same pixel for the same lightsource muliple times)
        device->CreateDepthStencilSurface(128 * 32, 128 * NCLUSTERSPERENTITY / 32, D3DFMT_D24S8, D3DMULTISAMPLE_NONE, 0, true, &prmBlendingDepthStencilBuffer, NULL);

#ifdef GENERATE_PATH_MAPS
        precompute();
#else
        loadPathMaps();
#endif  

        //create a texture for radion data (will be used for weight computations on the GPU)
        device->CreateTexture(2 * NRADIONS / 4096, 4096, 1, 0, D3DFMT_A32B32G32R32F, D3DPOOL_DEFAULT, &radionTexture, NULL);
        radionTexture->GetSurfaceLevel(0, &radionSurface);

        //fill texture with data
        uploadRadions();

        //create weights render target
        device->CreateTexture(4096, NRADIONS / 4096, 1, D3DUSAGE_RENDERTARGET, D3DFMT_R32F, D3DPOOL_DEFAULT, &weightsTexture, NULL);
        weightsTexture->GetSurfaceLevel(0, &weightsSurface);

        //create a sytem memory duplicate of the weights render target to be able to read back weights data
        device->CreateTexture(4096, NRADIONS / 4096, 1, 0, D3DFMT_R32F, D3DPOOL_SYSTEMMEM, &sysMemWeightsTexture, NULL);
        sysMemWeightsTexture->GetSurfaceLevel(0, &sysMemWeightsSurface);

        //create a vertex buffer to be able to visualize entry radion positions
        device->CreateVertexBuffer(NRADIONS * sizeof(float) * 6, D3DUSAGE_WRITEONLY, D3DFVF_XYZ | D3DFVF_NORMAL, D3DPOOL_DEFAULT, 
                &starterVertexBuffer, NULL);
        void* pData;
        starterVertexBuffer->Lock(0, 0, &pData, 0);
        fillRadionPosArray(pData);
        starterVertexBuffer->Unlock();

#ifdef GENERATE_PATH_MAPS
        savePathMaps();
#endif
}

PathMapEffect::~PathMapEffect(void)
{
        //release all resources allocated in the constructor
        starterVertexBuffer->Release();

        depthMapTexture->Release();
        depthMapDepthStencilBuffer->Release();
        fakeTexture->Release();
        fakeSurface->Release();

        sysMemWeightsTexture->Release();
        sysMemWeightsSurface->Release();
        weightsTexture->Release();
        weightsSurface->Release();
        radionTexture->Release();
        radionSurface->Release();

        prmBlendingDepthStencilBuffer->Release();
        delete kdtree;
        emptyTexture->Release();
        releaseTextures();
        releaseMeshes();
        releaseEntities();

        frameColorBuffer->Release();
        frameDepthStencilBuffer->Release();

        effect->Release();
        delete camera;
        delete lightCamera;
}

void PathMapEffect::loadMesh(LPCWSTR fileName)
{
        RenderMesh*             renderMesh = new RenderMesh();

        if(S_OK != 
                D3DXLoadMeshFromX(      fileName,
                                                        D3DXMESH_MANAGED,
                            device,
                                                        NULL,
                                                        &renderMesh->materialBuffer,
                                                        NULL,
                                                        (DWORD*)&renderMesh->nSubsets,
                                                        &renderMesh->mesh))
        {
                 MessageBox(NULL, fileName, L"Could not load mesh file!", MB_OK);
                 return;
        }
        //convert mesh into a vertex format using a position, a normal and a texture register
        //this is the only data our vertex shaders will require
        renderMesh->setVertexFormat(D3DFVF_XYZ |
                                                D3DFVF_NORMAL |
                                                D3DFVF_TEX1, device);

        //extract the D3DXMATERIAL* pointer from the generic D3DXBUFFER
        renderMesh->materials = (D3DXMATERIAL*)renderMesh->materialBuffer->GetBufferPointer();

        //for every material, load the specified texture
        renderMesh->textures = new LPDIRECT3DTEXTURE9[renderMesh->nSubsets];
        for(unsigned int i=0; i<renderMesh->nSubsets; i++)
        {
                wchar_t* wideTextureFileName = NULL;
                if(renderMesh->materials[i].pTextureFilename)
                {
                        wideTextureFileName = new wchar_t[MultiByteToWideChar(CP_ACP, 0, renderMesh->materials[i].pTextureFilename, -1, NULL, 0)];
                        MultiByteToWideChar(CP_ACP, 0, renderMesh->materials[i].pTextureFilename, -1, wideTextureFileName, 511);
                }
                renderMesh->textures[i] = loadTexture( wideTextureFileName );
                if(wideTextureFileName)
                        delete wideTextureFileName;
        }

        //lock the buffers and create ray tracing meshes from the data
        LPDIRECT3DVERTEXBUFFER9 vertexBuffer;
        renderMesh->mesh->GetVertexBuffer(&vertexBuffer);
        D3DVERTEX* vertexData;
        vertexBuffer->Lock(0,renderMesh->mesh->GetNumVertices()*renderMesh->mesh->GetNumBytesPerVertex(),(void**)&vertexData,0);
        LPDIRECT3DINDEXBUFFER9 indexBuffer;
        renderMesh->mesh->GetIndexBuffer(&indexBuffer);
        unsigned short* indexData;
        indexBuffer->Lock(0,renderMesh->mesh->GetNumFaces()*3*sizeof(unsigned short),(void**)&indexData,0);

        //create a ray-tracing Material for the texture
        if(renderMesh->textures[0])
        {
                D3DLOCKED_RECT lrect;
                HRESULT hr = renderMesh->textures[0]->LockRect(0, &lrect, NULL, D3DLOCK_READONLY);
                LPDIRECT3DSURFACE9 texsurf; 
                renderMesh->textures[0]->GetSurfaceLevel(0, &texsurf);
                D3DSURFACE_DESC desc;
                texsurf->GetDesc(&desc);
                renderMesh->rayTraceMaterial = new TexturedMaterial(lrect.pBits, lrect.Pitch, desc.Width);
                texsurf->Release();
                renderMesh->textures[0]->UnlockRect(0);
        }
        else
                renderMesh->rayTraceMaterial = new Material(Vector(1.0f, 1.0f, 1.0f), Vector(0.0f, 0.0f, 0.0f), 1.0);

        //create a TriangleMesh for ray-tracing
        renderMesh->rayTraceMesh = new TriangleMesh(renderMesh->rayTraceMaterial, indexData, renderMesh->mesh->GetNumFaces(),
                vertexData, renderMesh->mesh->GetNumVertices());

        vertexBuffer->Unlock();
        indexBuffer->Unlock();
        vertexBuffer->Release();
        indexBuffer->Release();

        //generate edge lines for seamless texture atlas rendering
        renderMesh->buildEdgeVertexBuffer(device);

        //store completed mesh
        renderMeshes.push_back(renderMesh);
}

LPDIRECT3DTEXTURE9 PathMapEffect::loadTexture(LPCWSTR fileName)
{
        if(fileName == NULL || wcscmp(fileName, L"") == 0)
                return NULL;
        wchar_t* mediaFileName = new wchar_t[wcslen(fileName) + 64];
        //we assume the x file was in folder .\media
        wcscpy(mediaFileName, L"media\\");
        wcscat(mediaFileName, fileName);
        LPDIRECT3DTEXTURE9 tex;

        D3DXIMAGE_INFO bobo;
        if(S_OK != 
                D3DXCreateTextureFromFileEx(device, mediaFileName, D3DX_DEFAULT, D3DX_DEFAULT, 1, 0,
                D3DFMT_FROM_FILE, D3DPOOL_MANAGED, D3DX_DEFAULT, D3DX_DEFAULT, 0, &bobo, NULL, &tex
                ))
        {
                 MessageBox(NULL, mediaFileName, L"Could not load texture file!", MB_OK);
                 return NULL;
        }
        
        materialTextures.push_back(tex);
        return tex;
}

void PathMapEffect::releaseTextures()
{
        std::vector<LPDIRECT3DTEXTURE9>::iterator i = materialTextures.begin();
        while(i != materialTextures.end())
        {
                (*i)->Release();
                i++;
        }
}

void PathMapEffect::releaseMeshes()
{
        std::vector<PathMapEffect::RenderMesh*>::iterator i = renderMeshes.begin();
        while(i != renderMeshes.end())
        {
                (*i)->mesh->Release();
                (*i)->materialBuffer->Release();
                delete (*i)->rayTraceMesh;
                delete (*i)->rayTraceMaterial;
                if((*i)->edgeVertexBuffer)
                        (*i)->edgeVertexBuffer->Release();
                delete *i;
                i++;
        }
}

void PathMapEffect::renderWithPRM()
{
        // first pass: render depth

        device->SetRenderTarget(0, fakeSurface);
        device->SetDepthStencilSurface(depthMapDepthStencilBuffer);
        device->Clear( 0, NULL, D3DCLEAR_ZBUFFER, D3DCOLOR_ARGB(255, 0, 0, 0), 1.0f, 0 );
        //no color writes
        device->SetRenderState(D3DRS_COLORWRITEENABLE, 0);
        //render back faces (avoids Z-fighting, no need for comparison bias)
        device->SetRenderState(D3DRS_CULLMODE, D3DCULL_CW);
        
        if( SUCCEEDED( device->BeginScene() ) )
        {
                effect->SetTechnique("Depth");
                UINT nPasses;
                effect->Begin(&nPasses, 0);
                effect->BeginPass(0);
                //loop through all entities, all mesh subsets, set proper transformations

                std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
                while(entityIterator != entities.end())
                {
                        effect->SetMatrix("modelToProjMatrix", 
                                &( entityIterator->modelWorldTransform * *lightCamera->GetViewMatrix() * *lightCamera->GetProjMatrix() ) );
                        effect->CommitChanges();
                        unsigned int nSubsets = entityIterator->renderMesh->nSubsets;
                        for(int i = 0; i< nSubsets; i++)
                        {
                                entityIterator->renderMesh->mesh->DrawSubset(i);
                        }
                        entityIterator++;
                }

                effect->EndPass();
                effect->End();
                device->EndScene();
        }

        //reset front face rendering, color writes
        device->SetRenderState(D3DRS_CULLMODE, D3DCULL_CCW);
        HRESULT hr = S_OK;
        hr = device->SetRenderState(D3DRS_COLORWRITEENABLE, D3DCOLORWRITEENABLE_ALPHA| D3DCOLORWRITEENABLE_BLUE| D3DCOLORWRITEENABLE_GREEN | D3DCOLORWRITEENABLE_RED);

        // 2. pass: compute weights

        // simple parallel computation for all pixels: no 3D, depth, shadows, only a full-texture quad
        hr = device->SetRenderState(D3DRS_ZENABLE, D3DZB_FALSE);
        hr = device->SetRenderState(D3DRS_STENCILENABLE, false);
        hr = device->SetRenderTarget(0,weightsSurface);
        hr = device->SetDepthStencilSurface(NULL);
        D3DXVECTOR3 lightPos = *lightCamera->GetEyePt();
        D3DXVECTOR3 lightDir = *lightCamera->GetWorldAhead();
        lightDir /= D3DXVec3Length( &lightDir );
        float ffl = parameters->Get(fLightScale) * 1.0f;
        D3DXVECTOR3 lightPower(ffl, ffl, ffl);

        device->Clear( 0, NULL, D3DCLEAR_TARGET,        D3DCOLOR_XRGB(0,0,0), 1.0f, 0 );

        if( SUCCEEDED( device->BeginScene() ) ){
                if( effect != NULL ){
                        D3DXHANDLE hTechnique = NULL;
                        hTechnique = effect->GetTechniqueByName((LPSTR)"ComputeWeights");
                        if(hTechnique==NULL){
                                return;
                        }
                        effect->SetTechnique( hTechnique );

                        UINT nPasses = 0;
                        effect->Begin( &nPasses, 0 );
                        for(UINT j=0;j<nPasses;j++){ 
                                effect->BeginPass(j);
                                hr = effect->SetTexture("radions", radionTexture);
                                static float opttme[9] = {0.5f, 0.0f, 0.0f,
                                                                                0.0f, -0.5f, 0.0f,
                                                                                0.5f + (0.5f / DEPTHMAPRES), 0.5f + (0.5f / DEPTHMAPRES), 1.0f};

                                //set global params
                                effect->SetFloatArray("occProjToTexMatrix", opttme, 9);
                                effect->SetMatrix("occWorldToProjMatrix", &(*lightCamera->GetViewMatrix() * *lightCamera->GetProjMatrix()));

                                effect->SetTexture("depthMap", depthMapTexture);
                                hr = effect->SetInt("nRadionColumns", NRADIONS / 4096);
                                hr = effect->SetFloatArray("lightPower", (float*)&lightPower, 3);
                                hr = effect->SetFloatArray("lightPos", (float*)&lightPos, 3);
                                hr = effect->SetFloatArray("lightDir", (float*)&lightDir, 3);
                                effect->CommitChanges();
                                drawFullScreenQuad(device);
                                effect->EndPass();
                        }
                        effect->End();
                }
                device->EndScene();
        }

        // read the weights back

        D3DLOCKED_RECT rData;
        hr = device->GetRenderTargetData(weightsSurface, sysMemWeightsSurface);
        hr =  sysMemWeightsSurface->LockRect(&rData, NULL, D3DLOCK_READONLY);
        float* allEntryWeights = (float*)rData.pBits;

        // average weights per cluster
        float sumWeightsAll = 0.0;
        unsigned int iRadion = 0;
        for(int iCluster=0; iCluster < NCLUSTERS; iCluster++)
        {
                weights[iCluster] = 0.0;
                float clusterrad = 0.0;
                for(int clusterSummer=0; clusterSummer < clusterLenghts[iCluster]; clusterSummer++, iRadion++)
                {
                        float radrad = bushStarters[iRadion].radiance.sum();
                        weights[iCluster] += allEntryWeights[iRadion] * radrad;
                        clusterrad += radrad;
                }
                weights[iCluster] /= clusterrad; //(double)clusterLenghts[iCluster];
//              if(iCluster != TEST_CLUST)
//                      weights[iCluster] = 0.0;
//              if(weights[iCluster] > 0.005)
//                      weights[iCluster] = 0.005;
                sumWeightsAll += weights[iCluster];
        }
        sysMemWeightsSurface->UnlockRect();

        // 3. pass: render scene using weights to combine PRM texture atlases

        
        device->SetRenderTarget(0, frameColorBuffer);
        device->SetDepthStencilSurface(frameDepthStencilBuffer);
        // use backface culling, depth test
        hr = device->SetRenderState(D3DRS_CULLMODE, D3DCULL_CCW);
        hr = device->SetRenderState(D3DRS_ZENABLE, D3DZB_TRUE);

        device->Clear( 0, NULL, D3DCLEAR_TARGET | D3DCLEAR_ZBUFFER,     D3DCOLOR_XRGB(0,0,0), 1.0f, 0 );
        
        if( SUCCEEDED( device->BeginScene() ) )
        {
                effect->SetTechnique("walk");
                UINT nPasses;
                effect->Begin(&nPasses, 0);
                effect->BeginPass(0);

                static float opttme[9] = {0.5f, 0.0f, 0.0f,
                                                                0.0f, -0.5f, 0.0f,
                                                                0.5f + (0.5f / DEPTHMAPRES), 0.5f + (0.5f / DEPTHMAPRES), 1.0f};

                //set global params
                effect->SetFloatArray("occProjToTexMatrix", opttme, 9);
                effect->SetMatrix("occWorldToProjMatrix", &(*lightCamera->GetViewMatrix() * *lightCamera->GetProjMatrix()));

                effect->SetTexture("depthMap", depthMapTexture);

                //loop through entities
                std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
                while(entityIterator != entities.end())
                {
                        //set entity params
                        // nearestIdx -> bushId
                        float weightsa[NCLUSTERSPERENTITY];
                        for(unsigned int y=0; y<NCLUSTERSPERENTITY; y++)
                                weightsa[y] = weights[entityIterator->nearClusterIndices[y]];
                        //## could sort weightsa, indicesa and use first few
                        hr = effect->SetFloatArray("weightsa", weightsa, NCLUSTERSPERENTITY);
                        INT* pr = (INT*)entityIterator->nearClusterIndices;
                        hr = effect->SetIntArray("indicesa", (INT*)entityIterator->nearClusterIndices, NCLUSTERSPERENTITY);
                        effect->SetMatrix("modelToWorldMatrix", &entityIterator->modelWorldTransform);
                        effect->SetMatrix("inverseTransposedModelToWorldMatrix", &entityIterator->inverseTransposedModelWorldTransform);
                        D3DXMATRIX t = entityIterator->modelWorldTransform;
                        if(parameters->Get(bLookFromLight))
                                effect->SetMatrix("modelToProjMatrix", &(entityIterator->modelWorldTransform * *lightCamera->GetViewMatrix()  * *lightCamera->GetProjMatrix()  ));
                        else
                                effect->SetMatrix("modelToProjMatrix", &(entityIterator->modelWorldTransform * *camera->GetViewMatrix()  * *camera->GetProjMatrix()  ));

                        effect->SetTexture("filteredAtlas", entityIterator->prmTexture);
                        hr = effect->SetFloatArray("lightPower", (float*)&lightPower, 3);
                        hr = effect->SetFloatArray("lightPos", (float*)&lightPos, 3);
                        hr = effect->SetFloatArray("lightDir", (float*)&lightDir, 3);

                        unsigned int nSubsets = entityIterator->renderMesh->nSubsets;
                        for(int i = 0; i< nSubsets; i++)
                        {
                                //set subset material texture
                                if(entityIterator->renderMesh->textures[i])
                                        effect->SetTexture("brdfMap", entityIterator->renderMesh->textures[i]);
                                else
                                        effect->SetTexture("brdfMap", emptyTexture);
                                effect->CommitChanges();
                                entityIterator->renderMesh->mesh->DrawSubset(i);
                        }
                        entityIterator++;
                }
                effect->EndPass();
                effect->End();
                device->EndScene();
        }
}

HRESULT PathMapEffect::RenderMesh::setVertexFormat(DWORD fvf, LPDIRECT3DDEVICE9 device)
{
        LPD3DXMESH tempMesh;

        //clone the mesh to the appropriate format
    HRESULT hr = mesh->CloneMeshFVF( mesh->GetOptions(), 
                                  fvf , 
                                  device, &tempMesh );
        if(hr == S_OK)
        {
                DWORD* padj = new DWORD[mesh->GetNumFaces() * 3];
                mesh->GenerateAdjacency(0.00001f, padj);
                delete padj;
                //forget he old mesh
                mesh->Release();
                //use the cloned mesh
                mesh = tempMesh;
        }
        return hr;
}

void PathMapEffect::renderFullScreen(float depth)
{
        float fLeftU = 0.0f, fTopV = 0.0f, fRightU = 1.0f, fBottomV = 1.0f;

    D3DSURFACE_DESC dtdsdRT;
    PDIRECT3DSURFACE9 pSurfRT;

    // Acquire render target width and height
    device->GetRenderTarget(0, &pSurfRT);
    pSurfRT->GetDesc(&dtdsdRT);
    pSurfRT->Release();

    // Ensure that we're directly mapping texels to pixels by offset by 0.5
    // For more info see the doc page titled "Directly Mapping Texels to Pixels"
    FLOAT fWidth5 = (FLOAT)dtdsdRT.Width - 0.5f;
    FLOAT fHeight5 = (FLOAT)dtdsdRT.Height - 0.5f;

    // Draw the quad
        struct D3DVERTEX{
                D3DXVECTOR3 pos;
                D3DXVECTOR2 tex0;
        };

    D3DVERTEX svQuad[4];

    svQuad[0].pos=D3DXVECTOR3(-1, 1, depth);
    svQuad[0].tex0 = D3DXVECTOR2(fLeftU, fTopV);

    svQuad[1].pos = D3DXVECTOR3(1, 1, depth);
    svQuad[1].tex0 = D3DXVECTOR2(fRightU, fTopV);

    svQuad[2].pos = D3DXVECTOR3(-1, -1, depth);
    svQuad[2].tex0 = D3DXVECTOR2(fLeftU, fBottomV);

    svQuad[3].pos = D3DXVECTOR3(1, -1, depth);
    svQuad[3].tex0 = D3DXVECTOR2(fRightU, fBottomV);
        
        device->SetRenderState(D3DRS_CULLMODE,D3DCULL_NONE);
        device->SetFVF(D3DFVF_XYZ | D3DFVF_TEX1 | D3DFVF_TEXCOORDSIZE2(0));
    device->DrawPrimitiveUP(D3DPT_TRIANGLESTRIP, 2, svQuad, sizeof(D3DVERTEX));
}

void PathMapEffect::move(float fElapsedTime)
{
        //apply some rotation to one of the entities
        if(parameters->Get(bTurbo))
        {
                camera->FrameMove(fElapsedTime * 10.0f);
                lightCamera->FrameMove(fElapsedTime * 10.0f);
        }
        else
        {
                camera->FrameMove(fElapsedTime);
                lightCamera->FrameMove(fElapsedTime);
        }
/*      D3DXMATRIX rot;
        D3DXMatrixRotationY(&rot, fElapsedTime * 0.1f);
        D3DXMatrixMultiply( &entities.at(0).modelWorldTransform, 
                                                &rot,
                                                &entities.at(0).modelWorldTransform);
        D3DXMatrixInverse(  &entities.at(0).inverseTransposedModelWorldTransform, NULL, 
                                                &entities.at(0).modelWorldTransform);
        D3DXMatrixTranspose( &entities.at(0).inverseTransposedModelWorldTransform, &entities.at(0).inverseTransposedModelWorldTransform);

        if(entities.size() > 1)
        {
                D3DXMatrixTranslation(&entities.at(1).modelWorldTransform, 0.0f, parameters->Get(fSecondModelHeight) * 5.0f , 0.0f);
                D3DXMatrixTranslation(&entities.at(1).inverseTransposedModelWorldTransform, 0.0f, -parameters->Get(fSecondModelHeight) * 5.0f , 0.0f);
        }*/
}

LRESULT PathMapEffect::handleMessage( HWND hWnd, UINT uMsg, WPARAM wParam, LPARAM lParam)

{
        if(parameters->Get(bMoveLight))
                return lightCamera->HandleMessages(hWnd, uMsg, wParam, lParam);
        else
                return camera->HandleMessages(hWnd, uMsg, wParam, lParam);
}

void PathMapEffect::addUiParameters(Parameters* parameters)
{
        PathMapEffect::parameters = parameters;
        parameters->Add( bLookFromLight, "Look from light", 1);
        parameters->Add( bMoveLight, "Move light", 1);
        parameters->Add( bDots, "Entry points", 1);
        parameters->Add( bTurbo, "Turbo", 1);
        parameters->Add( fLightScale, "Tone scale", 100, 'V', 'B', convert100);
}

void PathMapEffect::setUiParameterDefaults(Parameters* parameters)
{
}

Parameters*     PathMapEffect::parameters = NULL;

void PathMapEffect::render()
{
        if(method == Method::PRM)
                renderWithPRM();
        else
                showPRMTexture();
        if(parameters->Get(bDots))
        {
                float psf = 8.0f;
                HRESULT hr = device->SetRenderState(D3DRS_POINTSIZE, *((DWORD*)&psf));
                hr = device->SetRenderState(D3DRS_ZENABLE, D3DZB_TRUE);
                if( SUCCEEDED( device->BeginScene() ) )
                {
                        D3DXHANDLE hTechnique = NULL;
                        hTechnique = effect->GetTechniqueByName((LPSTR)"Dots");
                        effect->SetTechnique( hTechnique );
                        UINT nPasses = 0;
                        effect->Begin( &nPasses, 0 );
                        for(UINT j=0;j<nPasses;j++){
                                effect->BeginPass(j);
                                        hr = effect->SetTexture("depthMap",depthMapTexture);
                                        if(parameters->Get(bLookFromLight))
                                                effect->SetMatrix("worldToProjMatrix", &(*lightCamera->GetViewMatrix()  * *lightCamera->GetProjMatrix()  ));
                                        else
                                                effect->SetMatrix("worldToProjMatrix", &(*camera->GetViewMatrix()  * *camera->GetProjMatrix()  ));

                                                        static float opttme[9] = {0.5f, 0.0f, 0.0f,
                                                                0.0f, -0.5f, 0.0f,
                                                                0.5f + (0.5f / DEPTHMAPRES), 0.5f + (0.5f / DEPTHMAPRES), 1.0f};

                                        effect->SetFloatArray("occProjToTexMatrix", opttme, 9);
                                        D3DXVECTOR3 lightPos = *lightCamera->GetEyePt();
                                        hr = effect->SetFloatArray("lightPos", (float*)&lightPos, 3);
                                        effect->CommitChanges();
                                        hr = device->SetFVF(D3DFVF_XYZ |  D3DFVF_NORMAL);
                                        hr = device->SetStreamSource(0, starterVertexBuffer, 0, sizeof(float) * 6);
                                        UINT offset = 0;
                                        for(int iClusterDraw=0; iClusterDraw < NCLUSTERS; iClusterDraw++)
                                        {
                                                hr = effect->SetFloat("aClusterWeight", weights[iClusterDraw] * 10.0);
                                                effect->CommitChanges();
//                                              if(iClusterDraw == TEST_CLUST)
                                                device->DrawPrimitive(D3DPT_POINTLIST, offset, clusterLenghts[iClusterDraw]);
                                                offset += clusterLenghts[iClusterDraw];
                                        }
                                effect->EndPass();
                        }
                        effect->End();
                }
                device->EndScene();
        }
}

void PathMapEffect::Entity::createRayTraceEntity()
{
        rayTraceEntity = new Transformed(renderMesh->rayTraceMesh);
        rayTraceEntity->setTransforms(modelWorldTransform, inverseTransposedModelWorldTransform);
}

void PathMapEffect::releaseEntities()
{
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                delete entityIterator->rayTraceEntity;
                entityIterator->prmSurface->Release();
                entityIterator->prmTexture->Release();
                entityIterator++;
        }       
}

void PathMapEffect::sampleSphereDirection(Vector& outDir)
{
        do{
        outDir = Vector(-1.0 + 2.0 * (float)rand() / RAND_MAX, -1.0 + 2.0 * (float)rand() / RAND_MAX, -1.0 + 2.0 * (float)rand() / RAND_MAX);
        }while(outDir.norm2() > 1.0);
        outDir.normalize();
}

void PathMapEffect::sampleShootingDiffuseDirection(int depth, const Vector& normal, Vector& outDir)
{
        Vector u, v;
        if(fabsf(normal[0]) < 0.9f)
        {
                u.set(0.0f, -normal[2], normal[1]);
                u *= 1.0f / sqrtf(normal[2] * normal[2] + normal[1] * normal[1]);
        }
        else
        {
                u.set(normal[2], 0.0f, -normal[0]);
                u *= 1.0f / sqrtf(normal[2] * normal[2] + normal[0] * normal[0]);
        }
        v.setCrossProduct(normal, u);

        float phi = 2.0f * 3.14159265358979323846f * (float)rand() / RAND_MAX;
        float xi2 = (float)rand() / RAND_MAX;
        float cosPhi = cos(phi);
        float sinPhi = sin(phi);
        float cosTheta = sqrtf(xi2);
        float sinTheta = sqrtf(1.0f - xi2);

        outDir.setScaled(sinTheta * cosPhi, v);
        outDir.addScaled(sinTheta * sinPhi, u);
        outDir.addScaled(cosTheta, normal);
}

void PathMapEffect::createPRMTextures()
{
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                device->CreateTexture(128 * 32, 128 * NCLUSTERSPERENTITY / 32, 1, D3DUSAGE_RENDERTARGET,D3DFMT_A16B16G16R16F,D3DPOOL_DEFAULT,&entityIterator->prmTexture,NULL);
                entityIterator->prmTexture->GetSurfaceLevel(0,&entityIterator->prmSurface);
                entityIterator++;
        }       
}

void PathMapEffect::precompute()
{
        // generate entry points
        for(unsigned int cRad=0; cRad < NRADIONS; cRad++)
        {
                Radion starter;
                sampleSurfaceRadion(starter);
                //sampleSurfaceRadionUniform(starter);
                bushStarters.push_back(starter);
        }

        // sort entry radions into clusters
        Radion* entryArray = (Radion*)&*bushStarters.begin();
        clusterRadions(entryArray, bushStarters.size(), 0);
        clusterRadionsKMeans();

        // for every entity, find the most relevant clusters
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                entityIterator->findNearClusters(bushStarters, NCLUSTERS);
                entityIterator++;
        }

        // for every entry radion, shoot a bush of radions, and add their
        // contributions to all PRMs of all entities
        unsigned int iCluster = 0;
        std::vector<Radion> bushRadions;
        for(int iStarter = 0; iStarter < NRADIONS; )
        {
                bushRadions.clear();
                for(int iric=0; iric < clusterLenghts[iCluster]; iric++, iStarter++)
                        shootRadionBush(bushStarters[iStarter], bushRadions);
                std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
                while(entityIterator != entities.end())
                {
                        entityIterator->renderPRM(device, effect, bushRadions, iCluster, 1.0);
                        entityIterator++;
                }
                iCluster++;
        }
}

int compareX(const void* a, const void* b)
{
        if( ((const Radion*)a)->position.x < ((const Radion*)b)->position.x )
                return -1;
        else
                return 1;
}

int compareY(const void* a, const void* b)
{
        if( ((const Radion*)a)->position.y < ((const Radion*)b)->position.y )
                return -1;
        else
                return 1;
}

int compareZ(const void* a, const void* b)
{
        if( ((const Radion*)a)->position.z < ((const Radion*)b)->position.z )
                return -1;
        else
                return 1;
}

int compareDirX(const void* a, const void* b)
{
        if( ((const Radion*)a)->normal.x < ((const Radion*)b)->normal.x )
                return -1;
        else
                return 1;
}

int compareDirY(const void* a, const void* b)
{
        if( ((const Radion*)a)->normal.y < ((const Radion*)b)->normal.y )
                return -1;
        else
                return 1;
}

int compareDirZ(const void* a, const void* b)
{
        if( ((const Radion*)a)->normal.z < ((const Radion*)b)->normal.z )
                return -1;
        else
                return 1;
}


void PathMapEffect::clusterRadionsKMeans()
{
        Radion centroids[NCLUSTERS];
        std::vector<Radion> clusters[NCLUSTERS];
        //initial clusters: uniform length
        for(unsigned int i=0; i<NCLUSTERS; i++)
                clusterLenghts[i] = NRADIONS / NCLUSTERS;
        for(unsigned int iKMeans = 0; iKMeans < 10; iKMeans++)
        {
                //find centroids
                unsigned int iRadion = 0;
                for(unsigned int i=0; i<NCLUSTERS; i++)
                {
                        centroids[i].position = Vector::RGBBLACK;
                        centroids[i].normal = Vector::RGBBLACK;
                        for(unsigned int iRadionInCluster=0; iRadionInCluster < clusterLenghts[i]; iRadionInCluster++, iRadion++)
                        {
                                centroids[i].position += bushStarters.at(iRadion).position;
                                centroids[i].normal += bushStarters.at(iRadion).normal;
                        }
                        centroids[i].position *= 1.0f / (float)clusterLenghts[i];
                        centroids[i].normal.normalize();
                }
                
                for(unsigned int i=0; i<NCLUSTERS; i++)
                        clusters[i].clear();

                //sort radions to centroids
                iRadion = 0;
                for(unsigned int i=0; i<NCLUSTERS; i++)
                {
                        for(unsigned int iRadionInCluster=0; iRadionInCluster < clusterLenghts[i]; iRadionInCluster++, iRadion++)
                        {
                                unsigned int minimumDistanceClusterIndex = 0;
                                float minimumDistance = FLT_MAX;
                                for(unsigned int u=0; u<NCLUSTERS; u++)
                                {
                                        float dirDist = 2.0f - (bushStarters.at(iRadion).normal * centroids[u].normal);
                                        float dist = (bushStarters.at(iRadion).position - centroids[u].position).norm() * dirDist;
                                        if(dist < minimumDistance)
                                        {
                                                minimumDistanceClusterIndex = u;
                                                minimumDistance = dist;
                                        }
                                }
                                clusters[minimumDistanceClusterIndex].push_back( bushStarters.at(iRadion) );
                        }
                }
                //refill bushStarters, set cluster lengths
                iRadion = 0;
                for(unsigned int i=0; i<NCLUSTERS; i++)
                {
                        clusterLenghts[i] = clusters[i].size();
                        for(unsigned int iRadionInCluster=0; iRadionInCluster < clusterLenghts[i]; iRadionInCluster++, iRadion++)
                        {
                                bushStarters.at(iRadion) = clusters[i].at(iRadionInCluster);
                        }
                }
                //eliminate empty clusters
                for(unsigned int i=1; i<NCLUSTERS; i++)
                {
                        if(clusterLenghts[i] == 0)
                        {
                                clusterLenghts[i] = clusterLenghts[i-1] >> 1;
                                clusterLenghts[i-1] -= clusterLenghts[i-1] >> 1;
                        }
                }
                for(int i=NCLUSTERS - 1; i >= 0; i--)
                {
                        if(clusterLenghts[i] == 0)
                        {
                                clusterLenghts[i] = clusterLenghts[i+1] >> 1;
                                clusterLenghts[i+1] -= clusterLenghts[i+1] >> 1;
                        }
                }
        }
}

int PathMapEffect::clusterRadions(Radion* partition, int psize, char axis)
{
/*      if(psize < 2)
                return 1;
        if(axis == 0)
                qsort(partition, psize, sizeof(Radion), compareDirX);
        else if(axis == 1)
                qsort(partition, psize, sizeof(Radion), compareDirY);
        else if(axis == 2)
                qsort(partition, psize, sizeof(Radion), compareDirZ);
        else if(axis == 3)
                qsort(partition, psize, sizeof(Radion), compareX);
        else if(axis == 4)
                qsort(partition, psize, sizeof(Radion), compareY);
        else if(axis == 5)
                qsort(partition, psize, sizeof(Radion), compareZ);
        clusterRadions(partition, psize >> 1, (axis+1)%6);
        clusterRadions(partition + (psize >> 1), psize - (psize >> 1), (axis+1)%6);
        return 1;
/*/     if(psize < 2)
                return 1;
        if(axis == 0)
                qsort(partition, psize, sizeof(Radion), compareX);
        else if(axis == 1)
                qsort(partition, psize, sizeof(Radion), compareY);
        else if(axis == 2)
                qsort(partition, psize, sizeof(Radion), compareZ);
        clusterRadions(partition, psize >> 1, (axis+1)%3);
        clusterRadions(partition + (psize >> 1), psize - (psize >> 1), (axis+1)%3);
        return 1;//*/
}


void PathMapEffect::shootRadionBush(Radion& starter, std::vector<Radion>& bushRadions)
{
        bushRadions.push_back(starter);

        int nPhotonsShotForLight = 1;   //branching factor
        for(int iPhoton = 0; iPhoton < nPhotonsShotForLight; iPhoton++)
        {
                unsigned int depth = 0;
                ray.id = rayId++;
                ray.isShadowRay = false;
                ray.origin = starter.position;
                Vector power = starter.radiance;
                float prob = starter.probability;
                prob *= nPhotonsShotForLight;
                power *= 1.0f / nPhotonsShotForLight;
                sampleShootingDiffuseDirection(depth, starter.normal, ray.dir);

                for(;;depth++)
                {
                        kdtree->traverse(ray, hitRec, 0.0f, FLT_MAX);
                        if(hitRec.isIntersect)
                        {
                                Radion nr;
                                nr.position = hitRec.point;
                                nr.normal = hitRec.normal;
                                nr.radiance = power;
                                nr.radiance %= hitRec.material->getTextureDiffuseBrdf(hitRec.texUV);
                                nr.probability = prob;
                                bushRadions.push_back(nr);
                        }
                        else
                                break;
                        if(depth >= 3)
                                break;
                        float rrRandom = (float)rand() / RAND_MAX;
                        float diffuseAlbedo = hitRec.material->getTextureDiffuseAlbedo(hitRec.texUV);
                        float idealAlbedo = hitRec.material->getTextureIdealAlbedo(hitRec.texUV);
                        float refractiveAlbedo = hitRec.material->getRefractiveAlbedo();
                        if(rrRandom < diffuseAlbedo)
                        {
                                ray.id = rayId++;
                                ray.isShadowRay = false;
                                ray.origin = hitRec.point;
                                power %= hitRec.material->getTextureDiffuseBrdf(hitRec.texUV);
                                power *= 1.0f / diffuseAlbedo;
                                prob *= diffuseAlbedo;
                                sampleShootingDiffuseDirection(depth, hitRec.normal, ray.dir);
                        }
                        else
                        {
                                rrRandom -= diffuseAlbedo;
                                if(rrRandom < idealAlbedo)
                                {
                                        ray.dir.setIdealReflectedDirection(ray.dir, hitRec.normal);
                                        ray.id = rayId++;
                                        ray.isShadowRay = false;
                                        ray.origin = hitRec.point;
                                        power %= hitRec.material->getIdealBrdf();
                                        power *= 1.0f / idealAlbedo;
                                        prob *= idealAlbedo;
                                }
                                else
                                        break;
                        }
                }
        }
}

void PathMapEffect::sampleSurfaceRadion(Radion& starter)
{
        float randa = ((double)rand() / RAND_MAX) * entities.size();
//      float randa = ((double)rand() / RAND_MAX) * sumSurfaceArea;
        float uptonowSurfaceArea = 0.0;
        Entity* e = NULL;
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                uptonowSurfaceArea += 1.0f;
//              uptonowSurfaceArea += entityIterator->rayTraceEntity->getSurfaceArea();
                if(uptonowSurfaceArea >= randa)
                {
                        e = &*entityIterator;
                        break;
                }
                entityIterator++;
        }

        e->rayTraceEntity->sampleSurface(starter);

        float surfa = starter.radiance.z;
        starter.radiance = e->rayTraceEntity->getMaterial()->getTextureDiffuseBrdf(starter.radiance);
        starter.radiance *= (surfa / NRADIONS) * entities.size();
//      starter.radiance *= sumSurfaceArea / NRADIONS;
        starter.probability = (float)NRADIONS / (entities.size() * surfa );
}

void PathMapEffect::sampleSurfaceRadionUniform(Radion& starter)
{
        float randa = ((double)rand() / RAND_MAX) * sumSurfaceArea;
        float uptonowSurfaceArea = 0.0;
        Entity* e = NULL;
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                uptonowSurfaceArea += entityIterator->rayTraceEntity->getSurfaceArea();
                if(uptonowSurfaceArea >= randa)
                {
                        e = &*entityIterator;
                        break;
                }
                entityIterator++;
        }

        e->rayTraceEntity->sampleSurface(starter);

        float surfa = starter.radiance.z;
        starter.radiance = e->rayTraceEntity->getMaterial()->getTextureDiffuseBrdf(starter.radiance);
        starter.radiance *= sumSurfaceArea / NRADIONS;
        starter.probability = NRADIONS / sumSurfaceArea;
}


void PathMapEffect::Entity::renderPRM(LPDIRECT3DDEVICE9 device,
                                                        LPD3DXEFFECT effect,
                                                        std::vector<Radion>& bushRadions,
                                                        unsigned int bushId,
                                                        float scaleFactor
                                                        )
{
        HRESULT hr;     UINT nPasses=1;

        //# bushId -> nearestIdx
        int radinx = 0xffffff;
        for(int i=0; i<NCLUSTERSPERENTITY; i++)
                if(nearClusterIndices[i] == bushId)
                        radinx = i;
        if(radinx == 0xffffff)
                return;

        hr = device->SetRenderTarget(0, prmSurface);
        int tilex = radinx % 32;
        int tiley = radinx / 32;

        D3DVIEWPORT9 vp;
        vp.X = tilex * 128;                     vp.Y = tiley * 128;
        vp.Width = vp.Height = 128;
        vp.MinZ = 0.0; vp.MaxZ = 1.0;

        //clear everything
        device->SetRenderTarget(0, prmSurface);
        device->SetDepthStencilSurface(owner->prmBlendingDepthStencilBuffer);
        device->SetViewport(&vp);

        device->Clear( 0, NULL, D3DCLEAR_TARGET,
                                                D3DCOLOR_RGBA(0,0,0,0), 1.0f, 0 );

        unsigned int nRadions = bushRadions.size();
        for(int ir=0; ir < nRadions; ir++)
        {
                // 1. pass render depth map for virtual light source
                device->SetRenderTarget(0, owner->fakeSurface);
                device->SetDepthStencilSurface(owner->depthMapDepthStencilBuffer);

                device->SetRenderState(D3DRS_COLORWRITEENABLE, 0);

                device->Clear( 0, NULL, D3DCLEAR_ZBUFFER, D3DCOLOR_ARGB(255, 0, 0, 0), 1.0f, 0 );

                CFirstPersonCamera radionCamera;
                radionCamera.SetViewParams((D3DXVECTOR3*)&bushRadions[ir].position, (D3DXVECTOR3*)&(bushRadions[ir].position + bushRadions[ir].normal));
                radionCamera.SetProjParams( D3DX_PI/1.9, 1.0f, 0.1f, 300.0f );

                // render backfaces to avoid z-fighting
                device->SetRenderState(D3DRS_CULLMODE, D3DCULL_CW);
                if( SUCCEEDED( device->BeginScene() ) )
                {
                        effect->SetTechnique("Depth");
                        UINT nPasses;
                        effect->Begin(&nPasses, 0);
                        effect->BeginPass(0);
                        //loop through all entities, all mesh subsets, set proper transformations

                        std::vector<PathMapEffect::Entity>::iterator entityIterator = owner->entities.begin();
                        while(entityIterator != owner->entities.end())
                        {
                                effect->SetMatrix("modelToProjMatrix", 
                                        &( entityIterator->modelWorldTransform * *radionCamera.GetViewMatrix() * *radionCamera.GetProjMatrix() ) );
                                effect->CommitChanges();
                                unsigned int nSubsets = entityIterator->renderMesh->nSubsets;
                                for(int i = 0; i< nSubsets; i++)
                                {
                                        entityIterator->renderMesh->mesh->DrawSubset(i);
                                }
                                entityIterator++;
                        }

                        effect->EndPass();
                        effect->End();
                        device->EndScene();
                }

                device->SetRenderState(D3DRS_CULLMODE, D3DCULL_CCW);
                device->SetRenderState(D3DRS_COLORWRITEENABLE, D3DCOLORWRITEENABLE_RED | D3DCOLORWRITEENABLE_GREEN | D3DCOLORWRITEENABLE_BLUE | D3DCOLORWRITEENABLE_ALPHA);

                // blend-add contribution of virtual light source to PRM
                device->SetRenderTarget(0, prmSurface);
                device->SetDepthStencilSurface(owner->prmBlendingDepthStencilBuffer);
                device->SetViewport(&vp);
                device->SetRenderState(D3DRS_CULLMODE,  D3DCULL_NONE);
                device->SetRenderState(D3DRS_ALPHABLENDENABLE,  true);
                device->SetRenderState(D3DRS_SRCBLEND, D3DBLEND_ONE);
                device->SetRenderState(D3DRS_DESTBLEND, D3DBLEND_ONE);
                device->SetRenderState(D3DRS_BLENDOP, D3DBLENDOP_ADD);
                device->SetRenderState(D3DRS_ZENABLE, D3DZB_FALSE);

                // use stencil to avoid adding contribution twice
                HRESULT hr = device->Clear( 0, NULL, D3DCLEAR_STENCIL,
                                        D3DCOLOR_RGBA(0,0,0,0), 1.0f, 0 );
                if( SUCCEEDED( device->BeginScene() ) )
                {
                        Radion radion = bushRadions[ir];

                        if( effect != NULL ){
                                D3DXHANDLE hTechnique = effect->GetTechniqueByName((LPSTR)"BushToAtlas");
                                if(hTechnique==NULL){
                                        return;
                                }
                                effect->SetTechnique( hTechnique );
                                
                                effect->Begin( &nPasses, 0 );
                                for(UINT j=0;j<nPasses;j++){
                                        effect->BeginPass(j);
                                        Vector fakejake = bushRadions[ir].position;
                                        fakejake = bushRadions[ir].normal;
                                        fakejake = bushRadions[ir].radiance;
                                        hr=effect->SetFloatArray("lightPos", (float*)&bushRadions[ir].position, 3);
                                        hr=effect->SetFloatArray("lightDir", (float*)&bushRadions[ir].normal, 3);
                                        static float opttme[9] = {0.5f, 0.0f, 0.0f,
                                                                                        0.0f, -0.5f, 0.0f,
                                                                                        0.5f + (0.5f / DEPTHMAPRES), 0.5f + (0.5f / DEPTHMAPRES), 1.0f};
                                        effect->SetTexture("depthMap", owner->depthMapTexture);

                                        Vector bb0 = this->rayTraceEntity->bbox.minPoint;
                                        Vector bb1 = this->rayTraceEntity->bbox.maxPoint;
                                        //effect->SetFloat("cutNearness2", (bb1 - bb0).norm2() / 200.0f);
                                        float nne = 5.0f / bushRadions[ir].probability;
                                        float onne = (bb1 - bb0).norm2() / 200.0f;
                                        effect->SetFloat("cutNearness2", nne);
                                        //effect->SetFloat("cutNearness2", 0.0f);
                                        
                                        //set global params
                                        effect->SetFloatArray("occProjToTexMatrix", opttme, 9);
                                        effect->SetMatrix("occWorldToProjMatrix", &(*radionCamera.GetViewMatrix() * *radionCamera.GetProjMatrix()));
                                        effect->SetMatrix("modelToWorldMatrix", &(this->modelWorldTransform));
                                        effect->SetMatrix("inverseTransposedModelToWorldMatrix", &(this->inverseTransposedModelWorldTransform));
                                        Vector scaledPower = bushRadions[ir].radiance * scaleFactor;
                                        hr=effect->SetFloatArray("lightPower", (float*)&scaledPower, 3);
                                        
                                        hr = device->SetRenderState(D3DRS_STENCILENABLE, true);
                                        hr = device->SetRenderState(D3DRS_STENCILREF, 0xffffffff);
                                        hr = device->SetRenderState(D3DRS_STENCILFUNC, D3DCMP_ALWAYS);
                                        hr = device->SetRenderState(D3DRS_STENCILPASS, D3DSTENCILOP_REPLACE);

                                        effect->CommitChanges();

                                        renderMesh->mesh->DrawSubset(0);

                                        hr = device->SetRenderState(D3DRS_STENCILFUNC, D3DCMP_NOTEQUAL);
                                        scaledPower = bushRadions[ir].radiance * 1.0f ;
                                        hr=effect->SetFloatArray("lightPower", (float*)&scaledPower, 3);
                                        effect->CommitChanges();

                                        if(renderMesh->nEdges)
                                        {
                                                device->SetStreamSource(0, renderMesh->edgeVertexBuffer, 0, sizeof(D3DVERTEX));
                                                device->DrawPrimitive(D3DPT_LINELIST, 0, renderMesh->nEdges);
                                        }

                                        effect->EndPass();
                                }
                                effect->End();
                        }
                        device->EndScene();
                        device->SetRenderState(D3DRS_STENCILENABLE, false);
                        device->SetRenderState(D3DRS_CULLMODE,  D3DCULL_CCW);
                        device->SetRenderState(D3DRS_ALPHABLENDENABLE,  false);
                        device->SetRenderState(D3DRS_ZENABLE, D3DZB_TRUE);
                }
        }
}

void PathMapEffect::drawFullScreenQuad(LPDIRECT3DDEVICE9 device, float depth, float fLeftU, float fTopV, float fRightU, float fBottomV)
{
    D3DSURFACE_DESC dtdsdRT;
    PDIRECT3DSURFACE9 pSurfRT;

    // Acquire render target width and height
    device->GetRenderTarget(0, &pSurfRT);
    pSurfRT->GetDesc(&dtdsdRT);
    pSurfRT->Release();

    // Ensure that we're directly mapping texels to pixels by offset by 0.5
    // For more info see the doc page titled "Directly Mapping Texels to Pixels"
    FLOAT fWidth5 = (FLOAT)dtdsdRT.Width - 0.5f;
    FLOAT fHeight5 = (FLOAT)dtdsdRT.Height - 0.5f;

    // Draw the quad
    D3DVERTEX svQuad[4];

    svQuad[0].pos=D3DXVECTOR3(-1, 1, depth);
        svQuad[0].normal=D3DXVECTOR3(0,0,0);
    svQuad[0].tex0 = D3DXVECTOR2(fLeftU, fTopV);

    svQuad[1].pos = D3DXVECTOR3(1, 1, depth);
        svQuad[1].normal = D3DXVECTOR3(0,0,0);
    svQuad[1].tex0 = D3DXVECTOR2(fRightU, fTopV);

    svQuad[2].pos = D3DXVECTOR3(-1, -1, depth);
        svQuad[2].normal = D3DXVECTOR3(0,0,0);
    svQuad[2].tex0 = D3DXVECTOR2(fLeftU, fBottomV);

    svQuad[3].pos = D3DXVECTOR3(1, -1, depth);
        svQuad[3].normal = D3DXVECTOR3(0,0,0);
    svQuad[3].tex0 = D3DXVECTOR2(fRightU, fBottomV);
        
        device->SetRenderState(D3DRS_CULLMODE,D3DCULL_NONE);
        device->SetFVF(D3DFVF_XYZ | D3DFVF_NORMAL | D3DFVF_TEX1 | D3DFVF_TEXCOORDSIZE2(0));
    device->DrawPrimitiveUP(D3DPT_TRIANGLESTRIP, 2, svQuad, sizeof(D3DVERTEX));
}

void PathMapEffect::showPRMTexture()
{
        device->SetRenderTarget(0, frameColorBuffer);
        device->SetDepthStencilSurface(frameDepthStencilBuffer);

        device->Clear( 0, NULL, D3DCLEAR_TARGET | D3DCLEAR_ZBUFFER,     D3DCOLOR_XRGB(0,0,0), 1.0f, 0 );

        if( SUCCEEDED( device->BeginScene() ) ){
                if( effect != NULL ){
                        D3DXHANDLE hTechnique = effect->GetTechniqueByName((LPSTR)"ShowTex");
                        if(hTechnique==NULL){
                                return;
                        }
                        effect->SetTechnique( hTechnique );
                        UINT nPasses = 0;
                        effect->Begin( &nPasses, 0 );
                        for(UINT j=0;j<nPasses;j++){
                                effect->BeginPass(j);
                
                                //loop through entities
                                std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
                                //while(entityIterator != entities.end())
                                {
                                        //set entity params
                                        effect->SetMatrix("modelToWorldMatrix", &entityIterator->modelWorldTransform);
                                        effect->SetMatrix("inverseTransposedModelToWorldMatrix", &entityIterator->inverseTransposedModelWorldTransform);
                                        D3DXMATRIX t = entityIterator->modelWorldTransform;
                                        if(parameters->Get(bLookFromLight))
                                                effect->SetMatrix("modelToProjMatrix", &(entityIterator->modelWorldTransform * *lightCamera->GetViewMatrix()  * *lightCamera->GetProjMatrix()  ));
                                        else
                                                effect->SetMatrix("modelToProjMatrix", &(entityIterator->modelWorldTransform * *camera->GetViewMatrix()  * *camera->GetProjMatrix()  ));

                                        effect->SetTexture("filteredAtlas", entityIterator->prmTexture);

                                        unsigned int nSubsets = entityIterator->renderMesh->nSubsets;
                                        for(int i = 0; i< nSubsets; i++)
                                        {
                                                //set subset material texture
                                                effect->SetTexture("texToShow", entityIterator->prmTexture);
                                                effect->CommitChanges();
                                                drawFullScreenQuad(device);
//                                              effect->CommitChanges();
//                                              entityIterator->renderMesh->mesh->DrawSubset(i);
                                        }
                                        entityIterator++;
                                }
                                effect->EndPass();
                        }
                        effect->End();
                }
                device->EndScene();
        }

}

void PathMapEffect::RenderMesh::buildEdgeVertexBuffer(LPDIRECT3DDEVICE9 device)
{
        LPD3DXMESH pMesh = this->mesh;
        DWORD* pAdj = new DWORD[3 * pMesh->GetNumFaces()];
        for(int r=0; r<pMesh->GetNumFaces()*3; r++)
                pAdj[r] = 0xffff;
        pMesh->ConvertPointRepsToAdjacency(NULL, pAdj);

        LPDIRECT3DVERTEXBUFFER9 vertexBuffer;
        pMesh->GetVertexBuffer(&vertexBuffer);
        D3DVERTEX* vertexData;
        vertexBuffer->Lock(0,pMesh->GetNumVertices()*pMesh->GetNumBytesPerVertex(),(void**)&vertexData,0);
        LPDIRECT3DINDEXBUFFER9 indexBuffer;
        pMesh->GetIndexBuffer(&indexBuffer);
        unsigned short* indexData;
        indexBuffer->Lock(0,pMesh->GetNumFaces()*3*sizeof(unsigned short),(void**)&indexData,0);


        int countEdges = 0;
        for(int u=0; u<pMesh->GetNumFaces(); u++)
                for(int t=0; t<3; t++)
                        if(pAdj[u * 3 + t] >= pMesh->GetNumFaces())
                                countEdges++;

        if(countEdges == 0)
        {
                edgeVertexBuffer = NULL;
                nEdges = 0;
                vertexBuffer->Unlock();
                indexBuffer->Unlock();
                SAFE_RELEASE(vertexBuffer);
                SAFE_RELEASE(indexBuffer);
                return;
        }

        device->CreateVertexBuffer(2 * countEdges * pMesh->GetNumBytesPerVertex(),
                D3DUSAGE_WRITEONLY, pMesh->GetFVF(), D3DPOOL_DEFAULT, 
                &edgeVertexBuffer, NULL);
        D3DVERTEX* pevData;
        edgeVertexBuffer->Lock(0, 0, (void**)&pevData, 0);

        int iEdge = 0;
        for(int u=0; u<pMesh->GetNumFaces(); u++)
        {
                for(int t=0; t<3; t++)
                {
                        if(pAdj[u * 3 + t] >= pMesh->GetNumFaces())
                        {
                                unsigned short edge0 = indexData[u * 3 + t];
                                unsigned short edge1 = indexData[u * 3 + (t+1)%3];
                                unsigned short inner = indexData[u * 3 + (t+2)%3];

                                //push vertex edge0, edge1
                                pevData[iEdge * 2] = vertexData[edge0];
                                pevData[iEdge * 2 + 1] = vertexData[edge1];
                                
                                D3DXVECTOR2 ediff =  vertexData[edge1].tex0 - vertexData[edge0].tex0;
                                if(fabsf(ediff.x) > fabsf(ediff.y))
                                {
                                        if(vertexData[inner].tex0.y <= vertexData[edge0].tex0.y + 
                                                (vertexData[inner].tex0.x - vertexData[edge0].tex0.x) *
                                                (vertexData[edge1].tex0.y - vertexData[edge0].tex0.y) /
                                                (vertexData[edge1].tex0.x - vertexData[edge0].tex0.x))
                                        {
                                                pevData[iEdge * 2].tex0.y += 1.0/DBLATLASSIZE;
                                                pevData[iEdge * 2 + 1].tex0.y += 1.0/DBLATLASSIZE;
                                        }
                                        else
                                        {
                                                pevData[iEdge * 2].tex0.y -= 1.0/DBLATLASSIZE;
                                                pevData[iEdge * 2 + 1].tex0.y -= 1.0/DBLATLASSIZE;
                                        }
                                        D3DXVECTOR2 edgeDirUnitExt(1.0f, (vertexData[edge1].tex0.y - vertexData[edge0].tex0.y) /
                                                (vertexData[edge1].tex0.x - vertexData[edge0].tex0.x));
                                        edgeDirUnitExt *= 1.0/ATLASSIZE;
                                        if(vertexData[edge1].tex0.x > vertexData[edge0].tex0.x)
                                        {
                                                pevData[iEdge * 2].tex0 -= edgeDirUnitExt;
                                                pevData[iEdge * 2 + 1].tex0 += edgeDirUnitExt;
                                        }
                                        else
                                        {
                                                pevData[iEdge * 2].tex0 += edgeDirUnitExt;
                                                pevData[iEdge * 2 + 1].tex0 -= edgeDirUnitExt;
                                        }
                                }
                                else
                                {
                                        if(vertexData[inner].tex0.x <= vertexData[edge0].tex0.x + 
                                                (vertexData[inner].tex0.y - vertexData[edge0].tex0.y) *
                                                (vertexData[edge1].tex0.x - vertexData[edge0].tex0.x) /
                                                (vertexData[edge1].tex0.y - vertexData[edge0].tex0.y))
                                        {
                                                pevData[iEdge * 2].tex0.x += 1.0/DBLATLASSIZE;
                                                pevData[iEdge * 2 + 1].tex0.x += 1.0/DBLATLASSIZE;
                                        }
                                        else
                                        {
                                                pevData[iEdge * 2].tex0.x -= 1.0/DBLATLASSIZE;
                                                pevData[iEdge * 2 + 1].tex0.x -= 1.0/DBLATLASSIZE;
                                        }
                                        D3DXVECTOR2 edgeDirUnitExt((vertexData[edge1].tex0.x - vertexData[edge0].tex0.x) /
                                                (vertexData[edge1].tex0.y - vertexData[edge0].tex0.y), 1.0f);
                                        edgeDirUnitExt *= 1.0/ATLASSIZE;
                                        if(vertexData[edge1].tex0.y > vertexData[edge0].tex0.y)
                                        {
                                                pevData[iEdge * 2].tex0 -= edgeDirUnitExt;
                                                pevData[iEdge * 2 + 1].tex0 += edgeDirUnitExt;
                                        }
                                        else
                                        {
                                                pevData[iEdge * 2].tex0 += edgeDirUnitExt;
                                                pevData[iEdge * 2 + 1].tex0 -= edgeDirUnitExt;
                                        }
                                }

                                iEdge++;
                        }
                }
        }
        nEdges = iEdge;

        edgeVertexBuffer->Unlock();
        vertexBuffer->Unlock();
        indexBuffer->Unlock();
        SAFE_RELEASE(vertexBuffer);
        SAFE_RELEASE(indexBuffer);

        delete pAdj;
}

void PathMapEffect::uploadRadions()
{
        int nRadions = bushStarters.size();
        struct TMR{
                float pos[4];
                float dir[4];
//              float pow[4];
        };
        TMR* tm = new TMR[nRadions];
        for(int i=0; i<nRadions; i++)
        {
//*
                tm[i].pos[0] = bushStarters[i].position[0];
                tm[i].pos[1] = bushStarters[i].position[1];
                tm[i].pos[2] = bushStarters[i].position[2];
                tm[i].pos[3] = 1.0;
                tm[i].dir[0] = bushStarters[i].normal[0];
                tm[i].dir[1] = bushStarters[i].normal[1];
                tm[i].dir[2] = bushStarters[i].normal[2];
                tm[i].dir[3] = 1.0;
//              tm[i].pow[0] = bushStarters[i].radiance[0];
//              tm[i].pow[1] = bushStarters[i].radiance[1];
//              tm[i].pow[2] = bushStarters[i].radiance[2];
//              tm[i].pow[3] = 1.0;
/*/
                tm[i].pos[0] = 0.9;
                tm[i].pos[1] = 0.2;
                tm[i].pos[2] = 0.2;
                tm[i].pos[3] = 1.0;
                tm[i].dir[0] = 0.2;
                tm[i].dir[1] = 0.9;
                tm[i].dir[2] = 0.2;
                tm[i].dir[3] = 1.0;
                tm[i].pow[0] = 0.2;
                tm[i].pow[1] = 0.2;
                tm[i].pow[2] = 0.9;     //blue
                tm[i].pow[3] = 0.99;
//*/
        }

        RECT imgSize;
        imgSize.top = imgSize.left = 0;
        imgSize.bottom =  nRadions;
        imgSize.right = 2 * nRadions / 4096;
        D3DXLoadSurfaceFromMemory(radionSurface, NULL, NULL, tm, D3DFMT_A32B32G32R32F, 
                imgSize.right * 4 * sizeof(float), NULL, &imgSize, D3DX_FILTER_NONE, 0);

        delete [] tm;
}

void PathMapEffect::fillRadionPosArray(void* pData)
{
        int nRadions = NRADIONS;
        struct TMR{
                float pos[3];
                float rad[3];
        };
        TMR* tm = (TMR*)pData;
        for(int i=0; i<nRadions; i++)
        {
                tm[i].pos[0] = bushStarters[i].position[0];
                tm[i].pos[1] = bushStarters[i].position[1];
                tm[i].pos[2] = bushStarters[i].position[2];
                tm[i].rad[0] = bushStarters[i].radiance[0];
                tm[i].rad[1] = bushStarters[i].radiance[1];
                tm[i].rad[2] = bushStarters[i].radiance[2];
        }
}

const wchar_t* PathMapEffect::getWeightsString()
{
        static wchar_t ws[512];
        ws[0] = 0;
        char wan[64];
        float sumw = 0.0f;
        for(int i=0; i<33; i++)
        {
//              StrCpy(wan, L"%f ");
                if(i == 32)
                        sprintf(wan, "%.3f ", (double)sumw);
                else
                {
                        sumw += weights[i];
                        sprintf(wan, "%.3f ", (double)weights[i]);
                }
                wchar_t* wideValue = new wchar_t[MultiByteToWideChar(CP_ACP, 0, wan, -1, NULL, 0)];
                MultiByteToWideChar(CP_ACP, 0, wan, -1, wideValue, 511);
                StrCatW(ws, wideValue);
                delete wideValue;
        }
        return ws;
}

struct ClusterDist2{
        unsigned int index;
        double dist2;
        ClusterDist2(unsigned int index, double dist2) {this->index = index; this->dist2 = dist2;}
        bool operator<(const ClusterDist2& o) const {return dist2 > o.dist2;}
        bool operator>(const ClusterDist2& o) const {return dist2 < o.dist2;}
        bool operator<=(const ClusterDist2& o) const {return dist2 >= o.dist2;}
        bool operator>=(const ClusterDist2& o) const {return dist2 <= o.dist2;}
};

void PathMapEffect::Entity::findNearClusters(const std::vector<Radion>& starters, unsigned int nClusters)
{
        D3DXVECTOR3 cc = (D3DXVECTOR3&)rayTraceEntity->bbox.getCentre();

        std::priority_queue<ClusterDist2> nearestClusterDists;

        for(unsigned int iCluster = 0; iCluster < NCLUSTERS; iCluster++)
        {
                unsigned int iRadion = 0;
                unsigned int nRadionsPerCluster = owner->clusterLenghts[iCluster];
                double aClusterDist2 = 0.0;
                for(unsigned int iric=0; iric < nRadionsPerCluster; iric++, iRadion++)
                {
                        aClusterDist2 += (starters[iRadion++].position - cc).norm2();
                }
                nearestClusterDists.push( ClusterDist2(iCluster, aClusterDist2) );
        }
        for(unsigned int i=0; i<32; i++)
        {
                nearClusterIndices[i] = i;// nearestClusterDists.top().index;
                nearestClusterDists.pop();
        }
}



void PathMapEffect::savePathMaps()
{
        std::fstream entryPointFile("prm\\prmEntryPoints.dat", std::ios::out | std::ios::binary);

        for(unsigned int u=0; u < NRADIONS; u++)
        {
                entryPointFile << bushStarters.at(u);
        }
        
        for(unsigned int c=0; c < NCLUSTERS; c++)
        {
                entryPointFile.write((char*)&clusterLenghts[c], sizeof(unsigned int));
        }

        entryPointFile.flush();
        entryPointFile.close();

        unsigned int iEntity = 0;
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                wchar_t prmFileName[256];
                wsprintf(prmFileName, L"prm\\prm_%08d.hdr", iEntity);
                D3DXSaveSurfaceToFile(prmFileName, D3DXIFF_HDR, entityIterator->prmSurface, NULL, NULL);
                entityIterator++;
                iEntity++;
        }
}

void PathMapEffect::loadPathMaps()
{
        std::fstream entryPointFile("prm\\prmEntryPoints.dat", std::ios::in | std::ios::binary);

        bushStarters.clear();

        for(unsigned int u=0; u < NRADIONS; u++)
        {
                Radion ri;
                entryPointFile >> ri;
                bushStarters.push_back(ri);
        }
        
        for(unsigned int c=0; c < NCLUSTERS; c++)
        {
                entryPointFile.read((char*)&clusterLenghts[c], sizeof(unsigned int));           
        }

        entryPointFile.close();

        unsigned int iEntity = 0;
        std::vector<PathMapEffect::Entity>::iterator entityIterator = entities.begin();
        while(entityIterator != entities.end())
        {
                entityIterator->findNearClusters(bushStarters, NCLUSTERS);

                wchar_t prmFileName[256];
                wsprintf(prmFileName, L"prm\\prm_%08d.hdr", iEntity);
                D3DXLoadSurfaceFromFile(entityIterator->prmSurface, NULL, NULL, prmFileName, NULL, D3DX_DEFAULT, 0, NULL);
                entityIterator++;
                iEntity++;
        }
}