def solve_nonlinear(self, params, unknowns, resids):

        room_width = params['room_width']
        room_length = params['room_length']

        unknowns['shear_stress_ratio'] = 0.5 * ((0.5 * TOTAL_LOAD_PSI * room_width)  + BEAM_WEIGHT_LBS_PER_IN) * (CROSS_SECTIONAL_AREA_SQIN * YIELD_STRENGTH_PSI) / (room_length)
dataset/ETHPy150Open OpenMDAO/OpenMDAO/examples/beam_tutorial.py/ShearStress.solve_nonlinear/swapped operands of `/`