Paper presents results obtained when porting FEM 2D linear elastostatic local stiffness matrix calculations to Tesla architecture with OpenCL framework. Comparison with native NVIDIA CUDA implementations has been provided. I. MOTIVATION IN FINITE element simulations usually two computation stages most significantly impact performance of the whole process: