import fluidsim;
#pragma MSP stage(compute)
+float get_clamping(float depth, float v_right, float v_up, float v_left, float v_down)
+{
+ float total_out = max(v_right, 0.0)+max(v_up, 0.0)-min(v_left, 0.0)-min(v_down, 0.0);
+ return min(max_flow_fraction/max(total_out*delta_time, max_flow_fraction), depth/residual_depth);
+}
+
void main()
{
ivec2 size = imageSize(surface_out);