def METHOD_NAME(datafile):  # uses datafile pytest fixture
    mlcp = sn.MLCP()
    fname = datafile("diodeBridge_mlcp.dat")
    sn.mixedLinearComplementarity_newFromFilename(mlcp, fname)
    zsol = np.array(
        [
            9.85185185e-01,
            9.85185185e-01,
            -0.00000000e00,
            9.85185185e-04,
            0.00000000e00,
            0.00000000e00,
            9.85185185e-04,
        ]
    )
    return (mlcp, zsol)