CDensities.cpp

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <glut.h>
#include "lifevariables.h"
#include "CDensities.h"

// CDensities uses CVolVoxels
void CDensities::InitDensities_XY_plain_by_Slice(int iSliceNum)
{
    unsigned short * pD   = _pDensities;
    unsigned short dimX, dimY;
    int i,j;
  
    dimX = GVolVoxels.GetDimX();
    dimY = GVolVoxels.GetDimY();

    for ( j = 0; j < dimY; j++ )
        for ( i = 0; i < dimX; i++, pD++ )
        {
            (*pD) = GVolVoxels.GetVoxelValue(i,j,iSliceNum); 
        }
}

// CDensities uses CVolVoxels
void CDensities::InitDensities_XZ_plain_by_Slice(int iSliceNum)
{
    unsigned short * pD   = _pDensities;
    unsigned short dimX, dimZ;
        int i,j;
  
    dimX = GVolVoxels.GetDimX();
    dimZ = GVolVoxels.GetDimZ();

    for ( j = 0; j < dimZ; j++ ) 
    {
        for ( i = 0; i < dimX; i++, pD++ )
        {
            (*pD) = GVolVoxels.GetVoxelValue(i,iSliceNum,j);  
        }
    }
}


void CDensities::InitDensities_YZ_plain_by_Slice(int iSliceNum)
{
    unsigned short * pD   = _pDensities;
    unsigned short dimY, dimZ;
        int i,j;

    dimY = GVolVoxels.GetDimY();
    dimZ = GVolVoxels.GetDimZ();

    for ( j = 0; j < dimZ; j++ )
        for ( i = 0; i < dimY; i++, pD++ )
        {
            (*pD) = GVolVoxels.GetVoxelValue(iSliceNum,i,j);
        }
}

unsigned short CDensities::GetMaxValue()
{
    unsigned uNewVal = 0;
    unsigned uMaxCount = _dimX * _dimY;


    for ( int i = 0; i < uMaxCount; i++ )
    {
        if ( uNewVal < _pDensities[i] )
        {
             uNewVal = _pDensities[i];
        }
    }

    return uNewVal;
}

---