/*
original_author:  Inigo Quiles
description: Union operation of two SDFs 
use: <float> opUnion( in <float|float4> d1, in <float|float4> d2 [, <float> smooth_factor] ) 
*/

#ifndef FNC_OPUNION
#define FNC_OPUNION

float opUnion( float d1, float d2 ) { return min(d1, d2); }
float4  opUnion( float4 d1, float4 d2 ) { return (d1.a < d2.a) ? d1 : d2; }

// Soft union
float opUnion( float d1, float d2, float k ) {
    float h = saturate( 0.5 + 0.5*(d2-d1)/k );
    return lerp( d2, d1, h ) - k*h*(1.0-h); 
}

float4 opUnion( float4 d1, float4 d2, float k ) {
    float h = saturate( 0.5 + 0.5*(d2.a - d1.a)/k );
    return lerp( d2, d1, h ) - k*h*(1.0-h); 
}

#endif