int max_pts = chi("max_pts");
float max_rai = chf("max_rai");
int muscle_pts = npoints(1);
for(int i = 0; i<muscle_pts; i++)
{
float SpringMask = point(1, "SpringMask", i);
if (SpringMask > chf("SpringMask_Thr") )
{
vector muscle_pos = point(1, "P", i);
int bone_pts[] = nearpoints(2, muscle_pos, max_rai, max_pts);
foreach(int bone_pt; bone_pts)
{
int muscle_spring_pt = addpoint(0, muscle_pos);
// setpointattrib(0, "muscle_pt_index", muscle_spring_pt, i);
vector bone_pos = point(2, "P", bone_pt);
int bone_spring_pt = addpoint(0, bone_pos);
// setpointattrib(0, "bone_pt_index", bone_spring_pt, bone_pt);
int spring = addprim(0, "polyline", muscle_spring_pt, bone_spring_pt);
dict rel;
rel["muscle_pt_index"] = i;
rel["bone_pt_index"] = bone_pt;
setprimattrib(0, "extraSpring", spring, rel);
}
}
}