fixed a performance issue, neutral particles connected and disconnected parts are the same with Markus program. There is some difference in the calculation of the d-propagators.