レイマーチングの法線をdFdxとdFdyで導出する ref: https://qiita.com/gam0022/items/1878150981494fd66abe

dFdyによる法線の計算

vec3 getNormal(vec3 p) {
    vec3 dx = dFdx(p);
    vec3 dy = dFdy(p);
    return normalize(cross(dx, dy));
}

距離関数の勾配による法線の計算

const float EPS = 0.001;

vec3 getNormal(vec3 p) {
    return normalize(vec3(
        sceneDist(p + vec3(EPS, 0.0, 0.0)) - sceneDist(p + vec3(-EPS,  0.0,  0.0)),
        sceneDist(p + vec3(0.0, EPS, 0.0)) - sceneDist(p + vec3( 0.0, -EPS,  0.0)),
        sceneDist(p + vec3(0.0, 0.0, EPS)) - sceneDist(p + vec3( 0.0,  0.0, -EPS))
    ));
}