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 "../HMesh/AttributeVector.h"
#include "../LinAlg/Matrix.h"
#include "../LinAlg/Vector.h"

namespace  HMesh {
    
    class Harmonics
    {
        HMesh::Manifold& mani;
        HMesh::VertexAttributeVector<int> vtouched;
        
        int maximum_eigenvalue;
        
        bool is_initialized;
        
        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:
        
        /// Initial analysis of harmonics
        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);
        
        double compute_adf(HMesh::VertexAttributeVector<double>& adf, double t);
        
    };
    
}
#endif