roessler_analytical_fp_cdp Subroutine

public subroutine roessler_analytical_fp_cdp(fp1, fp2)

Arguments

Type IntentOptional Attributes Name
class(state_vector_cdp), intent(out) :: fp1
class(state_vector_cdp), intent(out) :: fp2

Source Code

    subroutine roessler_analytical_fp_cdp(fp1, fp2)
        class(state_vector_cdp), intent(out) :: fp1, fp2

        complex(dp) :: d

        d = sqrt(c_dp**2 - 4*one_im_cdp*a_dp*b_dp)

        fp1%z = ( c_dp - d)/(2*one_im_cdp*a_dp)
        fp1%x = one_im_cdp*a_dp*fp1%z
        fp1%y = -one_im_cdp*fp1%z

        fp2%z = ( c_dp + d)/(2*one_im_cdp*a_dp)
        fp2%x = one_im_cdp*a_dp*fp2%z
        fp2%y = -one_im_cdp*fp2%z
    end subroutine roessler_analytical_fp_cdp