void matrix_inverse(double T[4][4], double TIN[4][4]) {
  TIN[0][0] = T[0][0];
  TIN[0][1] = T[1][0];
  TIN[0][2] = T[2][0];
  TIN[0][3] = -(T[0][0] * T[0][3] + T[1][0] * T[1][3] + T[2][0] * T[2][3]);
  TIN[1][0] = T[0][1];
  TIN[1][1] = T[1][1];
  TIN[1][2] = T[2][1];
  TIN[1][3] = -(T[0][1] * T[0][3] + T[1][1] * T[1][3] + T[2][1] * T[2][3]);
  TIN[2][0] = T[0][2];
  TIN[2][1] = T[1][2];
  TIN[2][2] = T[2][2];
  TIN[2][3] = -(T[0][2] * T[0][3] + T[1][2] * T[1][3] + T[2][2] * T[2][3]);
  TIN[3][0] = 0;
  TIN[3][1] = 0;
  TIN[3][2] = 0;
  TIN[3][3] = 1;
}


x86 assembly (clang, -O0):
<matrix_inverse>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x38(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x58(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x18(%rax)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x20(%rax)
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x28(%rax)
mov    -0x8(%rbp),%rax
movsd  0x48(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x30(%rax)
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x38(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x48(%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x58(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x38(%rax)
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x40(%rax)
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x48(%rax)
mov    -0x8(%rbp),%rax
movsd  0x50(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x50(%rax)
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  0x38(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
movsd  0x50(%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x58(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x58(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x60(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x68(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x70(%rax)
mov    -0x10(%rbp),%rax
movsd  0xcf3(%rip),%xmm0
movsd  %xmm0,0x78(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrix_inverse>:
movsd  (%rdi),%xmm0
movsd  %xmm0,(%rsi)
movsd  0x20(%rdi),%xmm0
movsd  %xmm0,0x8(%rsi)
movsd  0x40(%rdi),%xmm0
movsd  %xmm0,0x10(%rsi)
movsd  (%rdi),%xmm0
movsd  0x20(%rdi),%xmm1
mulsd  0x38(%rdi),%xmm1
mulsd  0x18(%rdi),%xmm0
addsd  %xmm1,%xmm0
movsd  0x40(%rdi),%xmm1
mulsd  0x58(%rdi),%xmm1
addsd  %xmm0,%xmm1
movapd 0xeb7(%rip),%xmm0
xorpd  %xmm0,%xmm1
movlpd %xmm1,0x18(%rsi)
movsd  0x8(%rdi),%xmm1
movsd  %xmm1,0x20(%rsi)
movsd  0x28(%rdi),%xmm1
movsd  %xmm1,0x28(%rsi)
movsd  0x48(%rdi),%xmm1
movsd  %xmm1,0x30(%rsi)
movsd  0x8(%rdi),%xmm1
movsd  0x28(%rdi),%xmm2
mulsd  0x38(%rdi),%xmm2
mulsd  0x18(%rdi),%xmm1
movsd  0x48(%rdi),%xmm3
mulsd  0x58(%rdi),%xmm3
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm3
xorpd  %xmm0,%xmm3
movlpd %xmm3,0x38(%rsi)
movsd  0x10(%rdi),%xmm1
movsd  %xmm1,0x40(%rsi)
movsd  0x30(%rdi),%xmm1
movsd  %xmm1,0x48(%rsi)
movsd  0x50(%rdi),%xmm1
movsd  %xmm1,0x50(%rsi)
movsd  0x10(%rdi),%xmm1
movsd  0x30(%rdi),%xmm2
mulsd  0x38(%rdi),%xmm2
mulsd  0x18(%rdi),%xmm1
addsd  %xmm2,%xmm1
movsd  0x50(%rdi),%xmm2
mulsd  0x58(%rdi),%xmm2
addsd  %xmm1,%xmm2
xorpd  %xmm0,%xmm2
movlpd %xmm2,0x58(%rsi)
xorpd  %xmm0,%xmm0
movupd %xmm0,0x60(%rsi)
movq   $0x0,0x70(%rsi)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x78(%rsi)
ret

LLVM IR (-O0):
define void @matrix_inverse([4 x double]* %0, [4 x double]* %1) {
%3 = alloca [4 x double]*
%4 = alloca [4 x double]*
store [4 x double]* %0, [4 x double]** %3
store [4 x double]* %1, [4 x double]** %4
%5 = load [4 x double]*, [4 x double]** %3
%6 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0
%7 = getelementptr inbounds [4 x double], [4 x double]* %6, i64 0, i64 0
%8 = load double, double* %7
%9 = load [4 x double]*, [4 x double]** %4
%10 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 0
%11 = getelementptr inbounds [4 x double], [4 x double]* %10, i64 0, i64 0
store double %8, double* %11
%12 = load [4 x double]*, [4 x double]** %3
%13 = getelementptr inbounds [4 x double], [4 x double]* %12, i64 1
%14 = getelementptr inbounds [4 x double], [4 x double]* %13, i64 0, i64 0
%15 = load double, double* %14
%16 = load [4 x double]*, [4 x double]** %4
%17 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0
%18 = getelementptr inbounds [4 x double], [4 x double]* %17, i64 0, i64 1
store double %15, double* %18
%19 = load [4 x double]*, [4 x double]** %3
%20 = getelementptr inbounds [4 x double], [4 x double]* %19, i64 2
%21 = getelementptr inbounds [4 x double], [4 x double]* %20, i64 0, i64 0
%22 = load double, double* %21
%23 = load [4 x double]*, [4 x double]** %4
%24 = getelementptr inbounds [4 x double], [4 x double]* %23, i64 0
%25 = getelementptr inbounds [4 x double], [4 x double]* %24, i64 0, i64 2
store double %22, double* %25
%26 = load [4 x double]*, [4 x double]** %3
%27 = getelementptr inbounds [4 x double], [4 x double]* %26, i64 0
%28 = getelementptr inbounds [4 x double], [4 x double]* %27, i64 0, i64 0
%29 = load double, double* %28
%30 = load [4 x double]*, [4 x double]** %3
%31 = getelementptr inbounds [4 x double], [4 x double]* %30, i64 0
%32 = getelementptr inbounds [4 x double], [4 x double]* %31, i64 0, i64 3
%33 = load double, double* %32
%34 = load [4 x double]*, [4 x double]** %3
%35 = getelementptr inbounds [4 x double], [4 x double]* %34, i64 1
%36 = getelementptr inbounds [4 x double], [4 x double]* %35, i64 0, i64 0
%37 = load double, double* %36
%38 = load [4 x double]*, [4 x double]** %3
%39 = getelementptr inbounds [4 x double], [4 x double]* %38, i64 1
%40 = getelementptr inbounds [4 x double], [4 x double]* %39, i64 0, i64 3
%41 = load double, double* %40
%42 = fmul double %37, %41
%43 = call double @llvm.fmuladd.f64(double %29, double %33, double %42)
%44 = load [4 x double]*, [4 x double]** %3
%45 = getelementptr inbounds [4 x double], [4 x double]* %44, i64 2
%46 = getelementptr inbounds [4 x double], [4 x double]* %45, i64 0, i64 0
%47 = load double, double* %46
%48 = load [4 x double]*, [4 x double]** %3
%49 = getelementptr inbounds [4 x double], [4 x double]* %48, i64 2
%50 = getelementptr inbounds [4 x double], [4 x double]* %49, i64 0, i64 3
%51 = load double, double* %50
%52 = call double @llvm.fmuladd.f64(double %47, double %51, double %43)
%53 = fneg double %52
%54 = load [4 x double]*, [4 x double]** %4
%55 = getelementptr inbounds [4 x double], [4 x double]* %54, i64 0
%56 = getelementptr inbounds [4 x double], [4 x double]* %55, i64 0, i64 3
store double %53, double* %56
%57 = load [4 x double]*, [4 x double]** %3
%58 = getelementptr inbounds [4 x double], [4 x double]* %57, i64 0
%59 = getelementptr inbounds [4 x double], [4 x double]* %58, i64 0, i64 1
%60 = load double, double* %59
%61 = load [4 x double]*, [4 x double]** %4
%62 = getelementptr inbounds [4 x double], [4 x double]* %61, i64 1
%63 = getelementptr inbounds [4 x double], [4 x double]* %62, i64 0, i64 0
store double %60, double* %63
%64 = load [4 x double]*, [4 x double]** %3
%65 = getelementptr inbounds [4 x double], [4 x double]* %64, i64 1
%66 = getelementptr inbounds [4 x double], [4 x double]* %65, i64 0, i64 1
%67 = load double, double* %66
%68 = load [4 x double]*, [4 x double]** %4
%69 = getelementptr inbounds [4 x double], [4 x double]* %68, i64 1
%70 = getelementptr inbounds [4 x double], [4 x double]* %69, i64 0, i64 1
store double %67, double* %70
%71 = load [4 x double]*, [4 x double]** %3
%72 = getelementptr inbounds [4 x double], [4 x double]* %71, i64 2
%73 = getelementptr inbounds [4 x double], [4 x double]* %72, i64 0, i64 1
%74 = load double, double* %73
%75 = load [4 x double]*, [4 x double]** %4
%76 = getelementptr inbounds [4 x double], [4 x double]* %75, i64 1
%77 = getelementptr inbounds [4 x double], [4 x double]* %76, i64 0, i64 2
store double %74, double* %77
%78 = load [4 x double]*, [4 x double]** %3
%79 = getelementptr inbounds [4 x double], [4 x double]* %78, i64 0
%80 = getelementptr inbounds [4 x double], [4 x double]* %79, i64 0, i64 1
%81 = load double, double* %80
%82 = load [4 x double]*, [4 x double]** %3
%83 = getelementptr inbounds [4 x double], [4 x double]* %82, i64 0
%84 = getelementptr inbounds [4 x double], [4 x double]* %83, i64 0, i64 3
%85 = load double, double* %84
%86 = load [4 x double]*, [4 x double]** %3
%87 = getelementptr inbounds [4 x double], [4 x double]* %86, i64 1
%88 = getelementptr inbounds [4 x double], [4 x double]* %87, i64 0, i64 1
%89 = load double, double* %88
%90 = load [4 x double]*, [4 x double]** %3
%91 = getelementptr inbounds [4 x double], [4 x double]* %90, i64 1
%92 = getelementptr inbounds [4 x double], [4 x double]* %91, i64 0, i64 3
%93 = load double, double* %92
%94 = fmul double %89, %93
%95 = call double @llvm.fmuladd.f64(double %81, double %85, double %94)
%96 = load [4 x double]*, [4 x double]** %3
%97 = getelementptr inbounds [4 x double], [4 x double]* %96, i64 2
%98 = getelementptr inbounds [4 x double], [4 x double]* %97, i64 0, i64 1
%99 = load double, double* %98
%100 = load [4 x double]*, [4 x double]** %3
%101 = getelementptr inbounds [4 x double], [4 x double]* %100, i64 2
%102 = getelementptr inbounds [4 x double], [4 x double]* %101, i64 0, i64 3
%103 = load double, double* %102
%104 = call double @llvm.fmuladd.f64(double %99, double %103, double %95)
%105 = fneg double %104
%106 = load [4 x double]*, [4 x double]** %4
%107 = getelementptr inbounds [4 x double], [4 x double]* %106, i64 1
%108 = getelementptr inbounds [4 x double], [4 x double]* %107, i64 0, i64 3
store double %105, double* %108
%109 = load [4 x double]*, [4 x double]** %3
%110 = getelementptr inbounds [4 x double], [4 x double]* %109, i64 0
%111 = getelementptr inbounds [4 x double], [4 x double]* %110, i64 0, i64 2
%112 = load double, double* %111
%113 = load [4 x double]*, [4 x double]** %4
%114 = getelementptr inbounds [4 x double], [4 x double]* %113, i64 2
%115 = getelementptr inbounds [4 x double], [4 x double]* %114, i64 0, i64 0
store double %112, double* %115
%116 = load [4 x double]*, [4 x double]** %3
%117 = getelementptr inbounds [4 x double], [4 x double]* %116, i64 1
%118 = getelementptr inbounds [4 x double], [4 x double]* %117, i64 0, i64 2
%119 = load double, double* %118
%120 = load [4 x double]*, [4 x double]** %4
%121 = getelementptr inbounds [4 x double], [4 x double]* %120, i64 2
%122 = getelementptr inbounds [4 x double], [4 x double]* %121, i64 0, i64 1
store double %119, double* %122
%123 = load [4 x double]*, [4 x double]** %3
%124 = getelementptr inbounds [4 x double], [4 x double]* %123, i64 2
%125 = getelementptr inbounds [4 x double], [4 x double]* %124, i64 0, i64 2
%126 = load double, double* %125
%127 = load [4 x double]*, [4 x double]** %4
%128 = getelementptr inbounds [4 x double], [4 x double]* %127, i64 2
%129 = getelementptr inbounds [4 x double], [4 x double]* %128, i64 0, i64 2
store double %126, double* %129
%130 = load [4 x double]*, [4 x double]** %3
%131 = getelementptr inbounds [4 x double], [4 x double]* %130, i64 0
%132 = getelementptr inbounds [4 x double], [4 x double]* %131, i64 0, i64 2
%133 = load double, double* %132
%134 = load [4 x double]*, [4 x double]** %3
%135 = getelementptr inbounds [4 x double], [4 x double]* %134, i64 0
%136 = getelementptr inbounds [4 x double], [4 x double]* %135, i64 0, i64 3
%137 = load double, double* %136
%138 = load [4 x double]*, [4 x double]** %3
%139 = getelementptr inbounds [4 x double], [4 x double]* %138, i64 1
%140 = getelementptr inbounds [4 x double], [4 x double]* %139, i64 0, i64 2
%141 = load double, double* %140
%142 = load [4 x double]*, [4 x double]** %3
%143 = getelementptr inbounds [4 x double], [4 x double]* %142, i64 1
%144 = getelementptr inbounds [4 x double], [4 x double]* %143, i64 0, i64 3
%145 = load double, double* %144
%146 = fmul double %141, %145
%147 = call double @llvm.fmuladd.f64(double %133, double %137, double %146)
%148 = load [4 x double]*, [4 x double]** %3
%149 = getelementptr inbounds [4 x double], [4 x double]* %148, i64 2
%150 = getelementptr inbounds [4 x double], [4 x double]* %149, i64 0, i64 2
%151 = load double, double* %150
%152 = load [4 x double]*, [4 x double]** %3
%153 = getelementptr inbounds [4 x double], [4 x double]* %152, i64 2
%154 = getelementptr inbounds [4 x double], [4 x double]* %153, i64 0, i64 3
%155 = load double, double* %154
%156 = call double @llvm.fmuladd.f64(double %151, double %155, double %147)
%157 = fneg double %156
%158 = load [4 x double]*, [4 x double]** %4
%159 = getelementptr inbounds [4 x double], [4 x double]* %158, i64 2
%160 = getelementptr inbounds [4 x double], [4 x double]* %159, i64 0, i64 3
store double %157, double* %160
%161 = load [4 x double]*, [4 x double]** %4
%162 = getelementptr inbounds [4 x double], [4 x double]* %161, i64 3
%163 = getelementptr inbounds [4 x double], [4 x double]* %162, i64 0, i64 0
store double 0.000000e+00, double* %163
%164 = load [4 x double]*, [4 x double]** %4
%165 = getelementptr inbounds [4 x double], [4 x double]* %164, i64 3
%166 = getelementptr inbounds [4 x double], [4 x double]* %165, i64 0, i64 1
store double 0.000000e+00, double* %166
%167 = load [4 x double]*, [4 x double]** %4
%168 = getelementptr inbounds [4 x double], [4 x double]* %167, i64 3
%169 = getelementptr inbounds [4 x double], [4 x double]* %168, i64 0, i64 2
store double 0.000000e+00, double* %169
%170 = load [4 x double]*, [4 x double]** %4
%171 = getelementptr inbounds [4 x double], [4 x double]* %170, i64 3
%172 = getelementptr inbounds [4 x double], [4 x double]* %171, i64 0, i64 3
store double 1.000000e+00, double* %172
ret void
}

LLVM IR (-O2):
define void @matrix_inverse([4 x double]* nocapture %0, [4 x double]* nocapture writeonly %1) {
%3 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 0
%4 = load double, double* %3
%5 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 0, i64 0
store double %4, double* %5
%6 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 0
%7 = load double, double* %6
%8 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 0, i64 1
store double %7, double* %8
%9 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 0
%10 = load double, double* %9
%11 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 0, i64 2
store double %10, double* %11
%12 = load double, double* %3
%13 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 3
%14 = load double, double* %13
%15 = load double, double* %6
%16 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 3
%17 = load double, double* %16
%18 = fmul double %15, %17
%19 = tail call double @llvm.fmuladd.f64(double %12, double %14, double %18)
%20 = load double, double* %9
%21 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 3
%22 = load double, double* %21
%23 = tail call double @llvm.fmuladd.f64(double %20, double %22, double %19)
%24 = fneg double %23
%25 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 0, i64 3
store double %24, double* %25
%26 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 1
%27 = load double, double* %26
%28 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 1, i64 0
store double %27, double* %28
%29 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 1
%30 = load double, double* %29
%31 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 1, i64 1
store double %30, double* %31
%32 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 1
%33 = load double, double* %32
%34 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 1, i64 2
store double %33, double* %34
%35 = load double, double* %26
%36 = load double, double* %13
%37 = load double, double* %29
%38 = load double, double* %16
%39 = fmul double %37, %38
%40 = tail call double @llvm.fmuladd.f64(double %35, double %36, double %39)
%41 = load double, double* %32
%42 = load double, double* %21
%43 = tail call double @llvm.fmuladd.f64(double %41, double %42, double %40)
%44 = fneg double %43
%45 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 1, i64 3
store double %44, double* %45
%46 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 0, i64 2
%47 = load double, double* %46
%48 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 2, i64 0
store double %47, double* %48
%49 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 1, i64 2
%50 = load double, double* %49
%51 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 2, i64 1
store double %50, double* %51
%52 = getelementptr inbounds [4 x double], [4 x double]* %0, i64 2, i64 2
%53 = load double, double* %52
%54 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 2, i64 2
store double %53, double* %54
%55 = load double, double* %46
%56 = load double, double* %13
%57 = load double, double* %49
%58 = load double, double* %16
%59 = fmul double %57, %58
%60 = tail call double @llvm.fmuladd.f64(double %55, double %56, double %59)
%61 = load double, double* %52
%62 = load double, double* %21
%63 = tail call double @llvm.fmuladd.f64(double %61, double %62, double %60)
%64 = fneg double %63
%65 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 2, i64 3
store double %64, double* %65
%66 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 3, i64 0
%67 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 3, i64 3
%68 = bitcast double* %66 to i8*
store double 1.000000e+00, double* %67
ret void
}