From: bkraczek <bkraczek@ic...>  20081204 18:39:58

I'm new to libmesh, and I am not a great C++ programmer, so the answer to this may seem obvious. I am working from example 3 to implement a poisson solver with a RHS that includes the dot product of the weighting function with a vector that I set. So, working from line 351 of ex3.C, Fe(i) += JxW[qp]*fxy*phi[i][qp]; I want to change this to Fe(i) += JxW[qp]*dphi[i][qp]* gradGat ; where gradGat is the vector I desire to take the inner product with. How do I properly construct gradGat so I can get the desired dot product? Thanks, Brent 