nEvals = 200
phi = np.linspace(a, b, n)
y = np.zeros(n)
y2 = np.zeros(n)
mu_i = -0.6
mu_o = 0.5
alpha = 0.05
eta_i = 1.5c
eta_o = 1.0
eta = eta_i / eta_o
etaC = complex(eta, 0)
relerr = 10e-6
for i in range(n):
curPhi = phi[i]
y[i] = microfacetNoExp(mu_o, mu_i, etaC, alpha, curPhi)
quickTest = ll.microfacetFourierSeries(1.0, 1.0, etaC, alpha, n, 10e-4)
quickTest2 = microfacetFourierSeries(1.0, 1.0, etaC, alpha, n, 10e-4)
cOrig = ll.microfacetFourierSeries(mu_o, mu_i, etaC, alpha, n, 10e-4)
c = microfacetFourierSeries(mu_o, mu_i, etaC, alpha, n)
cnoExp = microfacetNoExpFourierSeries(mu_o, mu_i, etaC, alpha, nCoeffb, np.pi)
print("Our coeff: " + str(len(c)))
print("ll coeff: " + str(len(cOrig)))
plt.figure()
plt.plot(cOrig, 'b')
plt.plot(c, 'r--')
plt.savefig("Coeff.pdf")
def f(a):
return microfacetNoExp(mu_o, mu_i, etaC, alpha, a)
