running wrong

[jasonaebischerGIT]

when running lequ1_2233 = 1/160**2 in the warsaw basis from from 160 GeV to 60 GeV the output of lequ1_2233 is wrong by several orders of magnitude

[jackypheno]

@jasonaebischerGIT

Looks like its an issue related to the size of WCs, for lequ1_2233 I get sensible output for

'wc_in = Wilson({'lequ1_2233':1e-6}, 160, 'SMEFT', 'Warsaw')
wc = wc_in.match_run(60, 'SMEFT', 'Warsaw')
print(wc.dict['lequ1_2233'])'

(1.0333888075747854e-06+2.115327309804727e-18j)

But the moment we change the initial condition to 1e-5 it gives

(-1.498326492711292e-11-1.39516319299177e-22j)

The problem persists even when using the 'leading_log' approximation.

[DavidMStraub]

I would check what happens before flavour rotations. Perhaps it's not a bug but just an effect of flavour rotations (i.e. that the output basis is not the same as the input basis) like in #63.

[jackypheno]

@DavidMStraub Good idea! Let me check that.

[peterstangl]

@jackypheno did you actually find out if the results was just due to a change of the flavour basis or if there is actually an issue?

[peterstangl]

I changed the label since it is not clear that this is actually a bug.