def linearize(self, params, unknowns, resids):
        J = {}

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

        J['shear_stress_ratio','room_width'] = TOTAL_LOAD_PSI * room_length / (4.0 * YIELD_STRENGTH_PSI * CROSS_SECTIONAL_AREA_SQIN)
        J['shear_stress_ratio','room_length'] = (BEAM_WEIGHT_LBS_PER_IN + (TOTAL_LOAD_PSI * room_width / 2.0))/(2.0 * YIELD_STRENGTH_PSI * CROSS_SECTIONAL_AREA_SQIN)

        return J
dataset/ETHPy150Open OpenMDAO/OpenMDAO/examples/beam_tutorial.py/ShearStress.linearize/original