vec3d diag1 = {-1.0,-2.0,-3.0};
vec3f diag2 = {-1.0f,-2.0f,-3.0f};
mat3d@test1 = diag(diag1);
mat3f@test2 = diag(diag2);

vec4d diag3 = {-1.0,-2.0,-3.0,-4.0};
vec4f diag4 = {-1.0f,-2.0f,-3.0f,-4.0f};
mat4d@test3 = diag(diag3);
mat4f@test4 = diag(diag4);

mat3d diag5 = {-1.0,2.0,3.0, 4.0,-5.0,6.0, 7.0,8.0,-9.0};
mat3f diag6 = {-1.0f,2.0f,3.0f, 4.0f,-5.0f,6.0f, 7.0f,8.0f,-9.0f};
vec3d@test5 = diag(diag5);
vec3f@test6 = diag(diag6);

mat4d diag7 = {-1.0,2.0,3.0,4.0, 5.0,-6.0,7.0,8.0, 9.0,10.0,-11.0,12.0, 13.0,14.0,15.0,-16.0};
mat4f diag8 = {-1.0f,2.0f,3.0f,4.0f, 5.0f,-6.0f,7.0f,8.0f, 9.0f,10.0f,-11.0f,12.0f, 13.0f,14.0f,15.0f,-16.0f};

vec4d@test7 = diag(diag7);
vec4f@test8 = diag(diag8);
