src/expde/formulas/subdivIV.h

Go to the documentation of this file.

//    expde: expression templates for partial differential equations.
//    Copyright (C) 2001  Christoph Pflaum
//    This program is free software; you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation; either version 2 of the License, or
//    (at your option) any later version.
//
//    This program is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//    GNU General Public License for more details.
//
//                 SEE  Notice1.doc made by 
//                 LAWRENCE LIVERMORE NATIONAL LABORATORY
//

// ------------------------------------------------------------
//
// subdivIV.h
//
// ------------------------------------------------------------

#ifndef SUBDIVIV_H_
#define SUBDIVIV_H_

void Act_on_subdivision_partIV(Bo_description* Description, 
                                Calculator_FE* Act,
                                bool* corners, bool& stop) {
  int i;

  Edge_Corner_point poi_1,   poi_2,   poi_3,   poi_4;
  Edge_Corner_point poi_5,   poi_6,   poi_7,   poi_8,   poi_14;
  Edge_Corner_point poi_9,   poi_10,  poi_11,  poi_12,  poi_13;
  Edge_Corner_point poi_A,   poi_B,   poi_C,   poi_D;
  Edge_Corner_point poi_A2,  poi_B2,  poi_C2,  poi_D2;
  Edge_Corner_point poi_boA, poi_boB, poi_boC, poi_boD;

  dir_3D    dirA, dirB;
  dir_3D sur, suropp;
  bool later_stop, LA, LB, LC, LD, LcA, LcB, LcC, LcD, Lsp;
  int num_corners;
  double h_7, h_9, h_11, h_13;

  later_stop = false;

// stop = false;
for(i=0;i<6 && !stop;++i) {
  sur    = (dir_3D)i;
  if(i==Wdir || i==Ndir || i==Ddir) {   
    poi_A.corner = corner_sur(sur,SWdir);
    poi_B.corner = corner_sur(sur,SEdir);
    poi_C.corner = corner_sur(sur,NWdir);
    poi_D.corner = corner_sur(sur,NEdir);
  }
  else {
    poi_A.corner = corner_sur(sur,SWdir);
    poi_B.corner = corner_sur(sur,NWdir);
    poi_C.corner = corner_sur(sur,SEdir);
    poi_D.corner = corner_sur(sur,NEdir);
  }

  LcA = Description->corner_point(poi_A.corner);
  LcB = Description->corner_point(poi_B.corner);
  LcC = Description->corner_point(poi_C.corner);
  LcD = Description->corner_point(poi_D.corner);

  num_corners=0;
  if(LcA) ++num_corners;
  if(LcB) ++num_corners;
  if(LcC) ++num_corners;
  if(LcD) ++num_corners;
  if(num_corners == 3) {
    // Flaeche vom Typ d)
    poi_A.edge_point = corner_poi_typ;
    poi_B.edge_point = corner_poi_typ;
    poi_C.edge_point = corner_poi_typ;
    poi_D.edge_point = corner_poi_typ;
    
    poi_A2.edge_point = corner_poi_typ;
    poi_B2.edge_point = corner_poi_typ;
    poi_C2.edge_point = corner_poi_typ;
    poi_D2.edge_point = corner_poi_typ;
    
    poi_boA.edge_point = edge_poi_typ;
    poi_boB.edge_point = edge_poi_typ;
    poi_boC.edge_point = edge_poi_typ;
    poi_boD.edge_point = edge_poi_typ;
    suropp = opposite3D((dir_3D)i);
    if(i==Wdir || i==Ndir || i==Ddir) { 
      poi_A2.corner = corner_sur(suropp,SWdir);
      poi_B2.corner = corner_sur(suropp,SEdir);
      poi_C2.corner = corner_sur(suropp,NWdir);
      poi_D2.corner = corner_sur(suropp,NEdir);
      
      poi_boA.corner = poi_A.corner;   poi_boA.d = opposite3D((dir_3D)i);
      LA = Description->edge_point(poi_boA);
      poi_boB.corner = poi_B.corner;   poi_boB.d = opposite3D((dir_3D)i);
      LB = Description->edge_point(poi_boB);
      poi_boC.corner = poi_C.corner;   poi_boC.d = opposite3D((dir_3D)i);
      LC = Description->edge_point(poi_boC);
      poi_boD.corner = poi_D.corner;   poi_boD.d = opposite3D((dir_3D)i);
      LD = Description->edge_point(poi_boD);
    }
    else {
      poi_A2.corner = corner_sur(suropp,SWdir);
      poi_B2.corner = corner_sur(suropp,NWdir);
      poi_C2.corner = corner_sur(suropp,SEdir);
      poi_D2.corner = corner_sur(suropp,NEdir);
      
      poi_boA.corner = poi_A.corner;   poi_boA.d = opposite3D((dir_3D)i);
      LA = Description->edge_point(poi_boA);
      poi_boB.corner = poi_B.corner;   poi_boB.d = opposite3D((dir_3D)i);
      LB = Description->edge_point(poi_boB);
      poi_boC.corner = poi_C.corner;   poi_boC.d = opposite3D((dir_3D)i);
      LC = Description->edge_point(poi_boC);
      poi_boD.corner = poi_D.corner;   poi_boD.d = opposite3D((dir_3D)i);
      LD = Description->edge_point(poi_boD);
    }

    // Drehen:
    if(!LcC) {
      // D B  ->  // C D
      // C A  ->  // A B
      Lsp = LC;
      LC = LD;     LD = LB;
      LB = LA;     LA = Lsp;    

      poi_1 = poi_C;
      poi_C = poi_D;     poi_D = poi_B;
      poi_B = poi_A;     poi_A = poi_1;     
 
      poi_1 = poi_C2;
      poi_C2 = poi_D2;     poi_D2 = poi_B2;
      poi_B2 = poi_A2;     poi_A2 = poi_1;

      poi_1   = poi_boC;
      poi_boC = poi_boD;     poi_boD = poi_boB;
      poi_boB = poi_boA;     poi_boA = poi_1;     
    }
    if(!LcD) {
      // B A ->  // C D
      // D C ->  // A B
      Lsp = LD;
      LD = LA;      LA = Lsp;
      Lsp = LB;
      LB = LC;      LC = Lsp;

      poi_1 = poi_D;
      poi_D = poi_A;      poi_A = poi_1;
      poi_1 = poi_B;
      poi_B = poi_C;      poi_C = poi_1;

      poi_1  = poi_D2;
      poi_D2 = poi_A2;      poi_A2 = poi_1;
      poi_1  = poi_B2;
      poi_B2 = poi_C2;      poi_C2 = poi_1;

      poi_1   = poi_boD;
      poi_boD = poi_boA;      poi_boA = poi_1;
      poi_1   = poi_boB;
      poi_boB = poi_boC;      poi_boC = poi_1;
    }
    if(!LcB) {
      // A C  ->  // C D
      // B D  ->  // A B
      Lsp = LB;
      LB = LD;     LD = LC;
      LC = LA;     LA = Lsp;

      poi_1 = poi_B;
      poi_B = poi_D;     poi_D = poi_C;
      poi_C = poi_A;     poi_A = poi_1;

      poi_1  = poi_B2;
      poi_B2 = poi_D2;     poi_D2 = poi_C2;
      poi_C2 = poi_A2;     poi_A2 = poi_1;

      poi_1   = poi_boB;
      poi_boB = poi_boD;     poi_boD = poi_boC;
      poi_boC = poi_boA;     poi_boA = poi_1;
    }

    // Jetzt immer:
    // C D
    // A B

    poi_2.edge_point = corner_poi_typ;
    poi_3.edge_point = poi_4.edge_point = corner_poi_typ;
    poi_5.edge_point = poi_6.edge_point = edge_poi_typ;

    poi_2 = poi_C;        poi_3 = poi_D;
                          poi_4 = poi_B;

    poi_5.corner = poi_2.corner;   
    poi_5.d = direction_from_to(poi_C.corner,poi_A.corner);
    poi_6.corner = poi_4.corner;   
    poi_6.d = direction_from_to(poi_B.corner,poi_A.corner);

    //  Boundary cell : 2.3.
    // ----------------------
    poi_boA.corner = poi_B2.corner;   
    poi_boA.d = direction_from_to(poi_B.corner,poi_A.corner);
    LA = Description->edge_point(poi_boA);
    // Symmetrie a):
    if(LA && LD && LC && !LB &&
       corners[poi_B.corner] && corners[poi_C.corner] &&
       corners[poi_D.corner] && corners[poi_B2.corner]) {
      corners[poi_B.corner] = corners[poi_C.corner] =
        corners[poi_D.corner] = corners[poi_B2.corner] = false;
      if(later_stop) {
        stop=true;
        cout << "\n Achtung:" << endl;
      }
      later_stop=true;
      poi_7.edge_point = poi_9.edge_point = edge_poi_typ;
      poi_8.edge_point = corner_poi_typ;
      poi_10.edge_point = poi_11.edge_point = edge_poi_typ;
      // C D
      // A B
      poi_11 = poi_boC;   poi_10 = poi_boD;
      poi_8  = poi_B2;
      
      poi_7.corner = poi_8.corner;   
      poi_9.corner = poi_8.corner;   
      poi_7.d = direction_from_to(poi_B.corner,poi_A.corner);
      poi_9.d = direction_from_to(poi_B.corner,poi_D.corner);

      poi_12.edge_point = cell_poi_typ;

      
      if(poi_9.d==Wdir || poi_9.d==Edir)
        Act->Add_Bo_freedom(D3vector( 0.0,
                                     -0.25+0.5*y1DCoord(poi_4.corner),
                                     -0.25+0.5*z1DCoord(poi_4.corner)));
      if(poi_9.d==Ndir || poi_9.d==Sdir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_4.corner),
                                      0.0,
                                     -0.25+0.5*z1DCoord(poi_4.corner)));
      if(poi_9.d==Tdir || poi_9.d==Ddir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_4.corner),
                                     -0.25+0.5*y1DCoord(poi_4.corner),
                                      0.0));
      /*
      if(poi_9.d==Wdir || poi_9.d==Edir)
        Act->Add_Bo_freedom(D3vector(-0.05+0.1*x1DCoord(poi_4.corner),
                                     -0.25+0.5*y1DCoord(poi_4.corner),
                                     -0.25+0.5*z1DCoord(poi_4.corner)));
      if(poi_9.d==Ndir || poi_9.d==Sdir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_4.corner),
                                     -0.05+0.1*y1DCoord(poi_4.corner),
                                     -0.25+0.5*z1DCoord(poi_4.corner)));
      if(poi_9.d==Tdir || poi_9.d==Ddir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_4.corner),
                                     -0.25+0.5*y1DCoord(poi_4.corner),
                                     -0.05+0.1*z1DCoord(poi_4.corner)));
      */
      // D Flaeche
      Act->Tetraeder(poi_12,poi_5,poi_6,poi_3);
      Act->Tetraeder(poi_12,poi_2,poi_5,poi_3);
      Act->Tetraeder(poi_12,poi_6,poi_4,poi_3);

      // E Flaeche
      Act->Tetraeder(poi_12,poi_10,poi_4,poi_9);
      Act->Tetraeder(poi_12,poi_10,poi_3,poi_4);
      Act->Tetraeder(poi_12,poi_8,poi_9,poi_4);

      // W Flaeche
      Act->Tetraeder(poi_12,poi_11,poi_5,poi_2);

      // T Flaeche
      Act->Tetraeder(poi_12,poi_8,poi_7,poi_9);

      // N Flaeche
      // Anstelle: Act->Pyramid(poi_12,poi_10,poi_11,poi_2,poi_3);
      dirA = poi_10.d;
      dirB = direction_from_to(poi_3.corner,poi_2.corner);
      if(Calc_f_type(dirA,dirB)==an) {
        Act->Tetraeder(poi_12,poi_10,poi_2,poi_3);
        Act->Tetraeder(poi_12,poi_10,poi_11,poi_2);
      }
      else {
        Act->Tetraeder(poi_12,poi_11,poi_3,poi_10);
        Act->Tetraeder(poi_12,poi_11,poi_2,poi_3);
      }
      // S Flaeche
      // Anstelle: Act->Pyramid(poi_12,poi_7,poi_8,poi_4,poi_6);
      if(Calc_f_type(dirA,dirB)==an) {
        Act->Tetraeder(poi_12,poi_8,poi_4,poi_6);
        Act->Tetraeder(poi_12,poi_8,poi_6,poi_7);
      }
      else {
        Act->Tetraeder(poi_12,poi_8,poi_4,poi_7);
        Act->Tetraeder(poi_12,poi_7,poi_4,poi_6);
      }
      // Oberflaeche:
      Act->Tetraeder(poi_12,poi_7,poi_11,poi_9);
      Act->Tetraeder(poi_12,poi_9,poi_11,poi_10);
      // Anstelle: Act->Pyramid(poi_12,poi_7,poi_6,poi_5,poi_11);
      Act->Tetraeder(poi_12,poi_7,poi_5,poi_11);
      Act->Tetraeder(poi_12,poi_7,poi_6,poi_5);
    }
    poi_boA.corner = poi_C2.corner;   
    poi_boA.d = direction_from_to(poi_C.corner,poi_A.corner);
    LA = Description->edge_point(poi_boA);
    // Symmetrie b):
    if(LA && LD && !LC && LB && !stop &&
       corners[poi_B.corner] && corners[poi_C.corner] &&
       corners[poi_D.corner] && corners[poi_C2.corner]) {
      corners[poi_B.corner] = corners[poi_C.corner] =
        corners[poi_D.corner] = corners[poi_C2.corner] = false;
      if(later_stop) stop=true;
      later_stop=true;
      poi_7.edge_point = poi_9.edge_point = edge_poi_typ;
      poi_8.edge_point = corner_poi_typ;
      poi_10.edge_point = poi_11.edge_point = edge_poi_typ;
      // C D
      // A B
      poi_8 = poi_C2;     poi_10 = poi_boD;
                          poi_11 = poi_boB;
      
      poi_7.corner = poi_8.corner;   
      poi_9.corner = poi_8.corner;   
      poi_7.d = direction_from_to(poi_C.corner,poi_A.corner);
      poi_9.d = direction_from_to(poi_C.corner,poi_D.corner);

      poi_12.edge_point = cell_poi_typ;

      if(poi_9.d==Wdir || poi_9.d==Edir)
        Act->Add_Bo_freedom(D3vector( 0.0,
                                     -0.25+0.5*y1DCoord(poi_2.corner),
                                     -0.25+0.5*z1DCoord(poi_2.corner)));
      if(poi_9.d==Ndir || poi_9.d==Sdir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_2.corner),
                                      0.0,
                                     -0.25+0.5*z1DCoord(poi_2.corner)));
      if(poi_9.d==Tdir || poi_9.d==Ddir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_2.corner),
                                     -0.25+0.5*y1DCoord(poi_2.corner),
                                      0.0));
      /*
      if(poi_9.d==Wdir || poi_9.d==Edir)
        Act->Add_Bo_freedom(D3vector(-0.05+0.1*x1DCoord(poi_2.corner),
                                     -0.25+0.5*y1DCoord(poi_2.corner),
                                     -0.25+0.5*z1DCoord(poi_2.corner)));
      if(poi_9.d==Ndir || poi_9.d==Sdir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_2.corner),
                                     -0.05+0.1*y1DCoord(poi_2.corner),
                                     -0.25+0.5*z1DCoord(poi_2.corner)));
      if(poi_9.d==Tdir || poi_9.d==Ddir)
        Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_2.corner),
                                     -0.25+0.5*y1DCoord(poi_2.corner),
                                     -0.05+0.1*z1DCoord(poi_2.corner)));
      */    
      // D Flaeche
      Act->Tetraeder(poi_12,poi_5,poi_6,poi_3);
      Act->Tetraeder(poi_12,poi_2,poi_5,poi_3);
      Act->Tetraeder(poi_12,poi_6,poi_4,poi_3);

      // N Flaeche
      Act->Tetraeder(poi_12,poi_10,poi_2,poi_3);
      Act->Tetraeder(poi_12,poi_10,poi_9,poi_2);
      Act->Tetraeder(poi_12,poi_9,poi_8,poi_2);

      // S Flaeche
      Act->Tetraeder(poi_12,poi_11,poi_4,poi_6);

      // T Flaeche
      Act->Tetraeder(poi_12,poi_7,poi_8,poi_9);

      // E Flaeche
      // Anstelle: Act->Pyramid(poi_12,poi_10,poi_11,poi_2,poi_3);
      dirA = poi_10.d;
      dirB = direction_from_to(poi_3.corner,poi_4.corner);
      if(Calc_f_type(dirA,dirB)==an) {
        Act->Tetraeder(poi_12,poi_10,poi_3,poi_4);
        Act->Tetraeder(poi_12,poi_10,poi_4,poi_11);
      }
      else {
        Act->Tetraeder(poi_12,poi_11,poi_10,poi_3);
        Act->Tetraeder(poi_12,poi_11,poi_3,poi_4);
      }

      // W Flaeche
      // Anstelle: Act->Pyramid(poi_12,poi_7,poi_8,poi_4,poi_6);
      if(Calc_f_type(dirA,dirB)==an) {
        Act->Tetraeder(poi_12,poi_8,poi_5,poi_2);
        Act->Tetraeder(poi_12,poi_8,poi_7,poi_5);
      }
      else {
        Act->Tetraeder(poi_12,poi_7,poi_2,poi_8);
        Act->Tetraeder(poi_12,poi_7,poi_5,poi_2);
      }
               
      // Oberflaeche:
      Act->Tetraeder(poi_12,poi_11,poi_7,poi_9);
      Act->Tetraeder(poi_12,poi_11,poi_9,poi_10);
      // Anstelle: Act->Pyramid(poi_12,poi_7,poi_6,poi_5,poi_11);
      // falsch:
      //   Act->Tetraeder(poi_12,poi_11,poi_5,poi_7);
      //   Act->Tetraeder(poi_12,poi_11,poi_6,poi_5);
      Act->Tetraeder(poi_12,poi_11,poi_6,poi_7);
      Act->Tetraeder(poi_12,poi_7,poi_6,poi_5);
    }
    //  Boundary cell : 2.4.
    // ----------------------
    poi_boA.corner = poi_C2.corner;   
    poi_boA.d = direction_from_to(poi_C.corner,poi_A.corner);
    LA = Description->edge_point(poi_boA);
    if(LA && LD && !LC && !LB) {
      stop=true;
      poi_7.edge_point = poi_9.edge_point = edge_poi_typ;
      poi_8.edge_point = corner_poi_typ;
      poi_10.edge_point = edge_poi_typ;
      poi_12.edge_point = corner_poi_typ;
      poi_11.edge_point = poi_13.edge_point = edge_poi_typ;
      // C D
      // A B
      poi_12 = poi_C2;    poi_10 = poi_boD;
                          poi_8  = poi_B2;
      
      poi_7.corner = poi_8.corner;   
      poi_9.corner = poi_8.corner;   
      poi_7.d = direction_from_to(poi_B.corner,poi_A.corner);
      poi_9.d = direction_from_to(poi_B.corner,poi_D.corner);

      poi_11.corner = poi_12.corner;   
      poi_13.corner = poi_12.corner;   
      poi_11.d = direction_from_to(poi_C.corner,poi_D.corner);
      poi_13.d = direction_from_to(poi_C.corner,poi_A.corner);

      poi_14.edge_point = cell_poi_typ;
      /*  old version
      Act->Add_Bo_freedom(D3vector(-0.3+0.6*x1DCoord(poi_3.corner),
                                   -0.3+0.6*y1DCoord(poi_3.corner),
                                   -0.3+0.6*z1DCoord(poi_3.corner)));
      */
      
      Act->Add_Bo_freedom(D3vector(-0.1+0.2*x1DCoord(poi_3.corner),
                                   -0.1+0.2*y1DCoord(poi_3.corner),
                                   -0.1+0.2*z1DCoord(poi_3.corner)));
      
      /*
      Act->Add_Bo_freedom(D3vector(-0.25+0.5*x1DCoord(poi_3.corner),
                                   -0.25+0.5*y1DCoord(poi_3.corner),
                                   -0.25+0.5*z1DCoord(poi_3.corner)));
      */
      // D Flaeche
      Act->Tetraeder(poi_14,poi_5,poi_6,poi_3);
      Act->Tetraeder(poi_14,poi_2,poi_5,poi_3);
      Act->Tetraeder(poi_14,poi_6,poi_4,poi_3);

      // E Flaeche
      Act->Tetraeder(poi_14,poi_10,poi_4,poi_9);
      Act->Tetraeder(poi_14,poi_10,poi_3,poi_4);
      Act->Tetraeder(poi_14,poi_8,poi_9,poi_4);

      // N Flaeche
      Act->Tetraeder(poi_14,poi_10,poi_11,poi_2);
      Act->Tetraeder(poi_14,poi_10,poi_2,poi_3);
      Act->Tetraeder(poi_14,poi_11,poi_12,poi_2);

      // T Flaeche
      Act->Tetraeder(poi_14,poi_8,poi_7,poi_9);
      Act->Tetraeder(poi_14,poi_13,poi_12,poi_11);

      // S Flaeche
      // Anstelle: Act->Pyramid(poi_14,poi_8,poi_7,poi_6,poi_4);
      dirA = poi_10.d;
      dirB = poi_6.d;
      if(Calc_f_type(dirA,dirB)==an) {
        Act->Tetraeder(poi_14,poi_8,poi_4,poi_6);
        Act->Tetraeder(poi_14,poi_8,poi_6,poi_7);
      }
      else {
        Act->Tetraeder(poi_14,poi_7,poi_8,poi_4);
        Act->Tetraeder(poi_14,poi_7,poi_4,poi_6);
      }
      // W Flaeche
      // Anstelle: Act->Pyramid(poi_12,poi_13,poi_5,poi_2,poi_12);
      dirA = poi_10.d;
      dirB = poi_5.d;
      if(Calc_f_type(dirA,dirB)==an) {
        Act->Tetraeder(poi_14,poi_12,poi_5,poi_2);
        Act->Tetraeder(poi_14,poi_12,poi_13,poi_5);
      }
      else {
        Act->Tetraeder(poi_14,poi_12,poi_13,poi_2);
        Act->Tetraeder(poi_14,poi_13,poi_5,poi_2);
      }

      // Oberflaeche:
      Act->Tetraeder(poi_14,poi_9,poi_11,poi_10);
      // Anstelle: Act->Pyramid(poi_14,poi_7,poi_6,poi_5,poi_13);
      Act->Tetraeder(poi_14,poi_7,poi_5,poi_13);
      Act->Tetraeder(poi_14,poi_7,poi_6,poi_5);

      h_9 = Description->h(poi_9.corner, poi_9.d);
      h_7 = Description->h(poi_7.corner, poi_7.d);
      h_11 = Description->h(poi_11.corner, poi_11.d);
      h_13 = Description->h(poi_13.corner, poi_13.d);

      double H; 
      double angle_9;
      double angle_7;
      double angle_11;
      double angle_13;
      
      H = Description->Meshsize();
      
      angle_9  = calc_special_angle(h_7,h_9,H-h_11,H-h_9);
      angle_7  = calc_special_angle(H-h_13,H-h_7,h_9,h_7);
      angle_13 = calc_special_angle(h_11,h_13,H-h_7,H-h_13);
      angle_11  = calc_special_angle(H-h_9,H-h_11,h_13,h_11);
      
      /*
      if( (h_7-h_11) < (h_9-h_13) ) {
        h_min = (h_7-h_11);
        num_min = 9;
      }
      else {
        h_min = (h_9-h_13);
        num_min = 7;
      }
      if( h_min > (h_13-h_9) ) {
        h_min = (h_13-h_9);
        num_min = 11;
      }
      if( h_min > (h_11-h_7) ) {
        h_min = (h_11-h_7);
      }
      */      

      // Anstelle von Pyr 14, 7, 13, 11, 9
      //      if(num_min==9 || num_min==13) {
      if(parallel_version) 
        Tell_about_round_of_problem(MAX(angle_7,angle_11),
                                    MAX(angle_13,angle_9));
      if(MAX(angle_7,angle_11) > MAX(angle_13,angle_9)) {
        Act->Tetraeder(poi_14,poi_7,poi_13,poi_11);
        Act->Tetraeder(poi_14,poi_7,poi_11,poi_9);
      }
      else {
        Act->Tetraeder(poi_14,poi_7,poi_13,poi_9);
        Act->Tetraeder(poi_14,poi_9,poi_13,poi_11);
      }
    }
    
    //  Boundary cell : 2.6.
    // ----------------------
    poi_boA.corner = poi_C2.corner;   
    poi_boA.d = direction_from_to(poi_C.corner,poi_A.corner);
    LA = Description->edge_point(poi_boA);
    if(!LA && LD && !LC && !LB) {
      cout << "\n Fall 2.6 Randzelle. Noch nicht getestet! " << endl;

      stop=true;
      poi_7.edge_point = poi_9.edge_point = edge_poi_typ;
      poi_8.edge_point = poi_10.edge_point = corner_poi_typ;
      poi_13.edge_point = corner_poi_typ;
      poi_11.edge_point = poi_12.edge_point = edge_poi_typ;
      // C D
      // A B
      poi_13 = poi_C2;    poi_9 = poi_boD;
      poi_10  = poi_A2;   poi_8  = poi_B2;
      
      poi_7.corner = poi_10.corner;   
      poi_7.d = direction_from_to(poi_B2.corner,poi_B.corner);

      poi_11.corner = poi_8.corner;   
      poi_12.corner = poi_13.corner;   
      poi_11.d = direction_from_to(poi_B.corner,poi_D.corner);
      poi_12.d = direction_from_to(poi_C.corner,poi_D.corner);

      poi_14.edge_point = cell_poi_typ;
      Act->Add_Bo_freedom(D3vector(0.0,0.0,0.0));

      // D Flaeche
      Act->Tetraeder(poi_14,poi_5,poi_6,poi_3);
      Act->Tetraeder(poi_14,poi_2,poi_5,poi_3);
      Act->Tetraeder(poi_14,poi_6,poi_4,poi_3);

      // E Flaeche
      Act->Tetraeder(poi_14,poi_11,poi_9,poi_4);
      Act->Tetraeder(poi_14,poi_9,poi_3,poi_4);
      Act->Tetraeder(poi_14,poi_8,poi_11,poi_4);

      // T Flaeche
      Act->Tetraeder(poi_14,poi_12,poi_10,poi_13);
      Act->Tetraeder(poi_14,poi_12,poi_11,poi_10);
      Act->Tetraeder(poi_14,poi_11,poi_8,poi_10);

      // N Flaeche
      Act->Tetraeder(poi_14,poi_12,poi_2,poi_9);
      Act->Tetraeder(poi_14,poi_12,poi_13,poi_2);
      Act->Tetraeder(poi_14,poi_9,poi_2,poi_3);

      // S Flaeche
      Act->Tetraeder(poi_14,poi_8,poi_6,poi_7);
      Act->Tetraeder(poi_14,poi_8,poi_4,poi_6);
      Act->Tetraeder(poi_14,poi_8,poi_7,poi_10);

      // S Flaeche
      Act->Tetraeder(poi_14,poi_13,poi_7,poi_5);
      Act->Tetraeder(poi_14,poi_13,poi_5,poi_2);
      Act->Tetraeder(poi_14,poi_13,poi_10,poi_7);

      // Oberflaeche
      Act->Tetraeder(poi_14,poi_7,poi_6,poi_5);
      Act->Tetraeder(poi_14,poi_11,poi_12,poi_9);
    }
  }
}

}

#endif
   

---