; LLVM IR from InsightSoftwareConsortium/ITK :: Modules/ThirdParty/MINC/src/libminc/nifti (nifti1_io)

; Unoptimized IR (clang -disable-llvm-passes):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

%struct.mat44 = type { [4 x [4 x float]] }

; Function Attrs: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.start.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.end.p0(i64 immarg, ptr nocapture) #0

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #1

; Function Attrs: nounwind uwtable
define dso_local void @nifti_mat44_inverse(ptr dead_on_unwind noalias writable sret(%struct.mat44) align 4 %0, ptr noundef byval(%struct.mat44) align 8 %1) #2 {
  %3 = alloca double, align 8
  %4 = alloca double, align 8
  %5 = alloca double, align 8
  %6 = alloca double, align 8
  %7 = alloca double, align 8
  %8 = alloca double, align 8
  %9 = alloca double, align 8
  %10 = alloca double, align 8
  %11 = alloca double, align 8
  %12 = alloca double, align 8
  %13 = alloca double, align 8
  %14 = alloca double, align 8
  %15 = alloca double, align 8
  call void @llvm.lifetime.start.p0(i64 8, ptr %3) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %4) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %5) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %6) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %7) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %8) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %9) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %11) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %13) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #3
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #3
  %16 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %17 = getelementptr inbounds [4 x [4 x float]], ptr %16, i64 0, i64 0
  %18 = getelementptr inbounds [4 x float], ptr %17, i64 0, i64 0
  %19 = load float, ptr %18, align 8, !tbaa !5
  %20 = fpext float %19 to double
  store double %20, ptr %3, align 8, !tbaa !9
  %21 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %22 = getelementptr inbounds [4 x [4 x float]], ptr %21, i64 0, i64 0
  %23 = getelementptr inbounds [4 x float], ptr %22, i64 0, i64 1
  %24 = load float, ptr %23, align 4, !tbaa !5
  %25 = fpext float %24 to double
  store double %25, ptr %4, align 8, !tbaa !9
  %26 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %27 = getelementptr inbounds [4 x [4 x float]], ptr %26, i64 0, i64 0
  %28 = getelementptr inbounds [4 x float], ptr %27, i64 0, i64 2
  %29 = load float, ptr %28, align 8, !tbaa !5
  %30 = fpext float %29 to double
  store double %30, ptr %5, align 8, !tbaa !9
  %31 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %32 = getelementptr inbounds [4 x [4 x float]], ptr %31, i64 0, i64 1
  %33 = getelementptr inbounds [4 x float], ptr %32, i64 0, i64 0
  %34 = load float, ptr %33, align 8, !tbaa !5
  %35 = fpext float %34 to double
  store double %35, ptr %6, align 8, !tbaa !9
  %36 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %37 = getelementptr inbounds [4 x [4 x float]], ptr %36, i64 0, i64 1
  %38 = getelementptr inbounds [4 x float], ptr %37, i64 0, i64 1
  %39 = load float, ptr %38, align 4, !tbaa !5
  %40 = fpext float %39 to double
  store double %40, ptr %7, align 8, !tbaa !9
  %41 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %42 = getelementptr inbounds [4 x [4 x float]], ptr %41, i64 0, i64 1
  %43 = getelementptr inbounds [4 x float], ptr %42, i64 0, i64 2
  %44 = load float, ptr %43, align 8, !tbaa !5
  %45 = fpext float %44 to double
  store double %45, ptr %8, align 8, !tbaa !9
  %46 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %47 = getelementptr inbounds [4 x [4 x float]], ptr %46, i64 0, i64 2
  %48 = getelementptr inbounds [4 x float], ptr %47, i64 0, i64 0
  %49 = load float, ptr %48, align 8, !tbaa !5
  %50 = fpext float %49 to double
  store double %50, ptr %9, align 8, !tbaa !9
  %51 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %52 = getelementptr inbounds [4 x [4 x float]], ptr %51, i64 0, i64 2
  %53 = getelementptr inbounds [4 x float], ptr %52, i64 0, i64 1
  %54 = load float, ptr %53, align 4, !tbaa !5
  %55 = fpext float %54 to double
  store double %55, ptr %10, align 8, !tbaa !9
  %56 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %57 = getelementptr inbounds [4 x [4 x float]], ptr %56, i64 0, i64 2
  %58 = getelementptr inbounds [4 x float], ptr %57, i64 0, i64 2
  %59 = load float, ptr %58, align 8, !tbaa !5
  %60 = fpext float %59 to double
  store double %60, ptr %11, align 8, !tbaa !9
  %61 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %62 = getelementptr inbounds [4 x [4 x float]], ptr %61, i64 0, i64 0
  %63 = getelementptr inbounds [4 x float], ptr %62, i64 0, i64 3
  %64 = load float, ptr %63, align 4, !tbaa !5
  %65 = fpext float %64 to double
  store double %65, ptr %12, align 8, !tbaa !9
  %66 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %67 = getelementptr inbounds [4 x [4 x float]], ptr %66, i64 0, i64 1
  %68 = getelementptr inbounds [4 x float], ptr %67, i64 0, i64 3
  %69 = load float, ptr %68, align 4, !tbaa !5
  %70 = fpext float %69 to double
  store double %70, ptr %13, align 8, !tbaa !9
  %71 = getelementptr inbounds %struct.mat44, ptr %1, i32 0, i32 0
  %72 = getelementptr inbounds [4 x [4 x float]], ptr %71, i64 0, i64 2
  %73 = getelementptr inbounds [4 x float], ptr %72, i64 0, i64 3
  %74 = load float, ptr %73, align 4, !tbaa !5
  %75 = fpext float %74 to double
  store double %75, ptr %14, align 8, !tbaa !9
  %76 = load double, ptr %3, align 8, !tbaa !9
  %77 = load double, ptr %7, align 8, !tbaa !9
  %78 = fmul double %76, %77
  %79 = load double, ptr %11, align 8, !tbaa !9
  %80 = load double, ptr %3, align 8, !tbaa !9
  %81 = load double, ptr %10, align 8, !tbaa !9
  %82 = fmul double %80, %81
  %83 = load double, ptr %8, align 8, !tbaa !9
  %84 = fmul double %82, %83
  %85 = fneg double %84
  %86 = call double @llvm.fmuladd.f64(double %78, double %79, double %85)
  %87 = load double, ptr %6, align 8, !tbaa !9
  %88 = load double, ptr %4, align 8, !tbaa !9
  %89 = fmul double %87, %88
  %90 = load double, ptr %11, align 8, !tbaa !9
  %91 = fneg double %89
  %92 = call double @llvm.fmuladd.f64(double %91, double %90, double %86)
  %93 = load double, ptr %6, align 8, !tbaa !9
  %94 = load double, ptr %10, align 8, !tbaa !9
  %95 = fmul double %93, %94
  %96 = load double, ptr %5, align 8, !tbaa !9
  %97 = call double @llvm.fmuladd.f64(double %95, double %96, double %92)
  %98 = load double, ptr %9, align 8, !tbaa !9
  %99 = load double, ptr %4, align 8, !tbaa !9
  %100 = fmul double %98, %99
  %101 = load double, ptr %8, align 8, !tbaa !9
  %102 = call double @llvm.fmuladd.f64(double %100, double %101, double %97)
  %103 = load double, ptr %9, align 8, !tbaa !9
  %104 = load double, ptr %7, align 8, !tbaa !9
  %105 = fmul double %103, %104
  %106 = load double, ptr %5, align 8, !tbaa !9
  %107 = fneg double %105
  %108 = call double @llvm.fmuladd.f64(double %107, double %106, double %102)
  store double %108, ptr %15, align 8, !tbaa !9
  %109 = load double, ptr %15, align 8, !tbaa !9
  %110 = fpext double %109 to x86_fp80
  %111 = fcmp une x86_fp80 %110, 0xK00000000000000000000
  br i1 %111, label %112, label %117

112:                                              ; preds = %2
  %113 = load double, ptr %15, align 8, !tbaa !9
  %114 = fpext double %113 to x86_fp80
  %115 = fdiv x86_fp80 0xK3FFF8000000000000000, %114
  %116 = fptrunc x86_fp80 %115 to double
  store double %116, ptr %15, align 8, !tbaa !9
  br label %117

117:                                              ; preds = %112, %2
  %118 = load double, ptr %15, align 8, !tbaa !9
  %119 = load double, ptr %7, align 8, !tbaa !9
  %120 = load double, ptr %11, align 8, !tbaa !9
  %121 = load double, ptr %10, align 8, !tbaa !9
  %122 = load double, ptr %8, align 8, !tbaa !9
  %123 = fmul double %121, %122
  %124 = fneg double %123
  %125 = call double @llvm.fmuladd.f64(double %119, double %120, double %124)
  %126 = fmul double %118, %125
  %127 = fptrunc double %126 to float
  %128 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %129 = getelementptr inbounds [4 x [4 x float]], ptr %128, i64 0, i64 0
  %130 = getelementptr inbounds [4 x float], ptr %129, i64 0, i64 0
  store float %127, ptr %130, align 4, !tbaa !5
  %131 = load double, ptr %15, align 8, !tbaa !9
  %132 = load double, ptr %4, align 8, !tbaa !9
  %133 = fneg double %132
  %134 = load double, ptr %11, align 8, !tbaa !9
  %135 = load double, ptr %10, align 8, !tbaa !9
  %136 = load double, ptr %5, align 8, !tbaa !9
  %137 = fmul double %135, %136
  %138 = call double @llvm.fmuladd.f64(double %133, double %134, double %137)
  %139 = fmul double %131, %138
  %140 = fptrunc double %139 to float
  %141 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %142 = getelementptr inbounds [4 x [4 x float]], ptr %141, i64 0, i64 0
  %143 = getelementptr inbounds [4 x float], ptr %142, i64 0, i64 1
  store float %140, ptr %143, align 4, !tbaa !5
  %144 = load double, ptr %15, align 8, !tbaa !9
  %145 = load double, ptr %4, align 8, !tbaa !9
  %146 = load double, ptr %8, align 8, !tbaa !9
  %147 = load double, ptr %7, align 8, !tbaa !9
  %148 = load double, ptr %5, align 8, !tbaa !9
  %149 = fmul double %147, %148
  %150 = fneg double %149
  %151 = call double @llvm.fmuladd.f64(double %145, double %146, double %150)
  %152 = fmul double %144, %151
  %153 = fptrunc double %152 to float
  %154 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %155 = getelementptr inbounds [4 x [4 x float]], ptr %154, i64 0, i64 0
  %156 = getelementptr inbounds [4 x float], ptr %155, i64 0, i64 2
  store float %153, ptr %156, align 4, !tbaa !5
  %157 = load double, ptr %15, align 8, !tbaa !9
  %158 = load double, ptr %4, align 8, !tbaa !9
  %159 = fneg double %158
  %160 = load double, ptr %8, align 8, !tbaa !9
  %161 = fmul double %159, %160
  %162 = load double, ptr %14, align 8, !tbaa !9
  %163 = load double, ptr %4, align 8, !tbaa !9
  %164 = load double, ptr %13, align 8, !tbaa !9
  %165 = fmul double %163, %164
  %166 = load double, ptr %11, align 8, !tbaa !9
  %167 = fmul double %165, %166
  %168 = call double @llvm.fmuladd.f64(double %161, double %162, double %167)
  %169 = load double, ptr %7, align 8, !tbaa !9
  %170 = load double, ptr %5, align 8, !tbaa !9
  %171 = fmul double %169, %170
  %172 = load double, ptr %14, align 8, !tbaa !9
  %173 = call double @llvm.fmuladd.f64(double %171, double %172, double %168)
  %174 = load double, ptr %7, align 8, !tbaa !9
  %175 = load double, ptr %12, align 8, !tbaa !9
  %176 = fmul double %174, %175
  %177 = load double, ptr %11, align 8, !tbaa !9
  %178 = fneg double %176
  %179 = call double @llvm.fmuladd.f64(double %178, double %177, double %173)
  %180 = load double, ptr %10, align 8, !tbaa !9
  %181 = load double, ptr %5, align 8, !tbaa !9
  %182 = fmul double %180, %181
  %183 = load double, ptr %13, align 8, !tbaa !9
  %184 = fneg double %182
  %185 = call double @llvm.fmuladd.f64(double %184, double %183, double %179)
  %186 = load double, ptr %10, align 8, !tbaa !9
  %187 = load double, ptr %12, align 8, !tbaa !9
  %188 = fmul double %186, %187
  %189 = load double, ptr %8, align 8, !tbaa !9
  %190 = call double @llvm.fmuladd.f64(double %188, double %189, double %185)
  %191 = fmul double %157, %190
  %192 = fptrunc double %191 to float
  %193 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %194 = getelementptr inbounds [4 x [4 x float]], ptr %193, i64 0, i64 0
  %195 = getelementptr inbounds [4 x float], ptr %194, i64 0, i64 3
  store float %192, ptr %195, align 4, !tbaa !5
  %196 = load double, ptr %15, align 8, !tbaa !9
  %197 = load double, ptr %6, align 8, !tbaa !9
  %198 = fneg double %197
  %199 = load double, ptr %11, align 8, !tbaa !9
  %200 = load double, ptr %9, align 8, !tbaa !9
  %201 = load double, ptr %8, align 8, !tbaa !9
  %202 = fmul double %200, %201
  %203 = call double @llvm.fmuladd.f64(double %198, double %199, double %202)
  %204 = fmul double %196, %203
  %205 = fptrunc double %204 to float
  %206 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %207 = getelementptr inbounds [4 x [4 x float]], ptr %206, i64 0, i64 1
  %208 = getelementptr inbounds [4 x float], ptr %207, i64 0, i64 0
  store float %205, ptr %208, align 4, !tbaa !5
  %209 = load double, ptr %15, align 8, !tbaa !9
  %210 = load double, ptr %3, align 8, !tbaa !9
  %211 = load double, ptr %11, align 8, !tbaa !9
  %212 = load double, ptr %9, align 8, !tbaa !9
  %213 = load double, ptr %5, align 8, !tbaa !9
  %214 = fmul double %212, %213
  %215 = fneg double %214
  %216 = call double @llvm.fmuladd.f64(double %210, double %211, double %215)
  %217 = fmul double %209, %216
  %218 = fptrunc double %217 to float
  %219 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %220 = getelementptr inbounds [4 x [4 x float]], ptr %219, i64 0, i64 1
  %221 = getelementptr inbounds [4 x float], ptr %220, i64 0, i64 1
  store float %218, ptr %221, align 4, !tbaa !5
  %222 = load double, ptr %15, align 8, !tbaa !9
  %223 = load double, ptr %3, align 8, !tbaa !9
  %224 = fneg double %223
  %225 = load double, ptr %8, align 8, !tbaa !9
  %226 = load double, ptr %6, align 8, !tbaa !9
  %227 = load double, ptr %5, align 8, !tbaa !9
  %228 = fmul double %226, %227
  %229 = call double @llvm.fmuladd.f64(double %224, double %225, double %228)
  %230 = fmul double %222, %229
  %231 = fptrunc double %230 to float
  %232 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %233 = getelementptr inbounds [4 x [4 x float]], ptr %232, i64 0, i64 1
  %234 = getelementptr inbounds [4 x float], ptr %233, i64 0, i64 2
  store float %231, ptr %234, align 4, !tbaa !5
  %235 = load double, ptr %15, align 8, !tbaa !9
  %236 = load double, ptr %3, align 8, !tbaa !9
  %237 = load double, ptr %8, align 8, !tbaa !9
  %238 = fmul double %236, %237
  %239 = load double, ptr %14, align 8, !tbaa !9
  %240 = load double, ptr %3, align 8, !tbaa !9
  %241 = load double, ptr %13, align 8, !tbaa !9
  %242 = fmul double %240, %241
  %243 = load double, ptr %11, align 8, !tbaa !9
  %244 = fmul double %242, %243
  %245 = fneg double %244
  %246 = call double @llvm.fmuladd.f64(double %238, double %239, double %245)
  %247 = load double, ptr %6, align 8, !tbaa !9
  %248 = load double, ptr %5, align 8, !tbaa !9
  %249 = fmul double %247, %248
  %250 = load double, ptr %14, align 8, !tbaa !9
  %251 = fneg double %249
  %252 = call double @llvm.fmuladd.f64(double %251, double %250, double %246)
  %253 = load double, ptr %6, align 8, !tbaa !9
  %254 = load double, ptr %12, align 8, !tbaa !9
  %255 = fmul double %253, %254
  %256 = load double, ptr %11, align 8, !tbaa !9
  %257 = call double @llvm.fmuladd.f64(double %255, double %256, double %252)
  %258 = load double, ptr %9, align 8, !tbaa !9
  %259 = load double, ptr %5, align 8, !tbaa !9
  %260 = fmul double %258, %259
  %261 = load double, ptr %13, align 8, !tbaa !9
  %262 = call double @llvm.fmuladd.f64(double %260, double %261, double %257)
  %263 = load double, ptr %9, align 8, !tbaa !9
  %264 = load double, ptr %12, align 8, !tbaa !9
  %265 = fmul double %263, %264
  %266 = load double, ptr %8, align 8, !tbaa !9
  %267 = fneg double %265
  %268 = call double @llvm.fmuladd.f64(double %267, double %266, double %262)
  %269 = fmul double %235, %268
  %270 = fptrunc double %269 to float
  %271 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %272 = getelementptr inbounds [4 x [4 x float]], ptr %271, i64 0, i64 1
  %273 = getelementptr inbounds [4 x float], ptr %272, i64 0, i64 3
  store float %270, ptr %273, align 4, !tbaa !5
  %274 = load double, ptr %15, align 8, !tbaa !9
  %275 = load double, ptr %6, align 8, !tbaa !9
  %276 = load double, ptr %10, align 8, !tbaa !9
  %277 = load double, ptr %9, align 8, !tbaa !9
  %278 = load double, ptr %7, align 8, !tbaa !9
  %279 = fmul double %277, %278
  %280 = fneg double %279
  %281 = call double @llvm.fmuladd.f64(double %275, double %276, double %280)
  %282 = fmul double %274, %281
  %283 = fptrunc double %282 to float
  %284 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %285 = getelementptr inbounds [4 x [4 x float]], ptr %284, i64 0, i64 2
  %286 = getelementptr inbounds [4 x float], ptr %285, i64 0, i64 0
  store float %283, ptr %286, align 4, !tbaa !5
  %287 = load double, ptr %15, align 8, !tbaa !9
  %288 = load double, ptr %3, align 8, !tbaa !9
  %289 = fneg double %288
  %290 = load double, ptr %10, align 8, !tbaa !9
  %291 = load double, ptr %9, align 8, !tbaa !9
  %292 = load double, ptr %4, align 8, !tbaa !9
  %293 = fmul double %291, %292
  %294 = call double @llvm.fmuladd.f64(double %289, double %290, double %293)
  %295 = fmul double %287, %294
  %296 = fptrunc double %295 to float
  %297 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %298 = getelementptr inbounds [4 x [4 x float]], ptr %297, i64 0, i64 2
  %299 = getelementptr inbounds [4 x float], ptr %298, i64 0, i64 1
  store float %296, ptr %299, align 4, !tbaa !5
  %300 = load double, ptr %15, align 8, !tbaa !9
  %301 = load double, ptr %3, align 8, !tbaa !9
  %302 = load double, ptr %7, align 8, !tbaa !9
  %303 = load double, ptr %6, align 8, !tbaa !9
  %304 = load double, ptr %4, align 8, !tbaa !9
  %305 = fmul double %303, %304
  %306 = fneg double %305
  %307 = call double @llvm.fmuladd.f64(double %301, double %302, double %306)
  %308 = fmul double %300, %307
  %309 = fptrunc double %308 to float
  %310 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %311 = getelementptr inbounds [4 x [4 x float]], ptr %310, i64 0, i64 2
  %312 = getelementptr inbounds [4 x float], ptr %311, i64 0, i64 2
  store float %309, ptr %312, align 4, !tbaa !5
  %313 = load double, ptr %15, align 8, !tbaa !9
  %314 = load double, ptr %3, align 8, !tbaa !9
  %315 = fneg double %314
  %316 = load double, ptr %7, align 8, !tbaa !9
  %317 = fmul double %315, %316
  %318 = load double, ptr %14, align 8, !tbaa !9
  %319 = load double, ptr %3, align 8, !tbaa !9
  %320 = load double, ptr %10, align 8, !tbaa !9
  %321 = fmul double %319, %320
  %322 = load double, ptr %13, align 8, !tbaa !9
  %323 = fmul double %321, %322
  %324 = call double @llvm.fmuladd.f64(double %317, double %318, double %323)
  %325 = load double, ptr %6, align 8, !tbaa !9
  %326 = load double, ptr %4, align 8, !tbaa !9
  %327 = fmul double %325, %326
  %328 = load double, ptr %14, align 8, !tbaa !9
  %329 = call double @llvm.fmuladd.f64(double %327, double %328, double %324)
  %330 = load double, ptr %6, align 8, !tbaa !9
  %331 = load double, ptr %10, align 8, !tbaa !9
  %332 = fmul double %330, %331
  %333 = load double, ptr %12, align 8, !tbaa !9
  %334 = fneg double %332
  %335 = call double @llvm.fmuladd.f64(double %334, double %333, double %329)
  %336 = load double, ptr %9, align 8, !tbaa !9
  %337 = load double, ptr %4, align 8, !tbaa !9
  %338 = fmul double %336, %337
  %339 = load double, ptr %13, align 8, !tbaa !9
  %340 = fneg double %338
  %341 = call double @llvm.fmuladd.f64(double %340, double %339, double %335)
  %342 = load double, ptr %9, align 8, !tbaa !9
  %343 = load double, ptr %7, align 8, !tbaa !9
  %344 = fmul double %342, %343
  %345 = load double, ptr %12, align 8, !tbaa !9
  %346 = call double @llvm.fmuladd.f64(double %344, double %345, double %341)
  %347 = fmul double %313, %346
  %348 = fptrunc double %347 to float
  %349 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %350 = getelementptr inbounds [4 x [4 x float]], ptr %349, i64 0, i64 2
  %351 = getelementptr inbounds [4 x float], ptr %350, i64 0, i64 3
  store float %348, ptr %351, align 4, !tbaa !5
  %352 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %353 = getelementptr inbounds [4 x [4 x float]], ptr %352, i64 0, i64 3
  %354 = getelementptr inbounds [4 x float], ptr %353, i64 0, i64 2
  store float 0.000000e+00, ptr %354, align 4, !tbaa !5
  %355 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %356 = getelementptr inbounds [4 x [4 x float]], ptr %355, i64 0, i64 3
  %357 = getelementptr inbounds [4 x float], ptr %356, i64 0, i64 1
  store float 0.000000e+00, ptr %357, align 4, !tbaa !5
  %358 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %359 = getelementptr inbounds [4 x [4 x float]], ptr %358, i64 0, i64 3
  %360 = getelementptr inbounds [4 x float], ptr %359, i64 0, i64 0
  store float 0.000000e+00, ptr %360, align 4, !tbaa !5
  %361 = load double, ptr %15, align 8, !tbaa !9
  %362 = fpext double %361 to x86_fp80
  %363 = fcmp oeq x86_fp80 %362, 0xK00000000000000000000
  %364 = zext i1 %363 to i64
  %365 = select i1 %363, x86_fp80 0xK00000000000000000000, x86_fp80 0xK3FFF8000000000000000
  %366 = fptrunc x86_fp80 %365 to float
  %367 = getelementptr inbounds %struct.mat44, ptr %0, i32 0, i32 0
  %368 = getelementptr inbounds [4 x [4 x float]], ptr %367, i64 0, i64 3
  %369 = getelementptr inbounds [4 x float], ptr %368, i64 0, i64 3
  store float %366, ptr %369, align 4, !tbaa !5
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %13) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %11) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %9) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %8) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %7) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %6) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %5) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %4) #3
  call void @llvm.lifetime.end.p0(i64 8, ptr %3) #3
  ret void
}

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #2 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { nounwind }

!llvm.module.flags = !{!0, !1, !2, !3}
!llvm.ident = !{!4}

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"float", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"double", !7, i64 0}


; Optimized IR (opt -passes='default<O3>'):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

%struct.mat44 = type { [4 x [4 x float]] }

; Function Attrs: mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #0

; Function Attrs: mustprogress nofree norecurse nosync nounwind willreturn memory(argmem: readwrite) uwtable
define dso_local void @nifti_mat44_inverse(ptr dead_on_unwind noalias nocapture writable writeonly sret(%struct.mat44) align 4 %0, ptr nocapture noundef readonly byval(%struct.mat44) align 8 %1) local_unnamed_addr #1 {
  %3 = load float, ptr %1, align 8, !tbaa !5
  %4 = fpext float %3 to double
  %5 = getelementptr inbounds i8, ptr %1, i64 4
  %6 = load float, ptr %5, align 4, !tbaa !5
  %7 = fpext float %6 to double
  %8 = getelementptr inbounds i8, ptr %1, i64 8
  %9 = load float, ptr %8, align 8, !tbaa !5
  %10 = fpext float %9 to double
  %11 = getelementptr inbounds i8, ptr %1, i64 16
  %12 = load <4 x float>, ptr %11, align 8
  %13 = shufflevector <4 x float> %12, <4 x float> poison, <2 x i32> <i32 0, i32 poison>
  %14 = getelementptr inbounds i8, ptr %1, i64 20
  %15 = load <4 x float>, ptr %14, align 4
  %16 = shufflevector <4 x float> %15, <4 x float> poison, <2 x i32> <i32 0, i32 poison>
  %17 = getelementptr inbounds i8, ptr %1, i64 24
  %18 = load float, ptr %17, align 8, !tbaa !5
  %19 = getelementptr inbounds i8, ptr %1, i64 32
  %20 = load float, ptr %19, align 8, !tbaa !5
  %21 = insertelement <2 x float> %16, float %20, i64 1
  %22 = fpext <2 x float> %21 to <2 x double>
  %23 = getelementptr inbounds i8, ptr %1, i64 36
  %24 = load float, ptr %23, align 4, !tbaa !5
  %25 = insertelement <2 x float> %13, float %24, i64 1
  %26 = fpext <2 x float> %25 to <2 x double>
  %27 = getelementptr inbounds i8, ptr %1, i64 40
  %28 = load <4 x float>, ptr %27, align 8
  %29 = shufflevector <4 x float> %28, <4 x float> poison, <2 x i32> <i32 0, i32 poison>
  %30 = insertelement <2 x float> %29, float %18, i64 1
  %31 = fpext <2 x float> %30 to <2 x double>
  %32 = getelementptr inbounds i8, ptr %1, i64 12
  %33 = load float, ptr %32, align 4, !tbaa !5
  %34 = getelementptr inbounds i8, ptr %1, i64 28
  %35 = load float, ptr %34, align 4, !tbaa !5
  %36 = getelementptr inbounds i8, ptr %1, i64 44
  %37 = load float, ptr %36, align 4, !tbaa !5
  %38 = extractelement <2 x double> %22, i64 0
  %39 = fmul double %38, %4
  %40 = extractelement <2 x double> %26, i64 1
  %41 = fmul double %40, %4
  %42 = fneg double %41
  %43 = extractelement <2 x double> %31, i64 1
  %44 = fmul double %43, %42
  %45 = extractelement <2 x double> %31, i64 0
  %46 = tail call double @llvm.fmuladd.f64(double %39, double %45, double %44)
  %47 = extractelement <2 x double> %26, i64 0
  %48 = fmul double %47, %7
  %49 = fneg double %48
  %50 = tail call double @llvm.fmuladd.f64(double %49, double %45, double %46)
  %51 = fmul double %47, %40
  %52 = tail call double @llvm.fmuladd.f64(double %51, double %10, double %50)
  %53 = extractelement <2 x double> %22, i64 1
  %54 = fmul double %53, %7
  %55 = tail call double @llvm.fmuladd.f64(double %54, double %43, double %52)
  %56 = fmul double %38, %53
  %57 = fneg double %56
  %58 = tail call double @llvm.fmuladd.f64(double %57, double %10, double %55)
  %59 = fcmp une double %58, 0.000000e+00
  br i1 %59, label %60, label %64

60:                                               ; preds = %2
  %61 = fpext double %58 to x86_fp80
  %62 = fdiv x86_fp80 0xK3FFF8000000000000000, %61
  %63 = fptrunc x86_fp80 %62 to double
  br label %64

64:                                               ; preds = %60, %2
  %.0 = phi double [ %63, %60 ], [ %58, %2 ]
  %65 = shufflevector <2 x double> %31, <2 x double> poison, <4 x i32> <i32 0, i32 0, i32 1, i32 1>
  %66 = fpext float %37 to double
  %67 = fpext float %35 to double
  %68 = fpext float %33 to double
  %69 = fneg double %40
  %70 = fneg double %7
  %71 = shufflevector <2 x double> %31, <2 x double> %26, <2 x i32> <i32 1, i32 3>
  %72 = insertelement <2 x double> poison, double %69, i64 0
  %73 = insertelement <2 x double> %72, double %10, i64 1
  %74 = fmul <2 x double> %71, %73
  %75 = fmul double %38, %10
  %76 = fneg double %75
  %77 = fmul double %43, %70
  %78 = fmul double %7, %67
  %79 = fmul double %78, %45
  %80 = tail call double @llvm.fmuladd.f64(double %77, double %66, double %79)
  %81 = tail call double @llvm.fmuladd.f64(double %75, double %66, double %80)
  %82 = fneg double %38
  %83 = fmul double %82, %68
  %84 = tail call double @llvm.fmuladd.f64(double %83, double %45, double %81)
  %85 = extractelement <2 x double> %74, i64 1
  %86 = fneg double %85
  %87 = tail call double @llvm.fmuladd.f64(double %86, double %67, double %84)
  %88 = fmul double %40, %68
  %89 = insertelement <4 x double> %65, double %70, i64 1
  %90 = shufflevector <4 x double> %89, <4 x double> poison, <4 x i32> <i32 0, i32 1, i32 2, i32 2>
  %91 = shufflevector <2 x double> %22, <2 x double> %31, <4 x i32> <i32 0, i32 2, i32 poison, i32 poison>
  %92 = insertelement <4 x double> %91, double %7, i64 2
  %93 = insertelement <4 x double> %92, double %88, i64 3
  %94 = shufflevector <2 x double> %74, <2 x double> poison, <4 x i32> <i32 0, i32 1, i32 poison, i32 poison>
  %95 = insertelement <4 x double> %94, double %76, i64 2
  %96 = insertelement <4 x double> %95, double %87, i64 3
  %97 = tail call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %90, <4 x double> %93, <4 x double> %96)
  %98 = insertelement <4 x double> poison, double %.0, i64 0
  %99 = shufflevector <4 x double> %98, <4 x double> poison, <4 x i32> zeroinitializer
  %100 = fmul <4 x double> %97, %99
  %101 = fptrunc <4 x double> %100 to <4 x float>
  store <4 x float> %101, ptr %0, align 4, !tbaa !5
  %102 = fneg double %47
  %103 = fmul double %43, %53
  %104 = getelementptr inbounds i8, ptr %0, i64 16
  %105 = fmul double %53, %10
  %106 = fneg double %105
  %107 = fneg double %4
  %108 = fmul double %47, %10
  %109 = fmul double %43, %4
  %110 = fmul double %107, %67
  %111 = fmul double %110, %45
  %112 = tail call double @llvm.fmuladd.f64(double %109, double %66, double %111)
  %113 = fneg double %108
  %114 = tail call double @llvm.fmuladd.f64(double %113, double %66, double %112)
  %115 = fmul double %47, %68
  %116 = tail call double @llvm.fmuladd.f64(double %115, double %45, double %114)
  %117 = tail call double @llvm.fmuladd.f64(double %105, double %67, double %116)
  %118 = fneg double %53
  %119 = fmul double %118, %68
  %120 = insertelement <4 x double> poison, double %102, i64 0
  %121 = insertelement <4 x double> %120, double %4, i64 1
  %122 = insertelement <4 x double> %121, double %107, i64 2
  %123 = insertelement <4 x double> %122, double %119, i64 3
  %124 = insertelement <4 x double> poison, double %103, i64 0
  %125 = insertelement <4 x double> %124, double %106, i64 1
  %126 = insertelement <4 x double> %125, double %108, i64 2
  %127 = insertelement <4 x double> %126, double %117, i64 3
  %128 = tail call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %123, <4 x double> %65, <4 x double> %127)
  %129 = fmul <4 x double> %128, %99
  %130 = fptrunc <4 x double> %129 to <4 x float>
  store <4 x float> %130, ptr %104, align 4, !tbaa !5
  %131 = getelementptr inbounds i8, ptr %0, i64 32
  %132 = fmul double %38, %107
  %133 = fmul double %41, %67
  %134 = tail call double @llvm.fmuladd.f64(double %132, double %66, double %133)
  %135 = tail call double @llvm.fmuladd.f64(double %48, double %66, double %134)
  %136 = fneg double %51
  %137 = tail call double @llvm.fmuladd.f64(double %136, double %68, double %135)
  %138 = fneg double %54
  %139 = tail call double @llvm.fmuladd.f64(double %138, double %67, double %137)
  %140 = shufflevector <2 x double> %26, <2 x double> poison, <4 x i32> <i32 1, i32 poison, i32 poison, i32 poison>
  %141 = insertelement <4 x double> %140, double %107, i64 1
  %142 = insertelement <4 x double> %141, double %4, i64 2
  %143 = insertelement <4 x double> %142, double %56, i64 3
  %144 = shufflevector <2 x double> %26, <2 x double> %22, <4 x i32> <i32 0, i32 1, i32 2, i32 poison>
  %145 = insertelement <4 x double> %144, double %68, i64 3
  %146 = insertelement <4 x double> poison, double %57, i64 0
  %147 = insertelement <4 x double> %146, double %54, i64 1
  %148 = insertelement <4 x double> %147, double %49, i64 2
  %149 = insertelement <4 x double> %148, double %139, i64 3
  %150 = tail call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %143, <4 x double> %145, <4 x double> %149)
  %151 = fmul <4 x double> %150, %99
  %152 = fptrunc <4 x double> %151 to <4 x float>
  store <4 x float> %152, ptr %131, align 4, !tbaa !5
  %153 = getelementptr inbounds i8, ptr %0, i64 48
  %154 = getelementptr inbounds i8, ptr %0, i64 56
  store float 0.000000e+00, ptr %154, align 4, !tbaa !5
  store <2 x float> zeroinitializer, ptr %153, align 4, !tbaa !5
  %155 = fcmp oeq double %.0, 0.000000e+00
  %156 = select i1 %155, float 0.000000e+00, float 1.000000e+00
  %157 = getelementptr inbounds i8, ptr %0, i64 60
  store float %156, ptr %157, align 4, !tbaa !5
  ret void
}

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare <4 x double> @llvm.fmuladd.v4f64(<4 x double>, <4 x double>, <4 x double>) #2

attributes #0 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #1 = { mustprogress nofree norecurse nosync nounwind willreturn memory(argmem: readwrite) uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }

!llvm.module.flags = !{!0, !1, !2, !3}
!llvm.ident = !{!4}

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"float", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}


; Effective O3 passes: EarlyCSEPass,GlobalOptPass,InferFunctionAttrsPass,InstCombinePass,PostOrderFunctionAttrsPass,ReassociatePass,SLPVectorizerPass,SROAPass,TailCallElimPass,VectorCombinePass