I'm looking for the set of m*n matrices W that simultaneously solve:

diag(AW)=b and diag(WA)=c for a given n*m matrix A and vectors b and c. With diag() I mean the vector of diagonal elements.

I'm happy even for just a starting point, or a single solution.