CelluarAutomataGenerator.h

#pragma once
#include "Scene.h"
#include "CloudGenerator.h"
#include "CloudScene.h"

#include <stdlib.h>
#include <string>
#include <stdio.h>
#include <iostream>
#include <cstdlib>
#include <math.h>

namespace osgCloudyDay
{
        class CelluarAutomataGenerator : public osgCloudyDay::CloudGenerator
        {               
        public:
                CelluarAutomataGenerator(unsigned int width, unsigned int depth, unsigned int height, int type, osg::Vec4 color);
                ~CelluarAutomataGenerator(void);

                void Perform(void);
                void Update(void);

        protected:
                void Initialize();

                int m_type;

                void AddEllipsoide(osg::Vec3 radius, osg::Vec3 position, float lifetime);

                float GetCellDensity( int i, int j, int k );
                float GetPointDensity(osg::Vec3 Point); 

        private:
                void DefineInitValues();
                void Interpolate(float& t);

                inline int Index(int i, int j, int k);
                inline bool IsCellInSpace(int x, int y, int z);

                void UpdateEllisoides();
                bool IsOutside(osg::Vec3 pos);
                void Grow();
                void Extrinction();
                void InitCell( int x, int y, int z, float fProbSeed );
                void ShapeVolume();
                void SetupArrays();
                void CalcDensity();
                void SupplyVapor();
                void AdvectionWind();

                float* m_current_density;
                float** m_density2;

                unsigned int** m_hum;
                unsigned int** m_act;
                unsigned int** m_cld;
                unsigned int* m_fact;

                float* m_phum;
                float* m_pext;
                float* m_pact;  

                int m_depth;
                int m_width;
                int m_height;
                float f_depth;
                float f_width;
                float f_height;

                unsigned int t;
                unsigned int tdens;

                unsigned int m_elapsedSteps;

                osg::Vec3 bb_min;
                osg::Vec3 bb_max;

                osg::Vec3 windVector;
        

                osg::ref_ptr<osg::Vec3Array> m_positions;
                osg::ref_ptr<osg::Vec3Array> m_radius;
                osg::ref_ptr<osg::UIntArray> m_lifetimes;

                float m_time;

                osg::Vec3* randomPos;
                osg::Vec3 windDirection;
        };
}