- 01
 - 02
 - 03
 - 04
 - 05
 - 06
 - 07
 - 08
 - 09
 - 10
 - 11
 - 12
 - 13
 - 14
 - 15
 - 16
 - 17
 - 18
 - 19
 - 20
 - 21
 - 22
 - 23
 - 24
 - 25
 - 26
 - 27
 - 28
 - 29
 - 30
 
                        QVector<double*>*** ElemBoundCond = tbc->getElementsBC();
    QVector<double*>* preNodalVals = tbc->getPreNodalValues();
    QVector<double*>* nodalForces = tbc->getNodalForces();
    long elems = grid->getElemsNumber();
    long nodes = grid->getNodesNumber();
    long** Elements = grid->getElems();
    double** Nodes = grid->getNodes();
    long* ElemMaters = grid->getElemsMaterials();
    long NBW = grid->getBandWidth();
    emit signal_addLineToCalculationProtocol(">>> Solve heat problem ");
    long    NS[4],//+numbers of 4 nodes
            NCL = 1,//number or unknown values per node (1 for thermal analysis)
            ni, nj, nk,//face nodes
            gnN, dN;
    double  ESM[4][4],//+element stiffness matrix
            EF[4],//+element force matrix
            **C,//matrix conteins X, Y, Z coordinates, rows like |1 Xi Yi Zi|
            **CA,//conteins columns like |ai bi ci di|
            //*GF,//band matrix, pointers to Global Force and Global Sriffness Matrices in A
            V36,//+hexa and 36 element volume
            S_3,//face area div 3
            hS_3, hS_6, hS_12,//h * face area div 6 and div 12
            hST_3,//h * face area * T / 3
            qS_3;//q * face area / 3
    QVector<double*>* fcond;
    MaterialsDatabase* matdb = MaterialsDatabase::getInstance();
    double  Kxx, Kyy, Kzz;//thermal condactivities in X, Y and Z axes
                                 
        
достойно отдельного ГК