Ejemplo n.º 1
0
    u = 1.0 / np.exp(y)


def rmse(observedValues, predictedValues):

    error = np.sqrt((predictedValues - observedValues) ** 2).mean()

    return error


if __name__ == "__main__":

    app = QApplication(sys.argv)

    # Spin up a main instance similar to what will exist after GUI input
    main = StartupWindow(tableData, headers, colors, old_coefficients)
    main.schlumberger()
    main.compute()
    rectangleCoordinates = [
        ((0.46988146917194829, 1267.1859604732258), 3.9293328705180155, -916.78591585951676),
        ((4.3992143396899639, 350.4000446137091), 19.549059554293287, 475.15927793601833),
        ((25.627732384051342, 807.14947142848234), 36.227680367632061, 579.6626526304766),
    ]
    for rectangle in rectangleCoordinates:
        main.canvas.rectxy = rectangle
        main.newRectangle()

    voltageSpacing = [0.55, 0.95, 1.5, 2.5, 3.0, 4.5, 5.5, 9.0, 12.0, 20.0, 30.0, 70.0]
    apparentResistivity = [125.0, 110.0, 95.0, 40.0, 24.0, 15.0, 10.5, 8.0, 6.0, 6.5, 11.0, 25.0]

    inversionAnalysis(voltageSpacing, apparentResistivity, main.canvas.rectCoordinates, "schlumberger")