febiosoftware / FEBio

FEBio Suite Solver
https://febio.org/
MIT License
174 stars 65 forks source link

Faster normal projection #47

Closed spietz closed 1 year ago

spietz commented 2 years ago

I was able to reduce the overall solution time by a factor of 3 for structural contact analysis by modifying normal-projection and octree search to only look for intersections in octree nodes within the specified search radius. The result of the simulation was the same.

I noticed that the search radius seems to be compared to absolute values in the code, while it is stated in the documentation that it is relative to the diagonal of model's bounding box. Maybe double check this in the code and documentation?

SteveMaas1978 commented 1 year ago

@spietz

Thanks for contributing this! Sorry it took so long to get this merged.