PMNS_SNSI.cxx File Reference

#include "PMNS_SNSI.h"

Go to the source code of this file.

Functions
void	HermitianSquare (complexD(&A)[3][3])

Function Documentation

HermitianSquare()

void HermitianSquare

(

complexD(&)

A[3][3]

)

Square a Hermitian matrix

Parameters

A	- input matrix A

Definition at line 98 of file PMNS_SNSI.cxx.

{
  complexD tmp[3][3];
  for (int i = 0; i < 3; i++) {
    for (int j = i; j < 3; j++) {
      tmp[i][j] = A[i][j];
      if (i < j) tmp[j][i] = conj(A[i][j]);
      A[i][j] = 0;
    }
  }
 
  for (int i = 0; i < 3; i++) {
    for (int j = i; j < 3; j++) {
      for (int k = 0; k < 3; k++) { A[i][j] += tmp[i][k] * tmp[k][j]; }
    }
  }
}

Referenced by OscProb::PMNS_SNSI::UpdateHam().