Сгенерированный код тормозной схемы ТЭП70 (без второй кабины

// This file was generated by Trainz Air Scheme Compiler
// Created on 2012.12.10 19:21:25
// Scheme statistics:
// node count: 82
// link count: 109
//     unique: 104
//     to atm: 18
//  with name: 67
// matrix non-zeros: 215
//          fill-in: 47
// matrix density: 0.031975014872099940511600238
// object count: 23


// Place following lines in config.txt file to use this scheme:
// script   "tep70airscheme.gs"
// class    "tep70airscheme"

// INCLUDE ALL OBJECT SCRIPTS HERE
//include "someobject.gs"

include "SPAO_LimitValve.gs"
include "SPAO_EPK150.gs"
include "SPAO_BU367.gs"
include "SPAO_KM395.gs"
include "SPAO_KVT254.gs"
include "SPAO_BackValve.gs"
include "SPAO_Reductor.gs"
include "SPAO_VR292.gs"
include "SPAO_EVR305.gs"
include "SPAO_SwitchValve.gs"
include "SPAO_Relay.gs"
include "SPAO_VariableBrakeCylinder.gs"



include "splibrary.gs"
include "tep70airscheme.instance.gs"

final game class tep70airscheme isclass SPAirSchemeSpec
{
    public LIntMap nodeIds;
    public LIntMap linkIds;
    public LIntMap objectIds;
    public float[] linkNominalFlowRate;
    public int[] linkSrc;
    public int[] linkDst;
    public float[] alinkNominalFlowRate;
    public int[] alinkSrc;

    float _a00_00, _a00_11;
    float _a01_01, _a01_11;
    float _a02_02, _a02_28;
    float _a03_03, _a03_29;
    float _a04_04, _a04_13;
    float _a05_05, _a05_14;
    float _a06_06, _a06_39, _a06_81;
    float _a07_07, _a07_24;
    float _a08_08, _a08_23;
    float _a09_09, _a09_29;
    float _a10_10, _a10_20, _a10_80;
    float _a11_11;
    float _a12_12, _a12_28, _a12_35;
    float _a13_13, _a13_25;
    float _a14_14, _a14_26;
    float _a15_15, _a15_30, _a15_34;
    float _a16_16, _a16_31, _a16_36;
    float _a17_17, _a17_38, _a17_39;
    float _a18_18, _a18_27, _a18_33;
    float _a19_19, _a19_32, _a19_35;
    float _a20_20, _a20_27, _a20_80;
    float _a21_21, _a21_37, _a21_39;
    float _a22_22, _a22_23, _a22_39;
    float _a23_23, _a23_39;
    float _a24_24, _a24_32, _a24_38;
    float _a25_25, _a25_39;
    float _a26_26, _a26_34;
    float _a27_27, _a27_33, _a27_80;
    float _a28_28, _a28_29, _a28_35;
    float _a29_29, _a29_35;
    float _a30_30, _a30_34, _a30_36;
    float _a31_31, _a31_36, _a31_39;
    float _a32_32, _a32_35, _a32_38;
    float _a33_33, _a33_37, _a33_80;
    float _a34_34, _a34_36;
    float _a35_35, _a35_37, _a35_38;
    float _a36_36, _a36_39;
    float _a37_37, _a37_38, _a37_39, _a37_80;
    float _a38_38, _a38_39, _a38_80;
    float _a39_39, _a39_80, _a39_81;
    float _a40_40, _a40_74;
    float _a41_41, _a41_64;
    float _a42_42, _a42_64;
    float _a43_43, _a43_50;
    float _a44_44, _a44_62;
    float _a45_45, _a45_62;
    float _a46_46, _a46_60, _a46_81;
    float _a47_47, _a47_61;
    float _a48_48, _a48_74, _a48_77;
    float _a49_49, _a49_64, _a49_72;
    float _a50_50, _a50_63;
    float _a51_51, _a51_65, _a51_71;
    float _a52_52, _a52_73, _a52_78;
    float _a53_53, _a53_66, _a53_73;
    float _a54_54, _a54_67, _a54_79;
    float _a55_55, _a55_68, _a55_79;
    float _a56_56, _a56_69, _a56_78;
    float _a57_57, _a57_69, _a57_76;
    float _a58_58, _a58_70, _a58_75;
    float _a59_59, _a59_72, _a59_77;
    float _a60_60, _a60_76, _a60_81;
    float _a61_61, _a61_70, _a61_74;
    float _a62_62, _a62_78;
    float _a63_63, _a63_71;
    float _a64_64, _a64_72, _a64_75;
    float _a65_65, _a65_71, _a65_76;
    float _a66_66, _a66_73, _a66_79;
    float _a67_67, _a67_76, _a67_79;
    float _a68_68, _a68_77, _a68_79;
    float _a69_69, _a69_76, _a69_78;
    float _a70_70, _a70_74, _a70_75;
    float _a71_71, _a71_76;
    float _a72_72, _a72_75, _a72_77;
    float _a73_73, _a73_78, _a73_79, _a73_80;
    float _a74_74, _a74_75, _a74_77;
    float _a75_75, _a75_77;
    float _a76_76, _a76_78, _a76_79, _a76_81;
    float _a77_77, _a77_79;
    float _a78_78, _a78_79, _a78_80, _a78_81;
    float _a79_79, _a79_80, _a79_81;
    float _a80_80, _a80_81;
    float _a81_81;


    bool allowIdQuery = true;
    string currentObjectName;
    string currentObjectNameQualified;

    void SetCurrentObjectName(string name)
    {
        currentObjectName = name;
        currentObjectNameQualified = name + ".";
    }

    void ResetCurrentObjectName()
    {
        currentObjectName = null;
        currentObjectNameQualified = null;
    }

    public int GetNodeId(string name)
    {
        if (allowIdQuery) {
            if (!currentObjectName)
                return nodeIds.Get(name, SPAirScheme.InvalidNodeId);
            if (name and name.size() > 0)
                return nodeIds.Get(currentObjectNameQualified + name, SPAirScheme.InvalidNodeId);
            return nodeIds.Get(currentObjectName, SPAirScheme.InvalidNodeId);
        }
        return SPAirScheme.InvalidNodeId;
    }

    public int GetLinkId(string name)
    {
        if (allowIdQuery) {
            if (!currentObjectName)
                return nodeIds.Get(name, SPAirScheme.InvalidLinkId);
            if (name and name.size() > 0)
                return linkIds.Get(currentObjectNameQualified + name, SPAirScheme.InvalidLinkId);
            return linkIds.Get(currentObjectName, SPAirScheme.InvalidLinkId);
        }
        return SPAirScheme.InvalidLinkId;
    }

    public int GetObjectId(string name)
    {
        if (allowIdQuery) {
            if (!currentObjectName)
                return objectIds.Get(name, SPAirScheme.InvalidObjectId);
            if (name and name.size() > 0)
                return objectIds.Get(currentObjectNameQualified + name, SPAirScheme.InvalidObjectId);
            return objectIds.Get(currentObjectName, SPAirScheme.InvalidObjectId);
        }
        return SPAirScheme.InvalidObjectId;
    }

    tep70airscheme_instance CreateInstanceInternal(SPVehicle owner)
    {
        tep70airscheme_instance scheme = new tep70airscheme_instance();
        scheme.spec = me;
        scheme.owner = owner;

        float[] nodePressure = new float[82];
        float[] nodeVolume = new float[82];
        scheme.nodePressure = nodePressure;
        scheme.nodeVolume = nodeVolume;
        nodePressure[0] = 0;
        nodeVolume[0] = 0.3;
        nodePressure[1] = 0;
        nodeVolume[1] = 0.03;
        nodePressure[2] = 0;
        nodeVolume[2] = 0.1;
        nodePressure[3] = 0;
        nodeVolume[3] = 0.1;
        nodePressure[4] = 0;
        nodeVolume[4] = 1;
        nodePressure[5] = 0;
        nodeVolume[5] = 0.1;
        nodePressure[6] = 0;
        nodeVolume[6] = 0.05;
        nodePressure[7] = 0;
        nodeVolume[7] = 0.1;
        nodePressure[8] = 0;
        nodeVolume[8] = 0.1;
        nodePressure[9] = 0;
        nodeVolume[9] = 20;
        nodePressure[10] = 0;
        nodeVolume[10] = 20;
        nodePressure[11] = 0;
        nodeVolume[11] = 0.03;
        nodePressure[12] = 0;
        nodeVolume[12] = 0.1;
        nodePressure[13] = 0;
        nodeVolume[13] = 0.01;
        nodePressure[14] = 0;
        nodeVolume[14] = 0.1;
        nodePressure[15] = 0;
        nodeVolume[15] = 0.1;
        nodePressure[16] = 0;
        nodeVolume[16] = 0.1;
        nodePressure[17] = 0;
        nodeVolume[17] = 0.1;
        nodePressure[18] = 0;
        nodeVolume[18] = 0.1;
        nodePressure[19] = 0;
        nodeVolume[19] = 0.1;
        nodePressure[20] = 0;
        nodeVolume[20] = 0.1;
        nodePressure[21] = 0;
        nodeVolume[21] = 0.1;
        nodePressure[22] = 0;
        nodeVolume[22] = 500;
        nodePressure[23] = 0;
        nodeVolume[23] = 500;
        nodePressure[24] = 0;
        nodeVolume[24] = 20;
        nodePressure[25] = 0;
        nodeVolume[25] = 0.1;
        nodePressure[26] = 0;
        nodeVolume[26] = 0.1;
        nodePressure[27] = 0;
        nodeVolume[27] = 0.1;
        nodePressure[28] = 0;
        nodeVolume[28] = 0.1;
        nodePressure[29] = 0;
        nodeVolume[29] = 0.1;
        nodePressure[30] = 0;
        nodeVolume[30] = 0.05;
        nodePressure[31] = 0;
        nodeVolume[31] = 0.1;
        nodePressure[32] = 0;
        nodeVolume[32] = 0.1;
        nodePressure[33] = 0;
        nodeVolume[33] = 0.1;
        nodePressure[34] = 0;
        nodeVolume[34] = 0.5;
        nodePressure[35] = 0;
        nodeVolume[35] = 0.1;
        nodePressure[36] = 0;
        nodeVolume[36] = 78;
        nodePressure[37] = 0;
        nodeVolume[37] = 0.1;
        nodePressure[38] = 0;
        nodeVolume[38] = 0.1;
        nodePressure[39] = 0;
        nodeVolume[39] = 20;
        nodePressure[40] = 0;
        nodeVolume[40] = 1.5;
        nodePressure[41] = 0;
        nodeVolume[41] = 1;
        nodePressure[42] = 0;
        nodeVolume[42] = 0.1;
        nodePressure[43] = 0;
        nodeVolume[43] = 0.1;
        nodePressure[44] = 0;
        nodeVolume[44] = 0.05;
        nodePressure[45] = 0;
        nodeVolume[45] = 0.05;
        nodePressure[46] = 0;
        nodeVolume[46] = 0.1;
        nodePressure[47] = 0;
        nodeVolume[47] = 78;
        nodePressure[48] = 0;
        nodeVolume[48] = 0.1;
        nodePressure[49] = 0;
        nodeVolume[49] = 0.01;
        nodePressure[50] = 0;
        nodeVolume[50] = 0.1;
        nodePressure[51] = 0;
        nodeVolume[51] = 0.1;
        nodePressure[52] = 0;
        nodeVolume[52] = 0.1;
        nodePressure[53] = 0;
        nodeVolume[53] = 0.1;
        nodePressure[54] = 0;
        nodeVolume[54] = 0.1;
        nodePressure[55] = 0;
        nodeVolume[55] = 0.1;
        nodePressure[56] = 0;
        nodeVolume[56] = 0.1;
        nodePressure[57] = 0;
        nodeVolume[57] = 0.05;
        nodePressure[58] = 0;
        nodeVolume[58] = 0.4;
        nodePressure[59] = 0;
        nodeVolume[59] = 0.1;
        nodePressure[60] = 0;
        nodeVolume[60] = 0.1;
        nodePressure[61] = 0;
        nodeVolume[61] = 0.1;
        nodePressure[62] = 0;
        nodeVolume[62] = 5;
        nodePressure[63] = 0;
        nodeVolume[63] = 0.1;
        nodePressure[64] = 0;
        nodeVolume[64] = 0.1;
        nodePressure[65] = 0;
        nodeVolume[65] = 0.05;
        nodePressure[66] = 0;
        nodeVolume[66] = 0.1;
        nodePressure[67] = 0;
        nodeVolume[67] = 0.1;
        nodePressure[68] = 0;
        nodeVolume[68] = 20;
        nodePressure[69] = 0;
        nodeVolume[69] = 0.1;
        nodePressure[70] = 0;
        nodeVolume[70] = 0.1;
        nodePressure[71] = 0;
        nodeVolume[71] = 0.5;
        nodePressure[72] = 0;
        nodeVolume[72] = 0.4;
        nodePressure[73] = 0;
        nodeVolume[73] = 0.01;
        nodePressure[74] = 0;
        nodeVolume[74] = 0.1;
        nodePressure[75] = 0;
        nodeVolume[75] = 0.4;
        nodePressure[76] = 0;
        nodeVolume[76] = 78;
        nodePressure[77] = 0;
        nodeVolume[77] = 0.1;
        nodePressure[78] = 0;
        nodeVolume[78] = 0.01;
        nodePressure[79] = 0;
        nodeVolume[79] = 0.01;
        nodePressure[80] = 0;
        nodeVolume[80] = 0.1;
        nodePressure[81] = 0;
        nodeVolume[81] = 0.1;

        float[] linkFlowRate = new float[48];
        scheme.linkFlowRate = linkFlowRate;
        linkFlowRate[0] = 9.027353;
        linkFlowRate[1] = 0;
        linkFlowRate[2] = 0;
        linkFlowRate[3] = 0;
        linkFlowRate[4] = 0;
        linkFlowRate[5] = 0;
        linkFlowRate[6] = 2311.00244;
        linkFlowRate[7] = 0;
        linkFlowRate[8] = 0;
        linkFlowRate[9] = 0;
        linkFlowRate[10] = 0;
        linkFlowRate[11] = 0;
        linkFlowRate[12] = 0;
        linkFlowRate[13] = 0;
        linkFlowRate[14] = 0;
        linkFlowRate[15] = 0;
        linkFlowRate[16] = 0;
        linkFlowRate[17] = 0;
        linkFlowRate[18] = 0;
        linkFlowRate[19] = 0;
        linkFlowRate[20] = 1448.38867;
        linkFlowRate[21] = 0;
        linkFlowRate[22] = 1448.38867;
        linkFlowRate[23] = 0;
        linkFlowRate[24] = 0;
        linkFlowRate[25] = 0;
        linkFlowRate[26] = 0;
        linkFlowRate[27] = 0;
        linkFlowRate[28] = 0;
        linkFlowRate[29] = 0;
        linkFlowRate[30] = 0;
        linkFlowRate[31] = 0;
        linkFlowRate[32] = 0;
        linkFlowRate[33] = 0;
        linkFlowRate[34] = 0;
        linkFlowRate[35] = 159.242508;
        linkFlowRate[36] = 0;
        linkFlowRate[37] = 0;
        linkFlowRate[38] = 0;
        linkFlowRate[39] = 0;
        linkFlowRate[40] = 0;
        linkFlowRate[41] = 0;
        linkFlowRate[42] = 159.242508;
        linkFlowRate[43] = 159.242508;
        linkFlowRate[44] = 0;
        linkFlowRate[45] = 0;
        linkFlowRate[46] = 0;
        linkFlowRate[47] = 0;

        float[] alinkFlowRate = new float[19];
        scheme.alinkFlowRate = alinkFlowRate;
        alinkFlowRate[0] = 0;
        alinkFlowRate[1] = 0;
        alinkFlowRate[2] = 0;
        alinkFlowRate[3] = 0;
        alinkFlowRate[4] = 0;
        alinkFlowRate[5] = 0;
        alinkFlowRate[6] = 0;
        alinkFlowRate[7] = 0;
        alinkFlowRate[8] = 0;
        alinkFlowRate[9] = 0;
        alinkFlowRate[10] = 0;
        alinkFlowRate[11] = 0;
        alinkFlowRate[12] = 0;
        alinkFlowRate[13] = 6.780545;
        alinkFlowRate[14] = 0;
        alinkFlowRate[15] = 0;
        alinkFlowRate[16] = 401.2157;
        alinkFlowRate[17] = 0;
        alinkFlowRate[18] = 0;

        SPAirObject[] objects = new SPAirObject[23];
        scheme.objects = objects;
        objects[0] = new SPAO_LimitValve();
        objects[1] = new SPAO_EPK150();
        objects[2] = new SPAO_BU367();
        objects[3] = new SPAO_KM395();
        objects[4] = new SPAO_KVT254();
        objects[5] = new SPAO_BackValve();
        objects[6] = new SPAO_Reductor();
        objects[7] = new SPAO_BackValve();
        objects[8] = new SPAO_BackValve();
        objects[9] = new SPAO_VR292();
        objects[10] = new SPAO_EVR305();
        objects[11] = new SPAO_SwitchValve();
        objects[12] = new SPAO_SwitchValve();
        objects[13] = new SPAO_SwitchValve();
        objects[14] = new SPAO_Reductor();
        objects[15] = new SPAO_Relay();
        objects[16] = new SPAO_VariableBrakeCylinder();
        objects[17] = new SPAO_VariableBrakeCylinder();
        objects[18] = new SPAO_VariableBrakeCylinder();
        objects[19] = new SPAO_Relay();
        objects[20] = new SPAO_VariableBrakeCylinder();
        objects[21] = new SPAO_VariableBrakeCylinder();
        objects[22] = new SPAO_VariableBrakeCylinder();
        return scheme;
    }

    public SPAirScheme CreateInstance(SPVehicle owner)
    {
        SPAirScheme ret = CreateInstanceInternal(owner);
        return ret;
    }

    void InitInstanceInternal(tep70airscheme_instance scheme)
    {
        SPAirObject[] objects = scheme.objects;
        SetCurrentObjectName("kp3");
        objects[0].Init(scheme);
        SetCurrentObjectName("epka_1");
        objects[1].Init(scheme);
        SetCurrentObjectName("ubt_1");
        objects[2].Init(scheme);
        SetCurrentObjectName("k395_1");
        objects[3].Init(scheme);
        SetCurrentObjectName("k254_1");
        objects[4].Init(scheme);
        SetCurrentObjectName("ko3");
        objects[5].Init(scheme);
        SetCurrentObjectName("red1");
        objects[6].Init(scheme);
        SetCurrentObjectName("ko2");
        objects[7].Init(scheme);
        SetCurrentObjectName("ko4");
        objects[8].Init(scheme);
        SetCurrentObjectName("vr");
        objects[9].Init(scheme);
        SetCurrentObjectName("evr");
        objects[10].Init(scheme);
        SetCurrentObjectName("k3");
        objects[11].Init(scheme);
        SetCurrentObjectName("k4");
        objects[12].Init(scheme);
        SetCurrentObjectName("k5");
        objects[13].Init(scheme);
        SetCurrentObjectName("red3");
        objects[14].Init(scheme);
        SetCurrentObjectName("rd2");
        objects[15].Init(scheme);
        SetCurrentObjectName("tc_4_10");
        objects[16].Init(scheme);
        SetCurrentObjectName("tc_5_11");
        objects[17].Init(scheme);
        SetCurrentObjectName("tc_6_12");
        objects[18].Init(scheme);
        SetCurrentObjectName("rd1");
        objects[19].Init(scheme);
        SetCurrentObjectName("tc_3_9");
        objects[20].Init(scheme);
        SetCurrentObjectName("tc_2_8");
        objects[21].Init(scheme);
        SetCurrentObjectName("tc_1_7");
        objects[22].Init(scheme);
        ResetCurrentObjectName();

    }

    define float minimumPressureDelta = 1.0f;

    public void InitInstance(SPAirScheme scheme)
    {
        tep70airscheme_instance myScheme = cast<tep70airscheme_instance> scheme;
        if (!myScheme) {
            Exception("InitInstance: invalid scheme");
            return;
        }
        InitInstanceInternal(myScheme);

        Soup initialState = GetSchemeState("initial");
        if (initialState) {
            Soup initialProperties = initialState.GetNamedSoup("properties");
            scheme.SetProperties(initialProperties);
        }

    }

    public void UpdateInstance(tep70airscheme_instance scheme, float dt)
    {
        allowIdQuery = false;
        SPAirObject[] objects = scheme.objects;
        objects[0].Update(dt);
        objects[1].Update(dt);
        objects[2].Update(dt);
        objects[3].Update(dt);
        objects[4].Update(dt);
        objects[5].Update(dt);
        objects[6].Update(dt);
        objects[7].Update(dt);
        objects[8].Update(dt);
        objects[9].Update(dt);
        objects[10].Update(dt);
        objects[11].Update(dt);
        objects[12].Update(dt);
        objects[13].Update(dt);
        objects[14].Update(dt);
        objects[15].Update(dt);
        objects[16].Update(dt);
        objects[17].Update(dt);
        objects[18].Update(dt);
        objects[19].Update(dt);
        objects[20].Update(dt);
        objects[21].Update(dt);
        objects[22].Update(dt);
        allowIdQuery = true;

        float[] nodePressure = scheme.nodePressure;
        float[] nodeVolume = scheme.nodeVolume;
        float[] linkFlowRate = scheme.linkFlowRate;
        float[] alinkFlowRate = scheme.alinkFlowRate;
        float[] b = nodePressure;

        float pressure, volume, delta;
        _a00_00 = nodeVolume[0];
        _a01_01 = nodeVolume[1];
        pressure = nodePressure[2]; volume = nodeVolume[2]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a02_02 = volume + (Math.Sqrt(pressure) / pressure * dt) * (0.8124618);
        pressure = nodePressure[3]; volume = nodeVolume[3]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a03_03 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[0]);
        _a04_04 = nodeVolume[4];
        _a05_05 = nodeVolume[5];
        _a06_06 = nodeVolume[6];
        pressure = nodePressure[7]; volume = nodeVolume[7]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a07_07 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[1]);
        pressure = nodePressure[8]; volume = nodeVolume[8]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a08_08 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[2]);
        _a09_09 = nodeVolume[9];
        _a10_10 = nodeVolume[10];
        pressure = nodePressure[11]; volume = nodeVolume[11]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a11_11 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[3]);
        _a12_12 = nodeVolume[12];
        pressure = nodePressure[13]; volume = nodeVolume[13]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a13_13 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[4]);
        _a14_14 = nodeVolume[14];
        pressure = nodePressure[15]; volume = nodeVolume[15]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a15_15 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[5]);
        _a16_16 = nodeVolume[16];
        _a17_17 = nodeVolume[17];
        pressure = nodePressure[18]; volume = nodeVolume[18]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a18_18 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[6]);
        _a19_19 = nodeVolume[19];
        _a20_20 = nodeVolume[20];
        _a21_21 = nodeVolume[21];
        _a22_22 = nodeVolume[22];
        _a23_23 = nodeVolume[23];
        _a24_24 = nodeVolume[24];
        _a25_25 = nodeVolume[25];
        _a26_26 = nodeVolume[26];
        _a27_27 = nodeVolume[27];
        pressure = nodePressure[28]; volume = nodeVolume[28]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a28_28 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[7]);
        pressure = nodePressure[29]; volume = nodeVolume[29]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a29_29 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[8]);
        _a30_30 = nodeVolume[30];
        _a31_31 = nodeVolume[31];
        pressure = nodePressure[32]; volume = nodeVolume[32]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a32_32 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[9]);
        _a33_33 = nodeVolume[33];
        _a34_34 = nodeVolume[34];
        pressure = nodePressure[35]; volume = nodeVolume[35]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a35_35 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[10] + alinkFlowRate[11]);
        _a36_36 = nodeVolume[36];
        _a37_37 = nodeVolume[37];
        _a38_38 = nodeVolume[38];
        _a39_39 = nodeVolume[39];
        pressure = nodePressure[40]; volume = nodeVolume[40]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a40_40 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[12] + alinkFlowRate[13]);
        pressure = nodePressure[41]; volume = nodeVolume[41]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a41_41 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[14]);
        pressure = nodePressure[42]; volume = nodeVolume[42]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a42_42 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[15]);
        _a43_43 = nodeVolume[43];
        _a44_44 = nodeVolume[44];
        _a45_45 = nodeVolume[45];
        _a46_46 = nodeVolume[46];
        _a47_47 = nodeVolume[47];
        pressure = nodePressure[48]; volume = nodeVolume[48]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a48_48 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[16]);
        _a49_49 = nodeVolume[49];
        _a50_50 = nodeVolume[50];
        pressure = nodePressure[51]; volume = nodeVolume[51]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a51_51 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[17]);
        _a52_52 = nodeVolume[52];
        _a53_53 = nodeVolume[53];
        _a54_54 = nodeVolume[54];
        _a55_55 = nodeVolume[55];
        _a56_56 = nodeVolume[56];
        _a57_57 = nodeVolume[57];
        _a58_58 = nodeVolume[58];
        _a59_59 = nodeVolume[59];
        _a60_60 = nodeVolume[60];
        _a61_61 = nodeVolume[61];
        _a62_62 = nodeVolume[62];
        _a63_63 = nodeVolume[63];
        _a64_64 = nodeVolume[64];
        _a65_65 = nodeVolume[65];
        _a66_66 = nodeVolume[66];
        _a67_67 = nodeVolume[67];
        _a68_68 = nodeVolume[68];
        _a69_69 = nodeVolume[69];
        _a70_70 = nodeVolume[70];
        _a71_71 = nodeVolume[71];
        pressure = nodePressure[72]; volume = nodeVolume[72]; if(pressure < minimumPressureDelta) pressure = minimumPressureDelta;
        _a72_72 = volume + (Math.Sqrt(pressure) / pressure * dt) * (alinkFlowRate[18]);
        _a73_73 = nodeVolume[73];
        _a74_74 = nodeVolume[74];
        _a75_75 = nodeVolume[75];
        _a76_76 = nodeVolume[76];
        _a77_77 = nodeVolume[77];
        _a78_78 = nodeVolume[78];
        _a79_79 = nodeVolume[79];
        _a80_80 = nodeVolume[80];
        _a81_81 = nodeVolume[81];

        float flow;
        delta = nodePressure[0] - nodePressure[11];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2.56778049);
        _a11_11 = _a11_11 + flow; _a00_11 = -flow; _a00_00 = _a00_00 + flow;
        delta = nodePressure[11] - nodePressure[1];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[0]);
        _a01_01 = _a01_01 + flow; _a01_11 = -flow; _a11_11 = _a11_11 + flow;
        delta = nodePressure[2] - nodePressure[28];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[10]);
        _a28_28 = _a28_28 + flow; _a02_28 = -flow; _a02_02 = _a02_02 + flow;
        delta = nodePressure[3] - nodePressure[29];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[11]);
        _a29_29 = _a29_29 + flow; _a03_29 = -flow; _a03_03 = _a03_03 + flow;
        delta = nodePressure[4] - nodePressure[13];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (16.0486279);
        _a13_13 = _a13_13 + flow; _a04_13 = -flow; _a04_04 = _a04_04 + flow;
        delta = nodePressure[5] - nodePressure[14];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a14_14 = _a14_14 + flow; _a05_14 = -flow; _a05_05 = _a05_05 + flow;
        delta = nodePressure[6] - nodePressure[39];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[20]);
        _a39_39 = _a39_39 + flow; _a06_39 = -flow; _a06_06 = _a06_06 + flow;
        delta = nodePressure[81] - nodePressure[6];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[1]);
        _a06_06 = _a06_06 + flow; _a06_81 = -flow; _a81_81 = _a81_81 + flow;
        delta = nodePressure[7] - nodePressure[24];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[7]);
        _a24_24 = _a24_24 + flow; _a07_24 = -flow; _a07_07 = _a07_07 + flow;
        delta = nodePressure[8] - nodePressure[23];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a23_23 = _a23_23 + flow; _a08_23 = -flow; _a08_08 = _a08_08 + flow;
        delta = nodePressure[9] - nodePressure[29];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (678.054565);
        _a29_29 = _a29_29 + flow; _a09_29 = -flow; _a09_09 = _a09_09 + flow;
        delta = nodePressure[20] - nodePressure[10];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a10_10 = _a10_10 + flow; _a10_20 = -flow; _a20_20 = _a20_20 + flow;
        delta = nodePressure[80] - nodePressure[10];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (159.242508);
        _a10_10 = _a10_10 + flow; _a10_80 = -flow; _a80_80 = _a80_80 + flow;
        delta = nodePressure[28] - nodePressure[12];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[2] + linkFlowRate[3]);
        _a12_12 = _a12_12 + flow; _a12_28 = -flow; _a28_28 = _a28_28 + flow;
        delta = nodePressure[12] - nodePressure[35];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[17] + linkFlowRate[18]);
        _a35_35 = _a35_35 + flow; _a12_35 = -flow; _a12_12 = _a12_12 + flow;
        delta = nodePressure[13] - nodePressure[25];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[9]);
        _a25_25 = _a25_25 + flow; _a13_25 = -flow; _a13_13 = _a13_13 + flow;
        delta = nodePressure[14] - nodePressure[26];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a26_26 = _a26_26 + flow; _a14_26 = -flow; _a14_14 = _a14_14 + flow;
        delta = nodePressure[15] - nodePressure[30];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[13]);
        _a30_30 = _a30_30 + flow; _a15_30 = -flow; _a15_15 = _a15_15 + flow;
        delta = nodePressure[34] - nodePressure[15];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a15_15 = _a15_15 + flow; _a15_34 = -flow; _a34_34 = _a34_34 + flow;
        delta = nodePressure[16] - nodePressure[31];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[14]);
        _a31_31 = _a31_31 + flow; _a16_31 = -flow; _a16_16 = _a16_16 + flow;
        delta = nodePressure[36] - nodePressure[16];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a16_16 = _a16_16 + flow; _a16_36 = -flow; _a36_36 = _a36_36 + flow;
        delta = nodePressure[17] - nodePressure[38];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[19]);
        _a38_38 = _a38_38 + flow; _a17_38 = -flow; _a17_17 = _a17_17 + flow;
        delta = nodePressure[39] - nodePressure[17];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a17_17 = _a17_17 + flow; _a17_39 = -flow; _a39_39 = _a39_39 + flow;
        delta = nodePressure[18] - nodePressure[27];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a27_27 = _a27_27 + flow; _a18_27 = -flow; _a18_18 = _a18_18 + flow;
        delta = nodePressure[18] - nodePressure[33];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[16]);
        _a33_33 = _a33_33 + flow; _a18_33 = -flow; _a18_18 = _a18_18 + flow;
        delta = nodePressure[19] - nodePressure[32];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[15]);
        _a32_32 = _a32_32 + flow; _a19_32 = -flow; _a19_19 = _a19_19 + flow;
        delta = nodePressure[35] - nodePressure[19];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a19_19 = _a19_19 + flow; _a19_35 = -flow; _a35_35 = _a35_35 + flow;
        delta = nodePressure[27] - nodePressure[20];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[4]);
        _a20_20 = _a20_20 + flow; _a20_27 = -flow; _a27_27 = _a27_27 + flow;
        delta = nodePressure[37] - nodePressure[21];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[5]);
        _a21_21 = _a21_21 + flow; _a21_37 = -flow; _a37_37 = _a37_37 + flow;
        delta = nodePressure[21] - nodePressure[39];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a39_39 = _a39_39 + flow; _a21_39 = -flow; _a21_21 = _a21_21 + flow;
        delta = nodePressure[22] - nodePressure[23];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[6]);
        _a23_23 = _a23_23 + flow; _a22_23 = -flow; _a22_22 = _a22_22 + flow;
        delta = nodePressure[39] - nodePressure[22];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a22_22 = _a22_22 + flow; _a22_39 = -flow; _a39_39 = _a39_39 + flow;
        delta = nodePressure[32] - nodePressure[24];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a24_24 = _a24_24 + flow; _a24_32 = -flow; _a32_32 = _a32_32 + flow;
        delta = nodePressure[38] - nodePressure[24];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[8]);
        _a24_24 = _a24_24 + flow; _a24_38 = -flow; _a38_38 = _a38_38 + flow;
        delta = nodePressure[25] - nodePressure[39];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[21]);
        _a39_39 = _a39_39 + flow; _a25_39 = -flow; _a25_25 = _a25_25 + flow;
        delta = nodePressure[26] - nodePressure[34];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (100.303925);
        _a34_34 = _a34_34 + flow; _a26_34 = -flow; _a26_26 = _a26_26 + flow;
        delta = nodePressure[28] - nodePressure[29];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (10.271122);
        _a29_29 = _a29_29 + flow; _a28_29 = -flow; _a28_28 = _a28_28 + flow;
        delta = nodePressure[35] - nodePressure[29];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[12]);
        _a29_29 = _a29_29 + flow; _a29_35 = -flow; _a35_35 = _a35_35 + flow;
        delta = nodePressure[36] - nodePressure[30];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a30_30 = _a30_30 + flow; _a30_36 = -flow; _a36_36 = _a36_36 + flow;
        delta = nodePressure[31] - nodePressure[39];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[22]);
        _a39_39 = _a39_39 + flow; _a31_39 = -flow; _a31_31 = _a31_31 + flow;
        delta = nodePressure[33] - nodePressure[37];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a37_37 = _a37_37 + flow; _a33_37 = -flow; _a33_33 = _a33_33 + flow;
        delta = nodePressure[35] - nodePressure[37];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (2311.00244);
        _a37_37 = _a37_37 + flow; _a35_37 = -flow; _a35_35 = _a35_35 + flow;
        delta = nodePressure[40] - nodePressure[74];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[38]);
        _a74_74 = _a74_74 + flow; _a40_74 = -flow; _a40_40 = _a40_40 + flow;
        delta = nodePressure[41] - nodePressure[64];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[31]);
        _a64_64 = _a64_64 + flow; _a41_64 = -flow; _a41_41 = _a41_41 + flow;
        delta = nodePressure[64] - nodePressure[42];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (196.595688);
        _a42_42 = _a42_42 + flow; _a42_64 = -flow; _a64_64 = _a64_64 + flow;
        delta = nodePressure[43] - nodePressure[50];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a50_50 = _a50_50 + flow; _a43_50 = -flow; _a43_43 = _a43_43 + flow;
        delta = nodePressure[62] - nodePressure[44];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (159.242508);
        _a44_44 = _a44_44 + flow; _a44_62 = -flow; _a62_62 = _a62_62 + flow;
        delta = nodePressure[62] - nodePressure[45];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (159.242508);
        _a45_45 = _a45_45 + flow; _a45_62 = -flow; _a62_62 = _a62_62 + flow;
        delta = nodePressure[60] - nodePressure[46];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[23]);
        _a46_46 = _a46_46 + flow; _a46_60 = -flow; _a60_60 = _a60_60 + flow;
        delta = nodePressure[46] - nodePressure[81];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a81_81 = _a81_81 + flow; _a46_81 = -flow; _a46_46 = _a46_46 + flow;
        delta = nodePressure[61] - nodePressure[47];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a47_47 = _a47_47 + flow; _a47_61 = -flow; _a61_61 = _a61_61 + flow;
        delta = nodePressure[48] - nodePressure[74];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[39]);
        _a74_74 = _a74_74 + flow; _a48_74 = -flow; _a48_48 = _a48_48 + flow;
        delta = nodePressure[48] - nodePressure[77];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[44]);
        _a77_77 = _a77_77 + flow; _a48_77 = -flow; _a48_48 = _a48_48 + flow;
        delta = nodePressure[49] - nodePressure[64];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[32]);
        _a64_64 = _a64_64 + flow; _a49_64 = -flow; _a49_49 = _a49_49 + flow;
        delta = nodePressure[72] - nodePressure[49];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[24]);
        _a49_49 = _a49_49 + flow; _a49_72 = -flow; _a72_72 = _a72_72 + flow;
        delta = nodePressure[50] - nodePressure[63];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a63_63 = _a63_63 + flow; _a50_63 = -flow; _a50_50 = _a50_50 + flow;
        delta = nodePressure[51] - nodePressure[65];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[34]);
        _a65_65 = _a65_65 + flow; _a51_65 = -flow; _a51_51 = _a51_51 + flow;
        delta = nodePressure[71] - nodePressure[51];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a51_51 = _a51_51 + flow; _a51_71 = -flow; _a71_71 = _a71_71 + flow;
        delta = nodePressure[52] - nodePressure[73];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (159.242508);
        _a73_73 = _a73_73 + flow; _a52_73 = -flow; _a52_52 = _a52_52 + flow;
        delta = nodePressure[78] - nodePressure[52];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[25]);
        _a52_52 = _a52_52 + flow; _a52_78 = -flow; _a78_78 = _a78_78 + flow;
        delta = nodePressure[53] - nodePressure[66];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[35]);
        _a66_66 = _a66_66 + flow; _a53_66 = -flow; _a53_53 = _a53_53 + flow;
        delta = nodePressure[73] - nodePressure[53];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[26]);
        _a53_53 = _a53_53 + flow; _a53_73 = -flow; _a73_73 = _a73_73 + flow;
        delta = nodePressure[54] - nodePressure[67];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[36]);
        _a67_67 = _a67_67 + flow; _a54_67 = -flow; _a54_54 = _a54_54 + flow;
        delta = nodePressure[79] - nodePressure[54];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[27]);
        _a54_54 = _a54_54 + flow; _a54_79 = -flow; _a79_79 = _a79_79 + flow;
        delta = nodePressure[55] - nodePressure[68];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (159.242508);
        _a68_68 = _a68_68 + flow; _a55_68 = -flow; _a55_55 = _a55_55 + flow;
        delta = nodePressure[79] - nodePressure[55];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[28]);
        _a55_55 = _a55_55 + flow; _a55_79 = -flow; _a79_79 = _a79_79 + flow;
        delta = nodePressure[56] - nodePressure[69];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[37]);
        _a69_69 = _a69_69 + flow; _a56_69 = -flow; _a56_56 = _a56_56 + flow;
        delta = nodePressure[78] - nodePressure[56];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[29]);
        _a56_56 = _a56_56 + flow; _a56_78 = -flow; _a78_78 = _a78_78 + flow;
        delta = nodePressure[69] - nodePressure[57];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[30]);
        _a57_57 = _a57_57 + flow; _a57_69 = -flow; _a69_69 = _a69_69 + flow;
        delta = nodePressure[57] - nodePressure[76];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[42]);
        _a76_76 = _a76_76 + flow; _a57_76 = -flow; _a57_57 = _a57_57 + flow;
        delta = nodePressure[58] - nodePressure[70];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a70_70 = _a70_70 + flow; _a58_70 = -flow; _a58_58 = _a58_58 + flow;
        delta = nodePressure[58] - nodePressure[75];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (324.9847);
        _a75_75 = _a75_75 + flow; _a58_75 = -flow; _a58_58 = _a58_58 + flow;
        delta = nodePressure[59] - nodePressure[72];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a72_72 = _a72_72 + flow; _a59_72 = -flow; _a59_59 = _a59_59 + flow;
        delta = nodePressure[59] - nodePressure[77];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[45]);
        _a77_77 = _a77_77 + flow; _a59_77 = -flow; _a59_59 = _a59_59 + flow;
        delta = nodePressure[76] - nodePressure[60];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a60_60 = _a60_60 + flow; _a60_76 = -flow; _a76_76 = _a76_76 + flow;
        delta = nodePressure[70] - nodePressure[61];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a61_61 = _a61_61 + flow; _a61_70 = -flow; _a70_70 = _a70_70 + flow;
        delta = nodePressure[74] - nodePressure[61];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a61_61 = _a61_61 + flow; _a61_74 = -flow; _a74_74 = _a74_74 + flow;
        delta = nodePressure[62] - nodePressure[78];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (159.242508);
        _a78_78 = _a78_78 + flow; _a62_78 = -flow; _a62_62 = _a62_62 + flow;
        delta = nodePressure[63] - nodePressure[71];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (100.303925);
        _a71_71 = _a71_71 + flow; _a63_71 = -flow; _a63_63 = _a63_63 + flow;
        delta = nodePressure[75] - nodePressure[64];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[33]);
        _a64_64 = _a64_64 + flow; _a64_75 = -flow; _a75_75 = _a75_75 + flow;
        delta = nodePressure[76] - nodePressure[65];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a65_65 = _a65_65 + flow; _a65_76 = -flow; _a76_76 = _a76_76 + flow;
        delta = nodePressure[66] - nodePressure[79];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[46]);
        _a79_79 = _a79_79 + flow; _a66_79 = -flow; _a66_66 = _a66_66 + flow;
        delta = nodePressure[67] - nodePressure[76];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[43]);
        _a76_76 = _a76_76 + flow; _a67_76 = -flow; _a67_67 = _a67_67 + flow;
        delta = nodePressure[68] - nodePressure[77];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (1448.38867);
        _a77_77 = _a77_77 + flow; _a68_77 = -flow; _a68_68 = _a68_68 + flow;
        delta = nodePressure[72] - nodePressure[75];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[40] + linkFlowRate[41]);
        _a75_75 = _a75_75 + flow; _a72_75 = -flow; _a72_72 = _a72_72 + flow;
        delta = nodePressure[73] - nodePressure[80];
        if (delta < 0.0f) delta = -delta; if (delta < minimumPressureDelta) delta = minimumPressureDelta;
        flow = Math.Sqrt(delta) / delta * dt * (linkFlowRate[47]);
        _a80_80 = _a80_80 + flow; _a73_80 = -flow; _a73_73 = _a73_73 + flow;

        b[0] = nodePressure[0] * nodeVolume[0];
        b[1] = nodePressure[1] * nodeVolume[1];
        b[2] = nodePressure[2] * nodeVolume[2];
        b[3] = nodePressure[3] * nodeVolume[3];
        b[4] = nodePressure[4] * nodeVolume[4];
        b[5] = nodePressure[5] * nodeVolume[5];
        b[6] = nodePressure[6] * nodeVolume[6];
        b[7] = nodePressure[7] * nodeVolume[7];
        b[8] = nodePressure[8] * nodeVolume[8];
        b[9] = nodePressure[9] * nodeVolume[9];
        b[10] = nodePressure[10] * nodeVolume[10];
        b[11] = nodePressure[11] * nodeVolume[11];
        b[12] = nodePressure[12] * nodeVolume[12];
        b[13] = nodePressure[13] * nodeVolume[13];
        b[14] = nodePressure[14] * nodeVolume[14];
        b[15] = nodePressure[15] * nodeVolume[15];
        b[16] = nodePressure[16] * nodeVolume[16];
        b[17] = nodePressure[17] * nodeVolume[17];
        b[18] = nodePressure[18] * nodeVolume[18];
        b[19] = nodePressure[19] * nodeVolume[19];
        b[20] = nodePressure[20] * nodeVolume[20];
        b[21] = nodePressure[21] * nodeVolume[21];
        b[22] = nodePressure[22] * nodeVolume[22];
        b[23] = nodePressure[23] * nodeVolume[23];
        b[24] = nodePressure[24] * nodeVolume[24];
        b[25] = nodePressure[25] * nodeVolume[25];
        b[26] = nodePressure[26] * nodeVolume[26];
        b[27] = nodePressure[27] * nodeVolume[27];
        b[28] = nodePressure[28] * nodeVolume[28];
        b[29] = nodePressure[29] * nodeVolume[29];
        b[30] = nodePressure[30] * nodeVolume[30];
        b[31] = nodePressure[31] * nodeVolume[31];
        b[32] = nodePressure[32] * nodeVolume[32];
        b[33] = nodePressure[33] * nodeVolume[33];
        b[34] = nodePressure[34] * nodeVolume[34];
        b[35] = nodePressure[35] * nodeVolume[35];
        b[36] = nodePressure[36] * nodeVolume[36];
        b[37] = nodePressure[37] * nodeVolume[37];
        b[38] = nodePressure[38] * nodeVolume[38];
        b[39] = nodePressure[39] * nodeVolume[39];
        b[40] = nodePressure[40] * nodeVolume[40];
        b[41] = nodePressure[41] * nodeVolume[41];
        b[42] = nodePressure[42] * nodeVolume[42];
        b[43] = nodePressure[43] * nodeVolume[43];
        b[44] = nodePressure[44] * nodeVolume[44];
        b[45] = nodePressure[45] * nodeVolume[45];
        b[46] = nodePressure[46] * nodeVolume[46];
        b[47] = nodePressure[47] * nodeVolume[47];
        b[48] = nodePressure[48] * nodeVolume[48];
        b[49] = nodePressure[49] * nodeVolume[49];
        b[50] = nodePressure[50] * nodeVolume[50];
        b[51] = nodePressure[51] * nodeVolume[51];
        b[52] = nodePressure[52] * nodeVolume[52];
        b[53] = nodePressure[53] * nodeVolume[53];
        b[54] = nodePressure[54] * nodeVolume[54];
        b[55] = nodePressure[55] * nodeVolume[55];
        b[56] = nodePressure[56] * nodeVolume[56];
        b[57] = nodePressure[57] * nodeVolume[57];
        b[58] = nodePressure[58] * nodeVolume[58];
        b[59] = nodePressure[59] * nodeVolume[59];
        b[60] = nodePressure[60] * nodeVolume[60];
        b[61] = nodePressure[61] * nodeVolume[61];
        b[62] = nodePressure[62] * nodeVolume[62];
        b[63] = nodePressure[63] * nodeVolume[63];
        b[64] = nodePressure[64] * nodeVolume[64];
        b[65] = nodePressure[65] * nodeVolume[65];
        b[66] = nodePressure[66] * nodeVolume[66];
        b[67] = nodePressure[67] * nodeVolume[67];
        b[68] = nodePressure[68] * nodeVolume[68];
        b[69] = nodePressure[69] * nodeVolume[69];
        b[70] = nodePressure[70] * nodeVolume[70];
        b[71] = nodePressure[71] * nodeVolume[71];
        b[72] = nodePressure[72] * nodeVolume[72];
        b[73] = nodePressure[73] * nodeVolume[73];
        b[74] = nodePressure[74] * nodeVolume[74];
        b[75] = nodePressure[75] * nodeVolume[75];
        b[76] = nodePressure[76] * nodeVolume[76];
        b[77] = nodePressure[77] * nodeVolume[77];
        b[78] = nodePressure[78] * nodeVolume[78];
        b[79] = nodePressure[79] * nodeVolume[79];
        b[80] = nodePressure[80] * nodeVolume[80];
        b[81] = nodePressure[81] * nodeVolume[81];

        float basePiv, m;
        m = _a00_11 / _a00_00; // eliminating 00 current 11
        _a11_11 = _a11_11 - _a00_11 * m;
        b[11] = b[11] - b[0] * m;
        m = _a01_11 / _a01_01; // eliminating 01 current 11
        _a11_11 = _a11_11 - _a01_11 * m;
        b[11] = b[11] - b[1] * m;
        m = _a02_28 / _a02_02; // eliminating 02 current 28
        _a28_28 = _a28_28 - _a02_28 * m;
        b[28] = b[28] - b[2] * m;
        m = _a03_29 / _a03_03; // eliminating 03 current 29
        _a29_29 = _a29_29 - _a03_29 * m;
        b[29] = b[29] - b[3] * m;
        m = _a04_13 / _a04_04; // eliminating 04 current 13
        _a13_13 = _a13_13 - _a04_13 * m;
        b[13] = b[13] - b[4] * m;
        m = _a05_14 / _a05_05; // eliminating 05 current 14
        _a14_14 = _a14_14 - _a05_14 * m;
        b[14] = b[14] - b[5] * m;
        basePiv = _a06_06;
        m = _a06_39 / basePiv; // eliminating 06 current 39
        _a39_39 = _a39_39 - _a06_39 * m;
        _a39_81 = -_a06_81 * m;
        b[39] = b[39] - b[6] * m;
        m = _a06_81 / basePiv; // eliminating 06 current 81
        _a81_81 = _a81_81 - _a06_81 * m;
        b[81] = b[81] - b[6] * m;
        m = _a07_24 / _a07_07; // eliminating 07 current 24
        _a24_24 = _a24_24 - _a07_24 * m;
        b[24] = b[24] - b[7] * m;
        m = _a08_23 / _a08_08; // eliminating 08 current 23
        _a23_23 = _a23_23 - _a08_23 * m;
        b[23] = b[23] - b[8] * m;
        m = _a09_29 / _a09_09; // eliminating 09 current 29
        _a29_29 = _a29_29 - _a09_29 * m;
        b[29] = b[29] - b[9] * m;
        basePiv = _a10_10;
        m = _a10_20 / basePiv; // eliminating 10 current 20
        _a20_20 = _a20_20 - _a10_20 * m;
        _a20_80 = -_a10_80 * m;
        b[20] = b[20] - b[10] * m;
        m = _a10_80 / basePiv; // eliminating 10 current 80
        _a80_80 = _a80_80 - _a10_80 * m;
        b[80] = b[80] - b[10] * m;
        basePiv = _a12_12;
        m = _a12_28 / basePiv; // eliminating 12 current 28
        _a28_28 = _a28_28 - _a12_28 * m;
        _a28_35 = -_a12_35 * m;
        b[28] = b[28] - b[12] * m;
        m = _a12_35 / basePiv; // eliminating 12 current 35
        _a35_35 = _a35_35 - _a12_35 * m;
        b[35] = b[35] - b[12] * m;
        m = _a13_25 / _a13_13; // eliminating 13 current 25
        _a25_25 = _a25_25 - _a13_25 * m;
        b[25] = b[25] - b[13] * m;
        m = _a14_26 / _a14_14; // eliminating 14 current 26
        _a26_26 = _a26_26 - _a14_26 * m;
        b[26] = b[26] - b[14] * m;
        basePiv = _a15_15;
        m = _a15_30 / basePiv; // eliminating 15 current 30
        _a30_30 = _a30_30 - _a15_30 * m;
        _a30_34 = -_a15_34 * m;
        b[30] = b[30] - b[15] * m;
        m = _a15_34 / basePiv; // eliminating 15 current 34
        _a34_34 = _a34_34 - _a15_34 * m;
        b[34] = b[34] - b[15] * m;
        basePiv = _a16_16;
        m = _a16_31 / basePiv; // eliminating 16 current 31
        _a31_31 = _a31_31 - _a16_31 * m;
        _a31_36 = -_a16_36 * m;
        b[31] = b[31] - b[16] * m;
        m = _a16_36 / basePiv; // eliminating 16 current 36
        _a36_36 = _a36_36 - _a16_36 * m;
        b[36] = b[36] - b[16] * m;
        basePiv = _a17_17;
        m = _a17_38 / basePiv; // eliminating 17 current 38
        _a38_38 = _a38_38 - _a17_38 * m;
        _a38_39 = -_a17_39 * m;
        b[38] = b[38] - b[17] * m;
        m = _a17_39 / basePiv; // eliminating 17 current 39
        _a39_39 = _a39_39 - _a17_39 * m;
        b[39] = b[39] - b[17] * m;
        basePiv = _a18_18;
        m = _a18_27 / basePiv; // eliminating 18 current 27
        _a27_27 = _a27_27 - _a18_27 * m;
        _a27_33 = -_a18_33 * m;
        b[27] = b[27] - b[18] * m;
        m = _a18_33 / basePiv; // eliminating 18 current 33
        _a33_33 = _a33_33 - _a18_33 * m;
        b[33] = b[33] - b[18] * m;
        basePiv = _a19_19;
        m = _a19_32 / basePiv; // eliminating 19 current 32
        _a32_32 = _a32_32 - _a19_32 * m;
        _a32_35 = -_a19_35 * m;
        b[32] = b[32] - b[19] * m;
        m = _a19_35 / basePiv; // eliminating 19 current 35
        _a35_35 = _a35_35 - _a19_35 * m;
        b[35] = b[35] - b[19] * m;
        basePiv = _a20_20;
        m = _a20_27 / basePiv; // eliminating 20 current 27
        _a27_27 = _a27_27 - _a20_27 * m;
        _a27_80 = -_a20_80 * m;
        b[27] = b[27] - b[20] * m;
        m = _a20_80 / basePiv; // eliminating 20 current 80
        _a80_80 = _a80_80 - _a20_80 * m;
        b[80] = b[80] - b[20] * m;
        basePiv = _a21_21;
        m = _a21_37 / basePiv; // eliminating 21 current 37
        _a37_37 = _a37_37 - _a21_37 * m;
        _a37_39 = -_a21_39 * m;
        b[37] = b[37] - b[21] * m;
        m = _a21_39 / basePiv; // eliminating 21 current 39
        _a39_39 = _a39_39 - _a21_39 * m;
        b[39] = b[39] - b[21] * m;
        basePiv = _a22_22;
        m = _a22_23 / basePiv; // eliminating 22 current 23
        _a23_23 = _a23_23 - _a22_23 * m;
        _a23_39 = -_a22_39 * m;
        b[23] = b[23] - b[22] * m;
        m = _a22_39 / basePiv; // eliminating 22 current 39
        _a39_39 = _a39_39 - _a22_39 * m;
        b[39] = b[39] - b[22] * m;
        m = _a23_39 / _a23_23; // eliminating 23 current 39
        _a39_39 = _a39_39 - _a23_39 * m;
        b[39] = b[39] - b[23] * m;
        basePiv = _a24_24;
        m = _a24_32 / basePiv; // eliminating 24 current 32
        _a32_32 = _a32_32 - _a24_32 * m;
        _a32_38 = -_a24_38 * m;
        b[32] = b[32] - b[24] * m;
        m = _a24_38 / basePiv; // eliminating 24 current 38
        _a38_38 = _a38_38 - _a24_38 * m;
        b[38] = b[38] - b[24] * m;
        m = _a25_39 / _a25_25; // eliminating 25 current 39
        _a39_39 = _a39_39 - _a25_39 * m;
        b[39] = b[39] - b[25] * m;
        m = _a26_34 / _a26_26; // eliminating 26 current 34
        _a34_34 = _a34_34 - _a26_34 * m;
        b[34] = b[34] - b[26] * m;
        basePiv = _a27_27;
        m = _a27_33 / basePiv; // eliminating 27 current 33
        _a33_33 = _a33_33 - _a27_33 * m;
        _a33_80 = -_a27_80 * m;
        b[33] = b[33] - b[27] * m;
        m = _a27_80 / basePiv; // eliminating 27 current 80
        _a80_80 = _a80_80 - _a27_80 * m;
        b[80] = b[80] - b[27] * m;
        basePiv = _a28_28;
        m = _a28_29 / basePiv; // eliminating 28 current 29
        _a29_29 = _a29_29 - _a28_29 * m;
        _a29_35 = _a29_35 - _a28_35 * m;
        b[29] = b[29] - b[28] * m;
        m = _a28_35 / basePiv; // eliminating 28 current 35
        _a35_35 = _a35_35 - _a28_35 * m;
        b[35] = b[35] - b[28] * m;
        m = _a29_35 / _a29_29; // eliminating 29 current 35
        _a35_35 = _a35_35 - _a29_35 * m;
        b[35] = b[35] - b[29] * m;
        basePiv = _a30_30;
        m = _a30_34 / basePiv; // eliminating 30 current 34
        _a34_34 = _a34_34 - _a30_34 * m;
        _a34_36 = -_a30_36 * m;
        b[34] = b[34] - b[30] * m;
        m = _a30_36 / basePiv; // eliminating 30 current 36
        _a36_36 = _a36_36 - _a30_36 * m;
        b[36] = b[36] - b[30] * m;
        basePiv = _a31_31;
        m = _a31_36 / basePiv; // eliminating 31 current 36
        _a36_36 = _a36_36 - _a31_36 * m;
        _a36_39 = -_a31_39 * m;
        b[36] = b[36] - b[31] * m;
        m = _a31_39 / basePiv; // eliminating 31 current 39
        _a39_39 = _a39_39 - _a31_39 * m;
        b[39] = b[39] - b[31] * m;
        basePiv = _a32_32;
        m = _a32_35 / basePiv; // eliminating 32 current 35
        _a35_35 = _a35_35 - _a32_35 * m;
        _a35_38 = -_a32_38 * m;
        b[35] = b[35] - b[32] * m;
        m = _a32_38 / basePiv; // eliminating 32 current 38
        _a38_38 = _a38_38 - _a32_38 * m;
        b[38] = b[38] - b[32] * m;
        basePiv = _a33_33;
        m = _a33_37 / basePiv; // eliminating 33 current 37
        _a37_37 = _a37_37 - _a33_37 * m;
        _a37_80 = -_a33_80 * m;
        b[37] = b[37] - b[33] * m;
        m = _a33_80 / basePiv; // eliminating 33 current 80
        _a80_80 = _a80_80 - _a33_80 * m;
        b[80] = b[80] - b[33] * m;
        m = _a34_36 / _a34_34; // eliminating 34 current 36
        _a36_36 = _a36_36 - _a34_36 * m;
        b[36] = b[36] - b[34] * m;
        basePiv = _a35_35;
        m = _a35_37 / basePiv; // eliminating 35 current 37
        _a37_37 = _a37_37 - _a35_37 * m;
        _a37_38 = -_a35_38 * m;
        b[37] = b[37] - b[35] * m;
        m = _a35_38 / basePiv; // eliminating 35 current 38
        _a38_38 = _a38_38 - _a35_38 * m;
        b[38] = b[38] - b[35] * m;
        m = _a36_39 / _a36_36; // eliminating 36 current 39
        _a39_39 = _a39_39 - _a36_39 * m;
        b[39] = b[39] - b[36] * m;
        basePiv = _a37_37;
        m = _a37_38 / basePiv; // eliminating 37 current 38
        _a38_38 = _a38_38 - _a37_38 * m;
        _a38_39 = _a38_39 - _a37_39 * m;
        _a38_80 = -_a37_80 * m;
        b[38] = b[38] - b[37] * m;
        m = _a37_39 / basePiv; // eliminating 37 current 39
        _a39_39 = _a39_39 - _a37_39 * m;
        _a39_80 = -_a37_80 * m;
        b[39] = b[39] - b[37] * m;
        m = _a37_80 / basePiv; // eliminating 37 current 80
        _a80_80 = _a80_80 - _a37_80 * m;
        b[80] = b[80] - b[37] * m;
        basePiv = _a38_38;
        m = _a38_39 / basePiv; // eliminating 38 current 39
        _a39_39 = _a39_39 - _a38_39 * m;
        _a39_80 = _a39_80 - _a38_80 * m;
        b[39] = b[39] - b[38] * m;
        m = _a38_80 / basePiv; // eliminating 38 current 80
        _a80_80 = _a80_80 - _a38_80 * m;
        b[80] = b[80] - b[38] * m;
        basePiv = _a39_39;
        m = _a39_80 / basePiv; // eliminating 39 current 80
        _a80_80 = _a80_80 - _a39_80 * m;
        _a80_81 = -_a39_81 * m;
        b[80] = b[80] - b[39] * m;
        m = _a39_81 / basePiv; // eliminating 39 current 81
        _a81_81 = _a81_81 - _a39_81 * m;
        b[81] = b[81] - b[39] * m;
        m = _a40_74 / _a40_40; // eliminating 40 current 74
        _a74_74 = _a74_74 - _a40_74 * m;
        b[74] = b[74] - b[40] * m;
        m = _a41_64 / _a41_41; // eliminating 41 current 64
        _a64_64 = _a64_64 - _a41_64 * m;
        b[64] = b[64] - b[41] * m;
        m = _a42_64 / _a42_42; // eliminating 42 current 64
        _a64_64 = _a64_64 - _a42_64 * m;
        b[64] = b[64] - b[42] * m;
        m = _a43_50 / _a43_43; // eliminating 43 current 50
        _a50_50 = _a50_50 - _a43_50 * m;
        b[50] = b[50] - b[43] * m;
        m = _a44_62 / _a44_44; // eliminating 44 current 62
        _a62_62 = _a62_62 - _a44_62 * m;
        b[62] = b[62] - b[44] * m;
        m = _a45_62 / _a45_45; // eliminating 45 current 62
        _a62_62 = _a62_62 - _a45_62 * m;
        b[62] = b[62] - b[45] * m;
        basePiv = _a46_46;
        m = _a46_60 / basePiv; // eliminating 46 current 60
        _a60_60 = _a60_60 - _a46_60 * m;
        _a60_81 = -_a46_81 * m;
        b[60] = b[60] - b[46] * m;
        m = _a46_81 / basePiv; // eliminating 46 current 81
        _a81_81 = _a81_81 - _a46_81 * m;
        b[81] = b[81] - b[46] * m;
        m = _a47_61 / _a47_47; // eliminating 47 current 61
        _a61_61 = _a61_61 - _a47_61 * m;
        b[61] = b[61] - b[47] * m;
        basePiv = _a48_48;
        m = _a48_74 / basePiv; // eliminating 48 current 74
        _a74_74 = _a74_74 - _a48_74 * m;
        _a74_77 = -_a48_77 * m;
        b[74] = b[74] - b[48] * m;
        m = _a48_77 / basePiv; // eliminating 48 current 77
        _a77_77 = _a77_77 - _a48_77 * m;
        b[77] = b[77] - b[48] * m;
        basePiv = _a49_49;
        m = _a49_64 / basePiv; // eliminating 49 current 64
        _a64_64 = _a64_64 - _a49_64 * m;
        _a64_72 = -_a49_72 * m;
        b[64] = b[64] - b[49] * m;
        m = _a49_72 / basePiv; // eliminating 49 current 72
        _a72_72 = _a72_72 - _a49_72 * m;
        b[72] = b[72] - b[49] * m;
        m = _a50_63 / _a50_50; // eliminating 50 current 63
        _a63_63 = _a63_63 - _a50_63 * m;
        b[63] = b[63] - b[50] * m;
        basePiv = _a51_51;
        m = _a51_65 / basePiv; // eliminating 51 current 65
        _a65_65 = _a65_65 - _a51_65 * m;
        _a65_71 = -_a51_71 * m;
        b[65] = b[65] - b[51] * m;
        m = _a51_71 / basePiv; // eliminating 51 current 71
        _a71_71 = _a71_71 - _a51_71 * m;
        b[71] = b[71] - b[51] * m;
        basePiv = _a52_52;
        m = _a52_73 / basePiv; // eliminating 52 current 73
        _a73_73 = _a73_73 - _a52_73 * m;
        _a73_78 = -_a52_78 * m;
        b[73] = b[73] - b[52] * m;
        m = _a52_78 / basePiv; // eliminating 52 current 78
        _a78_78 = _a78_78 - _a52_78 * m;
        b[78] = b[78] - b[52] * m;
        basePiv = _a53_53;
        m = _a53_66 / basePiv; // eliminating 53 current 66
        _a66_66 = _a66_66 - _a53_66 * m;
        _a66_73 = -_a53_73 * m;
        b[66] = b[66] - b[53] * m;
        m = _a53_73 / basePiv; // eliminating 53 current 73
        _a73_73 = _a73_73 - _a53_73 * m;
        b[73] = b[73] - b[53] * m;
        basePiv = _a54_54;
        m = _a54_67 / basePiv; // eliminating 54 current 67
        _a67_67 = _a67_67 - _a54_67 * m;
        _a67_79 = -_a54_79 * m;
        b[67] = b[67] - b[54] * m;
        m = _a54_79 / basePiv; // eliminating 54 current 79
        _a79_79 = _a79_79 - _a54_79 * m;
        b[79] = b[79] - b[54] * m;
        basePiv = _a55_55;
        m = _a55_68 / basePiv; // eliminating 55 current 68
        _a68_68 = _a68_68 - _a55_68 * m;
        _a68_79 = -_a55_79 * m;
        b[68] = b[68] - b[55] * m;
        m = _a55_79 / basePiv; // eliminating 55 current 79
        _a79_79 = _a79_79 - _a55_79 * m;
        b[79] = b[79] - b[55] * m;
        basePiv = _a56_56;
        m = _a56_69 / basePiv; // eliminating 56 current 69
        _a69_69 = _a69_69 - _a56_69 * m;
        _a69_78 = -_a56_78 * m;
        b[69] = b[69] - b[56] * m;
        m = _a56_78 / basePiv; // eliminating 56 current 78
        _a78_78 = _a78_78 - _a56_78 * m;
        b[78] = b[78] - b[56] * m;
        basePiv = _a57_57;
        m = _a57_69 / basePiv; // eliminating 57 current 69
        _a69_69 = _a69_69 - _a57_69 * m;
        _a69_76 = -_a57_76 * m;
        b[69] = b[69] - b[57] * m;
        m = _a57_76 / basePiv; // eliminating 57 current 76
        _a76_76 = _a76_76 - _a57_76 * m;
        b[76] = b[76] - b[57] * m;
        basePiv = _a58_58;
        m = _a58_70 / basePiv; // eliminating 58 current 70
        _a70_70 = _a70_70 - _a58_70 * m;
        _a70_75 = -_a58_75 * m;
        b[70] = b[70] - b[58] * m;
        m = _a58_75 / basePiv; // eliminating 58 current 75
        _a75_75 = _a75_75 - _a58_75 * m;
        b[75] = b[75] - b[58] * m;
        basePiv = _a59_59;
        m = _a59_72 / basePiv; // eliminating 59 current 72
        _a72_72 = _a72_72 - _a59_72 * m;
        _a72_77 = -_a59_77 * m;
        b[72] = b[72] - b[59] * m;
        m = _a59_77 / basePiv; // eliminating 59 current 77
        _a77_77 = _a77_77 - _a59_77 * m;
        b[77] = b[77] - b[59] * m;
        basePiv = _a60_60;
        m = _a60_76 / basePiv; // eliminating 60 current 76
        _a76_76 = _a76_76 - _a60_76 * m;
        _a76_81 = -_a60_81 * m;
        b[76] = b[76] - b[60] * m;
        m = _a60_81 / basePiv; // eliminating 60 current 81
        _a81_81 = _a81_81 - _a60_81 * m;
        b[81] = b[81] - b[60] * m;
        basePiv = _a61_61;
        m = _a61_70 / basePiv; // eliminating 61 current 70
        _a70_70 = _a70_70 - _a61_70 * m;
        _a70_74 = -_a61_74 * m;
        b[70] = b[70] - b[61] * m;
        m = _a61_74 / basePiv; // eliminating 61 current 74
        _a74_74 = _a74_74 - _a61_74 * m;
        b[74] = b[74] - b[61] * m;
        m = _a62_78 / _a62_62; // eliminating 62 current 78
        _a78_78 = _a78_78 - _a62_78 * m;
        b[78] = b[78] - b[62] * m;
        m = _a63_71 / _a63_63; // eliminating 63 current 71
        _a71_71 = _a71_71 - _a63_71 * m;
        b[71] = b[71] - b[63] * m;
        basePiv = _a64_64;
        m = _a64_72 / basePiv; // eliminating 64 current 72
        _a72_72 = _a72_72 - _a64_72 * m;
        _a72_75 = _a72_75 - _a64_75 * m;
        b[72] = b[72] - b[64] * m;
        m = _a64_75 / basePiv; // eliminating 64 current 75
        _a75_75 = _a75_75 - _a64_75 * m;
        b[75] = b[75] - b[64] * m;
        basePiv = _a65_65;
        m = _a65_71 / basePiv; // eliminating 65 current 71
        _a71_71 = _a71_71 - _a65_71 * m;
        _a71_76 = -_a65_76 * m;
        b[71] = b[71] - b[65] * m;
        m = _a65_76 / basePiv; // eliminating 65 current 76
        _a76_76 = _a76_76 - _a65_76 * m;
        b[76] = b[76] - b[65] * m;
        basePiv = _a66_66;
        m = _a66_73 / basePiv; // eliminating 66 current 73
        _a73_73 = _a73_73 - _a66_73 * m;
        _a73_79 = -_a66_79 * m;
        b[73] = b[73] - b[66] * m;
        m = _a66_79 / basePiv; // eliminating 66 current 79
        _a79_79 = _a79_79 - _a66_79 * m;
        b[79] = b[79] - b[66] * m;
        basePiv = _a67_67;
        m = _a67_76 / basePiv; // eliminating 67 current 76
        _a76_76 = _a76_76 - _a67_76 * m;
        _a76_79 = -_a67_79 * m;
        b[76] = b[76] - b[67] * m;
        m = _a67_79 / basePiv; // eliminating 67 current 79
        _a79_79 = _a79_79 - _a67_79 * m;
        b[79] = b[79] - b[67] * m;
        basePiv = _a68_68;
        m = _a68_77 / basePiv; // eliminating 68 current 77
        _a77_77 = _a77_77 - _a68_77 * m;
        _a77_79 = -_a68_79 * m;
        b[77] = b[77] - b[68] * m;
        m = _a68_79 / basePiv; // eliminating 68 current 79
        _a79_79 = _a79_79 - _a68_79 * m;
        b[79] = b[79] - b[68] * m;
        basePiv = _a69_69;
        m = _a69_76 / basePiv; // eliminating 69 current 76
        _a76_76 = _a76_76 - _a69_76 * m;
        _a76_78 = -_a69_78 * m;
        b[76] = b[76] - b[69] * m;
        m = _a69_78 / basePiv; // eliminating 69 current 78
        _a78_78 = _a78_78 - _a69_78 * m;
        b[78] = b[78] - b[69] * m;
        basePiv = _a70_70;
        m = _a70_74 / basePiv; // eliminating 70 current 74
        _a74_74 = _a74_74 - _a70_74 * m;
        _a74_75 = -_a70_75 * m;
        b[74] = b[74] - b[70] * m;
        m = _a70_75 / basePiv; // eliminating 70 current 75
        _a75_75 = _a75_75 - _a70_75 * m;
        b[75] = b[75] - b[70] * m;
        m = _a71_76 / _a71_71; // eliminating 71 current 76
        _a76_76 = _a76_76 - _a71_76 * m;
        b[76] = b[76] - b[71] * m;
        basePiv = _a72_72;
        m = _a72_75 / basePiv; // eliminating 72 current 75
        _a75_75 = _a75_75 - _a72_75 * m;
        _a75_77 = -_a72_77 * m;
        b[75] = b[75] - b[72] * m;
        m = _a72_77 / basePiv; // eliminating 72 current 77
        _a77_77 = _a77_77 - _a72_77 * m;
        b[77] = b[77] - b[72] * m;
        basePiv = _a73_73;
        m = _a73_78 / basePiv; // eliminating 73 current 78
        _a78_78 = _a78_78 - _a73_78 * m;
        _a78_79 = -_a73_79 * m;
        _a78_80 = -_a73_80 * m;
        b[78] = b[78] - b[73] * m;
        m = _a73_79 / basePiv; // eliminating 73 current 79
        _a79_79 = _a79_79 - _a73_79 * m;
        _a79_80 = -_a73_80 * m;
        b[79] = b[79] - b[73] * m;
        m = _a73_80 / basePiv; // eliminating 73 current 80
        _a80_80 = _a80_80 - _a73_80 * m;
        b[80] = b[80] - b[73] * m;
        basePiv = _a74_74;
        m = _a74_75 / basePiv; // eliminating 74 current 75
        _a75_75 = _a75_75 - _a74_75 * m;
        _a75_77 = _a75_77 - _a74_77 * m;
        b[75] = b[75] - b[74] * m;
        m = _a74_77 / basePiv; // eliminating 74 current 77
        _a77_77 = _a77_77 - _a74_77 * m;
        b[77] = b[77] - b[74] * m;
        m = _a75_77 / _a75_75; // eliminating 75 current 77
        _a77_77 = _a77_77 - _a75_77 * m;
        b[77] = b[77] - b[75] * m;
        basePiv = _a76_76;
        m = _a76_78 / basePiv; // eliminating 76 current 78
        _a78_78 = _a78_78 - _a76_78 * m;
        _a78_79 = _a78_79 - _a76_79 * m;
        _a78_81 = -_a76_81 * m;
        b[78] = b[78] - b[76] * m;
        m = _a76_79 / basePiv; // eliminating 76 current 79
        _a79_79 = _a79_79 - _a76_79 * m;
        _a79_81 = -_a76_81 * m;
        b[79] = b[79] - b[76] * m;
        m = _a76_81 / basePiv; // eliminating 76 current 81
        _a81_81 = _a81_81 - _a76_81 * m;
        b[81] = b[81] - b[76] * m;
        m = _a77_79 / _a77_77; // eliminating 77 current 79
        _a79_79 = _a79_79 - _a77_79 * m;
        b[79] = b[79] - b[77] * m;
        basePiv = _a78_78;
        m = _a78_79 / basePiv; // eliminating 78 current 79
        _a79_79 = _a79_79 - _a78_79 * m;
        _a79_80 = _a79_80 - _a78_80 * m;
        _a79_81 = _a79_81 - _a78_81 * m;
        b[79] = b[79] - b[78] * m;
        m = _a78_80 / basePiv; // eliminating 78 current 80
        _a80_80 = _a80_80 - _a78_80 * m;
        _a80_81 = _a80_81 - _a78_81 * m;
        b[80] = b[80] - b[78] * m;
        m = _a78_81 / basePiv; // eliminating 78 current 81
        _a81_81 = _a81_81 - _a78_81 * m;
        b[81] = b[81] - b[78] * m;
        basePiv = _a79_79;
        m = _a79_80 / basePiv; // eliminating 79 current 80
        _a80_80 = _a80_80 - _a79_80 * m;
        _a80_81 = _a80_81 - _a79_81 * m;
        b[80] = b[80] - b[79] * m;
        m = _a79_81 / basePiv; // eliminating 79 current 81
        _a81_81 = _a81_81 - _a79_81 * m;
        b[81] = b[81] - b[79] * m;
        m = _a80_81 / _a80_80; // eliminating 80 current 81
        _a81_81 = _a81_81 - _a80_81 * m;
        b[81] = b[81] - b[80] * m;

        b[81] = b[81] / _a81_81;
        b[80] = (b[80] - b[81] * _a80_81) / _a80_80;
        b[79] = (b[79] - b[80] * _a79_80 - b[81] * _a79_81) / _a79_79;
        b[78] = (b[78] - b[79] * _a78_79 - b[80] * _a78_80 - b[81] * _a78_81) / _a78_78;
        b[77] = (b[77] - b[79] * _a77_79) / _a77_77;
        b[76] = (b[76] - b[78] * _a76_78 - b[79] * _a76_79 - b[81] * _a76_81) / _a76_76;
        b[75] = (b[75] - b[77] * _a75_77) / _a75_75;
        b[74] = (b[74] - b[75] * _a74_75 - b[77] * _a74_77) / _a74_74;
        b[73] = (b[73] - b[78] * _a73_78 - b[79] * _a73_79 - b[80] * _a73_80) / _a73_73;
        b[72] = (b[72] - b[75] * _a72_75 - b[77] * _a72_77) / _a72_72;
        b[71] = (b[71] - b[76] * _a71_76) / _a71_71;
        b[70] = (b[70] - b[74] * _a70_74 - b[75] * _a70_75) / _a70_70;
        b[69] = (b[69] - b[76] * _a69_76 - b[78] * _a69_78) / _a69_69;
        b[68] = (b[68] - b[77] * _a68_77 - b[79] * _a68_79) / _a68_68;
        b[67] = (b[67] - b[76] * _a67_76 - b[79] * _a67_79) / _a67_67;
        b[66] = (b[66] - b[73] * _a66_73 - b[79] * _a66_79) / _a66_66;
        b[65] = (b[65] - b[71] * _a65_71 - b[76] * _a65_76) / _a65_65;
        b[64] = (b[64] - b[72] * _a64_72 - b[75] * _a64_75) / _a64_64;
        b[63] = (b[63] - b[71] * _a63_71) / _a63_63;
        b[62] = (b[62] - b[78] * _a62_78) / _a62_62;
        b[61] = (b[61] - b[70] * _a61_70 - b[74] * _a61_74) / _a61_61;
        b[60] = (b[60] - b[76] * _a60_76 - b[81] * _a60_81) / _a60_60;
        b[59] = (b[59] - b[72] * _a59_72 - b[77] * _a59_77) / _a59_59;
        b[58] = (b[58] - b[70] * _a58_70 - b[75] * _a58_75) / _a58_58;
        b[57] = (b[57] - b[69] * _a57_69 - b[76] * _a57_76) / _a57_57;
        b[56] = (b[56] - b[69] * _a56_69 - b[78] * _a56_78) / _a56_56;
        b[55] = (b[55] - b[68] * _a55_68 - b[79] * _a55_79) / _a55_55;
        b[54] = (b[54] - b[67] * _a54_67 - b[79] * _a54_79) / _a54_54;
        b[53] = (b[53] - b[66] * _a53_66 - b[73] * _a53_73) / _a53_53;
        b[52] = (b[52] - b[73] * _a52_73 - b[78] * _a52_78) / _a52_52;
        b[51] = (b[51] - b[65] * _a51_65 - b[71] * _a51_71) / _a51_51;
        b[50] = (b[50] - b[63] * _a50_63) / _a50_50;
        b[49] = (b[49] - b[64] * _a49_64 - b[72] * _a49_72) / _a49_49;
        b[48] = (b[48] - b[74] * _a48_74 - b[77] * _a48_77) / _a48_48;
        b[47] = (b[47] - b[61] * _a47_61) / _a47_47;
        b[46] = (b[46] - b[60] * _a46_60 - b[81] * _a46_81) / _a46_46;
        b[45] = (b[45] - b[62] * _a45_62) / _a45_45;
        b[44] = (b[44] - b[62] * _a44_62) / _a44_44;
        b[43] = (b[43] - b[50] * _a43_50) / _a43_43;
        b[42] = (b[42] - b[64] * _a42_64) / _a42_42;
        b[41] = (b[41] - b[64] * _a41_64) / _a41_41;
        b[40] = (b[40] - b[74] * _a40_74) / _a40_40;
        b[39] = (b[39] - b[80] * _a39_80 - b[81] * _a39_81) / _a39_39;
        b[38] = (b[38] - b[39] * _a38_39 - b[80] * _a38_80) / _a38_38;
        b[37] = (b[37] - b[38] * _a37_38 - b[39] * _a37_39 - b[80] * _a37_80) / _a37_37;
        b[36] = (b[36] - b[39] * _a36_39) / _a36_36;
        b[35] = (b[35] - b[37] * _a35_37 - b[38] * _a35_38) / _a35_35;
        b[34] = (b[34] - b[36] * _a34_36) / _a34_34;
        b[33] = (b[33] - b[37] * _a33_37 - b[80] * _a33_80) / _a33_33;
        b[32] = (b[32] - b[35] * _a32_35 - b[38] * _a32_38) / _a32_32;
        b[31] = (b[31] - b[36] * _a31_36 - b[39] * _a31_39) / _a31_31;
        b[30] = (b[30] - b[34] * _a30_34 - b[36] * _a30_36) / _a30_30;
        b[29] = (b[29] - b[35] * _a29_35) / _a29_29;
        b[28] = (b[28] - b[29] * _a28_29 - b[35] * _a28_35) / _a28_28;
        b[27] = (b[27] - b[33] * _a27_33 - b[80] * _a27_80) / _a27_27;
        b[26] = (b[26] - b[34] * _a26_34) / _a26_26;
        b[25] = (b[25] - b[39] * _a25_39) / _a25_25;
        b[24] = (b[24] - b[32] * _a24_32 - b[38] * _a24_38) / _a24_24;
        b[23] = (b[23] - b[39] * _a23_39) / _a23_23;
        b[22] = (b[22] - b[23] * _a22_23 - b[39] * _a22_39) / _a22_22;
        b[21] = (b[21] - b[37] * _a21_37 - b[39] * _a21_39) / _a21_21;
        b[20] = (b[20] - b[27] * _a20_27 - b[80] * _a20_80) / _a20_20;
        b[19] = (b[19] - b[32] * _a19_32 - b[35] * _a19_35) / _a19_19;
        b[18] = (b[18] - b[27] * _a18_27 - b[33] * _a18_33) / _a18_18;
        b[17] = (b[17] - b[38] * _a17_38 - b[39] * _a17_39) / _a17_17;
        b[16] = (b[16] - b[31] * _a16_31 - b[36] * _a16_36) / _a16_16;
        b[15] = (b[15] - b[30] * _a15_30 - b[34] * _a15_34) / _a15_15;
        b[14] = (b[14] - b[26] * _a14_26) / _a14_14;
        b[13] = (b[13] - b[25] * _a13_25) / _a13_13;
        b[12] = (b[12] - b[28] * _a12_28 - b[35] * _a12_35) / _a12_12;
        b[11] = b[11] / _a11_11;
        b[10] = (b[10] - b[20] * _a10_20 - b[80] * _a10_80) / _a10_10;
        b[9] = (b[9] - b[29] * _a09_29) / _a09_09;
        b[8] = (b[8] - b[23] * _a08_23) / _a08_08;
        b[7] = (b[7] - b[24] * _a07_24) / _a07_07;
        b[6] = (b[6] - b[39] * _a06_39 - b[81] * _a06_81) / _a06_06;
        b[5] = (b[5] - b[14] * _a05_14) / _a05_05;
        b[4] = (b[4] - b[13] * _a04_13) / _a04_04;
        b[3] = (b[3] - b[29] * _a03_29) / _a03_03;
        b[2] = (b[2] - b[28] * _a02_28) / _a02_02;
        b[1] = (b[1] - b[11] * _a01_11) / _a01_01;
        b[0] = (b[0] - b[11] * _a00_11) / _a00_00;

    }

    void InitSpec()
    {
        nodeIds = new LIntMap(); nodeIds.Reserve(82);
        nodeIds.Set("k254_1.r", 0);
        nodeIds.Set("k254_1.vr", 1);
        nodeIds.Set("k395_1.s", 2);
        nodeIds.Set("k395_1.d", 3);
        nodeIds.Set("epka_1.kvv", 4);
        nodeIds.Set("tc_6_12", 5);
        nodeIds.Set("red1.in", 6);
        nodeIds.Set("epka_1.tm", 7);
        nodeIds.Set("kp3", 8);
        nodeIds.Set("ur", 9);
        nodeIds.Set("mvt", 10);
        nodeIds.Set("k254_1.b", 11);
        nodeIds.Set("k395_1.gr", 12);
        nodeIds.Set("epka_1.k", 13);
        nodeIds.Set("tc_5_11", 14);
        nodeIds.Set("rd2.out", 15);
        nodeIds.Set("ko4.out", 16);
        nodeIds.Set("ko3.out", 17);
        nodeIds.Set("k254_1.tc", 18);
        nodeIds.Set("ubt_1.k_tm", 19);
        nodeIds.Set("ubt_1.tc", 20);
        nodeIds.Set("ubt_1.gr", 21);
        nodeIds.Set("rs2", 22);
        nodeIds.Set("rs1", 23);
        nodeIds.Set("tm", 24);
        nodeIds.Set("epka_1.gr", 25);
        nodeIds.Set("tc_4_10", 26);
        nodeIds.Set("ubt_1.k_tc", 27);
        nodeIds.Set("k395_1.u", 28);
        nodeIds.Set("k395_1.ur", 29);
        nodeIds.Set("rd2.feed", 30);
        nodeIds.Set("ko4.in", 31);
        nodeIds.Set("ubt_1.tm", 32);
        nodeIds.Set("k254_1.gr", 33);
        nodeIds.Set("tc_1", 34);
        nodeIds.Set("k395_1.tm", 35);
        nodeIds.Set("rs4", 36);
        nodeIds.Set("ubt_1.k_gr", 37);
        nodeIds.Set("ko3.in", 38);
        nodeIds.Set("pm", 39);
        nodeIds.Set("evr.rk", 40);
        nodeIds.Set("vr.kdr", 41);
        nodeIds.Set("vr.tm", 42);
        nodeIds.Set("tc_1_7", 43);
        nodeIds.Set("rd1.set", 44);
        nodeIds.Set("rd2.set", 45);
        nodeIds.Set("ko2.in", 46);
        nodeIds.Set("rs5", 47);
        nodeIds.Set("evr.rd_tc", 48);
        nodeIds.Set("vr.sk", 49);
        nodeIds.Set("tc_2_8", 50);
        nodeIds.Set("rd1.out", 51);
        nodeIds.Set("k3.b", 52);
        nodeIds.Set("k5.a", 53);
        nodeIds.Set("k4.b", 54);
        nodeIds.Set("k4.a", 55);
        nodeIds.Set("k3.a", 56);
        nodeIds.Set("red3.in", 57);
        nodeIds.Set("vr.zr", 58);
        nodeIds.Set("evr.vr_tc", 59);
        nodeIds.Set("ko2.out", 60);
        nodeIds.Set("evr.zr", 61);
        nodeIds.Set("tc_set", 62);
        nodeIds.Set("tc_3_9", 63);
        nodeIds.Set("vr.mk", 64);
        nodeIds.Set("rd1.feed", 65);
        nodeIds.Set("vt1_kn12", 66);
        nodeIds.Set("kn10_vt3", 67);
        nodeIds.Set("rs8", 68);
        nodeIds.Set("red3.out", 69);
        nodeIds.Set("evr.vr_zr", 70);
        nodeIds.Set("tc_2", 71);
        nodeIds.Set("vr.tc", 72);
        nodeIds.Set("k5.out", 73);
        nodeIds.Set("evr.rd_zr", 74);
        nodeIds.Set("vr.zk", 75);
        nodeIds.Set("rs3", 76);
        nodeIds.Set("evr.tc", 77);
        nodeIds.Set("k3.out", 78);
        nodeIds.Set("k4.out", 79);
        nodeIds.Set("k5.b", 80);
        nodeIds.Set("red1.out", 81);

        linkIds = new LIntMap(); linkIds.Reserve(67);
        linkIds.Set("k254_1.vr_b", 0);
        linkIds.Set("k395_1.d_at", -1);
        linkIds.Set("red1.feed", 1);
        linkIds.Set("epka_1.tm_at", -2);
        linkIds.Set("kp3.dump", -3);
        linkIds.Set("k254_1.b_at", -4);
        linkIds.Set("k395_1.gr_u_red", 2);
        linkIds.Set("k395_1.gr_u", 3);
        linkIds.Set("epka_1.k_at", -5);
        linkIds.Set("rd2.dump", -6);
        linkIds.Set("k254_1.tc_at", -7);
        linkIds.Set("ubt_1.tc", 4);
        linkIds.Set("ubt_1.gr", 5);
        linkIds.Set("kn27", 6);
        linkIds.Set("tm_epka_1", 7);
        linkIds.Set("kn18", 8);
        linkIds.Set("epka_1.gr_k", 9);
        linkIds.Set("k395_1.u_s", 10);
        linkIds.Set("k395_1.u_at", -8);
        linkIds.Set("k395_1.ur_d", 11);
        linkIds.Set("k395_1.ur_at", -9);
        linkIds.Set("k395_1.ok", 12);
        linkIds.Set("rd2.feed", 13);
        linkIds.Set("ko4.in_out", 14);
        linkIds.Set("ubt_1.tm", 15);
        linkIds.Set("ubt_1.tm_at", -10);
        linkIds.Set("k254_1.gr_tc", 16);
        linkIds.Set("k395_1.tm_gr", 17);
        linkIds.Set("k395_1.u_tm_gr", 18);
        linkIds.Set("k395_1.u_tm_at", -11);
        linkIds.Set("k395_1.tm_at", -12);
        linkIds.Set("ko3.in_out", 19);
        linkIds.Set("kn14", 20);
        linkIds.Set("pm_epka_1", 21);
        linkIds.Set("kn6", 22);
        linkIds.Set("evr.leak_rk", -13);
        linkIds.Set("evr.vo", -14);
        linkIds.Set("vr.kdr_at", -15);
        linkIds.Set("vr.mk_at_sk", -16);
        linkIds.Set("ko2.in_out", 23);
        linkIds.Set("evr.rd_rel", -17);
        linkIds.Set("vr.sk_tc", 24);
        linkIds.Set("rd1.dump", -18);
        linkIds.Set("k3.b_out", 25);
        linkIds.Set("k5.a_out", 26);
        linkIds.Set("k4.b_out", 27);
        linkIds.Set("k4.a_out", 28);
        linkIds.Set("k3.a_out", 29);
        linkIds.Set("red3.feed", 30);
        linkIds.Set("vr.mk_kdr", 31);
        linkIds.Set("vr.mk_sk", 32);
        linkIds.Set("vr.mk_zk", 33);
        linkIds.Set("rd1.feed", 34);
        linkIds.Set("kn12", 35);
        linkIds.Set("vt3", 36);
        linkIds.Set("vt2", 37);
        linkIds.Set("vr.tc_at", -19);
        linkIds.Set("evr.vt", 38);
        linkIds.Set("evr.rd_feed", 39);
        linkIds.Set("vr.zk_tc_e", 40);
        linkIds.Set("vr.zk_tc", 41);
        linkIds.Set("kn11", 42);
        linkIds.Set("kn10", 43);
        linkIds.Set("evr.pk_evr", 44);
        linkIds.Set("evr.pk_vr", 45);
        linkIds.Set("vt1", 46);
        linkIds.Set("k5.b_out", 47);

        objectIds = new LIntMap(); objectIds.Reserve(23);
        objectIds.Set("kp3", 0);
        objectIds.Set("epka_1", 1);
        objectIds.Set("ubt_1", 2);
        objectIds.Set("k395_1", 3);
        objectIds.Set("k254_1", 4);
        objectIds.Set("ko3", 5);
        objectIds.Set("red1", 6);
        objectIds.Set("ko2", 7);
        objectIds.Set("ko4", 8);
        objectIds.Set("vr", 9);
        objectIds.Set("evr", 10);
        objectIds.Set("k3", 11);
        objectIds.Set("k4", 12);
        objectIds.Set("k5", 13);
        objectIds.Set("red3", 14);
        objectIds.Set("rd2", 15);
        objectIds.Set("tc_4_10", 16);
        objectIds.Set("tc_5_11", 17);
        objectIds.Set("tc_6_12", 18);
        objectIds.Set("rd1", 19);
        objectIds.Set("tc_3_9", 20);
        objectIds.Set("tc_2_8", 21);
        objectIds.Set("tc_1_7", 22);

        linkNominalFlowRate = new float[48];
        linkSrc = new int[48];
        linkDst = new int[48];
        linkNominalFlowRate[0] = 9.027353; linkSrc[0] = 1; linkDst[0] = 11;
        linkNominalFlowRate[1] = 1; linkSrc[1] = 6; linkDst[1] = 81;
        linkNominalFlowRate[2] = 36.1094131; linkSrc[2] = 12; linkDst[2] = 28;
        linkNominalFlowRate[3] = 100.303925; linkSrc[3] = 12; linkDst[3] = 28;
        linkNominalFlowRate[4] = 2507.59814; linkSrc[4] = 20; linkDst[4] = 27;
        linkNominalFlowRate[5] = 2507.59814; linkSrc[5] = 21; linkDst[5] = 37;
        linkNominalFlowRate[6] = 2311.00244; linkSrc[6] = 23; linkDst[6] = 22;
        linkNominalFlowRate[7] = 2311.00244; linkSrc[7] = 24; linkDst[7] = 7;
        linkNominalFlowRate[8] = 1448.38867; linkSrc[8] = 24; linkDst[8] = 38;
        linkNominalFlowRate[9] = 3.24984717; linkSrc[9] = 25; linkDst[9] = 13;
        linkNominalFlowRate[10] = 36.1094131; linkSrc[10] = 28; linkDst[10] = 2;
        linkNominalFlowRate[11] = 36.1094131; linkSrc[11] = 29; linkDst[11] = 3;
        linkNominalFlowRate[12] = 36.1094131; linkSrc[12] = 29; linkDst[12] = 35;
        linkNominalFlowRate[13] = 1; linkSrc[13] = 30; linkDst[13] = 15;
        linkNominalFlowRate[14] = 1; linkSrc[14] = 31; linkDst[14] = 16;
        linkNominalFlowRate[15] = 2507.59814; linkSrc[15] = 32; linkDst[15] = 19;
        linkNominalFlowRate[16] = 401.2157; linkSrc[16] = 33; linkDst[16] = 18;
        linkNominalFlowRate[17] = 1027.11218; linkSrc[17] = 35; linkDst[17] = 12;
        linkNominalFlowRate[18] = 324.9847; linkSrc[18] = 35; linkDst[18] = 12;
        linkNominalFlowRate[19] = 1; linkSrc[19] = 38; linkDst[19] = 17;
        linkNominalFlowRate[20] = 1448.38867; linkSrc[20] = 39; linkDst[20] = 6;
        linkNominalFlowRate[21] = 678.054565; linkSrc[21] = 39; linkDst[21] = 25;
        linkNominalFlowRate[22] = 1448.38867; linkSrc[22] = 39; linkDst[22] = 31;
        linkNominalFlowRate[23] = 1; linkSrc[23] = 46; linkDst[23] = 60;
        linkNominalFlowRate[24] = 144.437653; linkSrc[24] = 49; linkDst[24] = 72;
        linkNominalFlowRate[25] = 1; linkSrc[25] = 52; linkDst[25] = 78;
        linkNominalFlowRate[26] = 1; linkSrc[26] = 53; linkDst[26] = 73;
        linkNominalFlowRate[27] = 1; linkSrc[27] = 54; linkDst[27] = 79;
        linkNominalFlowRate[28] = 1; linkSrc[28] = 55; linkDst[28] = 79;
        linkNominalFlowRate[29] = 1; linkSrc[29] = 56; linkDst[29] = 78;
        linkNominalFlowRate[30] = 1; linkSrc[30] = 57; linkDst[30] = 69;
        linkNominalFlowRate[31] = 49.148922; linkSrc[31] = 64; linkDst[31] = 41;
        linkNominalFlowRate[32] = 2.56778049; linkSrc[32] = 64; linkDst[32] = 49;
        linkNominalFlowRate[33] = 1; linkSrc[33] = 64; linkDst[33] = 75;
        linkNominalFlowRate[34] = 1; linkSrc[34] = 65; linkDst[34] = 51;
        linkNominalFlowRate[35] = 159.242508; linkSrc[35] = 66; linkDst[35] = 53;
        linkNominalFlowRate[36] = 159.242508; linkSrc[36] = 67; linkDst[36] = 54;
        linkNominalFlowRate[37] = 159.242508; linkSrc[37] = 69; linkDst[37] = 56;
        linkNominalFlowRate[38] = 12.9993887; linkSrc[38] = 74; linkDst[38] = 40;
        linkNominalFlowRate[39] = 324.9847; linkSrc[39] = 74; linkDst[39] = 48;
        linkNominalFlowRate[40] = 1; linkSrc[40] = 75; linkDst[40] = 72;
        linkNominalFlowRate[41] = 64.19451; linkSrc[41] = 75; linkDst[41] = 72;
        linkNominalFlowRate[42] = 159.242508; linkSrc[42] = 76; linkDst[42] = 57;
        linkNominalFlowRate[43] = 159.242508; linkSrc[43] = 76; linkDst[43] = 67;
        linkNominalFlowRate[44] = 1448.38867; linkSrc[44] = 77; linkDst[44] = 48;
        linkNominalFlowRate[45] = 1448.38867; linkSrc[45] = 77; linkDst[45] = 59;
        linkNominalFlowRate[46] = 159.242508; linkSrc[46] = 79; linkDst[46] = 66;
        linkNominalFlowRate[47] = 1; linkSrc[47] = 80; linkDst[47] = 73;

        alinkNominalFlowRate = new float[19];
        alinkSrc = new int[19];
        alinkNominalFlowRate[0] = 36.1094131; alinkSrc[0] = 3;
        alinkNominalFlowRate[1] = 3610.94141; alinkSrc[1] = 7;
        alinkNominalFlowRate[2] = 1; alinkSrc[2] = 8;
        alinkNominalFlowRate[3] = 100.303925; alinkSrc[3] = 11;
        alinkNominalFlowRate[4] = 4.012157; alinkSrc[4] = 13;
        alinkNominalFlowRate[5] = 1; alinkSrc[5] = 15;
        alinkNominalFlowRate[6] = 88.62855; alinkSrc[6] = 18;
        alinkNominalFlowRate[7] = 100.303925; alinkSrc[7] = 28;
        alinkNominalFlowRate[8] = 1; alinkSrc[8] = 29;
        alinkNominalFlowRate[9] = 256.778046; alinkSrc[9] = 32;
        alinkNominalFlowRate[10] = 401.2157; alinkSrc[10] = 35;
        alinkNominalFlowRate[11] = 1027.11218; alinkSrc[11] = 35;
        alinkNominalFlowRate[12] = 1; alinkSrc[12] = 40;
        alinkNominalFlowRate[13] = 6.780545; alinkSrc[13] = 40;
        alinkNominalFlowRate[14] = 4.012157; alinkSrc[14] = 41;
        alinkNominalFlowRate[15] = 1027.11218; alinkSrc[15] = 42;
        alinkNominalFlowRate[16] = 401.2157; alinkSrc[16] = 48;
        alinkNominalFlowRate[17] = 1; alinkSrc[17] = 51;
        alinkNominalFlowRate[18] = 1; alinkSrc[18] = 72;

    }

    public Soup GetSchemePropertiesInternal(tep70airscheme_instance scheme)
    {
        Soup pressureSoup = Constructors.NewSoup();
        float[] nodePressure = scheme.nodePressure;
        pressureSoup.SetNamedTag("k254_1.r", nodePressure[0]);
        pressureSoup.SetNamedTag("k254_1.vr", nodePressure[1]);
        pressureSoup.SetNamedTag("k395_1.s", nodePressure[2]);
        pressureSoup.SetNamedTag("k395_1.d", nodePressure[3]);
        pressureSoup.SetNamedTag("epka_1.kvv", nodePressure[4]);
        pressureSoup.SetNamedTag("tc_6_12", nodePressure[5]);
        pressureSoup.SetNamedTag("red1.in", nodePressure[6]);
        pressureSoup.SetNamedTag("epka_1.tm", nodePressure[7]);
        pressureSoup.SetNamedTag("kp3", nodePressure[8]);
        pressureSoup.SetNamedTag("ur", nodePressure[9]);
        pressureSoup.SetNamedTag("mvt", nodePressure[10]);
        pressureSoup.SetNamedTag("k254_1.b", nodePressure[11]);
        pressureSoup.SetNamedTag("k395_1.gr", nodePressure[12]);
        pressureSoup.SetNamedTag("epka_1.k", nodePressure[13]);
        pressureSoup.SetNamedTag("tc_5_11", nodePressure[14]);
        pressureSoup.SetNamedTag("rd2.out", nodePressure[15]);
        pressureSoup.SetNamedTag("ko4.out", nodePressure[16]);
        pressureSoup.SetNamedTag("ko3.out", nodePressure[17]);
        pressureSoup.SetNamedTag("k254_1.tc", nodePressure[18]);
        pressureSoup.SetNamedTag("ubt_1.k_tm", nodePressure[19]);
        pressureSoup.SetNamedTag("ubt_1.tc", nodePressure[20]);
        pressureSoup.SetNamedTag("ubt_1.gr", nodePressure[21]);
        pressureSoup.SetNamedTag("rs2", nodePressure[22]);
        pressureSoup.SetNamedTag("rs1", nodePressure[23]);
        pressureSoup.SetNamedTag("tm", nodePressure[24]);
        pressureSoup.SetNamedTag("epka_1.gr", nodePressure[25]);
        pressureSoup.SetNamedTag("tc_4_10", nodePressure[26]);
        pressureSoup.SetNamedTag("ubt_1.k_tc", nodePressure[27]);
        pressureSoup.SetNamedTag("k395_1.u", nodePressure[28]);
        pressureSoup.SetNamedTag("k395_1.ur", nodePressure[29]);
        pressureSoup.SetNamedTag("rd2.feed", nodePressure[30]);
        pressureSoup.SetNamedTag("ko4.in", nodePressure[31]);
        pressureSoup.SetNamedTag("ubt_1.tm", nodePressure[32]);
        pressureSoup.SetNamedTag("k254_1.gr", nodePressure[33]);
        pressureSoup.SetNamedTag("tc_1", nodePressure[34]);
        pressureSoup.SetNamedTag("k395_1.tm", nodePressure[35]);
        pressureSoup.SetNamedTag("rs4", nodePressure[36]);
        pressureSoup.SetNamedTag("ubt_1.k_gr", nodePressure[37]);
        pressureSoup.SetNamedTag("ko3.in", nodePressure[38]);
        pressureSoup.SetNamedTag("pm", nodePressure[39]);
        pressureSoup.SetNamedTag("evr.rk", nodePressure[40]);
        pressureSoup.SetNamedTag("vr.kdr", nodePressure[41]);
        pressureSoup.SetNamedTag("vr.tm", nodePressure[42]);
        pressureSoup.SetNamedTag("tc_1_7", nodePressure[43]);
        pressureSoup.SetNamedTag("rd1.set", nodePressure[44]);
        pressureSoup.SetNamedTag("rd2.set", nodePressure[45]);
        pressureSoup.SetNamedTag("ko2.in", nodePressure[46]);
        pressureSoup.SetNamedTag("rs5", nodePressure[47]);
        pressureSoup.SetNamedTag("evr.rd_tc", nodePressure[48]);
        pressureSoup.SetNamedTag("vr.sk", nodePressure[49]);
        pressureSoup.SetNamedTag("tc_2_8", nodePressure[50]);
        pressureSoup.SetNamedTag("rd1.out", nodePressure[51]);
        pressureSoup.SetNamedTag("k3.b", nodePressure[52]);
        pressureSoup.SetNamedTag("k5.a", nodePressure[53]);
        pressureSoup.SetNamedTag("k4.b", nodePressure[54]);
        pressureSoup.SetNamedTag("k4.a", nodePressure[55]);
        pressureSoup.SetNamedTag("k3.a", nodePressure[56]);
        pressureSoup.SetNamedTag("red3.in", nodePressure[57]);
        pressureSoup.SetNamedTag("vr.zr", nodePressure[58]);
        pressureSoup.SetNamedTag("evr.vr_tc", nodePressure[59]);
        pressureSoup.SetNamedTag("ko2.out", nodePressure[60]);
        pressureSoup.SetNamedTag("evr.zr", nodePressure[61]);
        pressureSoup.SetNamedTag("tc_set", nodePressure[62]);
        pressureSoup.SetNamedTag("tc_3_9", nodePressure[63]);
        pressureSoup.SetNamedTag("vr.mk", nodePressure[64]);
        pressureSoup.SetNamedTag("rd1.feed", nodePressure[65]);
        pressureSoup.SetNamedTag("vt1_kn12", nodePressure[66]);
        pressureSoup.SetNamedTag("kn10_vt3", nodePressure[67]);
        pressureSoup.SetNamedTag("rs8", nodePressure[68]);
        pressureSoup.SetNamedTag("red3.out", nodePressure[69]);
        pressureSoup.SetNamedTag("evr.vr_zr", nodePressure[70]);
        pressureSoup.SetNamedTag("tc_2", nodePressure[71]);
        pressureSoup.SetNamedTag("vr.tc", nodePressure[72]);
        pressureSoup.SetNamedTag("k5.out", nodePressure[73]);
        pressureSoup.SetNamedTag("evr.rd_zr", nodePressure[74]);
        pressureSoup.SetNamedTag("vr.zk", nodePressure[75]);
        pressureSoup.SetNamedTag("rs3", nodePressure[76]);
        pressureSoup.SetNamedTag("evr.tc", nodePressure[77]);
        pressureSoup.SetNamedTag("k3.out", nodePressure[78]);
        pressureSoup.SetNamedTag("k4.out", nodePressure[79]);
        pressureSoup.SetNamedTag("k5.b", nodePressure[80]);
        pressureSoup.SetNamedTag("red1.out", nodePressure[81]);

        Soup volumeSoup = Constructors.NewSoup();
        float[] nodeVolume = scheme.nodeVolume;
        volumeSoup.SetNamedTag("k254_1.r", nodeVolume[0]);
        volumeSoup.SetNamedTag("k254_1.vr", nodeVolume[1]);
        volumeSoup.SetNamedTag("k395_1.s", nodeVolume[2]);
        volumeSoup.SetNamedTag("k395_1.d", nodeVolume[3]);
        volumeSoup.SetNamedTag("epka_1.kvv", nodeVolume[4]);
        volumeSoup.SetNamedTag("tc_6_12", nodeVolume[5]);
        volumeSoup.SetNamedTag("red1.in", nodeVolume[6]);
        volumeSoup.SetNamedTag("epka_1.tm", nodeVolume[7]);
        volumeSoup.SetNamedTag("kp3", nodeVolume[8]);
        volumeSoup.SetNamedTag("ur", nodeVolume[9]);
        volumeSoup.SetNamedTag("mvt", nodeVolume[10]);
        volumeSoup.SetNamedTag("k254_1.b", nodeVolume[11]);
        volumeSoup.SetNamedTag("k395_1.gr", nodeVolume[12]);
        volumeSoup.SetNamedTag("epka_1.k", nodeVolume[13]);
        volumeSoup.SetNamedTag("tc_5_11", nodeVolume[14]);
        volumeSoup.SetNamedTag("rd2.out", nodeVolume[15]);
        volumeSoup.SetNamedTag("ko4.out", nodeVolume[16]);
        volumeSoup.SetNamedTag("ko3.out", nodeVolume[17]);
        volumeSoup.SetNamedTag("k254_1.tc", nodeVolume[18]);
        volumeSoup.SetNamedTag("ubt_1.k_tm", nodeVolume[19]);
        volumeSoup.SetNamedTag("ubt_1.tc", nodeVolume[20]);
        volumeSoup.SetNamedTag("ubt_1.gr", nodeVolume[21]);
        volumeSoup.SetNamedTag("rs2", nodeVolume[22]);
        volumeSoup.SetNamedTag("rs1", nodeVolume[23]);
        volumeSoup.SetNamedTag("tm", nodeVolume[24]);
        volumeSoup.SetNamedTag("epka_1.gr", nodeVolume[25]);
        volumeSoup.SetNamedTag("tc_4_10", nodeVolume[26]);
        volumeSoup.SetNamedTag("ubt_1.k_tc", nodeVolume[27]);
        volumeSoup.SetNamedTag("k395_1.u", nodeVolume[28]);
        volumeSoup.SetNamedTag("k395_1.ur", nodeVolume[29]);
        volumeSoup.SetNamedTag("rd2.feed", nodeVolume[30]);
        volumeSoup.SetNamedTag("ko4.in", nodeVolume[31]);
        volumeSoup.SetNamedTag("ubt_1.tm", nodeVolume[32]);
        volumeSoup.SetNamedTag("k254_1.gr", nodeVolume[33]);
        volumeSoup.SetNamedTag("tc_1", nodeVolume[34]);
        volumeSoup.SetNamedTag("k395_1.tm", nodeVolume[35]);
        volumeSoup.SetNamedTag("rs4", nodeVolume[36]);
        volumeSoup.SetNamedTag("ubt_1.k_gr", nodeVolume[37]);
        volumeSoup.SetNamedTag("ko3.in", nodeVolume[38]);
        volumeSoup.SetNamedTag("pm", nodeVolume[39]);
        volumeSoup.SetNamedTag("evr.rk", nodeVolume[40]);
        volumeSoup.SetNamedTag("vr.kdr", nodeVolume[41]);
        volumeSoup.SetNamedTag("vr.tm", nodeVolume[42]);
        volumeSoup.SetNamedTag("tc_1_7", nodeVolume[43]);
        volumeSoup.SetNamedTag("rd1.set", nodeVolume[44]);
        volumeSoup.SetNamedTag("rd2.set", nodeVolume[45]);
        volumeSoup.SetNamedTag("ko2.in", nodeVolume[46]);
        volumeSoup.SetNamedTag("rs5", nodeVolume[47]);
        volumeSoup.SetNamedTag("evr.rd_tc", nodeVolume[48]);
        volumeSoup.SetNamedTag("vr.sk", nodeVolume[49]);
        volumeSoup.SetNamedTag("tc_2_8", nodeVolume[50]);
        volumeSoup.SetNamedTag("rd1.out", nodeVolume[51]);
        volumeSoup.SetNamedTag("k3.b", nodeVolume[52]);
        volumeSoup.SetNamedTag("k5.a", nodeVolume[53]);
        volumeSoup.SetNamedTag("k4.b", nodeVolume[54]);
        volumeSoup.SetNamedTag("k4.a", nodeVolume[55]);
        volumeSoup.SetNamedTag("k3.a", nodeVolume[56]);
        volumeSoup.SetNamedTag("red3.in", nodeVolume[57]);
        volumeSoup.SetNamedTag("vr.zr", nodeVolume[58]);
        volumeSoup.SetNamedTag("evr.vr_tc", nodeVolume[59]);
        volumeSoup.SetNamedTag("ko2.out", nodeVolume[60]);
        volumeSoup.SetNamedTag("evr.zr", nodeVolume[61]);
        volumeSoup.SetNamedTag("tc_set", nodeVolume[62]);
        volumeSoup.SetNamedTag("tc_3_9", nodeVolume[63]);
        volumeSoup.SetNamedTag("vr.mk", nodeVolume[64]);
        volumeSoup.SetNamedTag("rd1.feed", nodeVolume[65]);
        volumeSoup.SetNamedTag("vt1_kn12", nodeVolume[66]);
        volumeSoup.SetNamedTag("kn10_vt3", nodeVolume[67]);
        volumeSoup.SetNamedTag("rs8", nodeVolume[68]);
        volumeSoup.SetNamedTag("red3.out", nodeVolume[69]);
        volumeSoup.SetNamedTag("evr.vr_zr", nodeVolume[70]);
        volumeSoup.SetNamedTag("tc_2", nodeVolume[71]);
        volumeSoup.SetNamedTag("vr.tc", nodeVolume[72]);
        volumeSoup.SetNamedTag("k5.out", nodeVolume[73]);
        volumeSoup.SetNamedTag("evr.rd_zr", nodeVolume[74]);
        volumeSoup.SetNamedTag("vr.zk", nodeVolume[75]);
        volumeSoup.SetNamedTag("rs3", nodeVolume[76]);
        volumeSoup.SetNamedTag("evr.tc", nodeVolume[77]);
        volumeSoup.SetNamedTag("k3.out", nodeVolume[78]);
        volumeSoup.SetNamedTag("k4.out", nodeVolume[79]);
        volumeSoup.SetNamedTag("k5.b", nodeVolume[80]);
        volumeSoup.SetNamedTag("red1.out", nodeVolume[81]);

        Soup flowRateSoup = Constructors.NewSoup();
        float[] linkFlowRate = scheme.linkFlowRate;
        float[] alinkFlowRate = scheme.alinkFlowRate;
        flowRateSoup.SetNamedTag("k254_1.vr_b", linkFlowRate[0]);
        flowRateSoup.SetNamedTag("k395_1.d_at", alinkFlowRate[0]);
        flowRateSoup.SetNamedTag("red1.feed", linkFlowRate[1]);
        flowRateSoup.SetNamedTag("epka_1.tm_at", alinkFlowRate[1]);
        flowRateSoup.SetNamedTag("kp3.dump", alinkFlowRate[2]);
        flowRateSoup.SetNamedTag("k254_1.b_at", alinkFlowRate[3]);
        flowRateSoup.SetNamedTag("k395_1.gr_u_red", linkFlowRate[2]);
        flowRateSoup.SetNamedTag("k395_1.gr_u", linkFlowRate[3]);
        flowRateSoup.SetNamedTag("epka_1.k_at", alinkFlowRate[4]);
        flowRateSoup.SetNamedTag("rd2.dump", alinkFlowRate[5]);
        flowRateSoup.SetNamedTag("k254_1.tc_at", alinkFlowRate[6]);
        flowRateSoup.SetNamedTag("ubt_1.tc", linkFlowRate[4]);
        flowRateSoup.SetNamedTag("ubt_1.gr", linkFlowRate[5]);
        flowRateSoup.SetNamedTag("kn27", linkFlowRate[6]);
        flowRateSoup.SetNamedTag("tm_epka_1", linkFlowRate[7]);
        flowRateSoup.SetNamedTag("kn18", linkFlowRate[8]);
        flowRateSoup.SetNamedTag("epka_1.gr_k", linkFlowRate[9]);
        flowRateSoup.SetNamedTag("k395_1.u_s", linkFlowRate[10]);
        flowRateSoup.SetNamedTag("k395_1.u_at", alinkFlowRate[7]);
        flowRateSoup.SetNamedTag("k395_1.ur_d", linkFlowRate[11]);
        flowRateSoup.SetNamedTag("k395_1.ur_at", alinkFlowRate[8]);
        flowRateSoup.SetNamedTag("k395_1.ok", linkFlowRate[12]);
        flowRateSoup.SetNamedTag("rd2.feed", linkFlowRate[13]);
        flowRateSoup.SetNamedTag("ko4.in_out", linkFlowRate[14]);
        flowRateSoup.SetNamedTag("ubt_1.tm", linkFlowRate[15]);
        flowRateSoup.SetNamedTag("ubt_1.tm_at", alinkFlowRate[9]);
        flowRateSoup.SetNamedTag("k254_1.gr_tc", linkFlowRate[16]);
        flowRateSoup.SetNamedTag("k395_1.tm_gr", linkFlowRate[17]);
        flowRateSoup.SetNamedTag("k395_1.u_tm_gr", linkFlowRate[18]);
        flowRateSoup.SetNamedTag("k395_1.u_tm_at", alinkFlowRate[10]);
        flowRateSoup.SetNamedTag("k395_1.tm_at", alinkFlowRate[11]);
        flowRateSoup.SetNamedTag("ko3.in_out", linkFlowRate[19]);
        flowRateSoup.SetNamedTag("kn14", linkFlowRate[20]);
        flowRateSoup.SetNamedTag("pm_epka_1", linkFlowRate[21]);
        flowRateSoup.SetNamedTag("kn6", linkFlowRate[22]);
        flowRateSoup.SetNamedTag("evr.leak_rk", alinkFlowRate[12]);
        flowRateSoup.SetNamedTag("evr.vo", alinkFlowRate[13]);
        flowRateSoup.SetNamedTag("vr.kdr_at", alinkFlowRate[14]);
        flowRateSoup.SetNamedTag("vr.mk_at_sk", alinkFlowRate[15]);
        flowRateSoup.SetNamedTag("ko2.in_out", linkFlowRate[23]);
        flowRateSoup.SetNamedTag("evr.rd_rel", alinkFlowRate[16]);
        flowRateSoup.SetNamedTag("vr.sk_tc", linkFlowRate[24]);
        flowRateSoup.SetNamedTag("rd1.dump", alinkFlowRate[17]);
        flowRateSoup.SetNamedTag("k3.b_out", linkFlowRate[25]);
        flowRateSoup.SetNamedTag("k5.a_out", linkFlowRate[26]);
        flowRateSoup.SetNamedTag("k4.b_out", linkFlowRate[27]);
        flowRateSoup.SetNamedTag("k4.a_out", linkFlowRate[28]);
        flowRateSoup.SetNamedTag("k3.a_out", linkFlowRate[29]);
        flowRateSoup.SetNamedTag("red3.feed", linkFlowRate[30]);
        flowRateSoup.SetNamedTag("vr.mk_kdr", linkFlowRate[31]);
        flowRateSoup.SetNamedTag("vr.mk_sk", linkFlowRate[32]);
        flowRateSoup.SetNamedTag("vr.mk_zk", linkFlowRate[33]);
        flowRateSoup.SetNamedTag("rd1.feed", linkFlowRate[34]);
        flowRateSoup.SetNamedTag("kn12", linkFlowRate[35]);
        flowRateSoup.SetNamedTag("vt3", linkFlowRate[36]);
        flowRateSoup.SetNamedTag("vt2", linkFlowRate[37]);
        flowRateSoup.SetNamedTag("vr.tc_at", alinkFlowRate[18]);
        flowRateSoup.SetNamedTag("evr.vt", linkFlowRate[38]);
        flowRateSoup.SetNamedTag("evr.rd_feed", linkFlowRate[39]);
        flowRateSoup.SetNamedTag("vr.zk_tc_e", linkFlowRate[40]);
        flowRateSoup.SetNamedTag("vr.zk_tc", linkFlowRate[41]);
        flowRateSoup.SetNamedTag("kn11", linkFlowRate[42]);
        flowRateSoup.SetNamedTag("kn10", linkFlowRate[43]);
        flowRateSoup.SetNamedTag("evr.pk_evr", linkFlowRate[44]);
        flowRateSoup.SetNamedTag("evr.pk_vr", linkFlowRate[45]);
        flowRateSoup.SetNamedTag("vt1", linkFlowRate[46]);
        flowRateSoup.SetNamedTag("k5.b_out", linkFlowRate[47]);

        Soup objectsSoup = Constructors.NewSoup();
        SetCurrentObjectName("kp3");
        objectsSoup.SetNamedSoup("kp3", scheme.objects[0].GetProperties());
        SetCurrentObjectName("epka_1");
        objectsSoup.SetNamedSoup("epka_1", scheme.objects[1].GetProperties());
        SetCurrentObjectName("ubt_1");
        objectsSoup.SetNamedSoup("ubt_1", scheme.objects[2].GetProperties());
        SetCurrentObjectName("k395_1");
        objectsSoup.SetNamedSoup("k395_1", scheme.objects[3].GetProperties());
        SetCurrentObjectName("k254_1");
        objectsSoup.SetNamedSoup("k254_1", scheme.objects[4].GetProperties());
        SetCurrentObjectName("ko3");
        objectsSoup.SetNamedSoup("ko3", scheme.objects[5].GetProperties());
        SetCurrentObjectName("red1");
        objectsSoup.SetNamedSoup("red1", scheme.objects[6].GetProperties());
        SetCurrentObjectName("ko2");
        objectsSoup.SetNamedSoup("ko2", scheme.objects[7].GetProperties());
        SetCurrentObjectName("ko4");
        objectsSoup.SetNamedSoup("ko4", scheme.objects[8].GetProperties());
        SetCurrentObjectName("vr");
        objectsSoup.SetNamedSoup("vr", scheme.objects[9].GetProperties());
        SetCurrentObjectName("evr");
        objectsSoup.SetNamedSoup("evr", scheme.objects[10].GetProperties());
        SetCurrentObjectName("k3");
        objectsSoup.SetNamedSoup("k3", scheme.objects[11].GetProperties());
        SetCurrentObjectName("k4");
        objectsSoup.SetNamedSoup("k4", scheme.objects[12].GetProperties());
        SetCurrentObjectName("k5");
        objectsSoup.SetNamedSoup("k5", scheme.objects[13].GetProperties());
        SetCurrentObjectName("red3");
        objectsSoup.SetNamedSoup("red3", scheme.objects[14].GetProperties());
        SetCurrentObjectName("rd2");
        objectsSoup.SetNamedSoup("rd2", scheme.objects[15].GetProperties());
        SetCurrentObjectName("tc_4_10");
        objectsSoup.SetNamedSoup("tc_4_10", scheme.objects[16].GetProperties());
        SetCurrentObjectName("tc_5_11");
        objectsSoup.SetNamedSoup("tc_5_11", scheme.objects[17].GetProperties());
        SetCurrentObjectName("tc_6_12");
        objectsSoup.SetNamedSoup("tc_6_12", scheme.objects[18].GetProperties());
        SetCurrentObjectName("rd1");
        objectsSoup.SetNamedSoup("rd1", scheme.objects[19].GetProperties());
        SetCurrentObjectName("tc_3_9");
        objectsSoup.SetNamedSoup("tc_3_9", scheme.objects[20].GetProperties());
        SetCurrentObjectName("tc_2_8");
        objectsSoup.SetNamedSoup("tc_2_8", scheme.objects[21].GetProperties());
        SetCurrentObjectName("tc_1_7");
        objectsSoup.SetNamedSoup("tc_1_7", scheme.objects[22].GetProperties());
        ResetCurrentObjectName();
        Soup ret = Constructors.NewSoup();
        ret.SetNamedSoup("pressure", pressureSoup);
        ret.SetNamedSoup("volume", volumeSoup);
        ret.SetNamedSoup("flowrate", flowRateSoup);
        ret.SetNamedSoup("objects", objectsSoup);
        return ret;
    }

    public void SetSchemePropertiesInternal(tep70airscheme_instance scheme, Soup sp)
    {
        Soup pressureSoup = sp.GetNamedSoup("pressure");
        float[] p = scheme.nodePressure;
        p[0] = SPHelper.GetFloat(pressureSoup, "k254_1.r", p[0]);
        p[1] = SPHelper.GetFloat(pressureSoup, "k254_1.vr", p[1]);
        p[2] = SPHelper.GetFloat(pressureSoup, "k395_1.s", p[2]);
        p[3] = SPHelper.GetFloat(pressureSoup, "k395_1.d", p[3]);
        p[4] = SPHelper.GetFloat(pressureSoup, "epka_1.kvv", p[4]);
        p[5] = SPHelper.GetFloat(pressureSoup, "tc_6_12", p[5]);
        p[6] = SPHelper.GetFloat(pressureSoup, "red1.in", p[6]);
        p[7] = SPHelper.GetFloat(pressureSoup, "epka_1.tm", p[7]);
        p[8] = SPHelper.GetFloat(pressureSoup, "kp3", p[8]);
        p[9] = SPHelper.GetFloat(pressureSoup, "ur", p[9]);
        p[10] = SPHelper.GetFloat(pressureSoup, "mvt", p[10]);
        p[11] = SPHelper.GetFloat(pressureSoup, "k254_1.b", p[11]);
        p[12] = SPHelper.GetFloat(pressureSoup, "k395_1.gr", p[12]);
        p[13] = SPHelper.GetFloat(pressureSoup, "epka_1.k", p[13]);
        p[14] = SPHelper.GetFloat(pressureSoup, "tc_5_11", p[14]);
        p[15] = SPHelper.GetFloat(pressureSoup, "rd2.out", p[15]);
        p[16] = SPHelper.GetFloat(pressureSoup, "ko4.out", p[16]);
        p[17] = SPHelper.GetFloat(pressureSoup, "ko3.out", p[17]);
        p[18] = SPHelper.GetFloat(pressureSoup, "k254_1.tc", p[18]);
        p[19] = SPHelper.GetFloat(pressureSoup, "ubt_1.k_tm", p[19]);
        p[20] = SPHelper.GetFloat(pressureSoup, "ubt_1.tc", p[20]);
        p[21] = SPHelper.GetFloat(pressureSoup, "ubt_1.gr", p[21]);
        p[22] = SPHelper.GetFloat(pressureSoup, "rs2", p[22]);
        p[23] = SPHelper.GetFloat(pressureSoup, "rs1", p[23]);
        p[24] = SPHelper.GetFloat(pressureSoup, "tm", p[24]);
        p[25] = SPHelper.GetFloat(pressureSoup, "epka_1.gr", p[25]);
        p[26] = SPHelper.GetFloat(pressureSoup, "tc_4_10", p[26]);
        p[27] = SPHelper.GetFloat(pressureSoup, "ubt_1.k_tc", p[27]);
        p[28] = SPHelper.GetFloat(pressureSoup, "k395_1.u", p[28]);
        p[29] = SPHelper.GetFloat(pressureSoup, "k395_1.ur", p[29]);
        p[30] = SPHelper.GetFloat(pressureSoup, "rd2.feed", p[30]);
        p[31] = SPHelper.GetFloat(pressureSoup, "ko4.in", p[31]);
        p[32] = SPHelper.GetFloat(pressureSoup, "ubt_1.tm", p[32]);
        p[33] = SPHelper.GetFloat(pressureSoup, "k254_1.gr", p[33]);
        p[34] = SPHelper.GetFloat(pressureSoup, "tc_1", p[34]);
        p[35] = SPHelper.GetFloat(pressureSoup, "k395_1.tm", p[35]);
        p[36] = SPHelper.GetFloat(pressureSoup, "rs4", p[36]);
        p[37] = SPHelper.GetFloat(pressureSoup, "ubt_1.k_gr", p[37]);
        p[38] = SPHelper.GetFloat(pressureSoup, "ko3.in", p[38]);
        p[39] = SPHelper.GetFloat(pressureSoup, "pm", p[39]);
        p[40] = SPHelper.GetFloat(pressureSoup, "evr.rk", p[40]);
        p[41] = SPHelper.GetFloat(pressureSoup, "vr.kdr", p[41]);
        p[42] = SPHelper.GetFloat(pressureSoup, "vr.tm", p[42]);
        p[43] = SPHelper.GetFloat(pressureSoup, "tc_1_7", p[43]);
        p[44] = SPHelper.GetFloat(pressureSoup, "rd1.set", p[44]);
        p[45] = SPHelper.GetFloat(pressureSoup, "rd2.set", p[45]);
        p[46] = SPHelper.GetFloat(pressureSoup, "ko2.in", p[46]);
        p[47] = SPHelper.GetFloat(pressureSoup, "rs5", p[47]);
        p[48] = SPHelper.GetFloat(pressureSoup, "evr.rd_tc", p[48]);
        p[49] = SPHelper.GetFloat(pressureSoup, "vr.sk", p[49]);
        p[50] = SPHelper.GetFloat(pressureSoup, "tc_2_8", p[50]);
        p[51] = SPHelper.GetFloat(pressureSoup, "rd1.out", p[51]);
        p[52] = SPHelper.GetFloat(pressureSoup, "k3.b", p[52]);
        p[53] = SPHelper.GetFloat(pressureSoup, "k5.a", p[53]);
        p[54] = SPHelper.GetFloat(pressureSoup, "k4.b", p[54]);
        p[55] = SPHelper.GetFloat(pressureSoup, "k4.a", p[55]);
        p[56] = SPHelper.GetFloat(pressureSoup, "k3.a", p[56]);
        p[57] = SPHelper.GetFloat(pressureSoup, "red3.in", p[57]);
        p[58] = SPHelper.GetFloat(pressureSoup, "vr.zr", p[58]);
        p[59] = SPHelper.GetFloat(pressureSoup, "evr.vr_tc", p[59]);
        p[60] = SPHelper.GetFloat(pressureSoup, "ko2.out", p[60]);
        p[61] = SPHelper.GetFloat(pressureSoup, "evr.zr", p[61]);
        p[62] = SPHelper.GetFloat(pressureSoup, "tc_set", p[62]);
        p[63] = SPHelper.GetFloat(pressureSoup, "tc_3_9", p[63]);
        p[64] = SPHelper.GetFloat(pressureSoup, "vr.mk", p[64]);
        p[65] = SPHelper.GetFloat(pressureSoup, "rd1.feed", p[65]);
        p[66] = SPHelper.GetFloat(pressureSoup, "vt1_kn12", p[66]);
        p[67] = SPHelper.GetFloat(pressureSoup, "kn10_vt3", p[67]);
        p[68] = SPHelper.GetFloat(pressureSoup, "rs8", p[68]);
        p[69] = SPHelper.GetFloat(pressureSoup, "red3.out", p[69]);
        p[70] = SPHelper.GetFloat(pressureSoup, "evr.vr_zr", p[70]);
        p[71] = SPHelper.GetFloat(pressureSoup, "tc_2", p[71]);
        p[72] = SPHelper.GetFloat(pressureSoup, "vr.tc", p[72]);
        p[73] = SPHelper.GetFloat(pressureSoup, "k5.out", p[73]);
        p[74] = SPHelper.GetFloat(pressureSoup, "evr.rd_zr", p[74]);
        p[75] = SPHelper.GetFloat(pressureSoup, "vr.zk", p[75]);
        p[76] = SPHelper.GetFloat(pressureSoup, "rs3", p[76]);
        p[77] = SPHelper.GetFloat(pressureSoup, "evr.tc", p[77]);
        p[78] = SPHelper.GetFloat(pressureSoup, "k3.out", p[78]);
        p[79] = SPHelper.GetFloat(pressureSoup, "k4.out", p[79]);
        p[80] = SPHelper.GetFloat(pressureSoup, "k5.b", p[80]);
        p[81] = SPHelper.GetFloat(pressureSoup, "red1.out", p[81]);

        Soup volumeSoup = sp.GetNamedSoup("volume");
        float[] v = scheme.nodeVolume;
        v[0] = SPHelper.GetFloat(volumeSoup, "k254_1.r", v[0]);
        v[1] = SPHelper.GetFloat(volumeSoup, "k254_1.vr", v[1]);
        v[2] = SPHelper.GetFloat(volumeSoup, "k395_1.s", v[2]);
        v[3] = SPHelper.GetFloat(volumeSoup, "k395_1.d", v[3]);
        v[4] = SPHelper.GetFloat(volumeSoup, "epka_1.kvv", v[4]);
        v[5] = SPHelper.GetFloat(volumeSoup, "tc_6_12", v[5]);
        v[6] = SPHelper.GetFloat(volumeSoup, "red1.in", v[6]);
        v[7] = SPHelper.GetFloat(volumeSoup, "epka_1.tm", v[7]);
        v[8] = SPHelper.GetFloat(volumeSoup, "kp3", v[8]);
        v[9] = SPHelper.GetFloat(volumeSoup, "ur", v[9]);
        v[10] = SPHelper.GetFloat(volumeSoup, "mvt", v[10]);
        v[11] = SPHelper.GetFloat(volumeSoup, "k254_1.b", v[11]);
        v[12] = SPHelper.GetFloat(volumeSoup, "k395_1.gr", v[12]);
        v[13] = SPHelper.GetFloat(volumeSoup, "epka_1.k", v[13]);
        v[14] = SPHelper.GetFloat(volumeSoup, "tc_5_11", v[14]);
        v[15] = SPHelper.GetFloat(volumeSoup, "rd2.out", v[15]);
        v[16] = SPHelper.GetFloat(volumeSoup, "ko4.out", v[16]);
        v[17] = SPHelper.GetFloat(volumeSoup, "ko3.out", v[17]);
        v[18] = SPHelper.GetFloat(volumeSoup, "k254_1.tc", v[18]);
        v[19] = SPHelper.GetFloat(volumeSoup, "ubt_1.k_tm", v[19]);
        v[20] = SPHelper.GetFloat(volumeSoup, "ubt_1.tc", v[20]);
        v[21] = SPHelper.GetFloat(volumeSoup, "ubt_1.gr", v[21]);
        v[22] = SPHelper.GetFloat(volumeSoup, "rs2", v[22]);
        v[23] = SPHelper.GetFloat(volumeSoup, "rs1", v[23]);
        v[24] = SPHelper.GetFloat(volumeSoup, "tm", v[24]);
        v[25] = SPHelper.GetFloat(volumeSoup, "epka_1.gr", v[25]);
        v[26] = SPHelper.GetFloat(volumeSoup, "tc_4_10", v[26]);
        v[27] = SPHelper.GetFloat(volumeSoup, "ubt_1.k_tc", v[27]);
        v[28] = SPHelper.GetFloat(volumeSoup, "k395_1.u", v[28]);
        v[29] = SPHelper.GetFloat(volumeSoup, "k395_1.ur", v[29]);
        v[30] = SPHelper.GetFloat(volumeSoup, "rd2.feed", v[30]);
        v[31] = SPHelper.GetFloat(volumeSoup, "ko4.in", v[31]);
        v[32] = SPHelper.GetFloat(volumeSoup, "ubt_1.tm", v[32]);
        v[33] = SPHelper.GetFloat(volumeSoup, "k254_1.gr", v[33]);
        v[34] = SPHelper.GetFloat(volumeSoup, "tc_1", v[34]);
        v[35] = SPHelper.GetFloat(volumeSoup, "k395_1.tm", v[35]);
        v[36] = SPHelper.GetFloat(volumeSoup, "rs4", v[36]);
        v[37] = SPHelper.GetFloat(volumeSoup, "ubt_1.k_gr", v[37]);
        v[38] = SPHelper.GetFloat(volumeSoup, "ko3.in", v[38]);
        v[39] = SPHelper.GetFloat(volumeSoup, "pm", v[39]);
        v[40] = SPHelper.GetFloat(volumeSoup, "evr.rk", v[40]);
        v[41] = SPHelper.GetFloat(volumeSoup, "vr.kdr", v[41]);
        v[42] = SPHelper.GetFloat(volumeSoup, "vr.tm", v[42]);
        v[43] = SPHelper.GetFloat(volumeSoup, "tc_1_7", v[43]);
        v[44] = SPHelper.GetFloat(volumeSoup, "rd1.set", v[44]);
        v[45] = SPHelper.GetFloat(volumeSoup, "rd2.set", v[45]);
        v[46] = SPHelper.GetFloat(volumeSoup, "ko2.in", v[46]);
        v[47] = SPHelper.GetFloat(volumeSoup, "rs5", v[47]);
        v[48] = SPHelper.GetFloat(volumeSoup, "evr.rd_tc", v[48]);
        v[49] = SPHelper.GetFloat(volumeSoup, "vr.sk", v[49]);
        v[50] = SPHelper.GetFloat(volumeSoup, "tc_2_8", v[50]);
        v[51] = SPHelper.GetFloat(volumeSoup, "rd1.out", v[51]);
        v[52] = SPHelper.GetFloat(volumeSoup, "k3.b", v[52]);
        v[53] = SPHelper.GetFloat(volumeSoup, "k5.a", v[53]);
        v[54] = SPHelper.GetFloat(volumeSoup, "k4.b", v[54]);
        v[55] = SPHelper.GetFloat(volumeSoup, "k4.a", v[55]);
        v[56] = SPHelper.GetFloat(volumeSoup, "k3.a", v[56]);
        v[57] = SPHelper.GetFloat(volumeSoup, "red3.in", v[57]);
        v[58] = SPHelper.GetFloat(volumeSoup, "vr.zr", v[58]);
        v[59] = SPHelper.GetFloat(volumeSoup, "evr.vr_tc", v[59]);
        v[60] = SPHelper.GetFloat(volumeSoup, "ko2.out", v[60]);
        v[61] = SPHelper.GetFloat(volumeSoup, "evr.zr", v[61]);
        v[62] = SPHelper.GetFloat(volumeSoup, "tc_set", v[62]);
        v[63] = SPHelper.GetFloat(volumeSoup, "tc_3_9", v[63]);
        v[64] = SPHelper.GetFloat(volumeSoup, "vr.mk", v[64]);
        v[65] = SPHelper.GetFloat(volumeSoup, "rd1.feed", v[65]);
        v[66] = SPHelper.GetFloat(volumeSoup, "vt1_kn12", v[66]);
        v[67] = SPHelper.GetFloat(volumeSoup, "kn10_vt3", v[67]);
        v[68] = SPHelper.GetFloat(volumeSoup, "rs8", v[68]);
        v[69] = SPHelper.GetFloat(volumeSoup, "red3.out", v[69]);
        v[70] = SPHelper.GetFloat(volumeSoup, "evr.vr_zr", v[70]);
        v[71] = SPHelper.GetFloat(volumeSoup, "tc_2", v[71]);
        v[72] = SPHelper.GetFloat(volumeSoup, "vr.tc", v[72]);
        v[73] = SPHelper.GetFloat(volumeSoup, "k5.out", v[73]);
        v[74] = SPHelper.GetFloat(volumeSoup, "evr.rd_zr", v[74]);
        v[75] = SPHelper.GetFloat(volumeSoup, "vr.zk", v[75]);
        v[76] = SPHelper.GetFloat(volumeSoup, "rs3", v[76]);
        v[77] = SPHelper.GetFloat(volumeSoup, "evr.tc", v[77]);
        v[78] = SPHelper.GetFloat(volumeSoup, "k3.out", v[78]);
        v[79] = SPHelper.GetFloat(volumeSoup, "k4.out", v[79]);
        v[80] = SPHelper.GetFloat(volumeSoup, "k5.b", v[80]);
        v[81] = SPHelper.GetFloat(volumeSoup, "red1.out", v[81]);

        Soup flowrateSoup = sp.GetNamedSoup("flowrate");
        float[] afr = scheme.alinkFlowRate;
        float[] fr = scheme.linkFlowRate;
        fr[0] = SPHelper.GetFlowRate(flowrateSoup, "k254_1.vr_b", fr[0]);
        afr[0] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.d_at", afr[0]);
        fr[1] = SPHelper.GetFlowRate(flowrateSoup, "red1.feed", fr[1]);
        afr[1] = SPHelper.GetFlowRate(flowrateSoup, "epka_1.tm_at", afr[1]);
        afr[2] = SPHelper.GetFlowRate(flowrateSoup, "kp3.dump", afr[2]);
        afr[3] = SPHelper.GetFlowRate(flowrateSoup, "k254_1.b_at", afr[3]);
        fr[2] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.gr_u_red", fr[2]);
        fr[3] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.gr_u", fr[3]);
        afr[4] = SPHelper.GetFlowRate(flowrateSoup, "epka_1.k_at", afr[4]);
        afr[5] = SPHelper.GetFlowRate(flowrateSoup, "rd2.dump", afr[5]);
        afr[6] = SPHelper.GetFlowRate(flowrateSoup, "k254_1.tc_at", afr[6]);
        fr[4] = SPHelper.GetFlowRate(flowrateSoup, "ubt_1.tc", fr[4]);
        fr[5] = SPHelper.GetFlowRate(flowrateSoup, "ubt_1.gr", fr[5]);
        fr[6] = SPHelper.GetFlowRate(flowrateSoup, "kn27", fr[6]);
        fr[7] = SPHelper.GetFlowRate(flowrateSoup, "tm_epka_1", fr[7]);
        fr[8] = SPHelper.GetFlowRate(flowrateSoup, "kn18", fr[8]);
        fr[9] = SPHelper.GetFlowRate(flowrateSoup, "epka_1.gr_k", fr[9]);
        fr[10] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.u_s", fr[10]);
        afr[7] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.u_at", afr[7]);
        fr[11] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.ur_d", fr[11]);
        afr[8] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.ur_at", afr[8]);
        fr[12] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.ok", fr[12]);
        fr[13] = SPHelper.GetFlowRate(flowrateSoup, "rd2.feed", fr[13]);
        fr[14] = SPHelper.GetFlowRate(flowrateSoup, "ko4.in_out", fr[14]);
        fr[15] = SPHelper.GetFlowRate(flowrateSoup, "ubt_1.tm", fr[15]);
        afr[9] = SPHelper.GetFlowRate(flowrateSoup, "ubt_1.tm_at", afr[9]);
        fr[16] = SPHelper.GetFlowRate(flowrateSoup, "k254_1.gr_tc", fr[16]);
        fr[17] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.tm_gr", fr[17]);
        fr[18] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.u_tm_gr", fr[18]);
        afr[10] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.u_tm_at", afr[10]);
        afr[11] = SPHelper.GetFlowRate(flowrateSoup, "k395_1.tm_at", afr[11]);
        fr[19] = SPHelper.GetFlowRate(flowrateSoup, "ko3.in_out", fr[19]);
        fr[20] = SPHelper.GetFlowRate(flowrateSoup, "kn14", fr[20]);
        fr[21] = SPHelper.GetFlowRate(flowrateSoup, "pm_epka_1", fr[21]);
        fr[22] = SPHelper.GetFlowRate(flowrateSoup, "kn6", fr[22]);
        afr[12] = SPHelper.GetFlowRate(flowrateSoup, "evr.leak_rk", afr[12]);
        afr[13] = SPHelper.GetFlowRate(flowrateSoup, "evr.vo", afr[13]);
        afr[14] = SPHelper.GetFlowRate(flowrateSoup, "vr.kdr_at", afr[14]);
        afr[15] = SPHelper.GetFlowRate(flowrateSoup, "vr.mk_at_sk", afr[15]);
        fr[23] = SPHelper.GetFlowRate(flowrateSoup, "ko2.in_out", fr[23]);
        afr[16] = SPHelper.GetFlowRate(flowrateSoup, "evr.rd_rel", afr[16]);
        fr[24] = SPHelper.GetFlowRate(flowrateSoup, "vr.sk_tc", fr[24]);
        afr[17] = SPHelper.GetFlowRate(flowrateSoup, "rd1.dump", afr[17]);
        fr[25] = SPHelper.GetFlowRate(flowrateSoup, "k3.b_out", fr[25]);
        fr[26] = SPHelper.GetFlowRate(flowrateSoup, "k5.a_out", fr[26]);
        fr[27] = SPHelper.GetFlowRate(flowrateSoup, "k4.b_out", fr[27]);
        fr[28] = SPHelper.GetFlowRate(flowrateSoup, "k4.a_out", fr[28]);
        fr[29] = SPHelper.GetFlowRate(flowrateSoup, "k3.a_out", fr[29]);
        fr[30] = SPHelper.GetFlowRate(flowrateSoup, "red3.feed", fr[30]);
        fr[31] = SPHelper.GetFlowRate(flowrateSoup, "vr.mk_kdr", fr[31]);
        fr[32] = SPHelper.GetFlowRate(flowrateSoup, "vr.mk_sk", fr[32]);
        fr[33] = SPHelper.GetFlowRate(flowrateSoup, "vr.mk_zk", fr[33]);
        fr[34] = SPHelper.GetFlowRate(flowrateSoup, "rd1.feed", fr[34]);
        fr[35] = SPHelper.GetFlowRate(flowrateSoup, "kn12", fr[35]);
        fr[36] = SPHelper.GetFlowRate(flowrateSoup, "vt3", fr[36]);
        fr[37] = SPHelper.GetFlowRate(flowrateSoup, "vt2", fr[37]);
        afr[18] = SPHelper.GetFlowRate(flowrateSoup, "vr.tc_at", afr[18]);
        fr[38] = SPHelper.GetFlowRate(flowrateSoup, "evr.vt", fr[38]);
        fr[39] = SPHelper.GetFlowRate(flowrateSoup, "evr.rd_feed", fr[39]);
        fr[40] = SPHelper.GetFlowRate(flowrateSoup, "vr.zk_tc_e", fr[40]);
        fr[41] = SPHelper.GetFlowRate(flowrateSoup, "vr.zk_tc", fr[41]);
        fr[42] = SPHelper.GetFlowRate(flowrateSoup, "kn11", fr[42]);
        fr[43] = SPHelper.GetFlowRate(flowrateSoup, "kn10", fr[43]);
        fr[44] = SPHelper.GetFlowRate(flowrateSoup, "evr.pk_evr", fr[44]);
        fr[45] = SPHelper.GetFlowRate(flowrateSoup, "evr.pk_vr", fr[45]);
        fr[46] = SPHelper.GetFlowRate(flowrateSoup, "vt1", fr[46]);
        fr[47] = SPHelper.GetFlowRate(flowrateSoup, "k5.b_out", fr[47]);

        Soup objectsSoup = sp.GetNamedSoup("objects");
        SetCurrentObjectName("kp3");
        scheme.objects[0].SetProperties(objectsSoup.GetNamedSoup("kp3"));
        SetCurrentObjectName("epka_1");
        scheme.objects[1].SetProperties(objectsSoup.GetNamedSoup("epka_1"));
        SetCurrentObjectName("ubt_1");
        scheme.objects[2].SetProperties(objectsSoup.GetNamedSoup("ubt_1"));
        SetCurrentObjectName("k395_1");
        scheme.objects[3].SetProperties(objectsSoup.GetNamedSoup("k395_1"));
        SetCurrentObjectName("k254_1");
        scheme.objects[4].SetProperties(objectsSoup.GetNamedSoup("k254_1"));
        SetCurrentObjectName("ko3");
        scheme.objects[5].SetProperties(objectsSoup.GetNamedSoup("ko3"));
        SetCurrentObjectName("red1");
        scheme.objects[6].SetProperties(objectsSoup.GetNamedSoup("red1"));
        SetCurrentObjectName("ko2");
        scheme.objects[7].SetProperties(objectsSoup.GetNamedSoup("ko2"));
        SetCurrentObjectName("ko4");
        scheme.objects[8].SetProperties(objectsSoup.GetNamedSoup("ko4"));
        SetCurrentObjectName("vr");
        scheme.objects[9].SetProperties(objectsSoup.GetNamedSoup("vr"));
        SetCurrentObjectName("evr");
        scheme.objects[10].SetProperties(objectsSoup.GetNamedSoup("evr"));
        SetCurrentObjectName("k3");
        scheme.objects[11].SetProperties(objectsSoup.GetNamedSoup("k3"));
        SetCurrentObjectName("k4");
        scheme.objects[12].SetProperties(objectsSoup.GetNamedSoup("k4"));
        SetCurrentObjectName("k5");
        scheme.objects[13].SetProperties(objectsSoup.GetNamedSoup("k5"));
        SetCurrentObjectName("red3");
        scheme.objects[14].SetProperties(objectsSoup.GetNamedSoup("red3"));
        SetCurrentObjectName("rd2");
        scheme.objects[15].SetProperties(objectsSoup.GetNamedSoup("rd2"));
        SetCurrentObjectName("tc_4_10");
        scheme.objects[16].SetProperties(objectsSoup.GetNamedSoup("tc_4_10"));
        SetCurrentObjectName("tc_5_11");
        scheme.objects[17].SetProperties(objectsSoup.GetNamedSoup("tc_5_11"));
        SetCurrentObjectName("tc_6_12");
        scheme.objects[18].SetProperties(objectsSoup.GetNamedSoup("tc_6_12"));
        SetCurrentObjectName("rd1");
        scheme.objects[19].SetProperties(objectsSoup.GetNamedSoup("rd1"));
        SetCurrentObjectName("tc_3_9");
        scheme.objects[20].SetProperties(objectsSoup.GetNamedSoup("tc_3_9"));
        SetCurrentObjectName("tc_2_8");
        scheme.objects[21].SetProperties(objectsSoup.GetNamedSoup("tc_2_8"));
        SetCurrentObjectName("tc_1_7");
        scheme.objects[22].SetProperties(objectsSoup.GetNamedSoup("tc_1_7"));
        ResetCurrentObjectName();

    }

    public Soup GetProperties()
    {
        return inherited();
    }

    public void SetProperties(Soup sp)
    {
        inherited(sp);
    }

    public void Init(Asset asset)
    {
        inherited(asset);

        InitSpec();
    }
};