Skip to content
Snippets Groups Projects
Select Git revision
  • fbc87f9ab83428a3c3b28c56a84c0cf6b5ae93b2
  • trunk
  • RELEASE_6_5_DRIVEDB
  • RELEASE_6_6_DRIVEDB
  • RELEASE_7_0_DRIVEDB
  • RELEASE_7_2_DRIVEDB
  • RELEASE_7_3_DRIVEDB
  • RELEASE_6_0_DRIVEDB
  • RELEASE_6_1_DRIVEDB
  • RELEASE_6_2_DRIVEDB
  • RELEASE_6_3_DRIVEDB
  • RELEASE_6_4_DRIVEDB
  • tags/RELEASE_7_4
  • tags/RELEASE_7_3
  • RELEASE_5_41_DRIVEDB
  • RELEASE_5_42_DRIVEDB
  • RELEASE_5_43_DRIVEDB
  • tags/RELEASE_7_2
  • tags/RELEASE_7_1
  • tags/RELEASE_7_0
  • RELEASE_5_40_DRIVEDB
21 results

verifymsg

Blame
  • EoMs.jl 690 B
    __precompile__()
    module EoMs
    
        # This module contains the Equations of motion
        # of PPN two-body problem
        # Note: I assume G = c = 1
    
        export EoM_massless
    
        function EoM_massless(u, p, t)
    
            x, y, vx, vy = u
            M, β, γ = p
    
            r = √(x^2 + y^2)
            φ = -M / r
            v² = vx^2 + vy^2
            rv = x * vx + y * vy
    
            A = - (1 + 2 * (β + γ) * φ + γ * v²) * M / r^3
            B = M * (2 * γ + 2) * rv / r^3
    
            acc_x = A * x + B * vx
            acc_y = A * y + B * vy
    
            return [vx, vy, acc_x, acc_y]
    
        end
    
        function EoM_general(u, p, t)
    
            x1, y1, vx1, vy1, x2, y2, vx2, vy2 = u
            M1, M2, β, γ = p
    
        end
    
    end