The following line in `rotations3d.jl` appears to take up 90% or more of the basis generation time: ```Julia G = [ sum( coco_dot(CC[a][i], CC[b][i]) for i = 1:length(Mll) ) for a = 1:length(CC), b = 1:length(CC) ] ```