Commit: 1fff391816c6aa58b51693ce03495dc49e7b13f2 - fermisurfer (git)

fermisurfer Git

Commit MetaInfo

• Revision: 1fff391816c6aa58b51693ce03495dc49e7b13f2 (tree)
• Zeit: 2020-03-13 01:08:42
• Autor: Mitsuaki Kawamura <kawamitsuaki@gmai...>
• Commiter: Mitsuaki Kawamura

Log Message

Use global variables again.
Output 2D Fermi lines
Display reciplocal lattice vectors
Universal design (color: RGB, CMY, MCY)

Ändern Zusammenfassung

• modified: src/bz_lines.cpp
• modified: src/bz_lines.hpp
• modified: src/calc_nodeline.cpp
• modified: src/calc_nodeline.hpp
• modified: src/draw.cpp
• modified: src/equator.cpp
• modified: src/equator.hpp
• modified: src/fermi_patch.cpp
• modified: src/fermi_patch.hpp
• modified: src/fermisurfer.cpp
• modified: src/free_patch.cpp
• modified: src/free_patch.hpp
• modified: src/initialize.cpp
• modified: src/initialize.hpp
• modified: src/kumo.cpp
• modified: src/kumo.hpp
• modified: src/menu.cpp
• modified: src/menu.hpp
• modified: src/read_file.cpp
• modified: src/read_file.hpp
• modified: src/section.cpp
• modified: src/section.hpp
• modified: src/variable.hpp

Diff

--- a/src/bz_lines.cpp
+++ b/src/bz_lines.cpp
@@ -25,6 +25,7 @@ THE SOFTWARE.
  @brief Compute lines of BZ boundary
 */
 #include "basic_math.hpp"
+#include "variable.hpp"
 #include <wx/wx.h>
 #if defined(HAVE_CONFIG_H)
 #include <config.h>
@@ -107,14 +108,8 @@ static int bragg_vert(
 /**
  @brief Compute real number of Bragg plane at 1st BZ
 */
-static void check_bragg(
-  wxTextCtrl *terminal,
-  int nbzl,
-  GLfloat bzl[676][2][3],
-  GLfloat bragg[26][3],
-  GLfloat brnrm[26],
-  int *nbragg
-) {
+static void check_bragg()
+{
   int ibr, ibzl, ibzc;
   int ii, kk, bzflag, nbzcorner, nn;
   GLfloat thr = (GLfloat)0.0001, prod, bzc[676][3];
@@ -144,7 +139,7 @@ static void check_bragg(
   Then, compute real number Bragg plane of 1st BZ (::nbragg), 
   Re-order ::bragg and ::brnrm
   */
-  *nbragg = 0;
+  nbragg = 0;
   for (ibr = 0; ibr < 26; ibr++) {
     nn = 0;
 
@@ -153,31 +148,25 @@ static void check_bragg(
       if (fabsf(prod - brnrm[ibr]) < thr) nn += 1;
     }
     if (nn >= 3) {
-      for (kk = 0; kk < 3; kk++) bragg[*nbragg][kk] = bragg[ibr][kk];
-      brnrm[*nbragg] = brnrm[ibr];
-      *nbragg += 1;
+      for (kk = 0; kk < 3; kk++) bragg[nbragg][kk] = bragg[ibr][kk];
+      brnrm[nbragg] = brnrm[ibr];
+      nbragg += 1;
     }
   }
-  *terminal << wxString::Format(wxT("    Number of plane of 1st BZ : %d\n"), *nbragg);
+  *terminal << wxString::Format(wxT("    Number of plane of 1st BZ : %d\n"), nbragg);
 }/*static void check_bragg*/
 /**
  @brief Compute Brillouin zone boundariy lines
 
  Modify : ::nbzl, ::bzl
 */
-void bz_lines(
-  wxTextCtrl* terminal,
-  GLfloat bragg[26][3],
-  GLfloat brnrm[26],
-  int *nbzl,
-  GLfloat bzl[676][2][3],
-  int *nbragg
-) {
+void bz_lines()
+{
   /**/
   int ibr, jbr, nbr, i, j, lvert;
   GLfloat vert[2][3];
   /**/
-  *nbzl = 0;
+  nbzl = 0;
   /**/
   for (ibr = 0; ibr < 26; ++ibr) {
     for (jbr = 0; jbr < 26; ++jbr) {
@@ -191,12 +180,12 @@ void bz_lines(
       lvert = bragg_vert(bragg, brnrm, ibr, jbr, nbr, vert[1], vert[0]);
       if (lvert == 0) continue;
       /**/
-      for (i = 0; i < 2; ++i) for (j = 0; j < 3; ++j) bzl[*nbzl][i][j] = vert[i][j];
-      *nbzl = *nbzl + 1;
+      for (i = 0; i < 2; ++i) for (j = 0; j < 3; ++j) bzl[nbzl][i][j] = vert[i][j];
+      nbzl = nbzl + 1;
 
     }/*for (jbr = 0; jbr < 26; ++jbr)*/
   }/*for (ibr = 0; ibr < 26; ++ibr)*/
 
-  check_bragg(terminal, *nbzl, bzl, bragg, brnrm, nbragg);
+  check_bragg();
 
 }/*bz_lines*/

--- a/src/bz_lines.hpp
+++ b/src/bz_lines.hpp
@@ -22,4 +22,4 @@ OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
 
-void bz_lines(wxTextCtrl* terminal, GLfloat bragg[26][3], GLfloat brnrm[26], int* nbzl, GLfloat bzl[676][2][3], int* nbragg);
+void bz_lines();

--- a/src/calc_nodeline.cpp
+++ b/src/calc_nodeline.cpp
@@ -38,24 +38,15 @@ THE SOFTWARE.
 #include <cmath>
 #include <vector>
 #include "basic_math.hpp"
-
+#include "variable.hpp"
 /**
 @brief Compute node-line where \f$\Delta_{n k} = 0\f$
 
  Modify : ::nnl, ::kvnl, ::kvnl_rot
 
 */
-void calc_nodeline(
-int nb,
-int *ntri,
-GLfloat ****matp,
-GLfloat ****kvp,
-int nthreads,
-wxTextCtrl *terminal,
-int* nnl,
-GLfloat**** kvnl,
-GLfloat** kvnl_rot
-) {
+void calc_nodeline()
+{
   int ib, itri, i, j, ithread, nnl0;
   std::vector<std::vector<std::vector<std::vector<GLfloat>>>> kvnl_v;
 

--- a/src/calc_nodeline.hpp
+++ b/src/calc_nodeline.hpp
@@ -21,15 +21,4 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include <wx/wx.h>
-
-void calc_nodeline(int nb, int* ntri, GLfloat**** matp, GLfloat**** kvp, int nthreads,
-  wxTextCtrl* terminal, int* nnl, GLfloat**** kvnl, GLfloat** kvnl_rot);
+void calc_nodeline();

--- a/src/draw.cpp
+++ b/src/draw.cpp
@@ -53,8 +53,8 @@ static void draw_fermi() {
    First, rotate k-vector and normal vector
   */
 #pragma omp parallel default(none) \
-  shared(nb,draw_band,ntri,rot,nmlp,arw,nmlp_rot,kvp,kvp_rot,arw_rot,trans,side) \
-  private(ib)
+shared(nb,draw_band,ntri,rot,nmlp,arw,nmlp_rot,kvp,kvp_rot,arw_rot,trans,side) \
+private(ib)
   {
     int i, j, l, itri;
 
@@ -104,12 +104,11 @@ static void draw_fermi() {
   /*
   Arrow for 3D value
   */
-  if (blackback == 1) glColor3fv(white);
-  else glColor3fv(black);
   glNormal3f(0.0f, 0.0f, 1.0f);
   if (color_scale == 3) {
     for (ib = 0; ib < nb; ib++) {
       if (draw_band[ib] == 1) {
+        glColor3f(1.0 - clr[ib][0], 1.0 - clr[ib][1], 1.0 - clr[ib][2]);
         glVertexPointer(3, GL_FLOAT, 0, arw_rot[ib]);
         glDrawArrays(GL_LINES, 0, ntri[ib] * 3 * 2);
       }
@@ -214,10 +213,7 @@ static void draw_bz_lines() {
   /*
    Line color is oposit of BG color
   */
-  if (blackback == 1) 
-    for (i = 0; i<4; i++) linecolor[i] = white[i];
-  else
-    for (i = 0; i<4; i++) linecolor[i] = black[i];
+  for (i = 0; i < 4; i++) linecolor[i] = LineColor[i];
   /**/
   glLineWidth(linewidth);
   /*
@@ -301,20 +297,45 @@ static void draw_bz_lines() {
 static void draw_colorbar()
 {
   int i, j, k;
-  GLfloat mat2, vertices[366], colors[488];
-  /**/
+  GLfloat mat2, vertices[366], colors[488], vector[18], vector_color[24],
+    norm;
+
   glEnableClientState(GL_COLOR_ARRAY);
+  /*
+  Reciplocal lattice vectors
+  */
+  for (i = 0; i < 6; i++) {
+    vector[3 * i] = -1.2f;
+    vector[3 * i + 1] = -1.05f;
+    vector[3 * i + 2] = 0.0f;
+  }
+  for (i = 0; i < 3; i++) {
+    norm = sqrtf(bvec[i][0] * bvec[i][0]+ bvec[i][1] * bvec[i][1]+ bvec[i][2] * bvec[i][2]);
+    for (j = 0; j < 3; j++) {
+      vector[j + 6 * i] 
+        += (rot[j][0] * bvec[i][0]
+        + rot[j][1] * bvec[i][1]
+        + rot[j][2] * bvec[i][2]) * 0.2f / norm;
+    }
+    for (j = 0; j < 4; j++) {
+      vector_color[j + 8 * i] = BarColor[i * 2][j];
+      vector_color[j + 8 * i + 4] = BarColor[i * 2][j];
+    }
+  }
+  glNormal3f(0.0f, 0.0f, 1.0f);
+  glVertexPointer(3, GL_FLOAT, 0, vector);
+  glColorPointer(4, GL_FLOAT, 0, vector_color);
+  glDrawArrays(GL_LINES, 0, 6);
+  /*
+   Color bar/circle/cube
+  */
   if (color_scale == 1 || color_scale == 4) {
     for (i = 0; i < 5; i++) {
       for (j = 0; j < 2; j++) {
         vertices[0 + j * 3 + i * 6] = -1.0f + 0.5f*(GLfloat)i;
         vertices[1 + j * 3 + i * 6] = -1.0f - 0.1f*(GLfloat)j;
         vertices[2 + j * 3 + i * 6] = 0.0f;
-        if (     i == 0) for (k = 0; k < 4; k++) colors[k + 4 * j + 8 * i] = blue[k];
-        else if (i == 1) for (k = 0; k < 4; k++) colors[k + 4 * j + 8 * i] = cyan[k];
-        else if (i == 2) for (k = 0; k < 4; k++) colors[k + 4 * j + 8 * i] = green[k];
-        else if (i == 3) for (k = 0; k < 4; k++) colors[k + 4 * j + 8 * i] = yellow[k];
-        else if (i == 4) for (k = 0; k < 4; k++) colors[k + 4 * j + 8 * i] = red[k];
+        for (k = 0; k < 4; k++) colors[k + 4 * j + 8 * i] = BarColor[i][k];
       }
     }/*for (i = 0; i < 10; i++)*/
     glNormal3f(0.0f, 0.0f, 1.0f);
@@ -329,10 +350,7 @@ static void draw_colorbar()
     vertices[0] = 0.0f;
     vertices[1] = -1.0f;
     vertices[2] = 0.0f;
-    if (blackback == 1)
-      for (j = 0; j < 4; j++) colors[j] = wgray[j];
-    else
-      for (j = 0; j < 4; j++) colors[j] = bgray[j];
+    for (j = 0; j < 4; j++) colors[j] = 1.0f - BackGroundColor[j];
     /**/
     for (i = 0; i <= 60; i++) {
       /**/
@@ -408,8 +426,7 @@ static void draw_circles(
     vertices[2 + (i + 1) * 3] = 0.0f;
   }/*for (i = 0; i <= 20; i++)*/
   glNormal3f(0.0f, 0.0f, 1.0f);
-  if (blackback == 1) glColor3fv(white);
-  else glColor3fv(black);
+  glColor3fv(LineColor);
   glVertexPointer(3, GL_FLOAT, 0, vertices);
   glDrawArrays(GL_TRIANGLE_FAN, 0, 22);
   /**/
@@ -431,8 +448,7 @@ static void draw_fermi_line() {
     for (i = 0; i < 3; i++) vertices[i + 3 * ibzl] = bzl2d_proj[ibzl][i];
   }/*for (ibzl = 0; ibzl < nbzl2d; ++ibzl)*/
   glNormal3f(0.0f, 0.0f, 1.0f);
-  if (blackback == 1)glColor3fv(white);
-  else glColor3fv(black);
+  glColor3fv(LineColor);
   glVertexPointer(3, GL_FLOAT, 0, vertices);
   glDrawArrays(GL_LINE_LOOP, 0, nbzl2d);
   /*

--- a/src/equator.cpp
+++ b/src/equator.cpp
@@ -38,23 +38,13 @@ THE SOFTWARE.
 #include <cstdlib>
 #include <cmath>
 #include "basic_math.hpp"
+#include "variable.hpp"
 /**
 @brief Compute equator \f$\{\bf v}_{n k} \cdot {\bf k} = 0\f$
 
  Modify : ::nequator, ::kveq, ::kveq_rot
 */
-void equator(
-int nb,
-int nthreads,
-GLfloat eqvec[3],
-int *ntri,
-GLfloat ****kvp,
-GLfloat ****nmlp,
-wxTextCtrl* terminal,
-int *nequator,
-GLfloat ****kveq,
-GLfloat **kveq_rot
-) {
+void equator() {
   int ib, itri, i, j, ithread, nequator0;
   std::vector<std::vector<std::vector<std::vector<GLfloat>>>> kveq_v;
  

--- a/src/equator.hpp
+++ b/src/equator.hpp
@@ -21,14 +21,4 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include <wx/wx.h>
-void equator(int nb,  int nthreads,  GLfloat eqvec[3],  int* ntri,  GLfloat**** kvp,
-  GLfloat**** nmlp, wxTextCtrl* terminal, int* nequator, GLfloat**** kveq, GLfloat** kveq_rot);
+void equator();

--- a/src/fermi_patch.cpp
+++ b/src/fermi_patch.cpp
@@ -38,6 +38,7 @@ THE SOFTWARE.
 #include <cstdio>
 #include <cmath>
 #include "basic_math.hpp"
+#include "variable.hpp"
 /**
  @brief Store triangle patch
 
@@ -77,11 +78,6 @@ static void triangle(
   GLfloat mat1[3][3], //!<[in] The matrix element
   GLfloat kvec1[3][3], //!<[in] @f$k@f$-vector of corners
   GLfloat vf1[3][3], //!<[in] @f$v_f@f$-vector of corners
-  int fbz,
-  int nbragg,
-  GLfloat bragg[26][3],
-  GLfloat brnrm[26],
-  GLfloat brnrm_min,
   std::vector<std::vector<std::vector<GLfloat> > > &kvp_v,
   std::vector<std::vector<std::vector<GLfloat> > > &matp_v,
   std::vector<std::vector<std::vector<GLfloat> > > &nmlp_v
@@ -147,8 +143,7 @@ static void triangle(
             vf2[i][j] = vf1[sw[i]][j];
           }
         }
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
         return;
       }
       else if (brnrm[ibr] < prod[sw[1]]) {
@@ -168,8 +163,7 @@ static void triangle(
           vf2[1][i] = vf1[sw[0]][i] * a[0][1] + vf1[sw[1]][i] * a[1][0];
           vf2[2][i] = vf1[sw[0]][i] * a[0][2] + vf1[sw[2]][i] * a[2][0];
         }/*for (i = 0; i < 3; ++i)*/
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
 
         for (i = 0; i < 3; ++i) {
           kvec2[0][i] = kvec1[sw[0]][i] * a[0][1] + kvec1[sw[1]][i] * a[1][0];
@@ -186,8 +180,7 @@ static void triangle(
         }/*for (i = 0; i < 3; ++i)*/
         for (i = 0; i < 3; ++i) for (j = 0; j < 3; ++j)
           kvec2[i][j] += bshift * bragg[ibr][j];
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
 
         for (i = 0; i < 3; ++i) {
           kvec2[0][i] = kvec1[sw[2]][i];
@@ -204,8 +197,7 @@ static void triangle(
         }/*for (i = 0; i < 3; ++i)*/
         for (i = 0; i < 3; ++i) for (j = 0; j < 3; ++j)
           kvec2[i][j] += bshift * bragg[ibr][j];
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
         return;
       }
       else if (brnrm[ibr] < prod[sw[2]]) {
@@ -227,8 +219,7 @@ static void triangle(
         }/*for (i = 0; i < 3; ++i)*/
         for (i = 0; i < 3; ++i) for (j = 0; j < 3; ++j)
           kvec2[i][j] += bshift * bragg[ibr][j];
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
 
         for (i = 0; i < 3; ++i) {
           kvec2[0][i] = kvec1[sw[0]][i];
@@ -243,8 +234,7 @@ static void triangle(
           vf2[1][i] = vf1[sw[1]][i];
           vf2[2][i] = vf1[sw[0]][i] * a[0][2] + vf1[sw[2]][i] * a[2][0];
         }/*for (i = 0; i < 3; ++i)*/
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
 
         for (i = 0; i < 3; ++i) {
           kvec2[0][i] = kvec1[sw[1]][i] * a[1][2] + kvec1[sw[2]][i] * a[2][1];
@@ -259,8 +249,7 @@ static void triangle(
           vf2[1][i] = vf1[sw[1]][i];
           vf2[2][i] = vf1[sw[0]][i] * a[0][2] + vf1[sw[2]][i] * a[2][0];
         }/*for (i = 0; i < 3; ++i)*/
-        triangle(ibr + 1, mat2, kvec2, vf2, fbz, nbragg,
-          bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+        triangle(ibr + 1, mat2, kvec2, vf2, kvp_v, matp_v, nmlp_v);
         return;
       }
       else {
@@ -332,13 +321,6 @@ static void tetrahedron(
   GLfloat mat1[8][3], //!< [in] Matrix elements @f$\Delta_{n k}@f$
   GLfloat kvec1[8][3], //!< [in] @f$k@f$-vectors
   GLfloat vf1[8][3], //!< [in] @f$v_f@f$-vectors
-  int corner[6][4],
-  GLfloat bvec[3][3],
-  int fbz,
-  int nbragg,
-  GLfloat bragg[26][3],
-  GLfloat brnrm[26],
-  GLfloat brnrm_min,
   std::vector<std::vector<std::vector<GLfloat> > > &kvp_v,
   std::vector<std::vector<std::vector<GLfloat> > > &matp_v,
   std::vector<std::vector<std::vector<GLfloat> > > &nmlp_v
@@ -392,8 +374,7 @@ static void tetrahedron(
       }
       
       vol = a[1][0] * a[2][0] * a[3][0];
-      triangle(0, mat3, kvec3, vf3, fbz, nbragg,
-        bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+      triangle(0, mat3, kvec3, vf3, kvp_v, matp_v, nmlp_v);
     }
     else if (eig2[sw[1]] <= 0.0 && 0.0 < eig2[sw[2]]) {
       for (i = 0; i < 3; ++i) {
@@ -411,8 +392,7 @@ static void tetrahedron(
       }
       
       vol = a[1][2] * a[2][0] * a[3][0];
-      triangle(0, mat3, kvec3, vf3, fbz, nbragg,
-        bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+      triangle(0, mat3, kvec3, vf3, kvp_v, matp_v, nmlp_v);
       /**/
       for (i = 0; i < 3; ++i) {
         kvec3[0][i] = kvec2[sw[1]][i] * a[1][3] + kvec2[sw[3]][i] * a[3][1];
@@ -429,8 +409,7 @@ static void tetrahedron(
       }
 
       vol = a[1][3] * a[3][0] * a[2][1];
-      triangle(0, mat3, kvec3, vf3, fbz, nbragg,
-        bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+      triangle(0, mat3, kvec3, vf3, kvp_v, matp_v, nmlp_v);
     }
     else if (eig2[sw[2]] <= 0.0 && 0.0 < eig2[sw[3]]) {
       for (i = 0; i < 3; ++i) {
@@ -448,8 +427,7 @@ static void tetrahedron(
       }
 
       vol = a[0][3] * a[1][3] * a[2][3];
-      triangle(0, mat3, kvec3, vf3, fbz, nbragg,
-        bragg, brnrm, brnrm_min, kvp_v, matp_v, nmlp_v);
+      triangle(0, mat3, kvec3, vf3, kvp_v, matp_v, nmlp_v);
     }
     else {
     }
@@ -460,34 +438,7 @@ static void tetrahedron(
 
  Modify : ::ntri, nmlp, ::matp, ::kvp, ::clr, ::nmlp_rot, ::kvp_rot
 */
-void fermi_patch(
-  int nb,
-  int nthreads,
-  int ng[3],
-  int ng0[3],
-  int shiftk[3],
-  int fbz,
-  GLfloat ****eig,
-  GLfloat EF,
-  GLfloat *****mat,
-  GLfloat *****vf,
-  int corner[6][4],
-  GLfloat bvec[3][3],
-  int nbragg,
-  GLfloat bragg[26][3],
-  GLfloat brnrm[26],
-  GLfloat brnrm_min,
-  wxTextCtrl* terminal,
-  int *ntri,
-  GLfloat ****kvp,
-  GLfloat ****matp,
-  GLfloat **clr,
-  GLfloat *****arw,
-  GLfloat ****nmlp,
-  GLfloat **kvp_rot,
-  GLfloat **nmlp_rot,
-  GLfloat **arw_rot
-  )
+void fermi_patch()
 {
   int ntri0, ib, i0, i1, j0, start[3], last[3];
   int ithread;
@@ -610,8 +561,7 @@ private(j0,i0,i1,ithread)
               vf1[7][j] = vf[ib][ii0][ii1][ii2][j];
             }/*for (j = 0; j < 3; j++)*/
             /**/
-            tetrahedron(eig1, mat1, kvec1, vf1, corner,
-              bvec, fbz, nbragg, bragg, brnrm, brnrm_min, 
+            tetrahedron(eig1, mat1, kvec1, vf1, 
               kvp_v.at(ithread), matp_v.at(ithread), nmlp_v.at(ithread));
           }/*for (j0 = start[0]; j0 < ng[0]; ++j0)*/
         }/*for (j1 = start[1]; j1 < ng[1]; ++j1)*/

--- a/src/fermi_patch.hpp
+++ b/src/fermi_patch.hpp
@@ -21,17 +21,4 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include <wx/wx.h>
-
-void fermi_patch(int nb,  int nthreads,  int ng[3],  int ng0[3],  int shiftk[3],
-  int fbz, GLfloat**** eig, GLfloat EF, GLfloat***** mat, GLfloat***** vf, int corner[6][4], GLfloat bvec[3][3], int nbragg,
-  GLfloat bragg[26][3], GLfloat brnrm[26], GLfloat brnrm_min, wxTextCtrl* terminal, int* ntri, GLfloat**** kvp, GLfloat**** matp,
-  GLfloat** clr, GLfloat***** arw, GLfloat**** nmlp, GLfloat** kvp_rot, GLfloat** nmlp_rot, GLfloat** arw_rot);
+void fermi_patch();

--- a/src/fermisurfer.cpp
+++ b/src/fermisurfer.cpp
@@ -112,7 +112,6 @@ int interpol = 1;     //!< Ratio of interpolation
 /*
  Switch for some modes
 */
-int blackback = 1;   //!< Switch for black background
 int color_scale = 1; //!< Switch for full color scale mode
 int fbz = 1;         //!< Switch for 1st Brillouin zone mode
 int nodeline = 0;    //!< Switch for node lines
@@ -204,6 +203,10 @@ GLfloat  yellow[4] = { 1.0, 1.0, 0.0, 1.0 }; //!< Yellow color code
 GLfloat     red[4] = { 1.0, 0.0, 0.0, 1.0 }; //!< Red color code
 GLfloat   green[4] = { 0.0, 1.0, 0.0, 1.0 }; //!< Green color code
 GLfloat    blue[4] = { 0.0, 0.0, 1.0, 1.0 }; //!< Blue color code
+GLfloat BackGroundColor[4] = { 0.0, 0.0, 0.0, 1.0 };//!< BackGround color code
+GLfloat LineColor[4] = { 1.0, 1.0, 1.0, 1.0 };//!< Line color code
+GLfloat BarColor[5][4] = { { 0.0, 0.0, 1.0, 1.0 } , { 0.0, 1.0, 1.0, 1.0 } ,
+{ 0.0, 1.0, 0.0, 1.0 } ,{ 1.0, 1.0, 0.0, 1.0 } ,{ 1.0, 0.0, 0.0, 1.0 } };
 /*
  Others
 */
@@ -311,19 +314,17 @@ bool MyApp::OnInit()
   kveq = new GLfloat *** [nb];
   kveq_rot = new GLfloat *[nb];
   /**/
-  interpol_energy(avec, nb, interpol, ng0, ng,
-    terminal, eig, vf, mat, eig0, mat0);
-  init_corner(itet, corner);
-  bragg_vector(terminal, bvec, bragg, brnrm, &brnrm_min);
+  interpol_energy();
+  init_corner();
+  bragg_vector();
   myf->modify_band();
   /*
    Brillouin zone
   */
-  bz_lines(terminal, bragg, brnrm, &nbzl, bzl, &nbragg);
-  calc_2dbz(fbz, secvec, secscale, axis2d,
-    &nbzl2d, nbragg, bragg, brnrm, bzl2d, bzl2d_proj);
+  bz_lines();
+  calc_2dbz();
   /**/
-  max_and_min_bz(terminal, nb, ng0, eig0, mat0);
+  max_and_min_bz();
   /**/
   compute_patch_segment();
   /*

--- a/src/free_patch.cpp
+++ b/src/free_patch.cpp
@@ -44,36 +44,15 @@ THE SOFTWARE.
 #include <cstdio>
 #include "basic_math.hpp"
 #include "menu.hpp"
+#include "variable.hpp"
 /**
  @brief Free variables for patch before new patch is computed
 
  Free : ::nmlp, ::matp, ::clr, ::kvp, ::nmlp_rot, ::kvp_rot,
         ::kvnl, ::kvnl_rot, ::kv2d, ::clr2d
 */
-void free_patch(
-  int nb,
-  int refresh_patch,
-  int *ntri,
-  GLfloat ****kvp,
-  GLfloat **kvp_rot,
-  GLfloat ****matp,
-  GLfloat** clr,
-  GLfloat ****nmlp,
-  GLfloat **nmlp_rot,
-  GLfloat *****arw,
-  GLfloat **arw_rot,
-  int refresh_nodeline,
-  int *nnl,
-  GLfloat ****kvnl,
-  GLfloat **kvnl_rot,
-  int refresh_section,
-  GLfloat **kv2d,
-  GLfloat **clr2d,
-  int refresh_equator,
-  int *nequator,
-  GLfloat ****kveq,
-  GLfloat **kveq_rot
-) {
+void free_patch()
+{
   int ib, i0, i1, i2;
   /*
    Fermi patch
@@ -150,18 +129,7 @@ void free_patch(
 
  Modify : ::clr
 */
-void max_and_min(
-  int nb,
-  int nthreads,
-  wxTextCtrl *terminal,
-  MyFrame *myf,
-  int color_scale,
-  int *ntri,
-  GLfloat ****matp, 
-  GLfloat *patch_min,
-  GLfloat *patch_max,
-  GLfloat ****nmlp
-) 
+void max_and_min() 
 {
   int itri, ithread;
   GLfloat *max_th, *min_th;
@@ -201,11 +169,11 @@ shared(nb,ntri,matp,max_th,min_th) private(itri,ithread)
       }/*for (ib = 0; ib < nb; ib++)*/
     }/*End of parallel region*/
     /**/
-    *patch_max = max_th[0];
-    *patch_min = min_th[0];
+    patch_max = max_th[0];
+    patch_min = min_th[0];
     for (ithread = 1; ithread < nthreads; ithread++) {
-      if (max_th[ithread] > *patch_max) *patch_max = max_th[ithread];
-      if (min_th[ithread] < *patch_min) *patch_min = min_th[ithread];
+      if (max_th[ithread] > patch_max) patch_max = max_th[ithread];
+      if (min_th[ithread] < patch_min) patch_min = min_th[ithread];
     }
   }/*if (color_scale == 0 || color_scale == 4)*/
   else   if (color_scale == 2) {
@@ -232,11 +200,11 @@ shared(nb,ntri,matp,max_th,min_th) private(itri,ithread)
       }/*for (ib = 0; ib < nb; ib++)*/
     }/*End of parallel region*/
     /**/
-    *patch_min = min_th[0];
-    *patch_max = max_th[0];
+    patch_min = min_th[0];
+    patch_max = max_th[0];
     for (ithread = 1; ithread < nthreads; ithread++) {
-      if (max_th[ithread] < *patch_min) *patch_min = max_th[ithread];
-      if (max_th[ithread] > *patch_max) *patch_max = max_th[ithread];
+      if (max_th[ithread] < patch_min) patch_min = max_th[ithread];
+      if (max_th[ithread] > patch_max) patch_max = max_th[ithread];
     }
   }/*if (color_scale == 2)*/
   else   if (color_scale == 3) {
@@ -264,11 +232,11 @@ shared(nb,ntri,matp,min_th,max_th) private(itri,ithread)
       }/*for (ib = 0; ib < nb; ib++)*/
     }/*End of parallel region*/
     /**/
-    *patch_max = max_th[0];
-    *patch_min = min_th[0];
+    patch_max = max_th[0];
+    patch_min = min_th[0];
     for (ithread = 1; ithread < nthreads; ithread++) {
-      if (max_th[ithread] > *patch_max) *patch_max = max_th[ithread];
-      if (min_th[ithread] < *patch_min) *patch_min = min_th[ithread];
+      if (max_th[ithread] > patch_max) patch_max = max_th[ithread];
+      if (min_th[ithread] < patch_min) patch_min = min_th[ithread];
     }
   }/*if (color_scale == 3)*/
   else if (color_scale == 4 || color_scale == 7) {
@@ -297,19 +265,19 @@ shared(nb,ntri,nmlp,max_th,min_th) private(itri,ithread)
       }/*for (ib = 0; ib < nb; ib++)*/
     }/*End of parallel region*/
     /**/
-    *patch_max = max_th[0];
-    *patch_min = min_th[0];
+    patch_max = max_th[0];
+    patch_min = min_th[0];
     for (ithread = 1; ithread < nthreads; ithread++) {
-      if (max_th[ithread] > *patch_max) *patch_max = max_th[ithread];
-      if (min_th[ithread] < *patch_min) *patch_min = min_th[ithread];
+      if (max_th[ithread] > patch_max) patch_max = max_th[ithread];
+      if (min_th[ithread] < patch_min) patch_min = min_th[ithread];
     }    
   }/*if (color_scale == 5 || color_scale == 6)*/
 
   delete[] max_th;
   delete[] min_th;
 
-  myf->textbox_min->ChangeValue(wxString::Format(wxT("%f"), *patch_min));
-  myf->textbox_max->ChangeValue(wxString::Format(wxT("%f"), *patch_max));
+  myf->textbox_min->ChangeValue(wxString::Format(wxT("%f"), patch_min));
+  myf->textbox_max->ChangeValue(wxString::Format(wxT("%f"), patch_max));
 }/* max_and_min */
  /**
  @brief Compute Max. & Min. of matrix elements.
@@ -317,34 +285,14 @@ shared(nb,ntri,nmlp,max_th,min_th) private(itri,ithread)
 
  Modify : ::clr
  */
-void paint(
-  int nb,
-  int color_scale,
-  int blackback,
-  int *ntri,
-  GLfloat patch_min,
-  GLfloat patch_max,
-  GLfloat ****kvp,
-  GLfloat ****matp,
-  GLfloat **clr,
-  GLfloat ****nmlp,
-  GLfloat *****arw,
-  GLfloat red[4],
-  GLfloat green[4],
-  GLfloat blue[4],
-  GLfloat cyan[4],
-  GLfloat magenta[4],
-  GLfloat yellow[4],
-  GLfloat wgray[4],
-  GLfloat bgray[4]
-)
+void paint()
 {
   int itri, j;
   GLfloat origin[4];
 
   if (color_scale == 1) {
 #pragma omp parallel default(none) \
-shared(nb,ntri,matp,clr,cyan,blue,green,yellow,red,patch_max,patch_min) \
+shared(nb,ntri,matp,clr,BarColor,patch_max,patch_min) \
 private(itri, j)
     {
       int i, ib;
@@ -359,19 +307,23 @@ private(itri, j)
             mat2 = mat2 * 4.0f;
             /**/
             if (mat2 <= 1.0) {
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = cyan[j] * mat2 + blue[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[1][j] * mat2 + BarColor[0][j] * (1.0f - mat2);
             }
             else if (mat2 <= 2.0) {
               mat2 = mat2 - 1.0f;
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = green[j] * mat2 + cyan[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri]
+                = BarColor[2][j] * mat2 + BarColor[1][j] * (1.0f - mat2);
             }
             else if (mat2 <= 3.0) {
               mat2 = mat2 - 2.0f;
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = yellow[j] * mat2 + green[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri]
+                = BarColor[3][j] * mat2 + BarColor[2][j] * (1.0f - mat2);
             }
             else {
               mat2 = mat2 - 3.0f;
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = red[j] * mat2 + yellow[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[4][j] * mat2 + BarColor[3][j] * (1.0f - mat2);
             }
           }/*for (i = 0; i < 3; ++i)*/
         }/*for (itri = 0; itri < ntri[ib]; ++itri)*/
@@ -380,11 +332,10 @@ private(itri, j)
   }/*if (color_scale == 1 || color_scale == 2)*/
   else if (color_scale == 2) {
 
-    if (blackback == 1) for (j = 0; j < 4; ++j) origin[j] = wgray[j];
-    else for (j = 0; j < 4; ++j) origin[j] = bgray[j];
+    for (j = 0; j < 4; ++j) origin[j] = 1.0f - BackGroundColor[j];
 
 #pragma omp parallel default(none) \
-shared(nb,ntri,matp,clr,cyan,blue,green,yellow,red,magenta,bgray,wgray,blackback,patch_max,origin) \
+shared(nb,ntri,matp,clr,cyan,blue,green,yellow,red,magenta,bgray,wgray,patch_max,origin) \
 private(itri, j)
     {
       int i, ib;
@@ -441,7 +392,7 @@ private(itri, j)
   }/*if (color_scale == 2)*/
   else if (color_scale == 4) {
 #pragma omp parallel default(none) \
-shared(nb,ntri,nmlp,clr,cyan,blue,green,yellow,red,patch_max,patch_min) \
+shared(nb,ntri,nmlp,clr,BarColor,patch_max,patch_min) \
 private(itri, j)
     {
       int i, ib;
@@ -459,19 +410,23 @@ private(itri, j)
             mat2 = mat2 * 4.0f;
             /**/
             if (mat2 <= 1.0) {
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = cyan[j] * mat2 + blue[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri]
+                = BarColor[1][j] * mat2 + BarColor[0][j] * (1.0f - mat2);
             }
             else if (mat2 <= 2.0) {
               mat2 = mat2 - 1.0f;
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = green[j] * mat2 + cyan[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[2][j] * mat2 + BarColor[1][j] * (1.0f - mat2);
             }
             else if (mat2 <= 3.0) {
               mat2 = mat2 - 2.0f;
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = yellow[j] * mat2 + green[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri]
+                = BarColor[3][j] * mat2 + BarColor[2][j] * (1.0f - mat2);
             }
             else {
               mat2 = mat2 - 3.0f;
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = red[j] * mat2 + yellow[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[4][j] * mat2 + BarColor[3][j] * (1.0f - mat2);
             }
           }/*for (i = 0; i < 3; ++i)*/
         }/*for (itri = 0; itri < ntri[ib]; ++itri)*/
@@ -480,7 +435,7 @@ private(itri, j)
   }/*if (color_scale == 4)*/
   else if (color_scale == 3 || color_scale == 5) {
 #pragma omp parallel default(none) \
-shared(nb,ntri,matp,clr,cyan,blue,green,yellow,red,color_scale,kvp,arw,patch_max) \
+shared(nb,ntri,matp,clr,BarColor,color_scale,kvp,arw,patch_max) \
 private(itri, j)
     {
       int i, ib;
@@ -496,7 +451,8 @@ private(itri, j)
 #pragma omp for nowait
           for (itri = 0; itri < ntri[ib]; ++itri) {
             for (i = 0; i < 3; ++i) {
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = cyan[j] * mat2 + blue[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri]
+                = BarColor[1][j] * mat2 + BarColor[0][j] * (1.0f - mat2);
             }
           }
         }
@@ -505,7 +461,8 @@ private(itri, j)
 #pragma omp for nowait
           for (itri = 0; itri < ntri[ib]; ++itri) {
             for (i = 0; i < 3; ++i) {
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = green[j] * mat2 + cyan[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[2][j] * mat2 + BarColor[1][j] * (1.0f - mat2);
             }
           }
         }
@@ -514,7 +471,8 @@ private(itri, j)
 #pragma omp for nowait
           for (itri = 0; itri < ntri[ib]; ++itri) {
             for (i = 0; i < 3; ++i) {
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = yellow[j] * mat2 + green[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[3][j] * mat2 + BarColor[2][j] * (1.0f - mat2);
             }
           }
         }
@@ -523,7 +481,8 @@ private(itri, j)
 #pragma omp for nowait
           for (itri = 0; itri < ntri[ib]; ++itri) {
             for (i = 0; i < 3; ++i) {
-              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] = red[j] * mat2 + yellow[j] * (1.0f - mat2);
+              for (j = 0; j < 4; ++j) clr[ib][j + 4 * i + 12 * itri] 
+                = BarColor[4][j] * mat2 + BarColor[3][j] * (1.0f - mat2);
             }
           }
         }

--- a/src/free_patch.hpp
+++ b/src/free_patch.hpp
@@ -21,29 +21,7 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include "menu.hpp"
-void free_patch(int nb,  int refresh_patch,
-  int* ntri, GLfloat**** kvp, GLfloat** kvp_rot, GLfloat**** matp,
-  GLfloat** clr, GLfloat**** nmlp, GLfloat** nmlp_rot, GLfloat***** arw,
-  GLfloat** arw_rot, int refresh_nodeline, int* nnl, GLfloat**** kvnl,
-  GLfloat** kvnl_rot, int refresh_section, GLfloat** kv2d, GLfloat** clr2d,
-  int refresh_equator, int* nequator, GLfloat**** kveq, GLfloat** kveq_rot);
-void max_and_min(int nb, int nthreads, wxTextCtrl* terminal,
-  MyFrame* myf, int color_scale, int* ntri, GLfloat**** matp,
-  GLfloat *patch_min, GLfloat *patch_max, GLfloat**** nmlp);
-void paint(int nb,  int color_scale,  int blackback,
-  int* ntri,  GLfloat patch_min,  GLfloat patch_max,
-  GLfloat**** kvp,  GLfloat**** matp,  GLfloat** clr,
-  GLfloat**** nmlp,  GLfloat***** arw,
-  GLfloat red[4], GLfloat green[4], GLfloat blue[4],
-  GLfloat cyan[4], GLfloat magenta[4], GLfloat yellow[4],
-  GLfloat wgray[4], GLfloat bgray[4]);
+void free_patch();
+void max_and_min();
+void paint();
 

--- a/src/initialize.cpp
+++ b/src/initialize.cpp
@@ -33,16 +33,14 @@ THE SOFTWARE.
 #elif defined(HAVE_OPENGL_GL_H)
 #include <OpenGL/gl.h>
 #endif
-
+#include "variable.hpp"
 /**
  @brief Specify corners of tetrahedron
 
  Modify : ::corner
 */
-void init_corner(
-  int itet,
-  int corner[6][4]
-) {
+void init_corner()
+{
   int i, j;
   int corner1[16][6][4] = {
     /*
@@ -202,13 +200,8 @@ void init_corner(
 
  Modify : ::bragg, ::brnrm
 */
-void bragg_vector(
-  wxTextCtrl *terminal,
-  GLfloat bvec[3][3],
-  GLfloat bragg[26][3],
-  GLfloat brnrm[26],
-  GLfloat *brnrm_min
-) {
+void bragg_vector()
+{
   int i0, i1, i2, i, ibr;
   /**/
   ibr = 0;
@@ -241,23 +234,18 @@ void bragg_vector(
   /*
    Search min. of brnrm
   */
-  *brnrm_min = brnrm[0];
+  brnrm_min = brnrm[0];
   for (ibr = 1; ibr < 26; ibr++) {
-    if (*brnrm_min > brnrm[ibr]) *brnrm_min = brnrm[ibr];
+    if (brnrm_min > brnrm[ibr]) brnrm_min = brnrm[ibr];
   }
-  *terminal << wxString::Format(wxT("    Minimum Bragg norm : %f\n"), *brnrm_min);
+  *terminal << wxString::Format(wxT("    Minimum Bragg norm : %f\n"), brnrm_min);
 }/* bragg_vector */
 /**
  @brief Print max and minimum @f$\varepsilon_{n k}, \Delta_{n k}@f$
         in the whole Brillouine zone
 */
-void max_and_min_bz(
-  wxTextCtrl* terminal, 
-  int nb,
-  int ng0[3],
-  GLfloat**** eig0, 
-  GLfloat***** mat0
-) {
+void max_and_min_bz()
+{
   int ib, i0, i1, i2;
   GLfloat eigmin, eigmax, matmin, matmax;
   /**/

--- a/src/initialize.hpp
+++ b/src/initialize.hpp
@@ -21,16 +21,6 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#include <wx/wx.h>
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-
-void init_corner(int itet, int corner[6][4]);
-void bragg_vector(wxTextCtrl* terminal, GLfloat bvec[3][3], GLfloat bragg[26][3], GLfloat brnrm[26], GLfloat *brnrm_min);
-void max_and_min_bz(wxTextCtrl* terminal, int nb, int ng0[3], GLfloat ****eig0, GLfloat *****mat0);
+void init_corner();
+void bragg_vector();
+void max_and_min_bz();

--- a/src/kumo.cpp
+++ b/src/kumo.cpp
@@ -37,6 +37,7 @@ denser @f$k@f$-grid with French-curve (Kumo) interpolation
 #include <cstdlib>
 #include <cstdio>
 #include "basic_math.hpp"
+#include "variable.hpp"
 /**
  @brief Compute coefficient for the French-curve (Kumo) interpolation
  @f[
@@ -62,19 +63,8 @@ static void kumo_coef(
 
  Modify : ::eig, ::mat
 */
-void interpol_energy(
-  GLfloat avec[3][3],
-  int nb,
-  int interpol,
-  int ng0[3],
-  int ng[3],
-  wxTextCtrl *terminal,
-  GLfloat ****eig,
-  GLfloat *****vf,
-  GLfloat *****mat,
-  GLfloat**** eig0,
-  GLfloat***** mat0
-) {
+void interpol_energy() 
+{
   int ib, i0, i1, i2, ii;
 
   *terminal << wxT("    Interpolating ... ");

--- a/src/kumo.hpp
+++ b/src/kumo.hpp
@@ -21,15 +21,4 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include <wx/wx.h>
-
-void interpol_energy(GLfloat avec[3][3], int nb, int interpol, int ng0[3], int ng[3], wxTextCtrl* terminal,
-  GLfloat**** eig, GLfloat***** vf, GLfloat***** mat, GLfloat**** eig0, GLfloat***** mat0);
+void interpol_energy();

--- a/src/menu.cpp
+++ b/src/menu.cpp
@@ -52,16 +52,12 @@ THE SOFTWARE.
 
 void compute_patch_segment() {
   if (refresh_interpol == 1){
-    interpol_energy(avec, nb, interpol, ng0, ng,
-      terminal, eig, vf, mat, eig0, mat0);
+    interpol_energy();
     refresh_patch = 1;
     refresh_interpol = 0;
   }
   if (refresh_patch == 1) {
-    fermi_patch(nb, nthreads, ng, ng0, shiftk,
-      fbz, eig, EF, mat, vf, corner, bvec, nbragg,
-      bragg, brnrm, brnrm_min, terminal, ntri, kvp, matp,
-      clr, arw, nmlp, kvp_rot, nmlp_rot, arw_rot);
+    fermi_patch();
     refresh_color = 1;
     refresh_section = 1;
     refresh_equator = 1;
@@ -70,37 +66,28 @@ void compute_patch_segment() {
   }
   if (refresh_color == 1) {
     if (skip_minmax == 1) skip_minmax = 0;
-    else max_and_min(nb, nthreads, terminal, myf, color_scale, ntri, matp,
-      &patch_min, &patch_max, nmlp);
-    paint(nb, color_scale, blackback, ntri, patch_min, patch_max,
-      kvp, matp, clr, nmlp, arw, red, green, blue,
-      cyan, magenta, yellow, wgray, bgray);
+    else max_and_min();
+    paint();
     refresh_section = 1;
     refresh_color = 0;
   }
   if (refresh_nodeline == 1) {
-    calc_nodeline(nb, ntri, matp, kvp, nthreads, terminal, nnl, kvnl, kvnl_rot);
+    calc_nodeline();
     refresh_nodeline = 0;
   }
   if (refresh_section == 1) {
-    calc_2dbz(fbz, secvec, secscale, axis2d,
-      &nbzl2d, nbragg, bragg, brnrm, bzl2d, bzl2d_proj);
-    calc_section(fbz, nb, nthreads, secvec, secscale,
-      axis2d, ntri, kvp, clr, terminal, n2d, kv2d, clr2d);
+    calc_2dbz();
+    calc_section();
     refresh_section = 0;
   }
   if (refresh_equator == 1) {
-    equator(nb, nthreads, eqvec, ntri, kvp,
-      nmlp, terminal, nequator, kveq, kveq_rot);
+    equator();
     refresh_equator = 0;
   }
 }
 
 void refresh_patch_segment() {
-  free_patch(nb, refresh_patch, ntri, kvp, kvp_rot, matp,
-    clr, nmlp, nmlp_rot, arw, arw_rot, refresh_nodeline, nnl, kvnl,
-    kvnl_rot, refresh_section, kv2d, clr2d,
-    refresh_equator, nequator, kveq, kveq_rot);
+  free_patch();
   compute_patch_segment();
 }
 
@@ -138,7 +125,15 @@ enum
   itext_roty,
   itext_rotz,
   ibutton_rotate,
-  icheck_band
+  icheck_band,
+  itext_BackGroundR,
+  itext_BackGroundG,
+  itext_BackGroundB,
+  itext_LineColorR,
+  itext_LineColorG,
+  itext_LineColorB,
+  iradio_BarColor,
+  ibutton_section
 };
 
 void MyFrame::button_refresh(
@@ -150,32 +145,105 @@ void MyFrame::button_refresh(
 void MyFrame::button_compute(
   wxCommandEvent& event//!<[in] Selected menu
 ) {
-  free_patch(nb, refresh_patch, ntri, kvp, kvp_rot, matp,
-    clr, nmlp, nmlp_rot, arw, arw_rot, refresh_nodeline, nnl, kvnl,
-    kvnl_rot, refresh_section, kv2d, clr2d,
-    refresh_equator, nequator, kveq, kveq_rot);
+  free_patch();
   compute_patch_segment();
   Refresh(false);
 }
+void MyFrame::button_section(
+  wxCommandEvent& event//!<[in] Selected menu
+) {
+  int ib, i2d, i, j;
+  FILE *fp;
+  fp = fopen("fermi_line.dat", "w");
+  for (ib = 0; ib < nb; ib++) {
+    if (draw_band[ib] == 1) {
+      for (i2d = 0; i2d < n2d[ib]; i2d++) {
+        for (i = 0; i < 2; i++) {
+          fprintf(fp, "%15.5e %15.5e\n",
+            kv2d[ib][i * 3 + 6 * i2d], 
+            kv2d[ib][1 + i * 3 + 6 * i2d]);
+        }
+        fprintf(fp, "\n\n");
+      }
+    }/*if (draw_band[ib] == 1)*/
+  }/* for (ib = 0; ib < nb; ib++)*/
+  fclose(fp);
+  *terminal << wxT("  fermi_line.dat was written.\n");
+
+  fp = fopen("bz_line.dat", "w");
+  for (i2d = 0; i2d < nbzl2d; i2d++) {
+    fprintf(fp, "%15.5e %15.5e\n",
+      bzl2d_proj[i2d][0], bzl2d_proj[i2d][1]);
+  }
+  fprintf(fp, "%15.5e %15.5e\n",
+    bzl2d_proj[0][0], bzl2d_proj[0][1]);
+  fclose(fp);
+  *terminal << wxT("  bz_line.dat was written.\n");
+}
+/**
+@brief Change Line color color (::blackback)
+*/
+void MyFrame::textctrl_LineColor(
+  wxCommandEvent& event //!<[in] Selected menu
+)
+{
+  int ierr;
+  double dvalue;
+
+  if (event.GetId() == itext_LineColorR) {
+    if (event.GetString().ToDouble(&dvalue)) {
+      LineColor[0] = (GLfloat)dvalue;
+      Refresh(false);
+    }
+  }
+  else  if (event.GetId() == itext_LineColorG) {
+    if (event.GetString().ToDouble(&dvalue)) {
+      LineColor[1] = (GLfloat)dvalue;
+      Refresh(false);
+    }
+  }
+  else  if (event.GetId() == itext_LineColorB) {
+    if (event.GetString().ToDouble(&dvalue)) {
+      LineColor[2] = (GLfloat)dvalue;
+      Refresh(false);
+    }
+  }
+}
 /**
 @brief Change background color (::blackback)
 */
-void MyFrame::radio_background(
+void MyFrame::textctrl_BackGround(
   wxCommandEvent& event //!<[in] Selected menu
 )
 {
-  if (event.GetString().Cmp(wxT("Black")) == 0 && blackback != 1) {
-    glClearColor(0.0, 0.0, 0.0, 0.0);
-    blackback = 1;
-    // debug if (color_scale == 2 || color_scale == 3) paint();
+  int ierr;
+  double dvalue;
+
+  if (event.GetId() == itext_BackGroundR) {
+    if (event.GetString().ToDouble(&dvalue)) {
+      BackGroundColor[0] = (GLfloat)dvalue;
+      glClearColor(BackGroundColor[0], BackGroundColor[1],
+        BackGroundColor[2], BackGroundColor[3]);
+      Refresh(false);
+    }
+  }
+  else  if (event.GetId() == itext_BackGroundG) {
+    if (event.GetString().ToDouble(&dvalue)) {
+      BackGroundColor[1] = (GLfloat)dvalue;
+      glClearColor(BackGroundColor[0], BackGroundColor[1], 
+        BackGroundColor[2], BackGroundColor[3]);
+      Refresh(false);
+    }
   }
-  else if (event.GetString().Cmp(wxT("White")) == 0 && blackback != 0) {
-    glClearColor(1.0, 1.0, 1.0, 0.0);
-    blackback = 0;
-    // debug if (color_scale == 2 || color_scale == 3) paint();
+  else  if (event.GetId() == itext_BackGroundB) {
+    if (event.GetString().ToDouble(&dvalue)) {
+      BackGroundColor[2] = (GLfloat)dvalue;
+      glClearColor(BackGroundColor[0], BackGroundColor[1],
+        BackGroundColor[2], BackGroundColor[3]);
+      Refresh(false);
+    }
   }
-  Refresh(false);
-}/* bgcolor change*/
+}
  /**
  @brief Toggle the appearance of each band (::draw_band)
 */
@@ -208,6 +276,44 @@ void MyFrame::radio_brillouinzone(
   refresh_patch = 1;
 } /* menu_brillouinzone */
 /**
+ @brief Change Brillouin zone (::fbz)
+*/
+void MyFrame::radio_BarColor(
+  wxCommandEvent& event //!<[in] Selected menu
+)
+{
+  int ii;
+  if (event.GetString().Cmp(wxT("BGR")) == 0) {
+    for (ii = 0; ii < 4; ii++) {
+      BarColor[0][ii] = blue[ii];
+      BarColor[1][ii] = cyan[ii];
+      BarColor[2][ii] = green[ii];
+      BarColor[3][ii] = yellow[ii];
+      BarColor[4][ii] = red[ii];
+    }
+  }
+  else if (event.GetString().Cmp(wxT("CMY")) == 0) {
+    for (ii = 0; ii < 4; ii++) {
+      BarColor[0][ii] = cyan[ii];
+      BarColor[1][ii] = blue[ii];
+      BarColor[2][ii] = magenta[ii];
+      BarColor[3][ii] = red[ii];
+      BarColor[4][ii] = yellow[ii];
+    }
+  }
+  else if (event.GetString().Cmp(wxT("MCY")) == 0) {
+    for (ii = 0; ii < 4; ii++) {
+      BarColor[0][ii] = magenta[ii];
+      BarColor[1][ii] = blue[ii];
+      BarColor[2][ii] = cyan[ii];
+      BarColor[3][ii] = green[ii];
+      BarColor[4][ii] = yellow[ii];
+    }
+  }
+  paint();
+  Refresh(false);
+} /* menu_brillouinzone */
+/**
  @brief Toggle Colorbar (::lcolorbar)
 */
 void MyFrame::check_colorbar(
@@ -440,7 +546,7 @@ void MyFrame::radio_tetra(
 )
 {
   itet = wxAtoi(event.GetString()) - 1;
-  init_corner(itet, corner);
+  init_corner();
   refresh_patch = 1;
 }/*menu_tetra*/
  /**
@@ -659,7 +765,7 @@ wxT("8"), wxT("9"), wxT("10"), wxT("11"), wxT("12"), wxT("13"), wxT("14"),
   gbsizer->Add(new wxRadioBox(panel, iradio_tetra, wxT("Tetrahedron"),
     wxDefaultPosition, wxDefaultSize,
     WXSIZEOF(choices_tetra), choices_tetra,
-    4, wxRA_SPECIFY_COLS), wxGBPosition(3,2), wxGBSpan(4, 2));
+    4, wxRA_SPECIFY_COLS), wxGBPosition(3,2), wxGBSpan(5, 2));
 
   wxString choices_colorscale[] = { wxT("Input (1D)"), wxT("Input (2D)"),
     wxT("Input (3D)"), wxT("Fermi Velocity"), wxT("Band Index"),
@@ -670,91 +776,125 @@ wxT("8"), wxT("9"), wxT("10"), wxT("11"), wxT("12"), wxT("13"), wxT("14"),
     WXSIZEOF(choices_colorscale), choices_colorscale,
     1, wxRA_SPECIFY_COLS);
   Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radiovalue_colorscale, this, iradio_colorscale);
-  gbsizer->Add(radiobox_color, wxGBPosition(7, 0), wxGBSpan(2, 2));
+  gbsizer->Add(radiobox_color, wxGBPosition(7, 0), wxGBSpan(3, 2));
 
   wxString choices_bz[] = { wxT("First Brillouin zone"), wxT("Primitive Brillouin zone") };
   Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radio_brillouinzone, this, iradio_brillouinzone);
   gbsizer->Add(new wxRadioBox(panel, iradio_brillouinzone, wxT("Brillouin zone"),
     wxDefaultPosition, wxDefaultSize,
     WXSIZEOF(choices_bz), choices_bz,
-    1, wxRA_SPECIFY_COLS), wxGBPosition(7, 2), wxGBSpan(1, 2));
+    1, wxRA_SPECIFY_COLS), wxGBPosition(8, 2), wxGBSpan(1, 2));
 
   wxString choices_stereo[] = { wxT("None"), wxT("Parallel"), wxT("Cross") };
   Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radio_stereo, this, iradio_stereo);
   gbsizer->Add(new wxRadioBox(panel, iradio_stereo, wxT("Stereogram"),
     wxDefaultPosition, wxDefaultSize,
     WXSIZEOF(choices_stereo), choices_stereo,
-    1, wxRA_SPECIFY_COLS), wxGBPosition(8, 2), wxGBSpan(1, 1));
+    1, wxRA_SPECIFY_COLS), wxGBPosition(9, 2), wxGBSpan(1, 1));
 
   wxString choices_mouse[] = { wxT("Rotate"), wxT("Scale"), wxT("Translate") };
   Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radio_mouse, this, iradio_mouse);
   gbsizer->Add(new wxRadioBox(panel, iradio_mouse, wxT("Mouse Drag"),
     wxDefaultPosition, wxDefaultSize,
     WXSIZEOF(choices_mouse), choices_mouse,
-    1, wxRA_SPECIFY_COLS), wxGBPosition(8, 3), wxGBSpan(1, 1));
+    1, wxRA_SPECIFY_COLS), wxGBPosition(9, 3), wxGBSpan(1, 1));
+
+  gbsizer->Add(new wxStaticText(panel, wxID_ANY, wxT("Background (RGB) : ")),
+    wxGBPosition(10, 0), wxGBSpan(1, 1), wxALIGN_RIGHT);
+  Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_BackGround, this, itext_BackGroundR);
+  textbox_BackGroundR = new wxTextCtrl(panel, itext_BackGroundR, wxT(""));
+  gbsizer->Add(textbox_BackGroundR, wxGBPosition(10, 1), wxGBSpan(1, 1));
+  Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_BackGround, this, itext_BackGroundG);
+  textbox_BackGroundG = new wxTextCtrl(panel, itext_BackGroundG, wxT(""));
+  gbsizer->Add(textbox_BackGroundG, wxGBPosition(10, 2), wxGBSpan(1, 1));
+  Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_BackGround, this, itext_BackGroundB);
+  textbox_BackGroundB = new wxTextCtrl(panel, itext_BackGroundB, wxT(""));
+  gbsizer->Add(textbox_BackGroundB, wxGBPosition(10, 3), wxGBSpan(1, 1));
+  textbox_BackGroundR->ChangeValue(wxT("0"));
+  textbox_BackGroundG->ChangeValue(wxT("0"));
+  textbox_BackGroundB->ChangeValue(wxT("0"));
+
+  gbsizer->Add(new wxStaticText(panel, wxID_ANY, wxT("Line color (RGB) : ")),
+    wxGBPosition(11, 0), wxGBSpan(1, 1), wxALIGN_RIGHT);
+  Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_LineColor, this, itext_LineColorR);
+  textbox_LineColorR = new wxTextCtrl(panel, itext_LineColorR, wxT(""));
+  gbsizer->Add(textbox_LineColorR, wxGBPosition(11, 1), wxGBSpan(1, 1));
+  Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_LineColor, this, itext_LineColorG);
+  textbox_LineColorG = new wxTextCtrl(panel, itext_LineColorG, wxT(""));
+  gbsizer->Add(textbox_LineColorG, wxGBPosition(11, 2), wxGBSpan(1, 1));
+  Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_LineColor, this, itext_LineColorB);
+  textbox_LineColorB = new wxTextCtrl(panel, itext_LineColorB, wxT(""));
+  gbsizer->Add(textbox_LineColorB, wxGBPosition(11, 3), wxGBSpan(1, 1));
+  textbox_LineColorR->ChangeValue(wxT("1"));
+  textbox_LineColorG->ChangeValue(wxT("1"));
+  textbox_LineColorB->ChangeValue(wxT("1"));
 
   Bind(wxEVT_COMMAND_BUTTON_CLICKED, &MyFrame::textctrl_view, this, ibutton_rotate);
   gbsizer->Add(new wxButton(panel, ibutton_rotate, wxT("Rotate")), 
-    wxGBPosition(9, 0), wxGBSpan(1, 1));
+    wxGBPosition(12, 0), wxGBSpan(1, 1));
   Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_view, this, itext_rotx);
   textbox_rotatex = new wxTextCtrl(panel, itext_rotx, wxT(""));
-  gbsizer->Add(textbox_rotatex, wxGBPosition(9, 1), wxGBSpan(1, 1));
+  gbsizer->Add(textbox_rotatex, wxGBPosition(12, 1), wxGBSpan(1, 1));
   Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_view, this, itext_roty);
   textbox_rotatey = new wxTextCtrl(panel, itext_roty, wxT(""));
-  gbsizer->Add(textbox_rotatey, wxGBPosition(9, 2), wxGBSpan(1, 1));
+  gbsizer->Add(textbox_rotatey, wxGBPosition(12, 2), wxGBSpan(1, 1));
   Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_view, this, itext_rotz);
   textbox_rotatez = new wxTextCtrl(panel, itext_rotz, wxT(""));
-  gbsizer->Add(textbox_rotatez, wxGBPosition(9, 3), wxGBSpan(1, 1));
+  gbsizer->Add(textbox_rotatez, wxGBPosition(12, 3), wxGBSpan(1, 1));
   textbox_rotatex->ChangeValue(wxT("0"));
   textbox_rotatey->ChangeValue(wxT("0"));
   textbox_rotatez->ChangeValue(wxT("0"));
 
   gbsizer->Add(new wxStaticText(panel, wxID_ANY, wxT("Position : ")), 
-    wxGBPosition(10, 1), wxGBSpan(1, 1), wxALIGN_RIGHT);
+    wxGBPosition(13, 1), wxGBSpan(1, 1), wxALIGN_RIGHT);
   Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_view, this, itext_positionx);
   textbox_positionx = new wxTextCtrl(panel, itext_positionx, wxT(""));
-  gbsizer->Add(textbox_positionx, wxGBPosition(10, 2), wxGBSpan(1, 1));
+  gbsizer->Add(textbox_positionx, wxGBPosition(13, 2), wxGBSpan(1, 1));
   Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_view, this, itext_positiony);
   textbox_positiony = new wxTextCtrl(panel, itext_positiony, wxT(""));
-  gbsizer->Add(textbox_positiony, wxGBPosition(10, 3), wxGBSpan(1, 1));
+  gbsizer->Add(textbox_positiony, wxGBPosition(13, 3), wxGBSpan(1, 1));
   textbox_positionx->ChangeValue(wxT("0"));
   textbox_positiony->ChangeValue(wxT("0"));
 
   gbsizer->Add(new wxStaticText(panel, wxID_ANY, wxT("Scale : ")), 
-    wxGBPosition(11, 2), wxGBSpan(1, 1), wxALIGN_RIGHT);
+    wxGBPosition(14, 2), wxGBSpan(1, 1), wxALIGN_RIGHT);
   Bind(wxEVT_COMMAND_TEXT_UPDATED, &MyFrame::textctrl_view, this, itext_scale);
   textbox_scale = new wxTextCtrl(panel, itext_scale, wxT(""));
-  gbsizer->Add(textbox_scale, wxGBPosition(11, 3), wxGBSpan(1, 1));
+  gbsizer->Add(textbox_scale, wxGBPosition(14, 3), wxGBSpan(1, 1));
   textbox_scale->ChangeValue(wxString::Format(wxT("%f"), scl));
 
   Bind(wxEVT_COMMAND_CHECKBOX_CLICKED, &MyFrame::check_colorbar, this, icheck_colorbar);
   wxCheckBox* check = new wxCheckBox(panel, icheck_colorbar, wxT("Color bar"));
-  gbsizer->Add(check, wxGBPosition(11, 1), wxGBSpan(1, 1));
+  gbsizer->Add(check, wxGBPosition(14, 1), wxGBSpan(1, 1));
   check->SetValue(true);
   // debug fileMenu->Check(menu_colorbar_check, true);
 
   Bind(wxEVT_COMMAND_CHECKBOX_CLICKED, &MyFrame::checkvalue_equator, this, icheck_equator);
-  gbsizer->Add(new wxCheckBox(panel, icheck_equator, wxT("Equator")), wxGBPosition(12, 1), wxGBSpan(1, 1));
+  gbsizer->Add(new wxCheckBox(panel, icheck_equator, wxT("Equator")), wxGBPosition(15, 1), wxGBSpan(1, 1));
 
   Bind(wxEVT_COMMAND_CHECKBOX_CLICKED, &MyFrame::check_nodeline, this, icheck_nodeline);
-  gbsizer->Add(new wxCheckBox(panel, icheck_nodeline, wxT("Nodal line")), wxGBPosition(13, 1), wxGBSpan(1, 1));
+  gbsizer->Add(new wxCheckBox(panel, icheck_nodeline, wxT("Nodal line")), wxGBPosition(16, 1), wxGBSpan(1, 1));
 
   Bind(wxEVT_COMMAND_CHECKBOX_CLICKED, &MyFrame::radiovalue_section, this, icheck_section);
-  gbsizer->Add(new wxCheckBox(panel, icheck_section, wxT("Section")), wxGBPosition(14, 1), wxGBSpan(1, 1));
+  gbsizer->Add(new wxCheckBox(panel, icheck_section, wxT("Section")), wxGBPosition(17, 1), wxGBSpan(1, 1));
+
+  Bind(wxEVT_COMMAND_BUTTON_CLICKED, &MyFrame::button_section, this, ibutton_section);
+  gbsizer->Add(new wxButton(panel, ibutton_section, wxT("Section file")),
+    wxGBPosition(18, 1), wxGBSpan(1, 1));
 
   wxString choices_light[] = { wxT("Both"), wxT("Unoccupy"), wxT("Occupy") };
   Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radio_lighting, this, iradio_lighting);
   gbsizer->Add(new wxRadioBox(panel, iradio_lighting, wxT("Lighting"),
     wxDefaultPosition, wxDefaultSize,
     WXSIZEOF(choices_light), choices_light,
-    1, wxRA_SPECIFY_COLS), wxGBPosition(12, 2), wxGBSpan(3, 1));
+    1, wxRA_SPECIFY_COLS), wxGBPosition(15, 2), wxGBSpan(4, 1));
 
-  wxString choices_bg[] = { wxT("Black"), wxT("White") };
-  Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radio_background, this, iradio_background);
-  gbsizer->Add(new wxRadioBox(panel, iradio_background, wxT("Background color"),
+  wxString choices_BarColor[] = { wxT("BGR"), wxT("CMY"), wxT("MCY")};
+  Bind(wxEVT_COMMAND_RADIOBOX_SELECTED, &MyFrame::radio_BarColor, this, iradio_BarColor);
+  gbsizer->Add(new wxRadioBox(panel, iradio_BarColor, wxT("Bar Color"),
     wxDefaultPosition, wxDefaultSize,
-    WXSIZEOF(choices_bg), choices_bg,
-    1, wxRA_SPECIFY_COLS), wxGBPosition(12, 3), wxGBSpan(3, 1));
+    WXSIZEOF(choices_BarColor), choices_BarColor,
+    1, wxRA_SPECIFY_COLS), wxGBPosition(15, 3), wxGBSpan(4, 1));
 
   SetSizer(sizermain);
   SetAutoLayout(true);
@@ -790,7 +930,7 @@ void MyFrame::modify_band() {
     Bind(wxEVT_COMMAND_CHECKBOX_CLICKED, &MyFrame::check_band, this, icheck_band + ib);
     check[ib] = new wxCheckBox(panel, icheck_band + ib, 
       wxString::Format(wxT("Band %d"), ib));
-    gbsizer->Add(check[ib], wxGBPosition(10 + ib, 0), wxGBSpan(1, 1));
+    gbsizer->Add(check[ib], wxGBPosition(13 + ib, 0), wxGBSpan(1, 1));
     check[ib]->SetValue(true);
   }
   gbsizer->Layout();

--- a/src/menu.hpp
+++ b/src/menu.hpp
@@ -62,12 +62,19 @@ public:
   wxTextCtrl* textbox_min;
   wxTextCtrl* textbox_max;
   wxTextCtrl* textbox_linewidth;
+  wxTextCtrl* textbox_BackGroundR;
+  wxTextCtrl* textbox_BackGroundG;
+  wxTextCtrl* textbox_BackGroundB;
+  wxTextCtrl* textbox_LineColorR;
+  wxTextCtrl* textbox_LineColorG;
+  wxTextCtrl* textbox_LineColorB;
 
 private:
   void OnExit(wxCommandEvent& event);
   void check_band(wxCommandEvent& event);
   void radio_brillouinzone(wxCommandEvent& event);
-  void radio_background(wxCommandEvent& event);
+  void textctrl_BackGround(wxCommandEvent& event);
+  void textctrl_LineColor(wxCommandEvent& event);
   void check_colorbar(wxCommandEvent& event);
   void radiovalue_colorscale(wxCommandEvent& event);
   void checkvalue_equator(wxCommandEvent& event);
@@ -83,6 +90,8 @@ private:
   void textctrl_view(wxCommandEvent& event);
   void button_refresh(wxCommandEvent& event);
   void button_compute(wxCommandEvent& event);
+  void radio_BarColor(wxCommandEvent& event);
+  void button_section(wxCommandEvent& event);
   //wxDECLARE_EVENT_TABLE();
 };
 

--- a/src/read_file.cpp
+++ b/src/read_file.cpp
@@ -268,15 +268,8 @@ void read_batch()
       continue;
     }
 
-    if (strcmp(keyword, "background") == 0) {
-      ierr = fscanf(fp, "%s", value);
-      Text2Lower(value);
-      if (strcmp(value, "black") == 0) blackback = 1;
-      else if (strcmp(value, "white") == 0) blackback = 0;
-      else {
-        printf("Error! %s = %s", keyword, value);
-        exit(-1);
-      }
+    if (strcmp(keyword, "backgroundcolor") == 0) {
+      ierr = fscanf(fp, "%f%f%f", &BackGroundColor[0], &BackGroundColor[1], &BackGroundColor[2]);
     }
     else if (strcmp(keyword, "band") == 0) {
       for (ib = 0; ib < nb; ib++)
@@ -321,6 +314,9 @@ void read_batch()
     else if (strcmp(keyword, "interpol") == 0) {
       ierr = fscanf(fp, "%d", &interpol);
     }
+    else if (strcmp(keyword, "linecolor") == 0) {
+      ierr = fscanf(fp, "%f%f%f", &LineColor[0], &LineColor[1], &LineColor[2]);
+    }
     else if (strcmp(keyword, "linewidth") == 0) {
       ierr = fscanf(fp, "%f", &linewidth);
     }
@@ -388,8 +384,7 @@ void read_batch()
   }
   fclose(fp);
 
-  if (blackback == 1) glClearColor(0.0, 0.0, 0.0, 0.0);
-  else glClearColor(1.0, 1.0, 1.0, 0.0);
+  glClearColor(BackGroundColor[0], BackGroundColor[1], BackGroundColor[2], BackGroundColor[2]);
 
   thetax = 3.14159265f / 180.0f * thetax;
   thetay = 3.14159265f / 180.0f * thetay;

--- a/src/read_file.hpp
+++ b/src/read_file.hpp
@@ -21,16 +21,6 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include <wx/wx.h>
-
 int read_file();
 void read_batch();
 void read_bxsf();

--- a/src/section.cpp
+++ b/src/section.cpp
@@ -39,6 +39,7 @@ THE SOFTWARE.
 #include <cstdlib>
 #include <cmath>
 #include "basic_math.hpp"
+#include "variable.hpp"
 /**
  @brief Project 3D \f$k\f$-vector into 2D plane. 
 
@@ -177,18 +178,7 @@ int bragg_vert2d(
 
  Modify : ::nbzl2d, ::bzl2d_proj
 */
-void calc_2dbz(
-int fbz,
-GLfloat secvec[3],
-GLfloat secscale,
-GLfloat axis2d[2][3],
-int *nbzl2d,
-int nbragg,
-GLfloat bragg[26][3],
-GLfloat brnrm[26],
-GLfloat bzl2d[26][3],
-GLfloat bzl2d_proj[26][3]
-) {
+void calc_2dbz() {
   int jbr, nbr, i, j, lvert, ibzl;
   GLfloat vert[2][3], vec[26][2][3], prod, thr;
 
@@ -198,7 +188,7 @@ GLfloat bzl2d_proj[26][3]
   */
   set2daxis(secvec, axis2d);
 
-  *nbzl2d = 0;
+  nbzl2d = 0;
 
   for (jbr = 0; jbr < nbragg; ++jbr) {
     /**/
@@ -211,15 +201,15 @@ GLfloat bzl2d_proj[26][3]
     lvert = bragg_vert2d(nbragg, bragg, brnrm, secvec, secscale, jbr, nbr, vert[1], vert[0]);
     if (lvert == 0) continue;
     /**/
-    for (i = 0; i < 2; ++i) for (j = 0; j < 3; ++j) vec[*nbzl2d][i][j] = vert[i][j];
-    *nbzl2d += 1;
+    for (i = 0; i < 2; ++i) for (j = 0; j < 3; ++j) vec[nbzl2d][i][j] = vert[i][j];
+    nbzl2d += 1;
   }/*for (jbr = 0; jbr < nbragg; ++jbr)*/
   /*
    Order bz lines
   */
   for (i = 0; i < 3; i++) bzl2d[0][i] = vec[0][0][i];
   for (i = 0; i < 3; i++) bzl2d[1][i] = vec[0][1][i];
-  for (ibzl = 0; ibzl < *nbzl2d; ibzl++) {
+  for (ibzl = 0; ibzl < nbzl2d; ibzl++) {
 
     thr = 0.0f;
     for (j = 0; j < 2; j++) for (i = 0; i < 3; i++) thr += bzl2d[j][i] * bzl2d[j][i];
@@ -238,21 +228,21 @@ GLfloat bzl2d_proj[26][3]
       for (j = 0; j < 2; j++) for (i = 0; i < 3; i++) vec[ibzl][j][i] = 0.0;
   }/*for (ibzl = 1; ibzl < *nbzl2d; ibzl++)*/
 
-  for (jbr = 1; jbr < *nbzl2d - 1; jbr++) {
+  for (jbr = 1; jbr < nbzl2d - 1; jbr++) {
 
     thr = 0.0f;
     for (j = 0; j < 2; j++) for (i = 0; i < 3; i++) thr += bzl2d[jbr][i] * bzl2d[jbr][i];
     thr *= 0.001f;
 
     prod = 0.0;
-    for (ibzl = 0; ibzl < *nbzl2d; ibzl++) for (i = 0; i < 3; i++)
+    for (ibzl = 0; ibzl < nbzl2d; ibzl++) for (i = 0; i < 3; i++)
       prod += vec[ibzl][0][i] * vec[ibzl][0][i];
     if (prod < thr) {
-      *nbzl2d = jbr + 1;
+      nbzl2d = jbr + 1;
       break;
     }
 
-    for (ibzl = 1; ibzl < *nbzl2d; ibzl++) {
+    for (ibzl = 1; ibzl < nbzl2d; ibzl++) {
       prod = (bzl2d[jbr][0] - vec[ibzl][0][0]) * (bzl2d[jbr][0] - vec[ibzl][0][0])
            + (bzl2d[jbr][1] - vec[ibzl][0][1]) * (bzl2d[jbr][1] - vec[ibzl][0][1])
            + (bzl2d[jbr][2] - vec[ibzl][0][2]) * (bzl2d[jbr][2] - vec[ibzl][0][2]);
@@ -273,7 +263,7 @@ GLfloat bzl2d_proj[26][3]
   /*
    Project into 2D plane
   */
-  for (ibzl = 0; ibzl < *nbzl2d; ibzl++) {
+  for (ibzl = 0; ibzl < nbzl2d; ibzl++) {
     for (i = 0; i < 3; i++) bzl2d_proj[ibzl][i] = bzl2d[ibzl][i];
     proj_2d(axis2d, bzl2d_proj[ibzl]);
   }/*for (ibzl = 0; ibzl < *nbzl2d; ibzl++)*/
@@ -283,21 +273,7 @@ GLfloat bzl2d_proj[26][3]
 
  Modify : ::n2d, ::clr2d, ::kv2d
 */
-void calc_section(
-int fbz,
-int nb,
-int nthreads,
-GLfloat secvec[3],
-GLfloat secscale,
-GLfloat axis2d[2][3],
-int *ntri,
-GLfloat ****kvp,
-GLfloat ** clr,
-wxTextCtrl *terminal,
-int *n2d,
-GLfloat **kv2d,
-GLfloat **clr2d
-) {
+void calc_section() {
   int i, j, ib, itri, ithread, n2d0;
   std::vector<std::vector<std::vector<std::vector<GLfloat>>>> kv2d_v, clr2d_v;
 

--- a/src/section.hpp
+++ b/src/section.hpp
@@ -21,19 +21,5 @@ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
 THE SOFTWARE.
 */
-#if defined(HAVE_CONFIG_H)
-#include <config.h>
-#endif
-#if defined(HAVE_GL_GL_H)
-#include <GL/gl.h>
-#elif defined(HAVE_OPENGL_GL_H)
-#include <OpenGL/gl.h>
-#endif
-#include <wx/wx.h>
-
-void calc_2dbz(int fbz, GLfloat secvec[3], GLfloat secscale, GLfloat axis2d[2][3],
-  int *nbzl2d, int nbragg, GLfloat bragg[26][3], GLfloat brnrm[26], GLfloat bzl2d[26][3],
-  GLfloat bzl2d_proj[26][3]);
-void calc_section(int fbz, int nb, int nthreads, GLfloat secvec[3], GLfloat secscale,
-  GLfloat axis2d[2][3], int* ntri, GLfloat**** kvp, GLfloat** clr, wxTextCtrl* terminal,
-  int* n2d, GLfloat** kv2d, GLfloat** clr2d);
+void calc_2dbz();
+void calc_section();

--- a/src/variable.hpp
+++ b/src/variable.hpp
@@ -60,7 +60,6 @@ extern int interpol;     //!< Ratio of interpolation
 /*
  Switch for some modes
 */
-extern int blackback;   //!< Switch for black background
 extern int color_scale; //!< Switch for full color scale mode
 extern int fbz;         //!< Switch for 1st Brillouin zone mode
 extern int nodeline;    //!< Switch for node lines
@@ -150,6 +149,9 @@ extern GLfloat  yellow[4]; //!< Yellow color code
 extern GLfloat     red[4]; //!< Red color code
 extern GLfloat   green[4]; //!< Green color code
 extern GLfloat    blue[4]; //!< Blue color code
+extern GLfloat BackGroundColor[4];//!< Background color code
+extern GLfloat LineColor[4];//!< Line color code
+extern GLfloat BarColor[5][4];
 /*
  Others
 */