We investigate the weighted sum-rate (WSR) maximization linear precoder design for massive MIMO downlink and propose a unified matrix manifold optimization framework applicable to total power constraint (TPC), per-user power constraint (PUPC) and per-antenna power constraint (PAPC). Particularly, we prove that the precoders under TPC, PUPC and PAPC are on different Riemannian submanifolds, and transform the constrained problems in Euclidean space to unconstrained ones on manifolds. In accordance with this, we derive Riemannian ingredients including orthogonal projection, Riemannian gradient, Riemannian Hessian, retraction and vector transport, which are needed for precoder design in matrix manifold framework. Then, Riemannian design methods using Riemannian steepest descent, Riemannian conjugate gradient and Riemannian trust region are provided to design the WSR-maximization precoders under TPC, PUPC or PAPC. Riemannian methods are free of the inverse of large dimensional matrix, which is of great significance for practice. Complexity analysis and performance simulations demonstrate the advantages of the proposed precoder design.