Hi. I wanted to run blender script fourth_dimension.py but I'm getting error :KeyError: 'bpy_prop_collection[key]: key "Lamp" not found'. After changing name from "Lamp" to "Light" I'm getting another error. My guess is that's because I have a newer Blender version.
Can you provide version you used? Even better a requirements.txt file with all the libraries would be awesome. Thanks in advance!

To get this to compile I needed to move the $(CAIROFLAGS) to the end.

%.o: %.cpp $(DEPS)
$(CXX) $(CXXFLAGS) -c -o $@ $< $(CAIROFLAGS)

$(BINS): $(OBJ)
$(CXX) $(CXXFLAGS) -o $(BINS) $^ $(CAIROFLAGS)

I'm not sure what's causing it, I may have done something wrong, but on v2.78c the script is completely non-functional. The edges get placed somewhere far off and no rotation seems to occur at all. The only thing that does seem to happen is the light source slowly moves.

In MD.cpp line 319 it should be -= instead of += to make it equivalent to the commented out code above. It is totally my fault but I'll submit the issue to you anyways. I blame tiredness -_-
You also need to redo the j=0; after line 220 that ArnaudFiorini pointed out.

diff = distance(lens.origin, p);
if (diff > 0.001){
    index = 1.0 / diff;
} else {
    index = 1000;
}

Try these in the diff > 0.001 part of the if branch:

• index = exp(diff) / diff;
• index = 1.0 / sqrt(diff);
• index = 1.0 / cbrt(diff);
• index = erf(1.0 / diff);
• index = exp(1 / cbrt(diff)) * (1 / diff);
• index = 1.0 / ((diff) / (diff + 1));
• index = 1.0 / ((diff * (6.2 * diff + 0.5)) / (diff * (6.2 * diff + 1.7) + 0.06));
• index = sqrt(length(p));

Doesn't fit on a list:

index = 0.5 * (fabs(0.5 * diff) + sqrt(1 - (pow(fabs(fabs(diff) - 2) - 1, 2))) -
        1 / 112 * (3 * sqrt(33) - 7) * diff * diff + 3 * sqrt(1 - (diff / 7) * (diff / 7) ) - 3) *
        ((diff + 4) / fabs(diff + 4) - (diff - 4) / fabs(diff - 4)) - 3 * sqrt(1 - (diff / 7) * (diff / 7));
index = fabs(index);