def __init__(self, width, folding='Gauss'):
self.width = width
if folding == 'Gauss':
self.func = Gauss(width)
elif folding == 'Lorentz':
self.func = Lorentz(width)
elif folding == None:
self.func = None
else:
raise RuntimeError('unknown folding "' + folding + '"')
def __init__(self, width, folding='Gauss'):
self.width = width
if folding == 'Gauss':
self.func = Gauss(width)
elif folding == 'Lorentz':
self.func = Lorentz(width)
elif folding == 'ComplexGauss':
self.func = ComplexGauss(width)
elif folding == 'ComplexLorentz':
self.func = ComplexLorentz(width)
elif folding == 'RealLorentzPole':
self.func = LorentzPole(width, imag=False)
elif folding == 'ImaginaryLorentzPole':
self.func = LorentzPole(width, imag=True)
elif folding == 'Voigt':
self.func = Voigt(width)
elif folding is None:
self.func = None
else:
raise RuntimeError('unknown folding "' + folding + '"')
from math import exp, pi, sqrt
import numpy as np
from gpaw.gauss import Gauss, Lorentz
from gpaw.test import equal
from gpaw.utilities.folder import Folder
# Gauss and Lorentz functions
width = 0.5
x = 1.5
equal(Gauss(width).get(x),
exp(- x**2 / 2 / width**2) / sqrt(2 * pi) / width,
1.e-15)
equal(Lorentz(width).get(x),
width / (x**2 + width**2) / pi,
1.e-15)
# folder function
for name in ['Gauss', 'Lorentz']:
folder = Folder(width, name)
x = [0, 2]
y = [[2, 0, 1], [1, 1, 1]]
xl, yl = folder.fold(x, y, dx=.7)
# check first value
if name == 'Lorentz':
from math import exp, pi, sqrt
import numpy as np
from gpaw.gauss import Gauss, Lorentz
from gpaw.test import equal
from gpaw.utilities.folder import Folder, Voigt
# Gauss and Lorentz functions
width = 0.5
x = 1.5
equal(Gauss(width).get(x),
exp(- x**2 / 2 / width**2) / sqrt(2 * pi) / width,
1.e-15)
equal(Lorentz(width).get(x),
width / (x**2 + width**2) / pi,
1.e-15)
# folder function
for name in ['Gauss', 'Lorentz', 'Voigt']:
folder = Folder(width, name)
x = [0, 2]
y = [[2, 0, 1], [1, 1, 1]]
xl, yl = folder.fold(x, y, dx=.7)
# check first value
exec('func = {0}(width)'.format(name))
