I have a rather basic question about how PyBAMM handles discretizations.
From what I understand the gradient of a variable pybamm.grad(c), has shape N-1. Now, if this gradient is to be multiplied by a function parameter such as FunctionParameter("parameter X", {"c": c}}. Then I assume that this parameter will have shape N? Am I correct? If so, how does handle mutliplying the gradient with the function parameter?