8 #ifndef IMPSAXS_SOLVENT_ACCESSIBLE_SURFACE_H
9 #define IMPSAXS_SOLVENT_ACCESSIBLE_SURFACE_H
11 #include <IMP/saxs/saxs_config.h>
14 IMPSAXS_BEGIN_NAMESPACE
33 IMP::Floats get_solvent_accessibility(
const core::XYZRs& points,
34 float probe_radius = 1.8,
39 const float radius1,
const float radius2) {
40 float squared_radius_sum = (radius1+radius2)*(radius1+radius2);
42 algebra::get_squared_distance(sphere_center1, sphere_center2);
43 if(fabs(squared_radius_sum - squared_dist) < 0.0001)
45 if(squared_radius_sum > squared_dist)