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

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

        J['bending_stress_ratio','room_width'] = (room_length**2) * BEAM_HEIGHT_IN * (TOTAL_LOAD_PSI*room_width/2.0 + BEAM_WEIGHT_LBS_PER_IN) / (16.0 * I * YIELD_STRENGTH_PSI)
        J['bending_stress_ratio','room_length'] = (BEAM_WEIGHT_LBS_PER_IN + (TOTAL_LOAD_PSI*room_width/2.0)) * BEAM_HEIGHT_IN * room_length / (8.0 * I * YIELD_STRENGTH_PSI)

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