#define PROG_NAME "TestCornerFE"
#define PROG_DESC "???"
#define PROG_VERS "1.0"

#define TestCornerFE_C_COPYRIGHT \
  ""
  
#define PROG_INFO \
  "" \
  " "


#include <Map.h>
#define _GNU_SOURCE
#include <stdio.h>
#include <Octf.h>
#include <Triangulation.h>


#include <Octf.h>
//  Clock, Spin, Srot, Tors;


#include <Map.h>
//  Corner, CCorner;

TYPE
double Place_t = Octf_Place_t;
  
VAR
  Place_t ca.el[8];

PROCEDURE PrintCorner()

 {
  for (l = 1; l < 5; l++)
  {
    ??? a = Map.MakeWedge(l);
    {
      /*
      Octf.PrintPlace(stderr, a, 1);
      fprintf(stderr, "\n");
      */
      fprintf(stderr, "(" & Fmt.Int(l) & "," & Fmt.Int(l)
                  & "):" & "\n");
      fprintf(stderr, "Corner(a)== ");
      Octf.PrintPlace(stderr, Corner(a), 1); 

      fprintf(stderr, "CornerSpin(a)== ");
      Octf.PrintPlace(stderr, Corner(Spin(a)), 1); 

      fprintf(stderr, "CornerSrot(a)== ");
      Octf.PrintPlace(stderr, Corner(Srot(a)), 1); 

      fprintf(stderr, "CornerSpinSrot(a)== ");
      Octf.PrintPlace(stderr, Corner(Spin(Srot(a))), 1); 

      fprintf(stderr, "CornerClock(a)== ");
      Octf.PrintPlace(stderr, Corner(Clock(a)), 1); 

      fprintf(stderr, "CornerSpinClock(a)== ");
      Octf.PrintPlace(stderr, Corner(Spin(Clock(a))), 1);

      fprintf(stderr, "CornerTor(a)== ");
      Octf.PrintPlace(stderr, Corner(Tors(a)), 1); 

      fprintf(stderr, "CornerSpinTors(a)== ");
      Octf.PrintPlace(stderr, Corner(Spin(Tors(a))), 1); 
      
    }
  }
 } PrintCorner;


PROCEDURE PrintCCorner()

 {
  for (l = 1; l < 5; l++)
  {
    ??? a = Map.MakeWedge(l);
    {
      fprintf(stderr, "(" & Fmt.Int(l) & "," & Fmt.Int(l)
                  & "):" & "\n");
      fprintf(stderr, "CCorner(a)== ");
      Octf.PrintPlace(stderr, CCorner(a), 1); 

      fprintf(stderr, "CCornerSpin(a)== ");
      Octf.PrintPlace(stderr, CCorner(Spin(a)), 1); 

      fprintf(stderr, "CCornerSrot(a)== ");
      Octf.PrintPlace(stderr, CCorner(Srot(a)), 1); 

      fprintf(stderr, "CCornerSpinSrot(a)== ");
      Octf.PrintPlace(stderr, CCorner(Spin(Srot(a))), 1); 

      fprintf(stderr, "CCornerClock(a)== ");
      Octf.PrintPlace(stderr, CCorner(Clock(a)), 1); 

      fprintf(stderr, "CCornerSpinClock(a)== ");
      Octf.PrintPlace(stderr, CCorner(Spin(Clock(a))), 1); 

      fprintf(stderr, "CCornerTors(a)== ");
      Octf.PrintPlace(stderr, CCorner(Tors(a)), 1); 

      fprintf(stderr, "CCornerSpinTors(a)== ");
      Octf.PrintPlace(stderr, CCorner(Spin(Tors(a))), 1); 
      
    }
  }
 } PrintCCorner;

 void Teste()
 
  { for (i = 1; i < 5; i++){
      ca = MakeTetraTopo(i,i);
      fprintf(stderr, "(" & Fmt.Int(i) & "," & Fmt.Int(i)
                  & "):" & "\n");
      for (k = 0; k < 8; k++)
      {
        fprintf(stderr, "c[" & Fmt.Int(k) & "]== ");
        Octf.PrintPlace(stderr, ca.el[k],3);
      }
    }
 } Teste;
 

{
  PrintCorner();
  PrintCCorner(); 
  Teste();
} TestCorner.


  /****************************************************************************
(1,1):
Corne r(a)8:2:1
CornerSpi n(a)5:2:0
CornerSro t(a)0:0:0
CornerSpinSro t(a)9:0:1
CornerCloc k(a)5:0:1
CornerSpinCloc k(a)8:0:0
CornerTo r(a)9:2:0
CornerSpinTor s(a)0:2:1
(2,2):
Corne r(a)17:2:1
CornerSpi n(a)8:2:0
CornerSro t(a)0:0:0
CornerSpinSro t(a)27:0:1
CornerCloc k(a)23:0:1
CornerSpinCloc k(a)32:0:0
CornerTo r(a)33:2:0
CornerSpinTor s(a)3:2:1
(3,3):
Corne r(a)26:2:1
CornerSpi n(a)11:2:0
CornerSro t(a)0:0:0
CornerSpinSro t(a)57:0:1
CornerCloc k(a)53:0:1
CornerSpinCloc k(a)68:0:0
CornerTo r(a)69:2:0
CornerSpinTor s(a)6:2:1
(4,4):
Corne r(a)35:2:1
CornerSpi n(a)14:2:0
CornerSro t(a)0:0:0
CornerSpinSro t(a)99:0:1
CornerCloc k(a)95:0:1
CornerSpinCloc k(a)116:0:0
CornerTo r(a)117:2:0
CornerSpinTor s(a)9:2:1
(1,1):
CCorne r(a)5:2:1
CCornerSpi n(a)8:2:0
CCornerSro t(a)9:0:0
CCornerSpinSro t(a)0:0:1
CCornerCloc k(a)8:0:1
CCornerSpinCloc k(a)5:0:0
CCornerTor s(a)0:2:0
CCornerSpinTor s(a)9:2:1
(2,2):
CCorne r(a)8:2:1
CCornerSpi n(a)17:2:0
CCornerSro t(a)27:0:0
CCornerSpinSro t(a)0:0:1
CCornerCloc k(a)32:0:1
CCornerSpinCloc k(a)23:0:0
CCornerTor s(a)3:2:0
CCornerSpinTor s(a)33:2:1
(3,3):
CCorne r(a)11:2:1
CCornerSpi n(a)26:2:0
CCornerSro t(a)57:0:0
CCornerSpinSro t(a)0:0:1
CCornerCloc k(a)68:0:1
CCornerSpinCloc k(a)53:0:0
CCornerTor s(a)6:2:0
CCornerSpinTor s(a)69:2:1
(4,4):
CCorne r(a)14:2:1
CCornerSpi n(a)35:2:0
CCornerSro t(a)99:0:0
CCornerSpinSro t(a)0:0:1
CCornerCloc k(a)116:0:1
CCornerSpinCloc k(a)95:0:0
CCornerTor s(a)9:2:0
CCornerSpinTor s(a)117:2:1
(1,1):
c[0]==   8:2:1
c[1]==   5:2:0
c[2]==   0:0:0
c[3]==   9:0:1
c[4]==   5:0:1
c[5]==   8:0:0
c[6]==   9:2:0
c[7]==   0:2:1
(2,2):
c[0]==  17:2:1
c[1]==   8:2:0
c[2]==   0:0:0
c[3]==  27:0:1
c[4]==  23:0:1
c[5]==  32:0:0
c[6]==  33:2:0
c[7]==   3:2:1
(3,3):
c[0]==  26:2:1
c[1]==  11:2:0
c[2]==   0:0:0
c[3]==  57:0:1
c[4]==  53:0:1
c[5]==  68:0:0
c[6]==  69:2:0
c[7]==   6:2:1
(4,4):
c[0]==  35:2:1
c[1]==  14:2:0
c[2]==   0:0:0
c[3]==  99:0:1
c[4]==  95:0:1
c[5]== 116:0:0
c[6]== 117:2:0
c[7]==   9:2:1

***************************************************************************/