Large improvements by not calculating useless stuff. For instance don't advect solute charges for salt free systems...
No parallelization here, because it seems useless to me for now. PARALLEL DO could be done easily in advect, comp_rho, comp_j, propagate.
Also, I often tabulate function pbc to move it out of loops.
Large improvements by not calculating useless stuff. For instance don't advect solute charges for salt free systems...
No parallelization here, because it seems useless to me for now. PARALLEL DO could be done easily in advect, comp_rho, comp_j, propagate. Also, I often tabulate function pbc to move it out of loops.