Skip to content

Use Jordan matrices and generalized eigenvectors to implement the solution when single-scattering albedo = 1. #7

Description

@dhjx1996

If the single-scattering albedo equals $0$, the coefficient matrix will not be diagonalizable. It is still possible to express the coefficient matrix $A = G J G^{-1}$, but $J$ will be a Jordan matrix and we will need to use generalized eigenvectors to solve the system of ODEs. Currently, users may use single-scattering albedos that are close to $1$ to approximate conservative scattering and PythonicDISORT will warn if a single-scattering albedo is too close to $1$.

We would like PythonicDISORT to instead use the special case solution if the single-scattering albedo is too close to $1$. This may render the joblib parallelization obsolete.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

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