2
$${\Delta c^\star = k[c\rho + (grad c^\star)^2/2kc]}$$