Skip to content

Failed

gpkit.tests.from_paths.TestFiles.test_model_print_perf_py_mosek_conif (from gpkit.tests.from_paths.TestFiles-20200714225616)

Failing for the past 1 build (Since Failed #446 )
Took 0.92 sec.

Error Message

Invalid ConstraintSet element '[<gpkit.Climb object containing 2 top-level constraint(s) and 933 variable(s)>, <gpkit.Cruise object containing 2 top-level constraint(s) and 364 variable(s)>, <gpkit.Loiter object containing 2 top-level constraint(s) and 612 variable(s)>, <gpkit.Cruise object containing 2 top-level constraint(s) and 364 variable(s)>]' <class 'list'> was between Aircraft ========  Cost ----  1  Constraints -----------   Fuselage    FuelTank     Mission3.Aircraft.Fuselage.FuelTank.W >= f·W_{fuel-tot}     \mathcal{V}/Mission3.Aircraft.Fuselage.FuelTank.m_{fac} >= W_{fuel-tot}/\rho_{fuel}     FuselageSkin     m >= S_{wet}·\rho_{kevlar}·Mission3.Aircraft.Fuselage.FuselageSkin.t     Mission3.Aircraft.Fuselage.FuselageSkin.W >= m·Mission3.Aircraft.Fuselage.FuselageSkin.g     Mission3.Aircraft.Fuselage.FuselageSkin.t >= t_{min}     Mission3.Aircraft.Fuselage.FuselageSkin.I <= π·R³·Mission3.Aircraft.Fuselage.FuselageSkin.t     I_G >= m·(4·R² + 4·R·Mission3.Aircraft.Fuselage.FuselageSkin.t + Mission3.Aircraft.Fuselage.FuselageSkin.t²)     l_{body} = l_{body}     Mission3.Aircraft.Fuselage.FuselageSkin.E = Mission3.Aircraft.Fuselage.FuselageSkin.E    k_{body} = l_{body}/R    S_{wet} >= S_{body} + S_{nose} + S_{bulk}    S_{body} >= 2π·R·l_{body}    S_{nose}^1.6 >= (2π·R²)^1.6·(0.333 + 0.667·k_{nose}^1.6)    S_{bulk} >= R²·(0.0123·k_{bulk}² + 1.52·k_{bulk} + 0.502)    \mathcal{V}_{body} <= π·R²·l_{body}    Mission3.Aircraft.Fuselage.l <= 3·R·(k_{body}·k_{nose}·k_{bulk})^0.3333333333333333    Mission3.Aircraft.Fuselage.S >= π·R²    \mathcal{V}_{body} >= \mathcal{V}    Mission3.Aircraft.Fuselage.W/Mission3.Aircraft.Fuselage.m_{fac} >= Mission3.Aircraft.Fuselage.FuelTank.W + Mission3.Aircraft.Fuselage.FuselageSkin.W    Wing    Mission3.Aircraft.Wing.W/Mission3.Aircraft.Wing.mfac >= Mission3.Aircraft.Wing.WingSkin.W + Mission3.Aircraft.Wing.CapSpar.W + Mission3.Aircraft.Wing.WingCore.W     Planform     Mission3.Aircraft.Wing.Planform.b² = Mission3.Aircraft.Wing.Planform.S·Mission3.Aircraft.Wing.Planform.AR     cave[:] = cbave[:]·Mission3.Aircraft.Wing.Planform.S/Mission3.Aircraft.Wing.Planform.b     Mission3.Aircraft.Wing.Planform.croot = Mission3.Aircraft.Wing.Planform.S/Mission3.Aircraft.Wing.Planform.b·cbar[0]     Mission3.Aircraft.Wing.Planform.cmac = Mission3.Aircraft.Wing.Planform.croot·Mission3.Aircraft.Wing.Planform.cbarmac    WingSkin     Mission3.Aircraft.Wing.WingSkin.W >= CFRPFabric.rho·Mission3.Aircraft.Wing.Planform.S·2·Mission3.Aircraft.Wing.WingSkin.t·g     Mission3.Aircraft.Wing.WingSkin.t >= tmin     CFRPFabric.tau >= 1/Mission3.Aircraft.Wing.WingSkin.Jtbar/Mission3.Aircraft.Wing.Planform.croot²/Mission3.Aircraft.Wing.WingSkin.t·Mission3.Aircraft.Wing.WingSkin.Cmw·Mission3.Aircraft.Wing.Planform.S·Mission3.Aircraft.Wing.WingSkin.rhosl·Mission3.Aircraft.Wing.WingSkin.Vne²     CapSpar     I[:]/Mission3.Aircraft.Wing.CapSpar.mfac <= 2·w[:]·t[:]·(hin[:]/2)²     dm[:] >= (CFRPUD.rho·2·w[:]·t[:] + 2·tshear[:]·CFRPFabric.rho·(hin[:] + 2·t[:]) + FoamHD.rho·w[:]·hin[:])·Mission3.Aircraft.Wing.Planform.b/2·deta[:]     Mission3.Aircraft.Wing.CapSpar.W >= 2·dm[:].sum()·g     w[:] <= wlim·cave[:]     cave[:]·Mission3.Aircraft.Wing.Planform.tau >= hin[:] + 2·t[:]     Sy[:]·(hin[:]/2 + t[:]) <= I[:]     tshear[:] >= tmin     WingCore     Mission3.Aircraft.Wing.WingCore.W >= 2·(g·FoamHD.rho·Mission3.Aircraft.Wing.WingCore.Abar·cave[:]²·Mission3.Aircraft.Wing.Planform.b/2·deta[:]).sum()    mw·(1 + 2/Mission3.Aircraft.Wing.Planform.AR) >= 6.28    DF70    W_{DF70} <= Mission3.Aircraft.DF70.W/Mission3.Aircraft.DF70.m_{fac}    P_{sl-max} = P_{sl-max}    Empennage    HorizontalTail     Mission3.Aircraft.Empennage.HorizontalTail.W/Mission3.Aircraft.Empennage.HorizontalTail.mfac >= Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.W + Mission3.Aircraft.Empennage.HorizontalTail.WingCore.W      Planform      Mission3.Aircraft.Empennage.HorizontalTail.Planform.b² = Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·Mission3.Aircraft.Empennage.HorizontalTail.Planform.AR      cave[:] = cbave[:]·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S/Mission3.Aircraft.Empennage.HorizontalTail.Planform.b      Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot = Mission3.Aircraft.Empennage.HorizontalTail.Planform.S/Mission3.Aircraft.Empennage.HorizontalTail.Planform.b·cbar[0]      Mission3.Aircraft.Empennage.HorizontalTail.Planform.cmac = Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot·Mission3.Aircraft.Empennage.HorizontalTail.Planform.cbarmac     WingSkin      Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.W >= CFRPFabric.rho·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·2·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.t·g      Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.t >= tmin      CFRPFabric.tau >= 1/Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.Jtbar/Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot²/Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.t·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.Cmw·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.rhosl·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.Vne²      WingCore      Mission3.Aircraft.Empennage.HorizontalTail.WingCore.W >= 2·(g·FoamHD.rho·Mission3.Aircraft.Empennage.HorizontalTail.WingCore.Abar·cave[:]²·Mission3.Aircraft.Empennage.HorizontalTail.Planform.b/2·deta[:]).sum()     mh·(1 + 2/Mission3.Aircraft.Empennage.HorizontalTail.Planform.AR) <= 6.28     VerticalTail     Mission3.Aircraft.Empennage.VerticalTail.W/Mission3.Aircraft.Empennage.VerticalTail.mfac >= Mission3.Aircraft.Empennage.VerticalTail.WingSkin.W + Mission3.Aircraft.Empennage.VerticalTail.WingCore.W      Planform      Mission3.Aircraft.Empennage.VerticalTail.Planform.b² = Mission3.Aircraft.Empennage.VerticalTail.Planform.S·Mission3.Aircraft.Empennage.VerticalTail.Planform.AR      cave[:] = cbave[:]·Mission3.Aircraft.Empennage.VerticalTail.Planform.S/Mission3.Aircraft.Empennage.VerticalTail.Planform.b      Mission3.Aircraft.Empennage.VerticalTail.Planform.croot = Mission3.Aircraft.Empennage.VerticalTail.Planform.S/Mission3.Aircraft.Empennage.VerticalTail.Planform.b·cbar[0]      Mission3.Aircraft.Empennage.VerticalTail.Planform.cmac = Mission3.Aircraft.Empennage.VerticalTail.Planform.croot·Mission3.Aircraft.Empennage.VerticalTail.Planform.cbarmac     WingSkin      Mission3.Aircraft.Empennage.VerticalTail.WingSkin.W >= CFRPFabric.rho·Mission3.Aircraft.Empennage.VerticalTail.Planform.S·2·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.t·g      Mission3.Aircraft.Empennage.VerticalTail.WingSkin.t >= tmin      CFRPFabric.tau >= 1/Mission3.Aircraft.Empennage.VerticalTail.WingSkin.Jtbar/Mission3.Aircraft.Empennage.VerticalTail.Planform.croot²/Mission3.Aircraft.Empennage.VerticalTail.WingSkin.t·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.Cmw·Mission3.Aircraft.Empennage.VerticalTail.Planform.S·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.rhosl·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.Vne²      WingCore      Mission3.Aircraft.Empennage.VerticalTail.WingCore.W >= 2·(g·FoamHD.rho·Mission3.Aircraft.Empennage.VerticalTail.WingCore.Abar·cave[:]²·Mission3.Aircraft.Empennage.VerticalTail.Planform.b/2·deta[:]).sum()     TailBoom     I[:] <= π·t[:]·d[:]³/8     Sy[:] <= 2·I[:]/d[:]     dm[:] >= π·CFRPFabric.rho·d[:]·Mission3.Aircraft.Empennage.TailBoom.deta·t[:]·kfac·Mission3.Aircraft.Empennage.TailBoom.l     Mission3.Aircraft.Empennage.TailBoom.W/Mission3.Aircraft.Empennage.TailBoom.mfac >= g·dm[:].sum()     t[:] >= tmin     Mission3.Aircraft.Empennage.TailBoom.S = Mission3.Aircraft.Empennage.TailBoom.l·π·d[0]     Mission3.Aircraft.Empennage.TailBoom.b = 2·Mission3.Aircraft.Empennage.TailBoom.l    Mission3.Aircraft.Empennage.W/Mission3.Aircraft.Empennage.mfac >= Mission3.Aircraft.Empennage.HorizontalTail.W + Mission3.Aircraft.Empennage.VerticalTail.W + Mission3.Aircraft.Empennage.TailBoom.W    Mission3.Aircraft.Empennage.TailBoom.l >= lh    Mission3.Aircraft.Empennage.TailBoom.l >= lv    Pylon    Mission3.Aircraft.Pylon.S >= 2·Mission3.Aircraft.Pylon.l·Mission3.Aircraft.Pylon.h   W_{zfw} >= Mission3.Aircraft.Fuselage.W + Mission3.Aircraft.Wing.W + Mission3.Aircraft.DF70.W + Mission3.Aircraft.Empennage.W + Mission3.Aircraft.Pylon.W + W_{pay} + W_{avn}   Vh <= Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·lh/Mission3.Aircraft.Wing.Planform.S²·Mission3.Aircraft.Wing.Planform.b   Vv = Mission3.Aircraft.Empennage.VerticalTail.Planform.S·lv/Mission3.Aircraft.Wing.Planform.S/Mission3.Aircraft.Wing.Planform.b   Mission3.Aircraft.Wing.Planform.CLmax/mw <= Mission3.Aircraft.Empennage.HorizontalTail.Planform.CLmax/mh   w_{antenna} <= Mission3.Aircraft.Empennage.VerticalTail.Planform.croot·Mission3.Aircraft.Empennage.VerticalTail.Planform.lam   Mission3.Aircraft.Empennage.VerticalTail.Planform.b >= l_{antenna}   Mission3.Aircraft.Empennage.TailBoom.l >= lh + Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot   \mathcal{V}_{pay} <= (π/1.5)·k_{nose}·R³   \mathcal{V}_{body} >= \mathcal{V} + \mathcal{V}_{avn}   Mission3.Aircraft.DF70.h <= 2·R and AircraftLoading ===============  Cost ----  1  Constraints -----------   TailBoomBending    dx[:] = Mission3.Aircraft.Empennage.TailBoom.deta    qne[:]·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S <= Mission3.AircraftLoading.TailBoomBending.F    \bar{EI}[:] <= CFRPFabric.E·I[:]/Mission3.AircraftLoading.TailBoomBending.F/l²/2    Mr[:] >= \bar{M}[:-1]·Mission3.AircraftLoading.TailBoomBending.F·l    Mr[:]/Sy[:] <= CFRPFabric.sigma    Mission3.AircraftLoading.TailBoomBending.th = \theta[1]    Mission3.AircraftLoading.TailBoomBending.kappa >= \bar{\delta}[1]·Mission3.Aircraft.Empennage.HorizontalTail.Planform.CLmax·Mission3.AircraftLoading.TailBoomBending.Nsafety     Beam     \bar{M}[:-1] >= \bar{M}[1:] + 0.5·dx[:]·(\bar{S}[:-1] + \bar{S}[1:])     \bar{M}[1] >= Mission3.AircraftLoading.TailBoomBending.Beam.\bar{M}_{tip}     \theta[0] >= Mission3.AircraftLoading.TailBoomBending.Beam.\theta_{root}     \theta[1:] >= \theta[:-1] + 0.5·dx[:]·(\bar{M}[1:] + \bar{M}[:-1])/\bar{EI}[:]     \bar{\delta}[0] >= Mission3.AircraftLoading.TailBoomBending.Beam.\bar{\delta}_{root}     \bar{\delta}[1:] >= \bar{\delta}[:-1] + 0.5·dx[:]·(\theta[1:] + \theta[:-1])    TailBoomBending1    dx[:] = Mission3.Aircraft.Empennage.TailBoom.deta    qne[:]·Mission3.Aircraft.Empennage.VerticalTail.Planform.S <= Mission3.AircraftLoading.TailBoomBending1.F    \bar{EI}[:] <= CFRPFabric.E·I[:]/Mission3.AircraftLoading.TailBoomBending1.F/l²/2    Mr[:] >= \bar{M}[:-1]·Mission3.AircraftLoading.TailBoomBending1.F·l    Mr[:]/Sy[:] <= CFRPFabric.sigma    Mission3.AircraftLoading.TailBoomBending1.th = \theta[1]    Mission3.AircraftLoading.TailBoomBending1.kappa >= \bar{\delta}[1]·Mission3.Aircraft.Empennage.VerticalTail.Planform.CLmax·Mission3.AircraftLoading.TailBoomBending1.Nsafety     Beam     \bar{M}[:-1] >= \bar{M}[1:] + 0.5·dx[:]·(\bar{S}[:-1] + \bar{S}[1:])     \bar{M}[1] >= Mission3.AircraftLoading.TailBoomBending1.Beam.\bar{M}_{tip}     \theta[0] >= Mission3.AircraftLoading.TailBoomBending1.Beam.\theta_{root}     \theta[1:] >= \theta[:-1] + 0.5·dx[:]·(\bar{M}[1:] + \bar{M}[:-1])/\bar{EI}[:]     \bar{\delta}[0] >= Mission3.AircraftLoading.TailBoomBending1.Beam.\bar{\delta}_{root}     \bar{\delta}[1:] >= \bar{\delta}[:-1] + 0.5·dx[:]·(\theta[1:] + \theta[:-1])    SparLoading    S[:-1] >= S[1:] + 0.5·deta[:]·b/2·(q[:-1] + q[1:])    M[:-1] >= M[1:] + 0.5·deta[:]·b/2·(S[:-1] + S[1:])    N = Mission3.AircraftLoading.SparLoading.Nsafety·Nmax    q[:] >= N·W/b·cbar[:]    S[11] >= Stip    M[11] >= Mtip    th[0] >= throot    th[1:] >= th[:-1] + 0.5·deta[:]·b/2·(M[1:] + M[:-1])/CFRPUD.E/I[:]    w[0] >= wroot    w[1:] >= w[:-1] + 0.5·deta[:]·b/2·(th[1:] + th[:-1])    M[:-1]/Sy[:] <= CFRPUD.sigma    Mission3.AircraftLoading.SparLoading.kappa >= w[11]/(b/2)    FuselageLoading    FuselageSkinL     M_h >= Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.N_{max}·W_{cent}/4·l_{body}     Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.\sigma_{Kevlar} >= M_h·R/Mission3.Aircraft.Fuselage.FuselageSkin.I     Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.q >= W_{cent}·Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.N_{max}/l_{body}     \kappa·l_{body}/2 >= Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.q·(l_{body}/2)⁴/(8·Mission3.Aircraft.Fuselage.FuselageSkin.E·Mission3.Aircraft.Fuselage.FuselageSkin.I)     FuselageLanding     Mission3.AircraftLoading.FuselageLoading.FuselageLanding.F >= W_{cent}·Mission3.AircraftLoading.FuselageLoading.FuselageLanding.N_{max}     a >= Mission3.AircraftLoading.FuselageLoading.FuselageLanding.F/m     \dot{\omega} >= a/(l_{body}/2)     M_G >= I_G·\dot{\omega}     Mission3.AircraftLoading.FuselageLoading.FuselageLanding.\sigma_{Kevlar} >= M_G·R/Mission3.Aircraft.Fuselage.FuselageSkin.I    TailBoomFlexibility    Fne >= 1 + mh·Mission3.AircraftLoading.TailBoomBending.th    sph1·mw·Fne/mh/Vh + deda <= 1    sph2 <= Vh·CLhmin/Mission3.Aircraft.Wing.Planform.CLmax    deda >= mw·Mission3.Aircraft.Wing.Planform.S/b/4/π/lh    sph1 + sph2 >= SMcorr + CM/Mission3.Aircraft.Wing.Planform.CLmax.

Stacktrace

Traceback (most recent call last):
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/constraints/set.py", line 86, in __init__
    self._update(subconstraint)
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/constraints/set.py", line 114, in _update
    self.substitutions.update(constraint.substitutions)
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/keydict.py", line 157, in update
    self[k] = v
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/keydict.py", line 220, in __setitem__
    goodvals = ~np.isnan(value)
TypeError: ufunc 'isnan' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/tests/helpers.py", line 59, in test
    testfn(name, import_dict, path)(self)
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/tests/from_paths.py", line 48, in <lambda>
    lambda self: getattr(self, name)()))  # pylint:disable=undefined-variable
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/tests/from_paths.py", line 37, in test_fn
    mod.test()
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/jho/model/print_perf.py", line 194, in test
    M = Mission(DF70=True)
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/constraints/model.py", line 71, in __init__
    CostedConstraintSet.__init__(self, cost, constraints, substitutions)
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/constraints/costed.py", line 25, in __init__
    ConstraintSet.__init__(self, constraints, subs)
  File "/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/gpkit/constraints/set.py", line 88, in __init__
    raise badelement(self, i, constraint) from e
ValueError: Invalid ConstraintSet element '[<gpkit.Climb object containing 2 top-level constraint(s) and 933 variable(s)>, <gpkit.Cruise object containing 2 top-level constraint(s) and 364 variable(s)>, <gpkit.Loiter object containing 2 top-level constraint(s) and 612 variable(s)>, <gpkit.Cruise object containing 2 top-level constraint(s) and 364 variable(s)>]' <class 'list'> was between Aircraft
========

Cost
----
 1

Constraints
-----------
  Fuselage
   FuelTank
    Mission3.Aircraft.Fuselage.FuelTank.W >= f·W_{fuel-tot}
    \mathcal{V}/Mission3.Aircraft.Fuselage.FuelTank.m_{fac} >= W_{fuel-tot}/\rho_{fuel}

   FuselageSkin
    m >= S_{wet}·\rho_{kevlar}·Mission3.Aircraft.Fuselage.FuselageSkin.t
    Mission3.Aircraft.Fuselage.FuselageSkin.W >= m·Mission3.Aircraft.Fuselage.FuselageSkin.g
    Mission3.Aircraft.Fuselage.FuselageSkin.t >= t_{min}
    Mission3.Aircraft.Fuselage.FuselageSkin.I <= π·R³·Mission3.Aircraft.Fuselage.FuselageSkin.t
    I_G >= m·(4·R² + 4·R·Mission3.Aircraft.Fuselage.FuselageSkin.t + Mission3.Aircraft.Fuselage.FuselageSkin.t²)
    l_{body} = l_{body}
    Mission3.Aircraft.Fuselage.FuselageSkin.E = Mission3.Aircraft.Fuselage.FuselageSkin.E
   k_{body} = l_{body}/R
   S_{wet} >= S_{body} + S_{nose} + S_{bulk}
   S_{body} >= 2π·R·l_{body}
   S_{nose}^1.6 >= (2π·R²)^1.6·(0.333 + 0.667·k_{nose}^1.6)
   S_{bulk} >= R²·(0.0123·k_{bulk}² + 1.52·k_{bulk} + 0.502)
   \mathcal{V}_{body} <= π·R²·l_{body}
   Mission3.Aircraft.Fuselage.l <= 3·R·(k_{body}·k_{nose}·k_{bulk})^0.3333333333333333
   Mission3.Aircraft.Fuselage.S >= π·R²
   \mathcal{V}_{body} >= \mathcal{V}
   Mission3.Aircraft.Fuselage.W/Mission3.Aircraft.Fuselage.m_{fac} >= Mission3.Aircraft.Fuselage.FuelTank.W + Mission3.Aircraft.Fuselage.FuselageSkin.W

  Wing
   Mission3.Aircraft.Wing.W/Mission3.Aircraft.Wing.mfac >= Mission3.Aircraft.Wing.WingSkin.W + Mission3.Aircraft.Wing.CapSpar.W + Mission3.Aircraft.Wing.WingCore.W

   Planform
    Mission3.Aircraft.Wing.Planform.b² = Mission3.Aircraft.Wing.Planform.S·Mission3.Aircraft.Wing.Planform.AR
    cave[:] = cbave[:]·Mission3.Aircraft.Wing.Planform.S/Mission3.Aircraft.Wing.Planform.b
    Mission3.Aircraft.Wing.Planform.croot = Mission3.Aircraft.Wing.Planform.S/Mission3.Aircraft.Wing.Planform.b·cbar[0]
    Mission3.Aircraft.Wing.Planform.cmac = Mission3.Aircraft.Wing.Planform.croot·Mission3.Aircraft.Wing.Planform.cbarmac
   WingSkin
    Mission3.Aircraft.Wing.WingSkin.W >= CFRPFabric.rho·Mission3.Aircraft.Wing.Planform.S·2·Mission3.Aircraft.Wing.WingSkin.t·g
    Mission3.Aircraft.Wing.WingSkin.t >= tmin
    CFRPFabric.tau >= 1/Mission3.Aircraft.Wing.WingSkin.Jtbar/Mission3.Aircraft.Wing.Planform.croot²/Mission3.Aircraft.Wing.WingSkin.t·Mission3.Aircraft.Wing.WingSkin.Cmw·Mission3.Aircraft.Wing.Planform.S·Mission3.Aircraft.Wing.WingSkin.rhosl·Mission3.Aircraft.Wing.WingSkin.Vne²

   CapSpar
    I[:]/Mission3.Aircraft.Wing.CapSpar.mfac <= 2·w[:]·t[:]·(hin[:]/2)²
    dm[:] >= (CFRPUD.rho·2·w[:]·t[:] + 2·tshear[:]·CFRPFabric.rho·(hin[:] + 2·t[:]) + FoamHD.rho·w[:]·hin[:])·Mission3.Aircraft.Wing.Planform.b/2·deta[:]
    Mission3.Aircraft.Wing.CapSpar.W >= 2·dm[:].sum()·g
    w[:] <= wlim·cave[:]
    cave[:]·Mission3.Aircraft.Wing.Planform.tau >= hin[:] + 2·t[:]
    Sy[:]·(hin[:]/2 + t[:]) <= I[:]
    tshear[:] >= tmin

   WingCore
    Mission3.Aircraft.Wing.WingCore.W >= 2·(g·FoamHD.rho·Mission3.Aircraft.Wing.WingCore.Abar·cave[:]²·Mission3.Aircraft.Wing.Planform.b/2·deta[:]).sum()
   mw·(1 + 2/Mission3.Aircraft.Wing.Planform.AR) >= 6.28

  DF70
   W_{DF70} <= Mission3.Aircraft.DF70.W/Mission3.Aircraft.DF70.m_{fac}
   P_{sl-max} = P_{sl-max}

  Empennage
   HorizontalTail
    Mission3.Aircraft.Empennage.HorizontalTail.W/Mission3.Aircraft.Empennage.HorizontalTail.mfac >= Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.W + Mission3.Aircraft.Empennage.HorizontalTail.WingCore.W

    Planform
     Mission3.Aircraft.Empennage.HorizontalTail.Planform.b² = Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·Mission3.Aircraft.Empennage.HorizontalTail.Planform.AR
     cave[:] = cbave[:]·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S/Mission3.Aircraft.Empennage.HorizontalTail.Planform.b
     Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot = Mission3.Aircraft.Empennage.HorizontalTail.Planform.S/Mission3.Aircraft.Empennage.HorizontalTail.Planform.b·cbar[0]
     Mission3.Aircraft.Empennage.HorizontalTail.Planform.cmac = Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot·Mission3.Aircraft.Empennage.HorizontalTail.Planform.cbarmac
    WingSkin
     Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.W >= CFRPFabric.rho·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·2·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.t·g
     Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.t >= tmin
     CFRPFabric.tau >= 1/Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.Jtbar/Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot²/Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.t·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.Cmw·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.rhosl·Mission3.Aircraft.Empennage.HorizontalTail.WingSkin.Vne²

    WingCore
     Mission3.Aircraft.Empennage.HorizontalTail.WingCore.W >= 2·(g·FoamHD.rho·Mission3.Aircraft.Empennage.HorizontalTail.WingCore.Abar·cave[:]²·Mission3.Aircraft.Empennage.HorizontalTail.Planform.b/2·deta[:]).sum()
    mh·(1 + 2/Mission3.Aircraft.Empennage.HorizontalTail.Planform.AR) <= 6.28

   VerticalTail
    Mission3.Aircraft.Empennage.VerticalTail.W/Mission3.Aircraft.Empennage.VerticalTail.mfac >= Mission3.Aircraft.Empennage.VerticalTail.WingSkin.W + Mission3.Aircraft.Empennage.VerticalTail.WingCore.W

    Planform
     Mission3.Aircraft.Empennage.VerticalTail.Planform.b² = Mission3.Aircraft.Empennage.VerticalTail.Planform.S·Mission3.Aircraft.Empennage.VerticalTail.Planform.AR
     cave[:] = cbave[:]·Mission3.Aircraft.Empennage.VerticalTail.Planform.S/Mission3.Aircraft.Empennage.VerticalTail.Planform.b
     Mission3.Aircraft.Empennage.VerticalTail.Planform.croot = Mission3.Aircraft.Empennage.VerticalTail.Planform.S/Mission3.Aircraft.Empennage.VerticalTail.Planform.b·cbar[0]
     Mission3.Aircraft.Empennage.VerticalTail.Planform.cmac = Mission3.Aircraft.Empennage.VerticalTail.Planform.croot·Mission3.Aircraft.Empennage.VerticalTail.Planform.cbarmac
    WingSkin
     Mission3.Aircraft.Empennage.VerticalTail.WingSkin.W >= CFRPFabric.rho·Mission3.Aircraft.Empennage.VerticalTail.Planform.S·2·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.t·g
     Mission3.Aircraft.Empennage.VerticalTail.WingSkin.t >= tmin
     CFRPFabric.tau >= 1/Mission3.Aircraft.Empennage.VerticalTail.WingSkin.Jtbar/Mission3.Aircraft.Empennage.VerticalTail.Planform.croot²/Mission3.Aircraft.Empennage.VerticalTail.WingSkin.t·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.Cmw·Mission3.Aircraft.Empennage.VerticalTail.Planform.S·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.rhosl·Mission3.Aircraft.Empennage.VerticalTail.WingSkin.Vne²

    WingCore
     Mission3.Aircraft.Empennage.VerticalTail.WingCore.W >= 2·(g·FoamHD.rho·Mission3.Aircraft.Empennage.VerticalTail.WingCore.Abar·cave[:]²·Mission3.Aircraft.Empennage.VerticalTail.Planform.b/2·deta[:]).sum()

   TailBoom
    I[:] <= π·t[:]·d[:]³/8
    Sy[:] <= 2·I[:]/d[:]
    dm[:] >= π·CFRPFabric.rho·d[:]·Mission3.Aircraft.Empennage.TailBoom.deta·t[:]·kfac·Mission3.Aircraft.Empennage.TailBoom.l
    Mission3.Aircraft.Empennage.TailBoom.W/Mission3.Aircraft.Empennage.TailBoom.mfac >= g·dm[:].sum()
    t[:] >= tmin
    Mission3.Aircraft.Empennage.TailBoom.S = Mission3.Aircraft.Empennage.TailBoom.l·π·d[0]
    Mission3.Aircraft.Empennage.TailBoom.b = 2·Mission3.Aircraft.Empennage.TailBoom.l
   Mission3.Aircraft.Empennage.W/Mission3.Aircraft.Empennage.mfac >= Mission3.Aircraft.Empennage.HorizontalTail.W + Mission3.Aircraft.Empennage.VerticalTail.W + Mission3.Aircraft.Empennage.TailBoom.W
   Mission3.Aircraft.Empennage.TailBoom.l >= lh
   Mission3.Aircraft.Empennage.TailBoom.l >= lv

  Pylon
   Mission3.Aircraft.Pylon.S >= 2·Mission3.Aircraft.Pylon.l·Mission3.Aircraft.Pylon.h
  W_{zfw} >= Mission3.Aircraft.Fuselage.W + Mission3.Aircraft.Wing.W + Mission3.Aircraft.DF70.W + Mission3.Aircraft.Empennage.W + Mission3.Aircraft.Pylon.W + W_{pay} + W_{avn}
  Vh <= Mission3.Aircraft.Empennage.HorizontalTail.Planform.S·lh/Mission3.Aircraft.Wing.Planform.S²·Mission3.Aircraft.Wing.Planform.b
  Vv = Mission3.Aircraft.Empennage.VerticalTail.Planform.S·lv/Mission3.Aircraft.Wing.Planform.S/Mission3.Aircraft.Wing.Planform.b
  Mission3.Aircraft.Wing.Planform.CLmax/mw <= Mission3.Aircraft.Empennage.HorizontalTail.Planform.CLmax/mh
  w_{antenna} <= Mission3.Aircraft.Empennage.VerticalTail.Planform.croot·Mission3.Aircraft.Empennage.VerticalTail.Planform.lam
  Mission3.Aircraft.Empennage.VerticalTail.Planform.b >= l_{antenna}
  Mission3.Aircraft.Empennage.TailBoom.l >= lh + Mission3.Aircraft.Empennage.HorizontalTail.Planform.croot
  \mathcal{V}_{pay} <= (π/1.5)·k_{nose}·R³
  \mathcal{V}_{body} >= \mathcal{V} + \mathcal{V}_{avn}
  Mission3.Aircraft.DF70.h <= 2·R and AircraftLoading
===============

Cost
----
 1

Constraints
-----------
  TailBoomBending
   dx[:] = Mission3.Aircraft.Empennage.TailBoom.deta
   qne[:]·Mission3.Aircraft.Empennage.HorizontalTail.Planform.S <= Mission3.AircraftLoading.TailBoomBending.F
   \bar{EI}[:] <= CFRPFabric.E·I[:]/Mission3.AircraftLoading.TailBoomBending.F/l²/2
   Mr[:] >= \bar{M}[:-1]·Mission3.AircraftLoading.TailBoomBending.F·l
   Mr[:]/Sy[:] <= CFRPFabric.sigma
   Mission3.AircraftLoading.TailBoomBending.th = \theta[1]
   Mission3.AircraftLoading.TailBoomBending.kappa >= \bar{\delta}[1]·Mission3.Aircraft.Empennage.HorizontalTail.Planform.CLmax·Mission3.AircraftLoading.TailBoomBending.Nsafety

   Beam
    \bar{M}[:-1] >= \bar{M}[1:] + 0.5·dx[:]·(\bar{S}[:-1] + \bar{S}[1:])
    \bar{M}[1] >= Mission3.AircraftLoading.TailBoomBending.Beam.\bar{M}_{tip}
    \theta[0] >= Mission3.AircraftLoading.TailBoomBending.Beam.\theta_{root}
    \theta[1:] >= \theta[:-1] + 0.5·dx[:]·(\bar{M}[1:] + \bar{M}[:-1])/\bar{EI}[:]
    \bar{\delta}[0] >= Mission3.AircraftLoading.TailBoomBending.Beam.\bar{\delta}_{root}
    \bar{\delta}[1:] >= \bar{\delta}[:-1] + 0.5·dx[:]·(\theta[1:] + \theta[:-1])

  TailBoomBending1
   dx[:] = Mission3.Aircraft.Empennage.TailBoom.deta
   qne[:]·Mission3.Aircraft.Empennage.VerticalTail.Planform.S <= Mission3.AircraftLoading.TailBoomBending1.F
   \bar{EI}[:] <= CFRPFabric.E·I[:]/Mission3.AircraftLoading.TailBoomBending1.F/l²/2
   Mr[:] >= \bar{M}[:-1]·Mission3.AircraftLoading.TailBoomBending1.F·l
   Mr[:]/Sy[:] <= CFRPFabric.sigma
   Mission3.AircraftLoading.TailBoomBending1.th = \theta[1]
   Mission3.AircraftLoading.TailBoomBending1.kappa >= \bar{\delta}[1]·Mission3.Aircraft.Empennage.VerticalTail.Planform.CLmax·Mission3.AircraftLoading.TailBoomBending1.Nsafety

   Beam
    \bar{M}[:-1] >= \bar{M}[1:] + 0.5·dx[:]·(\bar{S}[:-1] + \bar{S}[1:])
    \bar{M}[1] >= Mission3.AircraftLoading.TailBoomBending1.Beam.\bar{M}_{tip}
    \theta[0] >= Mission3.AircraftLoading.TailBoomBending1.Beam.\theta_{root}
    \theta[1:] >= \theta[:-1] + 0.5·dx[:]·(\bar{M}[1:] + \bar{M}[:-1])/\bar{EI}[:]
    \bar{\delta}[0] >= Mission3.AircraftLoading.TailBoomBending1.Beam.\bar{\delta}_{root}
    \bar{\delta}[1:] >= \bar{\delta}[:-1] + 0.5·dx[:]·(\theta[1:] + \theta[:-1])

  SparLoading
   S[:-1] >= S[1:] + 0.5·deta[:]·b/2·(q[:-1] + q[1:])
   M[:-1] >= M[1:] + 0.5·deta[:]·b/2·(S[:-1] + S[1:])
   N = Mission3.AircraftLoading.SparLoading.Nsafety·Nmax
   q[:] >= N·W/b·cbar[:]
   S[11] >= Stip
   M[11] >= Mtip
   th[0] >= throot
   th[1:] >= th[:-1] + 0.5·deta[:]·b/2·(M[1:] + M[:-1])/CFRPUD.E/I[:]
   w[0] >= wroot
   w[1:] >= w[:-1] + 0.5·deta[:]·b/2·(th[1:] + th[:-1])
   M[:-1]/Sy[:] <= CFRPUD.sigma
   Mission3.AircraftLoading.SparLoading.kappa >= w[11]/(b/2)

  FuselageLoading
   FuselageSkinL
    M_h >= Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.N_{max}·W_{cent}/4·l_{body}
    Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.\sigma_{Kevlar} >= M_h·R/Mission3.Aircraft.Fuselage.FuselageSkin.I
    Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.q >= W_{cent}·Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.N_{max}/l_{body}
    \kappa·l_{body}/2 >= Mission3.AircraftLoading.FuselageLoading.FuselageSkinL.q·(l_{body}/2)⁴/(8·Mission3.Aircraft.Fuselage.FuselageSkin.E·Mission3.Aircraft.Fuselage.FuselageSkin.I)

   FuselageLanding
    Mission3.AircraftLoading.FuselageLoading.FuselageLanding.F >= W_{cent}·Mission3.AircraftLoading.FuselageLoading.FuselageLanding.N_{max}
    a >= Mission3.AircraftLoading.FuselageLoading.FuselageLanding.F/m
    \dot{\omega} >= a/(l_{body}/2)
    M_G >= I_G·\dot{\omega}
    Mission3.AircraftLoading.FuselageLoading.FuselageLanding.\sigma_{Kevlar} >= M_G·R/Mission3.Aircraft.Fuselage.FuselageSkin.I

  TailBoomFlexibility
   Fne >= 1 + mh·Mission3.AircraftLoading.TailBoomBending.th
   sph1·mw·Fne/mh/Vh + deda <= 1
   sph2 <= Vh·CLhmin/Mission3.Aircraft.Wing.Planform.CLmax
   deda >= mw·Mission3.Aircraft.Wing.Planform.S/b/4/π/lh
   sph1 + sph2 >= SMcorr + CM/Mission3.Aircraft.Wing.Planform.CLmax.
		

Standard Output

	

Standard Error

<string>:6: UnitStrippedWarning: The unit of the quantity is stripped when downcasting to ndarray.
/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/venv2_gpkit/lib/python3.7/site-packages/gpfit/fit_constraintset.py:46: VisibleDeprecationWarning: Creating an ndarray from ragged nested sequences (which is a list-or-tuple of lists-or-tuples-or ndarrays with different lengths or shapes) is deprecated. If you meant to do this, you must specify 'dtype=object' when creating the ndarray
  for k in range(fitdata["K"])]
/Users/jenkins/workspace/CE_gpkit_PR_research_models/buildnode/macys_VM/optimizer/mosek/venv2_gpkit/lib/python3.7/site-packages/numpy/core/_asarray.py:83: VisibleDeprecationWarning: Creating an ndarray from ragged nested sequences (which is a list-or-tuple of lists-or-tuples-or ndarrays with different lengths or shapes) is deprecated. If you meant to do this, you must specify 'dtype=object' when creating the ndarray
  return array(a, dtype, copy=False, order=order)