#ifndef EquaAngleEnergy_H
#define EquaAngleEnergy_H

/* Last edited on DATE TIME by USER */

#include <Energy.h>

typedef struct EquaAngleEnergy_data_t { } EquaAngleEnergy_data_t;
  /* The {Equalization Angle} energy measures the non-uniformity of diedral
    angles among consecutive walls incident to an internal @{edge->?} {e}. We say
    that an @{edge->?} {e} is internal if every wall incident to {e} is an inter-
    nal wall. Analosgously, we say that a wall {f} is internal if there exists
    exactely two tetrahedra incident to it. 

    This energy is computed only for @{edge->?}s that exist, are internal, have wall
    ring degree  (DRF) three or more, and have only existing walls and knots
    incident to it.

    The diedral angles not need projected. For two walls consecutives incidents
    to  same @{edge->?}, the ideal angle {I} is {2*Pi/n} where {n} is the wall ring
    degree correspondly to @{edge->?} {e}.

    Minimizing this energy tends to equalize the angles between consecutive
    walls incident to @{edge->?} {e}.

                        
                      f2  \    /f1
                           \  /
                            \/
                          e /\
                           /  \
                          /    \
                      f3         f4
                          
  */
 
Energy_t EquaAngleEnergy_new(void);
 
#endif