Beispiel #1
0
 def integrand(k):
     first = 6*(np.sin(k*r)-k*r*np.cos(k*r))*(((k*r)**2 -3)*np.sin(k*r)\
             +3*k*r*np.cos(k*r))/(k*r)**7 *linpwer.deltasq(k,1.)
     
     second = 6*(np.sin(r/k)-r/k*np.cos(r/k))*(((r/k)**2 -3)*np.sin(r/k)\
             +3*r/k*np.cos(r/k))/(r/k)**7 *linpwer.deltasq(1/k,1.)        
     
     return first + second
Beispiel #2
0
    def integrand(k):
        first = 9.*(np.sin(k*r)/(k*r)**3 -np.cos(k*r)/(k*r)**2)**2\
            *linpwer.deltasq(k,1)/k
        second =9*(np.sin(r/k)/(r/k)**3 -np.cos(r/k)/(r/k)**2)**2\
            *linpwer.deltasq(1/k,1)/k


        return math.sqrt(first + second)
Beispiel #3
0
#!/usr/bin/env python
import linpwer
import numpy as np
import math
import matplotlib.pyplot as plt

#This program calls the deltasq module to calculate the linear power spectrum

k=1000
x=np.linspace(.0001,k,10000)
delta=[]

for i in x[:]:
    delta.append(linpwer.deltasq(i,1)*10)

fig=plt.figure()
ax=fig.add_subplot(2,1,1)
ax.plot(x,delta)
ax.set_xscale('log')


#plt.plot(x,delta)
#plt.axes.set_xscale('log')
plt.show()