roessler_analytical_fp_rdp Subroutine

public subroutine roessler_analytical_fp_rdp(fp1, fp2)

Arguments

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

Source Code

    subroutine roessler_analytical_fp_rdp(fp1, fp2)
        class(state_vector_rdp), intent(out) :: fp1, fp2

        real(dp) :: d

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

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

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