Skip to content

cctk.quasiclassical rotation is broken #14

@corinwagen

Description

@corinwagen

Not sure why, but I get insane velocities when I set do_rotation to True in cctk.quasiclassical.get_quasiclassical_perturbation.

For instance, here's the velocities I get from running a 300 K quasiclassical initialization on test/static/methane_hpmodes.out:

[[  0.05273943  -0.06015768   0.08881368]
 [-10.454403     6.928813     3.9630766 ]
 [-22.108727    18.917265    -4.42791   ]
 [ 21.231005    -7.157976   -28.53397   ]
 [ 10.807802   -17.971182    27.94058   ]]

Without rotation, this is what I get:

[[ 0.          0.00822778  0.00959557]
 [ 0.04637501 -0.01431042 -0.02822099]
 [-0.01059805 -0.01377504 -0.03283611]
 [-0.01417696 -0.01534551 -0.00800313]
 [ 0.08252529 -0.05453904 -0.04518817]]

It's been a while since I've looked at this code — not sure what's wrong!

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions