high performance computing on graphics processing units: hgpu.org

We present a multi-frontal direct solver for two dimensional isogeometric finite element method simulations with NVIDIA CUDA and perform numerical experiments for linear, quadratic and cubic B-splines. We compare the computational cost O(Np^2) for 2D parallel shared memory implementation with the corresponding estimate O(N^1.5p^3) for a standard 2D sequential implementation. We conclude the presentation with observation that computational cost of the shared memory direct solver scales like p^2 when we increase the global continuity of the isogeometric solution, which is an adventage with respect to sequential isogeometric solver scalability of the order of p^3.