A Comparative Study on Different Parallel Solvers for Nonlinear Analysis of Complex Structures

Algorithm 1

extern “C” int VeksInit(UINT32 nElems, UINT32 nGlobalNodes, UINT32 aiElemConn[ ][8]);

// Funcition:       Initialize the solver

extern “C” int VeksSetK(double afK[ ][24][24]);

// Funcition:       Set element stiffness matrices

extern “C” int VeksSetDemF(double afF[ ][8][3]);

// Funcition:       Set element loads

extern “C” int VeksSetU(double afU[ ][3]);

// Funcition:       Set nodal displacements

extern “C” int VeksSetF(double afF[ ][3]);

// Funcition:       Set nodal loads

extern “C” int VeksFixU(char abU[ ][3], double afU[ ][3]);

// Funcition:       Set nodal displacements

extern “C” int VeksGetU(double afU[ ][3]);

// Funcition:       Get nodal displacements

extern “C” int VeksSetC(double afC[ ][8][8]);

// Funcition:       Set element conductivity matrices

extern “C” int VeksSetDemQ(double afQ[ ][8]);

// Funcition:       Set element thermal loads

extern “C” int VeksSetQ(double afQ[ ]);

// Funcition:       Set nodal thermal loads

extern “C” int VeksFixT(char abT[ ], double afT[ ]);

// Funcition:       Set nodal thermal constraints

extern “C” int VeksGetT(double afT[ ]);

// Funcition:       Get nodal temperatures

extern “C” __int64 VeksSolve(__int64 nMaxIter, double fTimeStep, double fUnbRatio, double fEps, unsigned int nPrint);

// Funcition:       Solve equations

// The procedures for GPU analysis

#include “Eks.h”

int nRet = VeksInit(nElems, nGlobalNodes, aiElemConn);

VeksSetK(afK);

VeksSetDemF(afDemF);

VeksFixU(abU, afU);

VeksSetC(afC); VeksSetDemQ(afDemQ);

VeksFixT(abT, afT);

VeksSolve(nSteps, fTimeStep, fUnbRatio, fEps, nStipePrint);

VeksGetU(afU);

VeksGetT(afT);