Closed adhishm closed 11 years ago
The function now browses all defects.
In addition, the safeguard to prevent division by zero has been placed within the Dislocation::stressFieldLocal
function, thus removing the need to check for distances (or pointer equality) in any functions calling the stressField function.
In the function
SlipPlane::calculateDislocationStresses
, currently only dislocations are browsed and their stress fields calculated. However, other defects may have stress fields too. Therefore, for each dislocation, the total stress field at its position should be calculated using contributions from all defects. Moreover, the classDefect
has a virtual functionDefect::stressField
precisely for this reason.