Numerical solution of PDEs with hybrid and heterogeneous computing models