Commit 5e5e9454 authored by Philipp Arras's avatar Philipp Arras
Browse files

Rewrite cepskernel in numpy language

parent cdb27fee
...@@ -28,11 +28,8 @@ from ..operators.symmetrizing_operator import SymmetrizingOperator ...@@ -28,11 +28,8 @@ from ..operators.symmetrizing_operator import SymmetrizingOperator
from ..sugar import makeOp from ..sugar import makeOp
def _ceps_kernel(dof_space, k, a, k0): def _ceps_kernel(k, a, k0):
res = 1. return (a/(1+np.sum((k.T/k0)**2, axis=1)))**2
for i in range(len(k0)):
res += (k[i]/k0[i])**2
return a**2/(res)**2
def CepstrumOperator(target, a, k0): def CepstrumOperator(target, a, k0):
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment