Program mix
*
      Implicit None
*
      Real su3(3,3)
      Real id_matrix_01(3,3),su_matrix_01(3,3)
      Real id_matrix_02(3,3),su_matrix_02(3,3)
      Real id_matrix_03(3,3),su_matrix_03(3,3)
      Real id_matrix_04(3,3),su_matrix_04(3,3)
      Real id_matrix_05(3,3),su_matrix_05(3,3)
      Real id_matrix_06(3,3),su_matrix_06(3,3)
      Real id_matrix_07(3,3),su_matrix_07(3,3)
      Real id_matrix_08(3,3),su_matrix_08(3,3)
*
*
*     (1,1)  (1,2) (1,3)
*     (2,1)  (2,2) (2,3)
*     (3,1)  (3,2) (3,3)
*
*     ________________________________________
*
*---Rotate from the ideal mixed states to the SU(3) states
*
      su3(1,1) = sqrt(2./3.)
      su3(1,2) = sqrt(1./3.)
      su3(1,3) = 0.00
      su3(2,1) = sqrt(1./3.)
      su3(2,2) = -sqrt(2./3.)
      su3(2,3) = 0.00
      su3(3,1) = 0.00
      su3(3,2) = 0.00
      su3(3,3) = 1.00
*
*
*
      id_matrix_01(1,1) = -0.91 
      id_matrix_01(1,2) = -0.07
      id_matrix_01(1,3) = 0.40
      id_matrix_01(2,1) = -0.41
      id_matrix_01(2,2) =  0.35
      id_matrix_01(2,3) = -0.84
      id_matrix_01(3,1) = 0.09
      id_matrix_01(3,2) = 0.93
      id_matrix_01(3,3) = 0.36
*
*---Weingarten
*
      id_matrix_02(1,1) = .819
      id_matrix_02(1,2) = .290
      id_matrix_02(1,3) = -.495
      id_matrix_02(2,1) = -0.399
      id_matrix_02(2,2) =  0.908
      id_matrix_02(2,3) = -0.128
      id_matrix_02(3,1) = 0.413
      id_matrix_02(3,2) = 0.302
      id_matrix_02(3,3) = 0.859
*
*---Cheng
*
      id_matrix_03(1,1) =  0.78
      id_matrix_03(1,2) =  0.51
      id_matrix_03(1,3) = -0.36
      id_matrix_03(2,1) = -0.54
      id_matrix_03(2,2) =  0.84
      id_matrix_03(2,3) =  0.03
      id_matrix_03(3,1) =  0.32
      id_matrix_03(3,2) =  0.18
      id_matrix_03(3,3) =  0.93
*
*---Close, Farrar and Li PRD55 (1997)
*
      id_matrix_04(1,1) =  0.86
      id_matrix_04(1,2) =  0.13
      id_matrix_04(1,3) = -0.50
      id_matrix_04(2,1) =  0.43
      id_matrix_04(2,2) = -0.61
      id_matrix_04(2,3) =  0.61
      id_matrix_04(3,1) =  0.22
      id_matrix_04(3,2) =  0.76
      id_matrix_04(3,3) =  0.60
*
*----Giacosa solution1
*
      id_matrix_05(1,1) =  0.86
      id_matrix_05(1,2) =  0.24
      id_matrix_05(1,3) =  0.45
      id_matrix_05(2,1) = -0.45
      id_matrix_05(2,2) = -0.06
      id_matrix_05(2,3) =  0.89
      id_matrix_05(3,1) = -0.24
      id_matrix_05(3,2) =  0.97
      id_matrix_05(3,3) = -0.06
*
*----Giacosa solution2
*
      id_matrix_06(1,1) =  0.81
      id_matrix_06(1,2) =  0.19
      id_matrix_06(1,3) =  0.54
      id_matrix_06(2,1) = -0.49
      id_matrix_06(2,2) =  0.72
      id_matrix_06(2,3) =  0.49
      id_matrix_06(3,1) = -0.30
      id_matrix_06(3,2) =  0.67
      id_matrix_06(3,3) = -0.68

*
*----Giacosa solution3
*
      id_matrix_07(1,1) =  0.79
      id_matrix_07(1,2) =  0.26
      id_matrix_07(1,3) =  0.56
      id_matrix_07(2,1) = -0.58
      id_matrix_07(2,2) =  0.02
      id_matrix_07(2,3) =  0.81
      id_matrix_07(3,1) = -0.20
      id_matrix_07(3,2) =  0.97
      id_matrix_07(3,3) = -0.16

*
*----Giacosa solution4
*
      id_matrix_08(1,1) =  0.82
      id_matrix_08(1,2) =  0.57
      id_matrix_08(1,3) = -0.07
      id_matrix_08(2,1) = -0.57
      id_matrix_08(2,2) =  0.82
      id_matrix_08(2,3) =  0.00
      id_matrix_08(3,1) = -0.06
      id_matrix_08(3,2) =  0.04
      id_matrix_08(3,3) = -0.99

*
      Call Rotate_it(su_matrix_08,id_matrix_08,su3)
*
      Call Print_out(id_matrix_08)
      Call Print_out(su_matrix_08)
*
      Call CheckIt(id_matrix_08)
      Call CheckIt(su_matrix_08)
*
      Stop
      End
*
      Subroutine Rotate_it(C,A,B)
      Real A(3,3),B(3,3),C(3,3)
*
*---Rotate the matrix
*
      C(1,1) = A(1,1) * B(1,1) + A(1,2) * B(1,2) 
      C(1,2) = A(1,1) * B(2,1) + A(1,2) * B(2,2)
      C(1,3) = A(1,3)
* 
      C(2,1) = A(2,1) * B(1,1) + A(2,2) * B(1,2)
      C(2,2) = A(2,1) * B(2,1) + A(2,2) * B(2,2)
      C(2,3) = A(2,3)
* 
      C(3,1) = A(3,1) * B(1,1) + A(3,2) * B(1,2) 
      C(3,2) = A(3,1) * B(2,1) + A(3,2) * B(2,2) 
      C(3,3) = A(3,3) 
*
      Return
      End
*
      Subroutine Print_out(matrix)
      Real matrix(3,3)
*
      Write(6,100) matrix(1,1),matrix(1,2),matrix(1,3)
      Write(6,100) matrix(2,1),matrix(2,2),matrix(2,3)
      Write(6,100) matrix(3,1),matrix(3,2),matrix(3,3)
      Write(6,101)
*
 100  Format(2X,F5.3,2X,F5.3,2X,F5.3)
 101  Format(1X,'_________________________________________')
      Return
      End
*
      Subroutine Checkit(a)
      Real a(3,3)
*
      Sum_1 = a(1,1)**2 + a(2,1)**2 + a(3,1)**2
      Sum_2 = a(1,2)**2 + a(2,2)**2 + a(3,2)**2
      Sum_3 = a(1,3)**2 + a(2,3)**2 + a(3,3)**2
*
      Write(6,100)Sum_1,Sum_2,Sum_3

      Sum_1 = a(1,1)**2 + a(1,2)**2 + a(1,3)**2
      Sum_2 = a(2,1)**2 + a(2,2)**2 + a(2,3)**2
      Sum_3 = a(3,1)**2 + a(3,2)**2 + a(3,3)**2
*
      Write(6,100)Sum_1,Sum_2,Sum_3


 100  Format(1X,F5.2,2X,F5.2,2X,F5.2)
*
      Return
      End