def METHOD_NAME(model: BDF, Kaa, Maa, TR, TRX):
    """static aero"""
    SKJ, D1JK, AJJi = pfaero(model)
    DJX = np.zeros((3, 3))
    KAA = np.zeros((3, 3))
    unused_KAAX = np.zeros((3, 3))
    KSALL = np.zeros((3, 3))
    KALX = np.zeros((3, 3))
    KRXA = np.zeros((3, 3))
    KLXA = np.zeros((3, 3))
    GKA = np.zeros((3, 3))
    GTKL = np.zeros((3, 3))
    QLL = np.zeros((3, 3))
    QKKS = np.zeros((3, 3))
    QKX = np.zeros((3, 3))
    WKK = np.zeros((3, 3))
    SRKT = np.zeros((3, 3))
    KALLA = np.zeros((3, 3))
    fi_q = np.zeros(3)
    u_l = np.zeros(3)
    u_x = np.zeros(3)
    PL = np.zeros(3)
    WGJ = np.zeros(3)
    FAJE = np.zeros(3)
    q_inf = 1.0
    KLL, KLR, unused_KRR, unused_KRL = partition_matrix(Kaa)
    MLL, MLR, MRR, MRL = partition_matrix(Maa)
    KLLi = np.linalg.inv(KLL)
    D = -KLLi @ KLR
    mr = D.T @ MLL @ D + MRL @ D + D.T @ MLR + MRR
    #ZZX = mr @ TR.T @ TRX
    MSLR = MLL @ D + MLR
    unused_MSRR = MRL @ D + D.T @ MSLR + MRR
    unused_M1RR = D.T @ MRL + MRR
    unused_M1RL = D.T @ MLL + MRL
    unused_RFSOP = TR.T @ TRX
    QAA = GKA.T @ QKKS @ GKA
    QAX = GKA.T @ QKX
    KAAA = -q_inf * QAA
    unused_KAAX = -q_inf * QAX
    unused_KSAA1 = KAA + KAAA
    # partition SUPORTs using KSAA1, KALX, KAAA
    #KLLA, KLRA, KRRA, KRLA = partition_matrix(KAAA)
    KALL = KLL - q_inf * QLL
    # TRX - boolean matrix that selects accelerations from the
    #       aerodynamic extra points
    # TR - transforms acceleratiosnf from aero ref point to supported
    #      DOFs form ALX
    KSALLi = np.linalg.inv(KSALL)
    unused_ALX = KSALLi @ KALX
    # WGJ - user input downwash vector
    # FAJE - user input pressure coefficient vector
    # PSA - external load vector
    # PZ - loads for trim calculation
    # IPZ - restrained elastic dimensional intercepts
    # IPZF2 - Unrestrained elastic dimensional intercepts
    # RINT - Rigid, unsplined dimensional intercepts
    # INTZ - Rigid, splined dimensional intercepts
    # HP0 - Perturbation in the support point deformations relative to mean axes due to external loads
    u_k = GTKL.T @ u_l
    # total downwash velocity
    w_j = D1JK @ u_k + DJX @ u_x + WGJ
    # pressure on aero elements
    unused_FFAJ = q_inf * AJJi @ w_j + q_inf * FAJE
    #  statibility
    QKX = WKK @ SKJ @ AJJi @ DJX
    unused_RSTAB = q_inf * SRKT.T @ QKX
    RINT = q_inf * SRKT.T @ (WKK @ SKJ @ AJJi @ w_j + SKJ @ fi_q)
    unused_KSAZX = D.T @ KLXA + KRXA
    unused_INTZ = GKA.T @ RINT
    mri = np.linalg.inv(mr)
    KALLi = np.linalg.inv(KALL)
    HP = mri @ (D.T @ MLL + MRL) @ KALLA @ ((MLL @ D + MLR) @ TR.T @ TRX + KALX)
    HP0 = -mri @ (D.T @ MLL + MRL) @ KALLi @ PL
    return HP, HP0