Skip to content
Snippets Groups Projects
Commit cd2aeece authored by Edoardo Pasetto's avatar Edoardo Pasetto
Browse files

Delete utility.py

parent 7a22b71d
No related branches found
No related tags found
No related merge requests found
# -*- coding: utf-8 -*-
"""
Created on Wed Mar 23 11:44:13 2022
@author: Utente
"""
import sympy as sp
import numpy as np
bw_500_val=0.00288
b1_bw_val=0.5
b1_bph_val=0.002
b1_bsed_val=0.02
ased_p_val=0.042
sx_val=0.011
sy_val=0.014
def get_aphy(A,B,C):
a,b,c=sp.symbols('a b c')
exp=a*(c**(1-b))
aphy=float(exp.subs([(a,A),(b,B),(c,C)]))
return aphy
def get_ased(X_val,wl_val):
X,w,ased_p,sx=sp.symbols('X w ased_p sx')
exp=X*ased_p*(sp.exp(-sx*(w-440)))
ased=float(exp.subs([(X,X_val),(w,wl_val),(ased_p,ased_p_val),(sx,sx_val)]))
return ased
def get_ayel(Y_val,wl_val):
Y,sy,w=sp.symbols('Y sy w')
exp=Y*(sp.exp(-sy*(w-440)))
ayel=float(exp.subs([(Y,Y_val),(sy,sy_val),(w,wl_val)]))
return ayel
####################################à
def get_bbw(wl_val):
b1_bw, bw_500, w =sp.symbols('b1_bw bw_500 w')
exp=b1_bw* bw_500* ((w/500)**(-1*(4/3)))
#print(exp)
bbw=float(exp.subs([(b1_bw,b1_bw_val),(bw_500,bw_500_val),(w,wl_val)]))
return bbw
def get_bbph(C_val, wl_val):
b1_bph, C, w= sp.symbols('b1_bph C w')
exp=0.3*b1_bph*(C**(0.62)) * (1+ 9*(sp.exp(-((w-685)**2/(2*100)))))
bbph=float(exp.subs([(b1_bph,b1_bph_val),(C,C_val),(w,wl_val)]))
return bbph
def get_bbsed(X_val, wl_val):
b1_bsed, X, w = sp.symbols('b1_bsed X w')
exp=b1_bsed* X * (((w)/(500))**(-1))
bbsed=float(exp.subs([(b1_bsed,b1_bsed_val),(X,X_val),(w,wl_val)]))
return bbsed
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment