saturatemat(mat, sat) float mat[4][4];
float sat;
{
  float mmat[4][4];
  float a, b, c, d, e, f, g, h, i;
  float rwgt, gwgt, bwgt;
  rwgt = (0.3086);
  gwgt = (0.6094);
  bwgt = (0.0820);
  a = (1.0 - sat) * rwgt + sat;
  b = (1.0 - sat) * rwgt;
  c = (1.0 - sat) * rwgt;
  d = (1.0 - sat) * gwgt;
  e = (1.0 - sat) * gwgt + sat;
  f = (1.0 - sat) * gwgt;
  g = (1.0 - sat) * bwgt;
  h = (1.0 - sat) * bwgt;
  i = (1.0 - sat) * bwgt + sat;
  mmat[0][0] = a;
  mmat[0][1] = b;
  mmat[0][2] = c;
  mmat[0][3] = 0.0;
  mmat[1][0] = d;
  mmat[1][1] = e;
  mmat[1][2] = f;
  mmat[1][3] = 0.0;
  mmat[2][0] = g;
  mmat[2][1] = h;
  mmat[2][2] = i;
  mmat[2][3] = 0.0;
  mmat[3][0] = 0.0;
  mmat[3][1] = 0.0;
  mmat[3][2] = 0.0;
  mmat[3][3] = 1.0;
  matrixmult(mmat, mat, mat);
}


x86 assembly (clang, -O0):
<saturatemat>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
cvtsd2ss %xmm0,%xmm0
mov    %rdi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
movss  0xedc(%rip),%xmm0
movss  %xmm0,-0x88(%rbp)
ff
movss  0xec8(%rip),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
movss  0xeb4(%rip),%xmm0
movss  %xmm0,-0x90(%rbp)
ff
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xea7(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x88(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movss  -0x14(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe6c(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x88(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x68(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe3e(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x88(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x6c(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe10(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x8c(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x70(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xde2(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x8c(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movss  -0x14(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x74(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xda7(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x8c(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x78(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd79(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x90(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x7c(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd4b(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x90(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x80(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd1d(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x90(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm2
movss  -0x14(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x84(%rbp)
ff
movss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x68(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  -0x6c(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x70(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x74(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x78(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x7c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x80(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x84(%rbp),%xmm0
ff
movss  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  0xc4b(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
lea    -0x60(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    $0x0,%al
call   1030 <matrixmult@plt>
mov    -0x4(%rbp),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<saturatemat>:
sub    $0x48,%rsp
cvtsd2ss %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xedc(%rip),%xmm1
subsd  %xmm0,%xmm1
movsd  0xed8(%rip),%xmm2
mulsd  %xmm1,%xmm2
cvtsd2ss %xmm2,%xmm3
addsd  %xmm0,%xmm2
cvtsd2ss %xmm2,%xmm2
movsd  0xec8(%rip),%xmm4
mulsd  %xmm1,%xmm4
cvtsd2ss %xmm4,%xmm5
addsd  %xmm0,%xmm4
cvtsd2ss %xmm4,%xmm4
mulsd  0xeb8(%rip),%xmm1
mov    %rdi,%rdx
cvtsd2ss %xmm1,%xmm6
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm2,(%rsp)
movss  %xmm3,0x4(%rsp)
movss  %xmm3,0x8(%rsp)
movl   $0x0,0xc(%rsp)
movss  %xmm5,0x10(%rsp)
movss  %xmm4,0x14(%rsp)
movss  %xmm5,0x18(%rsp)
movl   $0x0,0x1c(%rsp)
movss  %xmm6,0x20(%rsp)
movss  %xmm6,0x24(%rsp)
movss  %xmm0,0x28(%rsp)
xorps  %xmm0,%xmm0
movups %xmm0,0x2c(%rsp)
movl   $0x3f800000,0x3c(%rsp)
3f
mov    %rsp,%rdi
mov    %rdx,%rsi
xor    %eax,%eax
call   1030 <matrixmult@plt>
add    $0x48,%rsp
ret

LLVM IR (-O0):
define i32 @saturatemat([4 x float]* %0, double %1) {
%3 = alloca i32
%4 = alloca [4 x float]*
%5 = alloca float
%6 = alloca [4 x [4 x float]]
%7 = alloca float
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = fptrunc double %1 to float
store [4 x float]* %0, [4 x float]** %4
store float %19, float* %5
store float 0x3FD3C01A40000000, float* %16
store float 0x3FE3803460000000, float* %17
store float 0x3FB4FDF3C0000000, float* %18
%20 = load float, float* %5
%21 = fpext float %20 to double
%22 = fsub double 1.000000e+00, %21
%23 = load float, float* %16
%24 = fpext float %23 to double
%25 = load float, float* %5
%26 = fpext float %25 to double
%27 = call double @llvm.fmuladd.f64(double %22, double %24, double %26)
%28 = fptrunc double %27 to float
store float %28, float* %7
%29 = load float, float* %5
%30 = fpext float %29 to double
%31 = fsub double 1.000000e+00, %30
%32 = load float, float* %16
%33 = fpext float %32 to double
%34 = fmul double %31, %33
%35 = fptrunc double %34 to float
store float %35, float* %8
%36 = load float, float* %5
%37 = fpext float %36 to double
%38 = fsub double 1.000000e+00, %37
%39 = load float, float* %16
%40 = fpext float %39 to double
%41 = fmul double %38, %40
%42 = fptrunc double %41 to float
store float %42, float* %9
%43 = load float, float* %5
%44 = fpext float %43 to double
%45 = fsub double 1.000000e+00, %44
%46 = load float, float* %17
%47 = fpext float %46 to double
%48 = fmul double %45, %47
%49 = fptrunc double %48 to float
store float %49, float* %10
%50 = load float, float* %5
%51 = fpext float %50 to double
%52 = fsub double 1.000000e+00, %51
%53 = load float, float* %17
%54 = fpext float %53 to double
%55 = load float, float* %5
%56 = fpext float %55 to double
%57 = call double @llvm.fmuladd.f64(double %52, double %54, double %56)
%58 = fptrunc double %57 to float
store float %58, float* %11
%59 = load float, float* %5
%60 = fpext float %59 to double
%61 = fsub double 1.000000e+00, %60
%62 = load float, float* %17
%63 = fpext float %62 to double
%64 = fmul double %61, %63
%65 = fptrunc double %64 to float
store float %65, float* %12
%66 = load float, float* %5
%67 = fpext float %66 to double
%68 = fsub double 1.000000e+00, %67
%69 = load float, float* %18
%70 = fpext float %69 to double
%71 = fmul double %68, %70
%72 = fptrunc double %71 to float
store float %72, float* %13
%73 = load float, float* %5
%74 = fpext float %73 to double
%75 = fsub double 1.000000e+00, %74
%76 = load float, float* %18
%77 = fpext float %76 to double
%78 = fmul double %75, %77
%79 = fptrunc double %78 to float
store float %79, float* %14
%80 = load float, float* %5
%81 = fpext float %80 to double
%82 = fsub double 1.000000e+00, %81
%83 = load float, float* %18
%84 = fpext float %83 to double
%85 = load float, float* %5
%86 = fpext float %85 to double
%87 = call double @llvm.fmuladd.f64(double %82, double %84, double %86)
%88 = fptrunc double %87 to float
store float %88, float* %15
%89 = load float, float* %7
%90 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 0
%91 = getelementptr inbounds [4 x float], [4 x float]* %90, i64 0, i64 0
store float %89, float* %91
%92 = load float, float* %8
%93 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 0
%94 = getelementptr inbounds [4 x float], [4 x float]* %93, i64 0, i64 1
store float %92, float* %94
%95 = load float, float* %9
%96 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 0
%97 = getelementptr inbounds [4 x float], [4 x float]* %96, i64 0, i64 2
store float %95, float* %97
%98 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 0
%99 = getelementptr inbounds [4 x float], [4 x float]* %98, i64 0, i64 3
store float 0.000000e+00, float* %99
%100 = load float, float* %10
%101 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 1
%102 = getelementptr inbounds [4 x float], [4 x float]* %101, i64 0, i64 0
store float %100, float* %102
%103 = load float, float* %11
%104 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 1
%105 = getelementptr inbounds [4 x float], [4 x float]* %104, i64 0, i64 1
store float %103, float* %105
%106 = load float, float* %12
%107 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 1
%108 = getelementptr inbounds [4 x float], [4 x float]* %107, i64 0, i64 2
store float %106, float* %108
%109 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 1
%110 = getelementptr inbounds [4 x float], [4 x float]* %109, i64 0, i64 3
store float 0.000000e+00, float* %110
%111 = load float, float* %13
%112 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 2
%113 = getelementptr inbounds [4 x float], [4 x float]* %112, i64 0, i64 0
store float %111, float* %113
%114 = load float, float* %14
%115 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 2
%116 = getelementptr inbounds [4 x float], [4 x float]* %115, i64 0, i64 1
store float %114, float* %116
%117 = load float, float* %15
%118 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 2
%119 = getelementptr inbounds [4 x float], [4 x float]* %118, i64 0, i64 2
store float %117, float* %119
%120 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 2
%121 = getelementptr inbounds [4 x float], [4 x float]* %120, i64 0, i64 3
store float 0.000000e+00, float* %121
%122 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 3
%123 = getelementptr inbounds [4 x float], [4 x float]* %122, i64 0, i64 0
store float 0.000000e+00, float* %123
%124 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 3
%125 = getelementptr inbounds [4 x float], [4 x float]* %124, i64 0, i64 1
store float 0.000000e+00, float* %125
%126 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 3
%127 = getelementptr inbounds [4 x float], [4 x float]* %126, i64 0, i64 2
store float 0.000000e+00, float* %127
%128 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 3
%129 = getelementptr inbounds [4 x float], [4 x float]* %128, i64 0, i64 3
store float 1.000000e+00, float* %129
%130 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %6, i64 0, i64 0
%131 = load [4 x float]*, [4 x float]** %4
%132 = load [4 x float]*, [4 x float]** %4
%133 = call i32 ([4 x float]*, [4 x float]*, [4 x float]*, ...) bitcast (i32 (...)* @matrixmult to i32 ([4 x float]*, [4 x float]*, [4 x float]*, ...)*)([4 x float]* %130, [4 x float]* %131, [4 x float]* %132)
%134 = load i32, i32* %3
ret i32 %134
}

LLVM IR (-O2):
define i32 @saturatemat([4 x float]* %0, double %1) {
%3 = alloca [4 x [4 x float]]
%4 = fptrunc double %1 to float
%5 = bitcast [4 x [4 x float]]* %3 to i8*
%6 = fpext float %4 to double
%7 = fsub double 1.000000e+00, %6
%8 = tail call double @llvm.fmuladd.f64(double %7, double 0x3FD3C01A40000000, double %6)
%9 = fptrunc double %8 to float
%10 = fmul double %7, 0x3FD3C01A40000000
%11 = fptrunc double %10 to float
%12 = fmul double %7, 0x3FE3803460000000
%13 = fptrunc double %12 to float
%14 = tail call double @llvm.fmuladd.f64(double %7, double 0x3FE3803460000000, double %6)
%15 = fptrunc double %14 to float
%16 = fmul double %7, 0x3FB4FDF3C0000000
%17 = fptrunc double %16 to float
%18 = tail call double @llvm.fmuladd.f64(double %7, double 0x3FB4FDF3C0000000, double %6)
%19 = fptrunc double %18 to float
%20 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 0
%21 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 0, i64 0
store float %9, float* %21
%22 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 0, i64 1
store float %11, float* %22
%23 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 0, i64 2
store float %11, float* %23
%24 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 0, i64 3
store float 0.000000e+00, float* %24
%25 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 1, i64 0
store float %13, float* %25
%26 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 1, i64 1
store float %15, float* %26
%27 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 1, i64 2
store float %13, float* %27
%28 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 1, i64 3
store float 0.000000e+00, float* %28
%29 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 2, i64 0
store float %17, float* %29
%30 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 2, i64 1
store float %17, float* %30
%31 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 2, i64 2
store float %19, float* %31
%32 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 2, i64 3
%33 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %3, i64 0, i64 3, i64 3
%34 = bitcast float* %32 to i8*
store float 1.000000e+00, float* %33
%35 = call i32 ([4 x float]*, [4 x float]*, [4 x float]*, ...) bitcast (i32 (...)* @matrixmult to i32 ([4 x float]*, [4 x float]*, [4 x float]*, ...)*)([4 x float]* %20, [4 x float]* %0, [4 x float]* %0)
ret i32 undef
}