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>

#include <GLGraphics/ManifoldRenderer.h>

#include <GLGraphics/Console.h>

#define USE_SPARSE_MATRIX 0

class Harmonics

{

	HMesh::Manifold& mani;

    HMesh::VertexAttributeVector<int> vtouched;

	int maximum_eigenvalue;

    bool is_initialized;

    GLuint prog_P0;

	std::vector<CGLA::Vec3d> proj;

	LinAlg::CMatrix Q;

	LinAlg::CVector qnorm;

	LinAlg::CVector V;

	LinAlg::CVector S;

    GLGraphics::Console::variable<float> display_harmonics_time;

    GLGraphics::Console::variable<int> display_harmonics_diffuse;

    GLGraphics::Console::variable<int> display_harmonics_highlight;

	std::vector<float> max_eig_values;

	void make_laplace_operator();

	void make_laplace_operator_sparse();

public:

	/// Initial analysis of harmonics

	Harmonics(HMesh::Manifold& mani, GLGraphics::Console& cs);

	/// 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 GLGraphics::ManifoldRenderer

{

public:

    HarmonicsRenderer(Harmonics* h)

    {

        glNewList(display_list,GL_COMPILE);

        if(h) h->draw_adf();

        glEndList();

    }

};

#endif