Tesan Srdjan

float droop_amount = chf("droop_amount");     
float tension = chf("tension");            
int use_ramp = chi("use_ramp");              
float axis_blend = chf("axis_blend");      
int preserve_length = chi("preserve_length");


int prims[] = pointprims(0, @ptnum);
if (len(prims) == 0) return; 
int prim = prims[0]; 

int vtx = vertexindex(0, prim, @ptnum);
int count = primvertexcount(0, prim);


float t = float(vtx) / max(count - 1, 1);


if (use_ramp) {
    t = chramp("droop_ramp", t);
}


vector axis_dir = normalize(lerp({1,0,0}, {0,0,1}, axis_blend));


float arch = pow(t * (1.0 - t), tension);
float dy = -arch * droop_amount;

@P.y += dy;

/*
if (preserve_length) {
    float stretch = 1.0 + abs(dy) * 0.05;
    vector center = getbbox_center(0);
    vector dir = normalize(@P - center);
    @P = center + dir * length(@P - center) * stretch;
}
*/