Subversion Repositories gelsvn

/*
 *  harmonics.h
 *  GEL
 *
 *  Created by J. Andreas Bærentzen on 01/09/08.
 *  Copyright 2008 __MyCompanyName__. All rights reserved.
 *
 */

#ifndef __MESHEDIT_HARMONICS_H__
#define __MESHEDIT_HARMONICS_H__

#include <CGLA/Vec3d.h>
#include <HMesh/Manifold.h>
#include <LinAlg/Matrix.h>
#include <LinAlg/Vector.h>
#include "Renderer.h"

#define USE_SPARSE_MATRIX 1


class Harmonics
{
        HMesh::Manifold& mani;
        
        int maximum_eigenvalue;

        static bool is_initialized;
        static GLuint prog_P0;
        
        std::vector<CGLA::Vec3d> proj;
        
        LinAlg::CMatrix Q;
        LinAlg::CVector qnorm;
        LinAlg::CVector V;
        LinAlg::CVector S;
        std::vector<float> max_eig_values;

        void make_laplace_operator();
        void make_laplace_operator_sparse();

public:
        
        /// Initialize API (must be called first)
        static void init();
        
        /// Initial analysis of harmonics (must be called second)
        Harmonics(HMesh::Manifold& mani);
        
        /// Add a frequency to mesh reconstruction
        void add_frequency(int f, float scale = 1.0f);
        
        /// Reset the shape to use 0 eigenvalues
        void reset_shape();
        
        /// Do a partial reconstruct with an interval of eigenvalues
        void partial_reconstruct(int E0, int E1, float scale=1.0f);
        
        /// Parse keystrokes that would influence the interactive display
        void parse_key(unsigned char key);
        
        /// Draw with eigenvalues
        void draw_adf();
        
};


class HarmonicsRenderer: public ManifoldRenderer
        {
                
        public:
                HarmonicsRenderer(Harmonics* h)
                {
                        glNewList(display_list,GL_COMPILE);
                        if(h) h->draw_adf();
                        glEndList();
                }
        };


#endif