Skip to content

reformulate vol step size #28

@bernstei

Description

@bernstei

The volume changing step size is specified in terms of an absolute delta volume because you can't use a fractional change like vol *= (1 + rv) where rv is a random variable in a symmetric interval around 0, because a 5% compression isn't reversed by a 5% expansion. However, this is already essentially solved in the stretch step, which does vec *= exp(rv) with rv a random variable in an interval symmetric about 0. This is reversible, because a * exp(rv) * exp(-rv) = a, unlike a * (1 + rv) * (1 - rv) != a.

It might be cleaner to replace the current code with one that uses a fractional volume change step size, with an exp(rv) factir to ensure reversibility.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions