Commit 6e23c2ab authored by Marco Selig's avatar Marco Selig

hotfix in infer_power.

parent aa67fcf6
...@@ -502,7 +502,7 @@ def infer_power(m,domain=None,Sk=None,D=None,pindex=None,pundex=None,kindex=None ...@@ -502,7 +502,7 @@ def infer_power(m,domain=None,Sk=None,D=None,pindex=None,pundex=None,kindex=None
if(smoothness): if(smoothness):
if(not domain.discrete): if(not domain.discrete):
numerator = weight_power(domain,numerator,power=-1,pindex=pindex,pundex=pundex) numerator = weight_power(domain,numerator,power=-1,pindex=pindex,pundex=pundex) ## bare(!)
pk = numerator/denominator1 ## bare(!) pk = numerator/denominator1 ## bare(!)
## smoothness prior ## smoothness prior
...@@ -535,8 +535,8 @@ def infer_power(m,domain=None,Sk=None,D=None,pindex=None,pundex=None,kindex=None ...@@ -535,8 +535,8 @@ def infer_power(m,domain=None,Sk=None,D=None,pindex=None,pundex=None,kindex=None
else: else:
pk = numerator/denominator1 ## non-bare(!) pk = numerator/denominator1 ## non-bare(!)
## weight if ... ## weight if ...
if(not domain.discrete)and(not bare): if(not domain.discrete)and(bare):
pk = weight_power(domain,pk,power=1,pindex=pindex,pundex=pundex) ## non-bare(!) pk = weight_power(domain,pk,power=-1,pindex=pindex,pundex=pundex) ## bare(!)
return pk return pk
......
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