This note deals with two fully parallel methods for solving linear partial differentialalgebraic equations (PDAEs) of the form: Aut + B∆u = f(x, t) where A is a singular, symmetric and nonnegative matrix, while B is a symmetric positive define matrix. The stability and convergence of proposed methods are discussed. Some numerical experiments on high-performance computers are also reported.