https://github.com/JoeMcEwen/FAST-PT
Raw File
Tip revision: c9e831df408f3a9d83a9c34c7423616ae1de4179 authored by Jonathan Blazek on 24 October 2018, 14:13:07 UTC
Revert "ENH python packaging"
Tip revision: c9e831d
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