found bug in surface tension computation; started implementing grid structure; tested on density computation so far. it gives plausible results but they are a bit different than the ones without grid