Velocibox.fuse

--[[--/*

    Velocibox.fuse

    Based on https://www.shadertoy.com/view/lsdXD8 a WebGL
    shader from shadertoy.com converted to DCTL and embeddet
    into a Lua Fuse by nmbr73 (https://www.youtube.com/c/nmbr73)
    for use in DaFusion.

    Shader created by zackpudil

	V0.2 bugfixes for OpenCL and CUDA

*/--]]--

-- /* FUSE_COFIG::BEGIN
local FC_CATEGORY  = "AbstractShader"
local FC_SHORTCUT  = false        -- use shortcut (e.g. DTD for DancyTreeDoodle) if exists; otherwise the shadertoy ID is used
local FC_DEVEVELOP = true         -- set to 'true' to display reload/edit buttons in inspector
local FC_INFOBTN   = -1           -- show button to open the project page -1 on top, 0 hide, 1 on bottom
local FC_PREFIX    = "ST"         -- use 'ST' in the repo; to be replaced by 'BETA' for user downlaods and by "" if it goes officially into the reactor
local FC_SCPREFIX  = FC_PREFIX    -- prefix for shortcuts
local FC_SUBMENU   = "Shadertoys" -- you may want to prepend e.g. "Fuses\\" to not have the Shadertoys on your 'Add tool...' toplevel menu
-- FUSE_COFIG::END */

-- /*

local shadertoy_name       = "Velocibox"
local shadertoy_author     = "zackpudil"
local shadertoy_id         = "lsdXD8"
local shadertoy_license    = ""
local dctlfuse_code        = true
local dctlfuse_name        = "Velocibox"
local dctlfuse_author      = "nmbr73"
local dctlfuse_authorurl   = "https://www.youtube.com/c/nmbr73"
local dctlfuse_infourl     = "https://nmbr73.github.io/Shadertoys/"
local dctlfuse_versionNo   = 2
local dctlfuse_versionDate = ""



-- MANDATORY -----------------------------------------------------------------
local shadertoy_name       = "Velocibox"
local shadertoy_author     = "zackpudil"
local shadertoy_id         = "lsdXD8"
local shadertoy_license    = ""
local dctlfuse_category    = "Tunnel"
local dctlfuse_name        = "Velocibox"
local dctlfuse_author      = "nmbr73"
-- OPTIONAL ------------------------------------------------------------------
local dctlfuse_versionNo   = 2
local dctlfuse_versionDate = ""
local dctlfuse_authorurl   = "https://www.youtube.com/c/nmbr73"
local dctlfuse_authorlogo  = 'width="212" height="41" src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAANQAAAApCAYAAABN0gffAAABg2lDQ1BJQ0MgcHJvZmlsZQAAKJF9kT1Iw0AcxV9TRZFKh3YQEclQnSyIijhqFYpQIdQKrTqYXPoFTRqSFBdHwbXg4Mdi1cHFWVcHV0EQ/ABxc3NSdJES/5cUWsR4cNyPd/ced+8AoVFhmtU1Dmi6baaTCTGbWxV7XiEgjAhiGJaZZcxJUgq+4+seAb7exXmW/7k/R7+atxgQEIlnmWHaxBvE05u2wXmfOMpKskp8Tjxm0gWJH7muePzGueiywDOjZiY9TxwlFosdrHQwK5ka8RRxTNV0yheyHquctzhrlRpr3ZO/MJTXV5a5TnMISSxiCRJEKKihjApsxGnVSbGQpv2Ej3/Q9UvkUshVBiPHAqrQILt+8D/43a1VmJzwkkIJoPvFcT5GgJ5doFl3nO9jx2meAMFn4Epv+6sNYOaT9Hpbix0B4W3g4rqtKXvA5Q4w8GTIpuxKQZpCoQC8n9E35YDILdC35vXW2sfpA5ChrlI3wMEhMFqk7HWfd/d29vbvmVZ/P2tycqR3fEbRAAAABmJLR0QA/wD/AP+gvaeTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAAB3RJTUUH5QIPFAcPbC/jXQAAABl0RVh0Q29tbWVudABDcmVhdGVkIHdpdGggR0lNUFeBDhcAAA9tSURBVHja7Z15eFRVmsZ/d6lKbVlYQiAJmISQoCAIaABtFWVaHsOm3bboINJuOLajz0zPMz099rRiL04r0qLS6jgt2i4tERC3RLttlxmNyAAqi6wxBkKoEBKyp9Z77/xRRaXurZsNCES47z+hTp06yz3nPd/3veeci0AfsKBs9mRgiyH5i5Li0imcAiwom70FmGxInlJSXPoFFiwMAIjWI7BgwSKUBQsWoSxYONMhn4QyhGhsFY9gSXHpDuvxWjjbIJwEoeCkCxWnW/ywYMFy+SxYsAhlwYJFKAsWLHSBkyFKBAHjxuoOk7hoBHCnIdlbUlz6X9YwWDhT0CdRoouTCs8CSwxp/wNc3ot8x4hoLHMh8IpJXmM+SopLBWsYLZxJFmpAo2jq9LuPrQfA9DO1z7IsMzI7C1nu7J6qaVRWfoumaaelTcOHD+fR5csRRRFORhuEbtZOTUMIhyJ/NdCESH5NlkHoMrI5CjjiPvuBy/PzcnYAVFRWtRt/kJ+X4z7274rKqi1nHaGAPwHZwCWACmichO2CgYaCMWN47LHlpKamdvriwSBz5s6jpaX19LpBAT+eD95D9AeOu4xQdjYdU7+XQExBUbBVV+H4YhO2Lf+LeLACwe9Dc7pRR44mfMEl+M+/gGDeGDTZZix2sOGzC7AbPneHyWcjoUYCClANNABDo/0+s0glgCAICHGruCgODM1JCIdIeudlpMrjPMOsaQg3/muEUHEQfR0kv1GC4+XfgqLq62wGsXYv8qZ3cQgawfn30HL9zSiDh556l29B2Ww7MN5EfHCaZDcbNbuZBxCNreKxBZhrkvfcPsR1k02EDm/c5yuAKmAPsB74ITDM8va/YzB4jEIwSMqLz5C0/sme10ZNwL7+CdJqD9L00/tRUtJOeQw1HpOTCoDPJO/TJcWldxgm+dJovBKPVBOh4mPgz8ADhnTFTGxYUDbbzBE3tvNBYGnc50uBVqAccCcOjYXTwg+nA0TB1NQKwTCEgp1JziEEx03QuXvOrzaR9NYqPZnsLsLnXw52B9L2TxDajnR+L4jIn72Oq+gyWq+e3308ZokS3WINUBeNo84BpPgvPR4PQ4cM0Y1LR4ePuro6XC4X6elDEQQBJazgra0lFArp3KoE02y3M3x4hs7dCocVDh48iM1mY3hGBpIsoWnQ1tbG0aNHdaKBpmmkpCSTljYISRLRNI0jR+ppb2/vtt7uIEkyw4dnYLPJoEFDQwMtra3dlhff92MIBAJ4vbW6fOnp6bjdkVDj8OE6fD5fL2adh6aVpSgeT2Jbm5tI/d19iJUbowkyHXf8Ox0Xxq3PqoKz7HVQOjqf27BCWn72EIHRBSAIyA31pDz7OPLG9XG+lEjSB2/QfuUsVIcT4HogPrAKod/yWdhDT6acjYR6H5gDXBi1Vjp7f838edx66y26ibNz5y5279nN3DlzsNlssYkeCASoqKhg1arn+fKrreZR6qQLuP/+X+J0dnrHLS0tPLfqeW6/7TY8HnesLlXVaGho4P333+f5F17g8ssu44YbbyA3JwdRFBEEAQ1QwmFqamp44403efOtt/uk2k2fPo0lS5YwKC0tVq+iKHi9XkpLy1izdh2KoiT87u9mXsm9996DJHWuPwcOHOCGv7+JvJwc7rrrH5h4wUTcLlcsdnvyyZWse319L4IqESUlFSU5RZ8cCpK89mXEyg2xSCI461barpoHYmc7pKP1SFX74sydSnDmPPzjJ3UyI2sU7dcuIHVTKaid1k78ZhNE+5ufl7Omu2bm5+X8uYfvvzgbCRWIqjE/AgZFLZQQb1Fc0UlxDMOGpXPuuWNxu926gtxuN0VFRRQWFvLkyj/w9tvvJC6+sozL5dIRCmDRTQvJyEgM3TweNz/+8WLGjTuPcePGkZKSYtqJlJQURo8ezdixhTz8yKO9IpUgCNy0cCHZWVkJ3yUnJ5Obm0tBQQG/+e1DqKo+qLfZbLhcLh2hnE4nI7OzWLHi94wcOTLBwtnstq4bk+RES86OMsKe6HJpGu5PPyLp9T/EyKRlnEfrdQvRbPpyNYeTtrt/jhgIIjY3IdYdJDBuQnRY46y9zaQeR0q/uXvdEaqr0w9DgN3G2MhEGJAArzEuIrK5G4/90TY8axA5agxlBrvpg7GdxnodROTywVEZtMenmZmZSSAQYMOGDQSCQc4dO5aMjIzODqemcsftt7H1q60cqK7u8SE7HA7S09PZum0bDfUNnJNzDnm5ubEJKcsy06dPj1mPr3fupP5IPaNGjSIvLzfmPtrtdq6++mq2bdtO2bvv9cLVkxg1ahS7du3C6/WSmZlJQUFBrDybzcZVV32fXbt2sWbtul5NmHvv+UdGjRrVp0mmOt003vcIQswSCqguvSJtP1CJ6+lfgRIdakHCd+MSwiOyE8vzpNBx0cUxIqJpUZJ0kkn0+3Bt+AQUvVSvjPteZG+K2D6TURqfcszyRPeZJnf1vSmhFpTNHk9kr8ZInh0GpW+8SeFEhQajqPAsMMKQ1gxcaUi7sIsyb++BNPH1GE9fTIye6DimDknBo2G5ZmWLjV5K5YqisOLxJ3jnnVI0TSM3N4dljzxCdnZWnBUbxpw5s3nq6Wd6NalKS8v4/WOPoSgqbpeLZcseZsqUKYZFWuPtt99h2aPLCYfDeDwe/vOh31BUVKQj5zXXzOfd9/7SKyv1ySef8qtf/5pgMITdbufBBx/gihkzYmS22WwUzy5mdclrOmvUFUEnTZrU92VbFFHSR3Qrq7vL1iO0HIoNkZo1jo7pl/ZsTQQhlkdqOorsPYh8pA7Hh2XIG98ytEPCN2s+mj0pFir2h4WydzGp6SZ9IEPStVsAySkeiRNfeyRV1f79vPvuezE36JtvKvnwow+5edEinTuVPyY/Euf0MLHb29v543OrUKJ7Je0dHaxZs5aJEyfqTjY0NTWxdt06wuEwEBEtXl29WkcogEGDBpGWlkpjY1O39YbDYZ5b9TzBYChi5oNBSla/xrSpU3HFWYjMESPweDw9CgppaWnIsozf7ycQCCS4l36//7gGzF6xB/tbK+ncgRHw/2Bx3+RtQcD1yfu4nvg3kITEYZbt+G75Jb7J0/p18p0NMRSaShJQH3VZ5ejIdUmscCiE3++PrdiCINBQ35CQz+lwIEki4bDSff1RpU6WOy1AzaFDCUQMBoPUG+ppbmpG0zRdvOJyuXC73D0SCqChoV7vD9d6E0QISZIYPHgwNTU13Zbl8Xj4+uudvPjSS+zbuw/FEHc1Njb2eWyEYAD3G6tB61RFNXc6vqmXHEesI0RHVf87zZVG+0+X0THt8pi7ZxHqBCDItAFPATOATKDQKJ+fAloPjMXlBJrR2NjI7x5+mL1795209ti/3Yft41d0JAhdcS1K2uCTN/6hII5PP0ZNTsV3/iSdYni6CVUBhE1cLKNQ0WAiQKgm5XlN4iMniRvIXb2fQumVCy8LbcCXQCWRvYdCLPQZra2t7Ny5S+emnhBUFefmjXqLoqkELroYTep7HUpqKqEZCxAaapG++RLB3xZlaAfyx6+QsvE9pLvup21mMUS8D+M+kjc/LycmauXn5Uw5HkJpJoSoAe4j8QjRtcB/GNIuBcYa0oZgflXDmG8s5tcyjGLD5Ghaj9c3osemLos9ZL+aSeSo0XAip4ktHLd101AU5aQRSmxvxf6X1/R1OFwECscfl+n1XTwT37QZoKpIjUdJXr8a+1tPdXoHvgbcK39BcHQBwdGFdKfWnYiFEkwmeofZW4sWlM2eN9AHvaS4NLigbHasP5Ij5pvfFrWmbosaAwOOPTsR6r/VW5kJM1Edjt77r8d8WFGMxkcRsodHZNFy020M2rkZqWJjXKDahuvTDwnm5lsx1AniIyAXuA7r6v+AUIrs2zaDqo8gwhOmoiV1TyjXl5uQmpoQ62sR6r2gKrQsvhvVo98UVzzJhItm6AkFSLu/RFAUKiqrlpu4fU/m5+WsswjVM1YAP7dm8gARilQNefvmhJgqPGRwj+pe0gel2P/2Yic3h+QgXrc4gVBoQCDRyxeCMcn/RhL3S/cAp4xQHV0IEEZftJZevGfiFCObyOZyXTSeks/GiZyUZDf1nEKh0Clth9jUgFh7KIEAakpqj78NF47H/rfOLUXhSCVJe3cRzsjUT+6GI9jL/5ooXoweh9ZPd8Xk6EQ3KnJdqWcuk3irpKS49IcmeR84wbbdaRLDeXv521JgdHSQRE3DNni+c8LRN32biZwuPistlaZpTJs2FYchRmlvb6O29jCi2L93LqW2FqToXpXtmz0Ibd4EQokdPuzV+yOTMG1QwgFagMCESTiTsxBaD8XiJ9dLz6C63JHYSBSQ64/gefUFhNpdCa5m4MJLQOqf9VSOBvFLTL77yQCcE2YkXWqSNiJG/OjJFClZ9AGfRX1m9Wwh0dChQ3G53CQl2SkquoibFy2KnaA/RrLt27f353nRTuu4cxvJ/3RTRDuWAJvBSsgingdvjyznArQ9/CrtRdMTyglm5xCcdT1Ja1d0WrwDW0i57wbU4QVgdyLWbAU18cREePoP8J8/qd/6OCDcngVlszMBX0lxaWN/1aEFNBcROb0VaCMyrGf0G5NkWeaZp5+KrCuCgNPp1JEJoLa2lmf/+4/Hfdeqb4GTEPFx5G7cLZsY8SHizuglmjqJ1usWIdVUI2+IC3cEEfFwxTGKGWcAyvlX0fyTn6E6Y8euXjURJf76nSdUSXHpof6uQw1pLmB71Hp1kPiCjjMv8BcE3UtbjO7f7t27WfH4E1RXH/zO9U1JG0TTP/8Cz5hCkl5/DqG9zvwYiCCCexj+H91J26y5uhMY+Xk5/3K6LZTZyQbviTYieuLdGC3biNygPCmwp8tVRN4tcW58vyu//Zby8nJd3qqqKgRD0Lr/wAHKy8t15+/27tuHquoH8XBdHZ9t2ECS3R4Xo3QkvDClubmZ8vJy3SZpY2MjAYMq1djYlNC+5uZm2tr1b7hqamzi888/T7jDZYZDh7xs3bqV/9u0mZaWFtM81dXVlJeX69p9+PBhpBOIPUJZ59Bx/wu99gtCWSO7J1VKGs033IL8/bkkVexBrtiNtH8PQlsTatpQ1Mw8QgXnERwzlvCgIXAKXlojRCf0UhOrsbQPZPhTF4qecct7cRcbxmb3TszwYG/a2QVBd5QUlwaxYKE/3ey+kKerxZ/eX/+wn6BruLSX+az/m8rCaYF1YsCCBYtQFiwMYJdvAMDMRTO7xmHBwhlPqB0kavnBLmKmHV3EPIutobBwJuD/AbGgUv9INd91AAAAAElFTkSuQmCC"'

-- // ------------------------------------------------------------------------
-- // Registry declaration
-- // ------------------------------------------------------------------------
FuRegisterClass(
  -- >>> SCHNIPP::FUREGISTERCLASS
  (FC_PREFIX==nil and "ST_" or (FC_PREFIX=="" and "" or FC_PREFIX.."_"))..dctlfuse_name, CT_SourceTool, { REGS_Category = (FC_SUBMENU==nil and "Shadertoys (dev)" or (FC_SUBMENU=="" and "Fuses" or FC_SUBMENU))..(FC_CATEGORY==nil and "\\"..dctlfuse_category or ( FC_CATEGORY ~= "" and "\\"..FC_CATEGORY or "")), REGS_OpIconString = (FC_SCPREFIX==nil and "ST-" or (SC_SCPREFIX=="" and "" or FC_SCPREFIX.."-"))..shadertoy_id,
  REGS_OpDescription = "Shadertoy '"..shadertoy_name.."' (ID: "..shadertoy_id..") created by "..shadertoy_author.." and ported by "..dctlfuse_author..(shadertoy_license == "" and ". Copyright "..shadertoy_author.." (CC BY-NC-SA 3.0)" or ". "..shadertoy_license)..". This port is by no means meant to take advantage of anyone or to do anyone wrong : Contact us on Discord (https://discord.gg/75FUn4N4pv) and/or GitHub (https://github.com/nmbr73/Shadertoys) if you see your rights abused or your intellectual property violated by this work.",
  REG_Fuse_NoEdit = not(FC_DEVELOP==nil or FC_DEVELOP==true), REG_Fuse_NoReload = not(FC_DEVELOP==nil or FC_DEVELOP==true), REG_Fuse_TilePic = FC_TITLEPIC, REGS_Company = dctlfuse_company==nil and dctlfuse_author or dctlfuse_company, REGS_URL = dctlfuse_authorurl==nil and "https://nmbr73.github.io/Shadertoys/" or dctlfuse_authorurl,
  REG_Version = FC_VERSIONNO~=nil and FC_VERSIONNO or (dctlfuse_versionNo==nil and 1 or dctlfuse_versionNo),
  -- <<< SCHNAPP::FUREGISTERCLASS

  REG_NoObjMatCtrls      = true,
  REG_NoMotionBlurCtrls  = true,

  REG_Source_GlobalCtrls = false,
  REG_Source_SizeCtrls   = true,
  REG_Source_AspectCtrls = true,
  REG_Source_DepthCtrls  = true,
  REG_OpNoMask           = true,
  REG_Fuse_NoEdit        = not(dctlfuse_code),
  REG_Fuse_NoReload      = not(dctlfuse_code),

  REG_TimeVariant        = true,
  })

-- */

dctlfuse_versionNo =  dctlfuse_versionNo/10

-- // ------------------------------------------------------------------------
-- // DCTL kernel parameters
-- // ------------------------------------------------------------------------

ST_VelociboxParams =
[[
  float wr, wg, wb; // wall
  float kr, kg, kb; // kite
  float jr, jg, jb; // joist

  float freq;
  float exp;
  int   width;
  int   height;
  float iTime;
  int   compOrder;
]]



-- // ------------------------------------------------------------------------
-- DCTL kernel implementation
-- // ------------------------------------------------------------------------

ST_VelociboxKernel =
    [[


#if defined(DEVICE_IS_METAL)

  typedef float3x3 mat3;

  __DEVICE__ inline float  fract1f(float f) {return fract(f);}
  __DEVICE__ inline float3 reflect_f3(float3 I, float3 N) {return reflect(I,N);}

  #define swixy(V) (V).xy
  #define swizy(V) (V).zy
  #define swixz(V) (V).xy

  #define swixyy(V) (V).xyy
  #define swiyxy(V) (V).yxy
  #define swiyyx(V) (V).yyx

#else

  __DEVICE__ inline float  fract1f(float f) {return f - _floor(f);}
  __DEVICE__ inline float3 reflect_f3(float3 I, float3 N) {return I - 2.0f * dot(N, I) * N;}

  #define swixy(V) to_float2((V).x,(V).y)
  #define swizy(V) to_float2((V).z,(V).y)
  #define swixz(V) to_float2((V).x,(V).z)

  #define swixyy(V) to_float3((V).x,(V).y,(V).y)
  #define swiyxy(V) to_float3((V).y,(V).x,(V).y)
  #define swiyyx(V) to_float3((V).y,(V).y,(V).x)

#endif

#if defined(DEVICE_IS_CUDA)
__DEVICE__ float2 abs_f2(float2 a) {return (to_float2(_fabs(a.x), _fabs(a.y)));}
__DEVICE__ float3 abs_f3(float3 a) {return (to_float3(_fabs(a.x), _fabs(a.y),_fabs(a.z)));}
#define lpowf _powf
__DEVICE__ float3 pow_f3(float3 a, float3 b) {float3 r; r.x = lpowf(a.x,b.x); r.y = lpowf(a.y,b.y); r.z = lpowf(a.z,b.z); return r;}
#else
__DEVICE__ float2 abs_f2(float2 a) {return (_fabs(a));}
__DEVICE__ float3 abs_f3(float3 a) {return (_fabs(a));}
__DEVICE__ float3 pow_f3(float3 a, float3 b) {return (_powf(a,b));}
#endif

//############################################################################


__DEVICE__ float hash(float n)
{
  return fract1f(_sinf(n)*4578.543123f);
}

__DEVICE__ float2 path(float z, float iTime)
{
  float2 a = to_float2_s(0.0f);
  float2 b = to_float2(2.0f*_cosf(z*0.3f), 0.0f);

  if(iTime >= 0.0f && iTime <= 4.0f)
  {
      if(iTime <= 3.0f)
        return a;

      return _mix(a, b, fmod(iTime, 1.0f));

  } else
  {
    return b;
  }
}

__DEVICE__ float gTime(float iTime)
{
  float s;

  if (iTime <= 3.0f)
    s = 7.0f;
  else
    s = 15.0f;

  return s*iTime;
}

__DEVICE__ float len(float3 p, float l)
{
  p = pow_f3(abs_f3(p), to_float3_s(l));
  return _powf(p.x + p.y + p.z, 1.0f/l);
}

__DEVICE__ float dSegment(float3 p, float3 a, float3 b, float r)
{
  float3 pa = p - a;
  float3 ba = b - a;

  float h = _clampf(dot(pa, ba)/dot(ba, ba), 0.0f, 1.0f);

  return len(pa - ba*h, 8.0f) - r;
}

__DEVICE__ float dObstacles(float3 p, float diff, float iTime)
{
  // p.xy -= path(p.z);
  float2 xy= path(p.z,iTime); p.x -= xy.y; p.y-=p.x;
  float c = floor((p.z + diff)/(diff*2.0f));
  float rnd = hash(c);
  p.z = fmod(p.z + diff, diff*2.0f) - diff;

  // p.xy = p.yx;

  float o;
  if(iTime >= 0.0f && iTime <= 4.0)
  {
      if(iTime <= 3.0f)
        o = 4.0f;
      else
        o = _mix(4.0f, 0.0f, fmod(iTime, 1.0f));
  } else
  {
    o = 0.0;
  }

  p.y += path(p.z,iTime).x - o - (rnd < 0.5f ? 0.0f : 2.0f);

  return dSegment(p, to_float3(5.5f, 1.0f, 0.0f), to_float3(-5.5f, 1.0f, 0.0f), 0.25f);
}

__DEVICE__ float2 opU(float2 a, float2 b)
{
  return a.x < b.x ? a : b;
}

__DEVICE__ float3 boxPos(float3 p, float diff, float o, float iTime)
{
  float3 bp = to_float3(0.0f, -1.0f, gTime(iTime));
  float c = floor((p.z + diff)/diff*2.0f);
  float rnd = hash(c);

  bp.y = 1.8f*_cosf(_sinf(iTime + o)*3.0f + iTime + o);

  return bp;
}

__DEVICE__ float2 map(float3 p, float iTime)
{
  // the tunnel distance estimate was taken from Shane's (https://www.shadertoy.com/view/MlXSWX)
  float2 tun = abs_f2( swixy(p) - path(p.z, iTime))*to_float2(0.4f, 0.4f);
  float2 t = to_float2(1.0f - _fmaxf(tun.x, tun.y), 0.0f);

  float3 bp = boxPos(p, 2.5f, 0.0f, iTime);
  float3 bp2 = boxPos(p, 2.5f, 0.1f, iTime);
  bp2.z += 0.5f;

  float2 xy;
  xy=path(bp.z ,iTime); bp.x  += xy.x; bp.y  += xy.y; //bp.xy += path(bp.z);
  xy=path(bp2.z,iTime); bp2.x += xy.x; bp2.y += xy.y; //bp2.xy += path(bp2.z);

  float2 s = to_float2(_fminf(len(p - bp2, 1.0f) - 0.3f, dSegment(p, bp, bp2, 0.04f)), 1.0f);
  float2 o = to_float2(dObstacles(p, 2.5f,iTime), 2.0f);

  return opU(t, opU(s, o));
}

__DEVICE__ float2 intersect(float3 ro, float3 rd, float iTime)
{
  float td = 0.0f;
  float mid = -1.0f;
  float tmax = 50.0f;

  for(int i = 0; i < 256; i++)
  {
      float2 s = map(ro + rd*td, iTime);
      td += s.x;
      mid = s.y;
      if(_fabs(s.x) < 0.005f || td >= tmax)
        break;
  }

  if(td >= tmax)
    mid = -1.;

  return to_float2(td, mid);
}


__DEVICE__ float3 normal(float3 p, float iTime)
{
  float2 h = to_float2(0.001f, 0.0f);
  float3 n = to_float3(
      map(p + swixyy(h), iTime).x - map(p - swixyy(h), iTime).x,
      map(p + swiyxy(h), iTime).x - map(p - swiyxy(h), iTime).x,
      map(p + swiyyx(h), iTime).x - map(p - swiyyx(h), iTime).x
      );

  return normalize(n);
}

// iq's ambient occlusion.
__DEVICE__ float ao(float3 p, float3 n, float iTime)
{
  float r = 0.0f;
  float w = 1.0f;
  float d;

  for(float i = 1.0f; i < 6.6f; i+=1)
  {
      d = i/6.0f;
      r += w*(d - map(p + n*d, iTime).x);
      w *= 0.5f;
  }

  return 1.0f - _clampf(r, 0.0f, 1.0f);
}

__DEVICE__ float3 lighting(float3 p, float3 lp, float3 rd,float iTime)
{
  float3 l = lp - p;
  float dist = _fmaxf(length(l), 0.01f);
  float atten = _fminf(1.0f/(1.0f + dist*0.5f), 0.2f);
  l /= dist;

  float3 n = normal(p,iTime);
  float3 r = reflect_f3(-l, n);

  float dif = _clampf(dot(l, n), 0.0f, 1.0f);
  float spe = _powf(_clampf(dot(r, -rd), 0.0f, 1.0f), 8.0f);
  float fre = _powf(_clampf(1.0f + dot(n, rd), 0.0f, 1.0f), 2.0f);
  float dom = smoothstep(-1.0f, 1.0f, r.y);

  float3 lin = to_float3_s(0.2f);
  lin += 1.0f*dif*to_float3(1.0f, 0.97f, 0.85f);
  lin += 2.5f*spe*to_float3(1.0f, 0.97f, 0.85f)*dif;
  lin += 2.5f*fre*to_float3_s(1.0f);
  lin += 0.5f*dom*to_float3_s(1.0f);

  return lin*atten*ao(p, n,iTime);
}



__KERNEL__ void ST_VelociboxKernel(
    __CONSTANTREF__ ST_VelociboxParams*   params,
    __TEXTURE2D__                         iChannel0,
    __TEXTURE2D_WRITE__                   dst
    )
{

  DEFINE_KERNEL_ITERATORS_XY(x, y);

  if (x >= params->width || y >= params->height)
    return;

  // ----- Standard parameters
  float2 fragCoord   = to_float2(x, y);
  float2 iResolution = to_float2(params->width, params->height);
  float  iTime       = params->iTime * params->freq;
  float4 fragColor   = to_float4_s(0.0f);


  float3 wallColor  = to_float3(params->wr,params->wg,params->wb);
  float3 kiteColor  = to_float3(params->kr,params->kg,params->kb);
  float3 joistColor = to_float3(params->jr,params->jg,params->jb);
  float exp = params->exp;
  //##########################################################################

  float2 uv = -1.0f + 2.0f* ( swixy(fragCoord) / swixy(iResolution));
  uv.x *= iResolution.x/iResolution.y;

  float3 ro = to_float3(0.0f, 0.0f, -2.0f);
  float3 rd = normalize(to_float3_aw(uv, 2.0f));

  // -----
  // render
  ro.z += gTime(iTime) - 2.0f;

  float3 la = ro + to_float3(0.0f, 0.0f, 2.0f);

  { // camera

    float3 rl = to_float3(_sinf(iTime), _cosf(iTime), 0.0f);
    float3 f = normalize(la - ro);
    float3 r = cross(rl, f);
    float3 u = cross(f, r);

    float x = r.x*rd.x + u.x*rd.y + f.x*rd.z;
    float y = r.y*rd.x + u.y*rd.y + f.y*rd.z;
    float z = r.z*rd.x + u.z*rd.y + f.z*rd.z;

    rd = to_float3(x,y,z);

  }

  float3 lp1 = ro + to_float3(0.0f, 0.0f, 1.0f);
  float3 lp2 = ro + to_float3(0.0f, 0.0f, 5.0f);

  float2 xy;
  xy=path(la.z,iTime); la.x += xy.x; la.y+=xy.y; //la.xy += path(la.z);
  xy=path(ro.z,iTime); ro.x += xy.x; ro.y+=xy.y; //ro.xy += path(ro.z);

  float3 col = to_float3_s(0.5f);

  float2 i = intersect(ro, rd, iTime);
  float3 p = ro + rd*i.x;

  if (i.y == 0.0f)
  {
      float2 guv = normal(p,iTime).y == 0.0f ? swizy(p) : swixz(p) - path(p.z, iTime);
      col = 0.45f +
      wallColor // to_float3_s(1.0f) // wall
          *smoothstep(-0.05f, 0.05f,_fabs(fract1f(guv.x) - 0.5f)*2.0f)
          *smoothstep(-0.05f, 0.05f,_fabs(fract1f(guv.y) - 0.5f)*2.0f);
  }
  if (i.y == 1.0f) col = kiteColor*0.5f;  //to_float3(1.0f, 0.1f, 1.0f)*0.5f;  // kite
  if (i.y == 2.0f) col = joistColor;      //to_float3(0.0f, 0.3f, 1.0f);       // joist

  if(i.y != -1.0f) col *= lighting(p, lp1, rd,iTime)+lighting(p, lp2, rd,iTime);

  //col = pow_f3(col, to_float3_s(0.454545f));
  col = pow_f3(col, to_float3_s(exp));


  fragColor = to_float4_aw(col, 1.0f);


  //##########################################################################

  _tex2DVec4Write(dst, x, y, fragColor);

}

]]

-- /*

-- // ------------------------------------------------------------------------
-- // Create
-- // ------------------------------------------------------------------------

function Create()
--------------------------------------------
StandardShaderFuseTopControls();
--------------------------------------------
  ----- In/Out
  InImage1 = self:AddInput("Image", "Image", {
    LINKID_DataType = "Image",
    LINK_Main = 1,
	LINK_Visible = false,
    INP_Required = false
  })



  OutImage = self:AddOutput("Output", "Output", {
    LINKID_DataType = "Image",
    LINK_Main = 1,
  })

  ----- Inspector Panel Controls



  self:BeginControlNest("Kite Color", "KColor", false, {})
  InKR = self:AddInput("Red",   "KRed",   { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 1.0, ICS_Name = "Color", IC_ControlGroup = 1, IC_ControlID = 0, })
	InKG = self:AddInput("Green", "KGreen", { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 0.1,                     IC_ControlGroup = 1, IC_ControlID = 1, })
	InKB = self:AddInput("Blue",  "KBlue",  { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 1.0,                     IC_ControlGroup = 1, IC_ControlID = 2, })
  self:EndControlNest()

  self:BeginControlNest("Joist Color", "JColor", false, {})
  InJR = self:AddInput("Red",   "JRed",   { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 0.0, ICS_Name = "Color", IC_ControlGroup = 2, IC_ControlID = 0, })
	InJG = self:AddInput("Green", "JGreen", { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 0.3,                     IC_ControlGroup = 2, IC_ControlID = 1, })
	InJB = self:AddInput("Blue",  "JBlue",  { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 1.0,                     IC_ControlGroup = 2, IC_ControlID = 2, })
  self:EndControlNest()

  self:BeginControlNest("Wall Color", "WColor", false, {})
  InWR = self:AddInput("Red",   "WRed",   { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 1.0, ICS_Name = "Color", IC_ControlGroup = 3, IC_ControlID = 0, })
	InWG = self:AddInput("Green", "WGreen", { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 1.0,                     IC_ControlGroup = 3, IC_ControlID = 1, })
	InWB = self:AddInput("Blue",  "WBlue",  { LINKID_DataType = "Number", INPID_InputControl = "ColorControl", INP_MinScale = 0.0, INP_MaxScale = 1.0, INP_Default  = 1.0,                     IC_ControlGroup = 3, IC_ControlID = 2, })
  self:EndControlNest()


  InFreq = self:AddInput("Frequency", "Frequency", {
    LINKID_DataType    = "Number",
    INPID_InputControl = "SliderControl",
    INP_Default        = 1.0,
    INP_MinScale 	     = 0.0,
    INP_MaxScale 	     = 5.0,
  })

  InExp = self:AddInput("Darkness", "Exposure", {
    LINKID_DataType    = "Number",
    INPID_InputControl = "SliderControl",
    INP_Default        = 0.454545,
    INP_MinScale 	     = 0.01,
    INP_MaxScale 	     = 10.0,
  })


  ------ Information

  if FC_INFOBTN~=0 then
    InfoButton = self:AddInput("Info ...", "InfoButton", {
      LINKID_DataType = "Text",
      INPID_InputControl = "ButtonControl",
      INP_DoNotifyChanged = true,
      INP_External = false,
      IC_ControlPage = FC_INFOBTN,
    })
  end

--------------------------------------------
StandardShaderFuseBottomControls();
--------------------------------------------

end



-- // ------------------------------------------------------------------------
-- // Process
-- // ------------------------------------------------------------------------

function Process(req)

    local framerate = self.Comp:GetPrefs("Comp.FrameFormat.Rate") -- get the frame rate of the comp set in the preferences

    --This creates an image for us to work on.
    local imgattrs = {
        IMG_Document = self.Comp,
        { IMG_Channel = "Red", },
        { IMG_Channel = "Green", },
        { IMG_Channel = "Blue", },
        { IMG_Channel = "Alpha", },
        IMG_Width = Width,
        IMG_Height = Height,
        IMG_XScale = XAspect,
        IMG_YScale = YAspect,
        IMAT_OriginalWidth = realwidth,
        IMAT_OriginalHeight = realheight,
        IMG_Quality = not req:IsQuick(),
        IMG_MotionBlurQuality = not req:IsNoMotionBlur(),
        }

    if not req:IsStampOnly() then
        imgattrs.IMG_ProxyScale = 1
    end

    if SourceDepth ~= 0 then
        imgattrs.IMG_Depth = SourceDepth
    end


    -- Extern Texture or create a new one
    if (InImage1:GetValue(req) ~= nil) then
        src1 = InImage1:GetValue(req)
    else
        src1 = Image(imgattrs)
        local p = Pixel({R=0,G=0,B=0,A=0}) -- Initial black Image to avoid see random picture from Memory
        src1:Fill(p)
    end

    local dst = Image {IMG_Like = src1, IMG_DeferAlloc = true} -- create an Output Image


    -- Interface for running DCTL-Code

    node = DVIPComputeNode(req, "ST_VelociboxKernel", ST_VelociboxKernel, "ST_VelociboxParams", ST_VelociboxParams)


    if not pcall(function ()  -- necessary to avoid memory leakage
        params.freq         = InFreq:GetValue(req).Value
        params.exp         = InExp:GetValue(req).Value
        params.wr            = InWR:GetValue(req).Value
        params.wg            = InWG:GetValue(req).Value
        params.wb            = InWB:GetValue(req).Value
        params.kr            = InKR:GetValue(req).Value
        params.kg            = InKG:GetValue(req).Value
        params.kb            = InKB:GetValue(req).Value
        params.jr            = InJR:GetValue(req).Value
        params.jg            = InJG:GetValue(req).Value
        params.jb            = InJB:GetValue(req).Value
        params.compOrder    = src1:IsMask() and 1 or 15
        params.width        = src1.DataWindow:Width()
        params.height       = src1.DataWindow:Height()
        params.iTime        = req.Time / framerate
        node:SetParamBlock(params)
    end) then
        params              = node:GetParamBlock(ST_VelociboxParams)
        params.freq         = InFreq:GetValue(req).Value
        params.exp         = InExp:GetValue(req).Value
        params.wr            = InWR:GetValue(req).Value
        params.wg            = InWG:GetValue(req).Value
        params.wb            = InWB:GetValue(req).Value
        params.kr            = InKR:GetValue(req).Value
        params.kg            = InKG:GetValue(req).Value
        params.kb            = InKB:GetValue(req).Value
        params.jr            = InJR:GetValue(req).Value
        params.jg            = InJG:GetValue(req).Value
        params.jb            = InJB:GetValue(req).Value
        params.compOrder    = src1:IsMask() and 1 or 15
        params.width        = src1.DataWindow:Width()
        params.height       = src1.DataWindow:Height()
        params.iTime        = req.Time / framerate
        node:SetParamBlock(params)
    end


    node:AddSampler("RowSampler", TEX_FILTER_MODE_LINEAR,TEX_ADDRESS_MODE_BORDER, TEX_NORMALIZED_COORDS_TRUE)

    node:AddInput("iChannel0", src1)
    node:AddOutput("dst",      dst)

    local success = node:RunSession(req)
    if not success then
        dst = nil
        dump(node:GetErrorLog()) -- Errormessages from DCTL-Compiler
    end
    OutImage:Set(req, dst)
end



-- // ------------------------------------------------------------------------
-- // Callback
-- // ------------------------------------------------------------------------

function NotifyChanged(inp, param, time)
	if (param ~= nil) then
		if (param.Value == 1) then

		end
	end
end

-- */

-- /* ====================== DO NOT TOUCH OR APPEND ANY CODE HERE ===========================================

-- >>> SCHNIPP::SCHADERFUSECONTROLS.v2

function ShaderFuseControls_AuthImg() if FC_AUTHIMG ~= nil and FC_AUTHIMG ~= ""  then self:AddInput('<p align="'..FC_AUTHIMGALIGN..'"><a href="https://github.com/nmbr73/Shadertoys"><img '..FC_AUTHIMG..' /></a></p>',"DctlFuseLogoLabel", {ICS_ControlPage = (FC_AUTHIMGPOS == 2 and 'Info' or nil),IC_ControlPage=(FC_AUTHIMGPOS == 2 and 1 or FC_AUTHIMGPOS),LINKID_DataType = "Text",INPID_InputControl = "LabelControl",LBLC_MultiLine = true,IC_NoLabel = true,IC_NoReset = true,INP_External = false,INP_Passive = true,}) end end
function ShaderFuseControls_InfoBtn() self:AddInput(dctlfuse_name.." Info ...", "DctlFuseInfoButton", {ICS_ControlPage = (FC_INFOBTNPOS == 2 and 'Info' or nil),IC_ControlPage = (FC_INFOBTNPOS == 2 and 1 or FC_INFOBTNPOS),LINKID_DataType = "Text",INPID_InputControl = "ButtonControl",INP_DoNotifyChanged = false,INP_External = false,BTNCS_Execute = 'bmd.openurl("'..(dctlfuse_infourl~=nil and dctlfuse_infourl or 'https://nmbr73.github.io/Shadertoys/Shaders/'..dctlfuse_category..'/'..dctlfuse_name..'.html')..'")'}) end

function ShaderFuseControls_InfoTxt()
  self:AddInput('<p align="'..FC_INFOTXTALIGN ..'">'
    ..'Shadertoy <a href="https://www.shadertoy.com/view/'..shadertoy_id..'" style="color:white; text-decoration:none; ">'..shadertoy_name..'</a> by <a href="https://www.shadertoy.com/user/'..shadertoy_author..'" style="color:yellow; text-decoration:none; ">'..shadertoy_author..'</a><br />'
    ..'<span style="color:#ff6060; ">'..(shadertoy_license == "" and '&copy; '..shadertoy_author..' (CC BY-NC-SA 3.0)' or shadertoy_license)..'</span><br />'
    ..'DCTLified and DaFused by <a href="'..(dctlfuse_authorurl==nil and "https://nmbr73.github.io/Shadertoys/" or dctlfuse_authorurl)..'" style="color:yellow; text-decoration:none; ">'..dctlfuse_author..'</a><br />'
    ..'<span style="color:#4060ff; ">'..(FC_VERSIONNO==nil and '' or 'Version '..FC_VERSIONNO)..(FC_VERSIONNO~=nil and FC_VERSIONDATE~=nil and ' - ' or '')..(FC_VERSIONDATE==nil and '' or FC_VERSIONDATE)..(FC_VERSIONNO~=nil and FC_VERSIONDATE~=nil and '<br />' or '')..'</span></p>'
    ,"DctlFuseInfoLabel", {ICS_ControlPage =(FC_INFOTXTPOS == 2 and 'Info' or nil),IC_ControlPage = (FC_INFOTXTPOS == 2 and 1 or FC_INFOTXTPOS),LINKID_DataType="Text",INPID_InputControl="LabelControl",LBLC_MultiLine=true,IC_NoLabel=true,IC_NoReset=true,INP_External=false,INP_Passive=true})
end


function StandardShaderFuseTopControls()
  if FC_ITEMORDER         == nil then FC_ITEMORDER       = {ShaderFuseControls_InfoTxt,ShaderFuseControls_InfoBtn,ShaderFuseControls_AuthImg} end
  if FC_AUTHBASEDLAYOUT   == nil then FC_AUTHBASEDLAYOUT = true                  end
  if FC_INFOBTNPOS        == nil then FC_INFOBTNPOS      = 2                     end
  if FC_INFOTXTPOS        == nil then FC_INFOTXTPOS      = 2                     end
  if FC_INFOTXTALIGN      == nil then FC_INFOTXTALIGN    = 'center'              end
  if FC_AUTHIMGPOS        == nil then FC_AUTHIMGPOS      = 2                     end
  if FC_AUTHIMGALIGN      == nil then FC_AUTHIMGALIGN    = 'center'              end
  if FC_VERSIONNO         == nil then FC_VERSIONNO       = dctlfuse_versionNo    end
  if FC_VERSIONDATE       == nil then FC_VERSIONDATE     = dctlfuse_versionDate  end
  if FC_AUTHIMG           == nil then FC_AUTHIMG         = ''                    end
  if dctlfuse_authorlogo  ~= nil then FC_AUTHIMG         = dctlfuse_authorlogo   end

  if FC_AUTHBASEDLAYOUT then
    if dctlfuse_author=='JiPi' then
      FC_ITEMORDER       = {ShaderFuseControls_AuthImg,ShaderFuseControls_InfoBtn,ShaderFuseControls_InfoTxt}
      FC_INFOBTNPOS      = 1
      FC_AUTHIMGPOS      = -1
      FC_AUTHIMGALIGN    = 'center'
      FC_INFOTXTPOS      = 1
      FC_INFOTXTALIGN    = 'center'
    elseif dctlfuse_author=='nmbr73' then
      FC_ITEMORDER       = {ShaderFuseControls_InfoBtn,ShaderFuseControls_InfoTxt,ShaderFuseControls_AuthImg}
      FC_INFOBTNPOS      = 1
      FC_AUTHIMGPOS      = 1
      FC_AUTHIMGALIGN    = 'right'
      FC_INFOTXTPOS      = 1
      FC_INFOTXTALIGN    = 'left'
    end
  end

end

function StandardShaderFuseBottomControls()

  if FC_INFOBTNPOS == 1 or FC_INFOTXTPOS == 1 or FC_AUTHIMGPOS==1 then
    self:AddInput( '<br />',"DctlFuseSeparatorLabel",{LINKID_DataType="Text",INPID_InputControl="LabelControl",LBLC_MultiLine=true,IC_NoLabel=true,IC_NoReset=true,INP_External=false,INP_Passive=true})
  end

  FC_ITEMORDER[1]();
  FC_ITEMORDER[2]();
  FC_ITEMORDER[3]();

end
-- <<< SCHNAPP::SCHADERFUSECONTROLS
-- */