I tried to add other objects, or render after animation ends, but give weird results also.

You have a Z-scale of 0 for all objects. Setting to 1 for all solves the issue.

@Sergey Sharybin (sergey), there's still something weird going on here, rendering the .blend as is, results in the blank canvas as in the report, setting scale of the red outline to != 0 results in this:

I forget to mention that some of the objects are imported from Adobe Illustrator as svg.

The problem (or at least one problem) here is that transform_quick_inverse(), which is called from qbvh_instance_motion_push() during traversal, fails if the z-scale is 0 and returns a null matrix, which then produces weird values for the rest of the traversal. Using transform_inverse() solves that problem, but 1. will most likely slow down rendering and 2. it seems that there still are issues. Generally, non-invertible transforms seem to produce various problems with motion blur...

In the Blender transform code pseudoinverse_m4_m4 is used to handle this kind of degenerate matrix, but doing an SVD decomposition here would be too slow and complex.

The comments in transform_quick_inverse mention inverting the individual components of the matrix on construction. The translation and rotation components would work reliably, it's the degenerate 3x3 scale matrix that we need to handle then. We could add a small epsilon to the diagonal of this 3x3 matrix if it is degenerate (same trick as in transform_inverse).

It's still possible then that the interpolation between two non-degenerate matrices produces a degenerate matrix, but that would be rare, and it would only happen for a very specific time sample so I don't think you'd notice.

By the way, if the object has no scale or uniform scale (which is almost always the case, shearing is rare), then composing the interpolated matrix and its inverse could be optimized quite a bit, no matrix-matrix multiplication or matrix inversion necessary then.

Actually, inverting the individual components isn't necessary at all, you can just add an epsilon in transform_decompose if the scale matrix is degenerate. It would be a nice optimization but not needed to fix this bug.

There's not only motion blur issue here, there is also "regular" instance push/pop issue involved here. It is caused by difference in the scale factor applied on the intersection distance in instance push and pop functions. Here's a quick patch for this issue:

Not really sure it could be less creepy, but on the other hand number of math operations seems to be the same, so it should be all the same performance.

In order to solve motion blur issue the same type of changes is to be applied to motion pop function. In order to motion blur to work correctly fix to quick matrix inverse is still to be applied tho. Think something like this will work:

That det == 0.0 check could be done in transform_decompose ahead of time, but it won't really matter, so these fixes definitely look fine to me.

@Brecht Van Lommel (brecht), yeah we can do some preliminary checks there and maybe even move something out of the render kernel. Didn't have time to look into details -- crunchtime with PhD thesis pre-defense event at uni. Hopefully will have more time on thursday.