https://github.com/JoeMcEwen/FAST-PT
Raw File
Tip revision: 137418e116f80d982427f13f4aed9273503952a8 authored by Jonathan Blazek on 17 October 2018, 13:18:51 UTC
Added LICENSE
Tip revision: 137418e
OV.py
import numpy as np
from J_table import J_table 
import sys

# Ordering is \alpha, \beta, l_1, l_2, l, A coeficient 
l_mat= np.array([[-2,0,0,0,0,2./3],\
          [-2,0,0,2,0,-2./3],\
          [0,-2,0,0,0,-2./3],\
          [0,-2,0,2,0,2./3]], dtype=float)
def OV():
	table=np.zeros(10,dtype=float)
	for i in range(l_mat.shape[0]):
		x=J_table(l_mat[i])
		table=np.row_stack((table,x))
		
	return table[1:,:] 
	
back to top