; LLVM IR from stan-dev/math :: lib/sundials_6.1.1/src/sundials (sundials_iterative)

; 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"

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

; Function Attrs: nounwind
declare double @sqrt(double noundef) #1

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

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

; Function Attrs: nounwind uwtable
define dso_local i32 @SUNQRfact(i32 noundef %0, ptr noundef %1, ptr noundef %2, i32 noundef %3) #3 {
  %5 = alloca i32, align 4
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca i32, align 4
  %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 i32, align 4
  %15 = alloca i32, align 4
  %16 = alloca i32, align 4
  %17 = alloca i32, align 4
  %18 = alloca i32, align 4
  %19 = alloca i32, align 4
  store i32 %0, ptr %5, align 4, !tbaa !5
  store ptr %1, ptr %6, align 8, !tbaa !9
  store ptr %2, ptr %7, align 8, !tbaa !9
  store i32 %3, ptr %8, align 4, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %9) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %11) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %13) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %14) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %15) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %16) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %17) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %18) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %19) #4
  store i32 0, ptr %19, align 4, !tbaa !5
  %20 = load i32, ptr %8, align 4, !tbaa !5
  switch i32 %20, label %208 [
    i32 0, label %21
  ]

21:                                               ; preds = %4
  store i32 0, ptr %19, align 4, !tbaa !5
  store i32 0, ptr %16, align 4, !tbaa !5
  br label %22

22:                                               ; preds = %204, %21
  %23 = load i32, ptr %16, align 4, !tbaa !5
  %24 = load i32, ptr %5, align 4, !tbaa !5
  %25 = icmp slt i32 %23, %24
  br i1 %25, label %26, label %207

26:                                               ; preds = %22
  store i32 0, ptr %15, align 4, !tbaa !5
  br label %27

27:                                               ; preds = %95, %26
  %28 = load i32, ptr %15, align 4, !tbaa !5
  %29 = load i32, ptr %16, align 4, !tbaa !5
  %30 = sub nsw i32 %29, 1
  %31 = icmp slt i32 %28, %30
  br i1 %31, label %32, label %98

32:                                               ; preds = %27
  %33 = load i32, ptr %15, align 4, !tbaa !5
  %34 = mul nsw i32 2, %33
  store i32 %34, ptr %14, align 4, !tbaa !5
  %35 = load ptr, ptr %6, align 8, !tbaa !9
  %36 = load i32, ptr %15, align 4, !tbaa !5
  %37 = sext i32 %36 to i64
  %38 = getelementptr inbounds ptr, ptr %35, i64 %37
  %39 = load ptr, ptr %38, align 8, !tbaa !9
  %40 = load i32, ptr %16, align 4, !tbaa !5
  %41 = sext i32 %40 to i64
  %42 = getelementptr inbounds double, ptr %39, i64 %41
  %43 = load double, ptr %42, align 8, !tbaa !11
  store double %43, ptr %11, align 8, !tbaa !11
  %44 = load ptr, ptr %6, align 8, !tbaa !9
  %45 = load i32, ptr %15, align 4, !tbaa !5
  %46 = add nsw i32 %45, 1
  %47 = sext i32 %46 to i64
  %48 = getelementptr inbounds ptr, ptr %44, i64 %47
  %49 = load ptr, ptr %48, align 8, !tbaa !9
  %50 = load i32, ptr %16, align 4, !tbaa !5
  %51 = sext i32 %50 to i64
  %52 = getelementptr inbounds double, ptr %49, i64 %51
  %53 = load double, ptr %52, align 8, !tbaa !11
  store double %53, ptr %12, align 8, !tbaa !11
  %54 = load ptr, ptr %7, align 8, !tbaa !9
  %55 = load i32, ptr %14, align 4, !tbaa !5
  %56 = sext i32 %55 to i64
  %57 = getelementptr inbounds double, ptr %54, i64 %56
  %58 = load double, ptr %57, align 8, !tbaa !11
  store double %58, ptr %9, align 8, !tbaa !11
  %59 = load ptr, ptr %7, align 8, !tbaa !9
  %60 = load i32, ptr %14, align 4, !tbaa !5
  %61 = add nsw i32 %60, 1
  %62 = sext i32 %61 to i64
  %63 = getelementptr inbounds double, ptr %59, i64 %62
  %64 = load double, ptr %63, align 8, !tbaa !11
  store double %64, ptr %10, align 8, !tbaa !11
  %65 = load double, ptr %9, align 8, !tbaa !11
  %66 = load double, ptr %11, align 8, !tbaa !11
  %67 = load double, ptr %10, align 8, !tbaa !11
  %68 = load double, ptr %12, align 8, !tbaa !11
  %69 = fmul double %67, %68
  %70 = fneg double %69
  %71 = call double @llvm.fmuladd.f64(double %65, double %66, double %70)
  %72 = load ptr, ptr %6, align 8, !tbaa !9
  %73 = load i32, ptr %15, align 4, !tbaa !5
  %74 = sext i32 %73 to i64
  %75 = getelementptr inbounds ptr, ptr %72, i64 %74
  %76 = load ptr, ptr %75, align 8, !tbaa !9
  %77 = load i32, ptr %16, align 4, !tbaa !5
  %78 = sext i32 %77 to i64
  %79 = getelementptr inbounds double, ptr %76, i64 %78
  store double %71, ptr %79, align 8, !tbaa !11
  %80 = load double, ptr %10, align 8, !tbaa !11
  %81 = load double, ptr %11, align 8, !tbaa !11
  %82 = load double, ptr %9, align 8, !tbaa !11
  %83 = load double, ptr %12, align 8, !tbaa !11
  %84 = fmul double %82, %83
  %85 = call double @llvm.fmuladd.f64(double %80, double %81, double %84)
  %86 = load ptr, ptr %6, align 8, !tbaa !9
  %87 = load i32, ptr %15, align 4, !tbaa !5
  %88 = add nsw i32 %87, 1
  %89 = sext i32 %88 to i64
  %90 = getelementptr inbounds ptr, ptr %86, i64 %89
  %91 = load ptr, ptr %90, align 8, !tbaa !9
  %92 = load i32, ptr %16, align 4, !tbaa !5
  %93 = sext i32 %92 to i64
  %94 = getelementptr inbounds double, ptr %91, i64 %93
  store double %85, ptr %94, align 8, !tbaa !11
  br label %95

95:                                               ; preds = %32
  %96 = load i32, ptr %15, align 4, !tbaa !5
  %97 = add nsw i32 %96, 1
  store i32 %97, ptr %15, align 4, !tbaa !5
  br label %27, !llvm.loop !13

98:                                               ; preds = %27
  %99 = load i32, ptr %16, align 4, !tbaa !5
  %100 = mul nsw i32 2, %99
  store i32 %100, ptr %17, align 4, !tbaa !5
  %101 = load ptr, ptr %6, align 8, !tbaa !9
  %102 = load i32, ptr %16, align 4, !tbaa !5
  %103 = sext i32 %102 to i64
  %104 = getelementptr inbounds ptr, ptr %101, i64 %103
  %105 = load ptr, ptr %104, align 8, !tbaa !9
  %106 = load i32, ptr %16, align 4, !tbaa !5
  %107 = sext i32 %106 to i64
  %108 = getelementptr inbounds double, ptr %105, i64 %107
  %109 = load double, ptr %108, align 8, !tbaa !11
  store double %109, ptr %11, align 8, !tbaa !11
  %110 = load ptr, ptr %6, align 8, !tbaa !9
  %111 = load i32, ptr %16, align 4, !tbaa !5
  %112 = add nsw i32 %111, 1
  %113 = sext i32 %112 to i64
  %114 = getelementptr inbounds ptr, ptr %110, i64 %113
  %115 = load ptr, ptr %114, align 8, !tbaa !9
  %116 = load i32, ptr %16, align 4, !tbaa !5
  %117 = sext i32 %116 to i64
  %118 = getelementptr inbounds double, ptr %115, i64 %117
  %119 = load double, ptr %118, align 8, !tbaa !11
  store double %119, ptr %12, align 8, !tbaa !11
  %120 = load double, ptr %12, align 8, !tbaa !11
  %121 = fcmp oeq double %120, 0.000000e+00
  br i1 %121, label %122, label %123

122:                                              ; preds = %98
  store double 1.000000e+00, ptr %9, align 8, !tbaa !11
  store double 0.000000e+00, ptr %10, align 8, !tbaa !11
  br label %172

123:                                              ; preds = %98
  %124 = load double, ptr %12, align 8, !tbaa !11
  %125 = call double @llvm.fabs.f64(double %124)
  %126 = load double, ptr %11, align 8, !tbaa !11
  %127 = call double @llvm.fabs.f64(double %126)
  %128 = fcmp oge double %125, %127
  br i1 %128, label %129, label %150

129:                                              ; preds = %123
  %130 = load double, ptr %11, align 8, !tbaa !11
  %131 = load double, ptr %12, align 8, !tbaa !11
  %132 = fdiv double %130, %131
  store double %132, ptr %13, align 8, !tbaa !11
  %133 = load double, ptr %13, align 8, !tbaa !11
  %134 = load double, ptr %13, align 8, !tbaa !11
  %135 = call double @llvm.fmuladd.f64(double %133, double %134, double 1.000000e+00)
  %136 = fcmp ole double %135, 0.000000e+00
  br i1 %136, label %137, label %138

137:                                              ; preds = %129
  br label %143

138:                                              ; preds = %129
  %139 = load double, ptr %13, align 8, !tbaa !11
  %140 = load double, ptr %13, align 8, !tbaa !11
  %141 = call double @llvm.fmuladd.f64(double %139, double %140, double 1.000000e+00)
  %142 = call double @sqrt(double noundef %141) #4
  br label %143

143:                                              ; preds = %138, %137
  %144 = phi double [ 0.000000e+00, %137 ], [ %142, %138 ]
  %145 = fdiv double -1.000000e+00, %144
  store double %145, ptr %10, align 8, !tbaa !11
  %146 = load double, ptr %10, align 8, !tbaa !11
  %147 = fneg double %146
  %148 = load double, ptr %13, align 8, !tbaa !11
  %149 = fmul double %147, %148
  store double %149, ptr %9, align 8, !tbaa !11
  br label %171

150:                                              ; preds = %123
  %151 = load double, ptr %12, align 8, !tbaa !11
  %152 = load double, ptr %11, align 8, !tbaa !11
  %153 = fdiv double %151, %152
  store double %153, ptr %13, align 8, !tbaa !11
  %154 = load double, ptr %13, align 8, !tbaa !11
  %155 = load double, ptr %13, align 8, !tbaa !11
  %156 = call double @llvm.fmuladd.f64(double %154, double %155, double 1.000000e+00)
  %157 = fcmp ole double %156, 0.000000e+00
  br i1 %157, label %158, label %159

158:                                              ; preds = %150
  br label %164

159:                                              ; preds = %150
  %160 = load double, ptr %13, align 8, !tbaa !11
  %161 = load double, ptr %13, align 8, !tbaa !11
  %162 = call double @llvm.fmuladd.f64(double %160, double %161, double 1.000000e+00)
  %163 = call double @sqrt(double noundef %162) #4
  br label %164

164:                                              ; preds = %159, %158
  %165 = phi double [ 0.000000e+00, %158 ], [ %163, %159 ]
  %166 = fdiv double 1.000000e+00, %165
  store double %166, ptr %9, align 8, !tbaa !11
  %167 = load double, ptr %9, align 8, !tbaa !11
  %168 = fneg double %167
  %169 = load double, ptr %13, align 8, !tbaa !11
  %170 = fmul double %168, %169
  store double %170, ptr %10, align 8, !tbaa !11
  br label %171

171:                                              ; preds = %164, %143
  br label %172

172:                                              ; preds = %171, %122
  %173 = load double, ptr %9, align 8, !tbaa !11
  %174 = load ptr, ptr %7, align 8, !tbaa !9
  %175 = load i32, ptr %17, align 4, !tbaa !5
  %176 = sext i32 %175 to i64
  %177 = getelementptr inbounds double, ptr %174, i64 %176
  store double %173, ptr %177, align 8, !tbaa !11
  %178 = load double, ptr %10, align 8, !tbaa !11
  %179 = load ptr, ptr %7, align 8, !tbaa !9
  %180 = load i32, ptr %17, align 4, !tbaa !5
  %181 = add nsw i32 %180, 1
  %182 = sext i32 %181 to i64
  %183 = getelementptr inbounds double, ptr %179, i64 %182
  store double %178, ptr %183, align 8, !tbaa !11
  %184 = load double, ptr %9, align 8, !tbaa !11
  %185 = load double, ptr %11, align 8, !tbaa !11
  %186 = load double, ptr %10, align 8, !tbaa !11
  %187 = load double, ptr %12, align 8, !tbaa !11
  %188 = fmul double %186, %187
  %189 = fneg double %188
  %190 = call double @llvm.fmuladd.f64(double %184, double %185, double %189)
  %191 = load ptr, ptr %6, align 8, !tbaa !9
  %192 = load i32, ptr %16, align 4, !tbaa !5
  %193 = sext i32 %192 to i64
  %194 = getelementptr inbounds ptr, ptr %191, i64 %193
  %195 = load ptr, ptr %194, align 8, !tbaa !9
  %196 = load i32, ptr %16, align 4, !tbaa !5
  %197 = sext i32 %196 to i64
  %198 = getelementptr inbounds double, ptr %195, i64 %197
  store double %190, ptr %198, align 8, !tbaa !11
  %199 = fcmp oeq double %190, 0.000000e+00
  br i1 %199, label %200, label %203

200:                                              ; preds = %172
  %201 = load i32, ptr %16, align 4, !tbaa !5
  %202 = add nsw i32 %201, 1
  store i32 %202, ptr %19, align 4, !tbaa !5
  br label %203

203:                                              ; preds = %200, %172
  br label %204

204:                                              ; preds = %203
  %205 = load i32, ptr %16, align 4, !tbaa !5
  %206 = add nsw i32 %205, 1
  store i32 %206, ptr %16, align 4, !tbaa !5
  br label %22, !llvm.loop !15

207:                                              ; preds = %22
  br label %385

208:                                              ; preds = %4
  %209 = load i32, ptr %5, align 4, !tbaa !5
  %210 = sub nsw i32 %209, 1
  store i32 %210, ptr %18, align 4, !tbaa !5
  store i32 0, ptr %19, align 4, !tbaa !5
  store i32 0, ptr %16, align 4, !tbaa !5
  br label %211

211:                                              ; preds = %278, %208
  %212 = load i32, ptr %16, align 4, !tbaa !5
  %213 = load i32, ptr %18, align 4, !tbaa !5
  %214 = icmp slt i32 %212, %213
  br i1 %214, label %215, label %281

215:                                              ; preds = %211
  %216 = load i32, ptr %16, align 4, !tbaa !5
  %217 = mul nsw i32 2, %216
  store i32 %217, ptr %14, align 4, !tbaa !5
  %218 = load ptr, ptr %6, align 8, !tbaa !9
  %219 = load i32, ptr %16, align 4, !tbaa !5
  %220 = sext i32 %219 to i64
  %221 = getelementptr inbounds ptr, ptr %218, i64 %220
  %222 = load ptr, ptr %221, align 8, !tbaa !9
  %223 = load i32, ptr %18, align 4, !tbaa !5
  %224 = sext i32 %223 to i64
  %225 = getelementptr inbounds double, ptr %222, i64 %224
  %226 = load double, ptr %225, align 8, !tbaa !11
  store double %226, ptr %11, align 8, !tbaa !11
  %227 = load ptr, ptr %6, align 8, !tbaa !9
  %228 = load i32, ptr %16, align 4, !tbaa !5
  %229 = add nsw i32 %228, 1
  %230 = sext i32 %229 to i64
  %231 = getelementptr inbounds ptr, ptr %227, i64 %230
  %232 = load ptr, ptr %231, align 8, !tbaa !9
  %233 = load i32, ptr %18, align 4, !tbaa !5
  %234 = sext i32 %233 to i64
  %235 = getelementptr inbounds double, ptr %232, i64 %234
  %236 = load double, ptr %235, align 8, !tbaa !11
  store double %236, ptr %12, align 8, !tbaa !11
  %237 = load ptr, ptr %7, align 8, !tbaa !9
  %238 = load i32, ptr %14, align 4, !tbaa !5
  %239 = sext i32 %238 to i64
  %240 = getelementptr inbounds double, ptr %237, i64 %239
  %241 = load double, ptr %240, align 8, !tbaa !11
  store double %241, ptr %9, align 8, !tbaa !11
  %242 = load ptr, ptr %7, align 8, !tbaa !9
  %243 = load i32, ptr %14, align 4, !tbaa !5
  %244 = add nsw i32 %243, 1
  %245 = sext i32 %244 to i64
  %246 = getelementptr inbounds double, ptr %242, i64 %245
  %247 = load double, ptr %246, align 8, !tbaa !11
  store double %247, ptr %10, align 8, !tbaa !11
  %248 = load double, ptr %9, align 8, !tbaa !11
  %249 = load double, ptr %11, align 8, !tbaa !11
  %250 = load double, ptr %10, align 8, !tbaa !11
  %251 = load double, ptr %12, align 8, !tbaa !11
  %252 = fmul double %250, %251
  %253 = fneg double %252
  %254 = call double @llvm.fmuladd.f64(double %248, double %249, double %253)
  %255 = load ptr, ptr %6, align 8, !tbaa !9
  %256 = load i32, ptr %16, align 4, !tbaa !5
  %257 = sext i32 %256 to i64
  %258 = getelementptr inbounds ptr, ptr %255, i64 %257
  %259 = load ptr, ptr %258, align 8, !tbaa !9
  %260 = load i32, ptr %18, align 4, !tbaa !5
  %261 = sext i32 %260 to i64
  %262 = getelementptr inbounds double, ptr %259, i64 %261
  store double %254, ptr %262, align 8, !tbaa !11
  %263 = load double, ptr %10, align 8, !tbaa !11
  %264 = load double, ptr %11, align 8, !tbaa !11
  %265 = load double, ptr %9, align 8, !tbaa !11
  %266 = load double, ptr %12, align 8, !tbaa !11
  %267 = fmul double %265, %266
  %268 = call double @llvm.fmuladd.f64(double %263, double %264, double %267)
  %269 = load ptr, ptr %6, align 8, !tbaa !9
  %270 = load i32, ptr %16, align 4, !tbaa !5
  %271 = add nsw i32 %270, 1
  %272 = sext i32 %271 to i64
  %273 = getelementptr inbounds ptr, ptr %269, i64 %272
  %274 = load ptr, ptr %273, align 8, !tbaa !9
  %275 = load i32, ptr %18, align 4, !tbaa !5
  %276 = sext i32 %275 to i64
  %277 = getelementptr inbounds double, ptr %274, i64 %276
  store double %268, ptr %277, align 8, !tbaa !11
  br label %278

278:                                              ; preds = %215
  %279 = load i32, ptr %16, align 4, !tbaa !5
  %280 = add nsw i32 %279, 1
  store i32 %280, ptr %16, align 4, !tbaa !5
  br label %211, !llvm.loop !16

281:                                              ; preds = %211
  %282 = load ptr, ptr %6, align 8, !tbaa !9
  %283 = load i32, ptr %18, align 4, !tbaa !5
  %284 = sext i32 %283 to i64
  %285 = getelementptr inbounds ptr, ptr %282, i64 %284
  %286 = load ptr, ptr %285, align 8, !tbaa !9
  %287 = load i32, ptr %18, align 4, !tbaa !5
  %288 = sext i32 %287 to i64
  %289 = getelementptr inbounds double, ptr %286, i64 %288
  %290 = load double, ptr %289, align 8, !tbaa !11
  store double %290, ptr %11, align 8, !tbaa !11
  %291 = load ptr, ptr %6, align 8, !tbaa !9
  %292 = load i32, ptr %5, align 4, !tbaa !5
  %293 = sext i32 %292 to i64
  %294 = getelementptr inbounds ptr, ptr %291, i64 %293
  %295 = load ptr, ptr %294, align 8, !tbaa !9
  %296 = load i32, ptr %18, align 4, !tbaa !5
  %297 = sext i32 %296 to i64
  %298 = getelementptr inbounds double, ptr %295, i64 %297
  %299 = load double, ptr %298, align 8, !tbaa !11
  store double %299, ptr %12, align 8, !tbaa !11
  %300 = load double, ptr %12, align 8, !tbaa !11
  %301 = fcmp oeq double %300, 0.000000e+00
  br i1 %301, label %302, label %303

302:                                              ; preds = %281
  store double 1.000000e+00, ptr %9, align 8, !tbaa !11
  store double 0.000000e+00, ptr %10, align 8, !tbaa !11
  br label %352

303:                                              ; preds = %281
  %304 = load double, ptr %12, align 8, !tbaa !11
  %305 = call double @llvm.fabs.f64(double %304)
  %306 = load double, ptr %11, align 8, !tbaa !11
  %307 = call double @llvm.fabs.f64(double %306)
  %308 = fcmp oge double %305, %307
  br i1 %308, label %309, label %330

309:                                              ; preds = %303
  %310 = load double, ptr %11, align 8, !tbaa !11
  %311 = load double, ptr %12, align 8, !tbaa !11
  %312 = fdiv double %310, %311
  store double %312, ptr %13, align 8, !tbaa !11
  %313 = load double, ptr %13, align 8, !tbaa !11
  %314 = load double, ptr %13, align 8, !tbaa !11
  %315 = call double @llvm.fmuladd.f64(double %313, double %314, double 1.000000e+00)
  %316 = fcmp ole double %315, 0.000000e+00
  br i1 %316, label %317, label %318

317:                                              ; preds = %309
  br label %323

318:                                              ; preds = %309
  %319 = load double, ptr %13, align 8, !tbaa !11
  %320 = load double, ptr %13, align 8, !tbaa !11
  %321 = call double @llvm.fmuladd.f64(double %319, double %320, double 1.000000e+00)
  %322 = call double @sqrt(double noundef %321) #4
  br label %323

323:                                              ; preds = %318, %317
  %324 = phi double [ 0.000000e+00, %317 ], [ %322, %318 ]
  %325 = fdiv double -1.000000e+00, %324
  store double %325, ptr %10, align 8, !tbaa !11
  %326 = load double, ptr %10, align 8, !tbaa !11
  %327 = fneg double %326
  %328 = load double, ptr %13, align 8, !tbaa !11
  %329 = fmul double %327, %328
  store double %329, ptr %9, align 8, !tbaa !11
  br label %351

330:                                              ; preds = %303
  %331 = load double, ptr %12, align 8, !tbaa !11
  %332 = load double, ptr %11, align 8, !tbaa !11
  %333 = fdiv double %331, %332
  store double %333, ptr %13, align 8, !tbaa !11
  %334 = load double, ptr %13, align 8, !tbaa !11
  %335 = load double, ptr %13, align 8, !tbaa !11
  %336 = call double @llvm.fmuladd.f64(double %334, double %335, double 1.000000e+00)
  %337 = fcmp ole double %336, 0.000000e+00
  br i1 %337, label %338, label %339

338:                                              ; preds = %330
  br label %344

339:                                              ; preds = %330
  %340 = load double, ptr %13, align 8, !tbaa !11
  %341 = load double, ptr %13, align 8, !tbaa !11
  %342 = call double @llvm.fmuladd.f64(double %340, double %341, double 1.000000e+00)
  %343 = call double @sqrt(double noundef %342) #4
  br label %344

344:                                              ; preds = %339, %338
  %345 = phi double [ 0.000000e+00, %338 ], [ %343, %339 ]
  %346 = fdiv double 1.000000e+00, %345
  store double %346, ptr %9, align 8, !tbaa !11
  %347 = load double, ptr %9, align 8, !tbaa !11
  %348 = fneg double %347
  %349 = load double, ptr %13, align 8, !tbaa !11
  %350 = fmul double %348, %349
  store double %350, ptr %10, align 8, !tbaa !11
  br label %351

351:                                              ; preds = %344, %323
  br label %352

352:                                              ; preds = %351, %302
  %353 = load i32, ptr %18, align 4, !tbaa !5
  %354 = mul nsw i32 2, %353
  store i32 %354, ptr %17, align 4, !tbaa !5
  %355 = load double, ptr %9, align 8, !tbaa !11
  %356 = load ptr, ptr %7, align 8, !tbaa !9
  %357 = load i32, ptr %17, align 4, !tbaa !5
  %358 = sext i32 %357 to i64
  %359 = getelementptr inbounds double, ptr %356, i64 %358
  store double %355, ptr %359, align 8, !tbaa !11
  %360 = load double, ptr %10, align 8, !tbaa !11
  %361 = load ptr, ptr %7, align 8, !tbaa !9
  %362 = load i32, ptr %17, align 4, !tbaa !5
  %363 = add nsw i32 %362, 1
  %364 = sext i32 %363 to i64
  %365 = getelementptr inbounds double, ptr %361, i64 %364
  store double %360, ptr %365, align 8, !tbaa !11
  %366 = load double, ptr %9, align 8, !tbaa !11
  %367 = load double, ptr %11, align 8, !tbaa !11
  %368 = load double, ptr %10, align 8, !tbaa !11
  %369 = load double, ptr %12, align 8, !tbaa !11
  %370 = fmul double %368, %369
  %371 = fneg double %370
  %372 = call double @llvm.fmuladd.f64(double %366, double %367, double %371)
  %373 = load ptr, ptr %6, align 8, !tbaa !9
  %374 = load i32, ptr %18, align 4, !tbaa !5
  %375 = sext i32 %374 to i64
  %376 = getelementptr inbounds ptr, ptr %373, i64 %375
  %377 = load ptr, ptr %376, align 8, !tbaa !9
  %378 = load i32, ptr %18, align 4, !tbaa !5
  %379 = sext i32 %378 to i64
  %380 = getelementptr inbounds double, ptr %377, i64 %379
  store double %372, ptr %380, align 8, !tbaa !11
  %381 = fcmp oeq double %372, 0.000000e+00
  br i1 %381, label %382, label %384

382:                                              ; preds = %352
  %383 = load i32, ptr %5, align 4, !tbaa !5
  store i32 %383, ptr %19, align 4, !tbaa !5
  br label %384

384:                                              ; preds = %382, %352
  br label %385

385:                                              ; preds = %384, %207
  %386 = load i32, ptr %19, align 4, !tbaa !5
  call void @llvm.lifetime.end.p0(i64 4, ptr %19) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %18) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %17) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %16) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %15) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %14) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %13) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %11) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %9) #4
  ret i32 %386
}

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

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nounwind "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) }
attributes #3 = { 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 #4 = { 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 = !{!"int", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"any pointer", !7, i64 0}
!11 = !{!12, !12, i64 0}
!12 = !{!"double", !7, i64 0}
!13 = distinct !{!13, !14}
!14 = !{!"llvm.loop.mustprogress"}
!15 = distinct !{!15, !14}
!16 = distinct !{!16, !14}


; 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"

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

; Function Attrs: nofree norecurse nosync nounwind memory(readwrite, inaccessiblemem: none) uwtable
define dso_local i32 @SUNQRfact(i32 noundef %0, ptr nocapture noundef readonly %1, ptr nocapture noundef %2, i32 noundef %3) local_unnamed_addr #1 {
  %cond = icmp eq i32 %3, 0
  br i1 %cond, label %.preheader159, label %65

.preheader159:                                    ; preds = %4
  %5 = icmp sgt i32 %0, 0
  br i1 %5, label %.preheader.preheader, label %.loopexit

.preheader.preheader:                             ; preds = %.preheader159
  %wide.trip.count177 = zext nneg i32 %0 to i64
  br label %.preheader

.preheader:                                       ; preds = %.preheader.preheader, %56
  %indvars.iv174 = phi i64 [ 0, %.preheader.preheader ], [ %indvars.iv.next175, %56 ]
  %indvars.iv170 = phi i32 [ -1, %.preheader.preheader ], [ %indvars.iv.next171, %56 ]
  %.0165 = phi i32 [ 0, %.preheader.preheader ], [ %.1, %56 ]
  %6 = icmp ugt i64 %indvars.iv174, 1
  br i1 %6, label %.lr.ph162, label %._crit_edge163

.lr.ph162:                                        ; preds = %.preheader
  %wide.trip.count172 = zext i32 %indvars.iv170 to i64
  %.pre180 = load ptr, ptr %1, align 8, !tbaa !5
  %.phi.trans.insert181 = getelementptr inbounds double, ptr %.pre180, i64 %indvars.iv174
  %.pre182 = load double, ptr %.phi.trans.insert181, align 8, !tbaa !9
  br label %7

7:                                                ; preds = %.lr.ph162, %7
  %8 = phi double [ %.pre182, %.lr.ph162 ], [ %25, %7 ]
  %9 = phi ptr [ %.pre180, %.lr.ph162 ], [ %13, %7 ]
  %indvars.iv167 = phi i64 [ 0, %.lr.ph162 ], [ %indvars.iv.next168, %7 ]
  %10 = shl nuw nsw i64 %indvars.iv167, 1
  %11 = getelementptr inbounds double, ptr %9, i64 %indvars.iv174
  %indvars.iv.next168 = add nuw nsw i64 %indvars.iv167, 1
  %12 = getelementptr inbounds ptr, ptr %1, i64 %indvars.iv.next168
  %13 = load ptr, ptr %12, align 8, !tbaa !5
  %14 = getelementptr inbounds double, ptr %13, i64 %indvars.iv174
  %15 = load double, ptr %14, align 8, !tbaa !9
  %16 = getelementptr inbounds double, ptr %2, i64 %10
  %17 = load double, ptr %16, align 8, !tbaa !9
  %18 = or disjoint i64 %10, 1
  %19 = getelementptr inbounds double, ptr %2, i64 %18
  %20 = load double, ptr %19, align 8, !tbaa !9
  %21 = fneg double %20
  %22 = fmul double %15, %21
  %23 = tail call double @llvm.fmuladd.f64(double %17, double %8, double %22)
  store double %23, ptr %11, align 8, !tbaa !9
  %24 = fmul double %15, %17
  %25 = tail call double @llvm.fmuladd.f64(double %20, double %8, double %24)
  store double %25, ptr %14, align 8, !tbaa !9
  %exitcond173.not = icmp eq i64 %indvars.iv.next168, %wide.trip.count172
  br i1 %exitcond173.not, label %._crit_edge163, label %7, !llvm.loop !11

._crit_edge163:                                   ; preds = %7, %.preheader
  %26 = shl nuw nsw i64 %indvars.iv174, 1
  %27 = getelementptr inbounds ptr, ptr %1, i64 %indvars.iv174
  %28 = load ptr, ptr %27, align 8, !tbaa !5
  %29 = getelementptr inbounds double, ptr %28, i64 %indvars.iv174
  %30 = load double, ptr %29, align 8, !tbaa !9
  %indvars.iv.next175 = add nuw nsw i64 %indvars.iv174, 1
  %31 = getelementptr inbounds ptr, ptr %1, i64 %indvars.iv.next175
  %32 = load ptr, ptr %31, align 8, !tbaa !5
  %33 = getelementptr inbounds double, ptr %32, i64 %indvars.iv174
  %34 = load double, ptr %33, align 8, !tbaa !9
  %35 = fcmp oeq double %34, 0.000000e+00
  br i1 %35, label %56, label %36

36:                                               ; preds = %._crit_edge163
  %37 = tail call double @llvm.fabs.f64(double %34)
  %38 = tail call double @llvm.fabs.f64(double %30)
  %39 = fcmp ult double %37, %38
  br i1 %39, label %48, label %40

40:                                               ; preds = %36
  %41 = fdiv double %30, %34
  %42 = tail call double @llvm.fmuladd.f64(double %41, double %41, double 1.000000e+00)
  %43 = fcmp ugt double %42, 0.000000e+00
  %sqrt = tail call double @llvm.sqrt.f64(double %42)
  %44 = fdiv double -1.000000e+00, %sqrt
  %45 = select i1 %43, double %44, double 0xFFF0000000000000
  %46 = fneg double %45
  %47 = fmul double %41, %46
  br label %56

48:                                               ; preds = %36
  %49 = fdiv double %34, %30
  %50 = tail call double @llvm.fmuladd.f64(double %49, double %49, double 1.000000e+00)
  %51 = fcmp ugt double %50, 0.000000e+00
  %sqrt156 = tail call double @llvm.sqrt.f64(double %50)
  %52 = fdiv double 1.000000e+00, %sqrt156
  %53 = select i1 %51, double %52, double 0x7FF0000000000000
  %54 = fneg double %53
  %55 = fmul double %49, %54
  br label %56

56:                                               ; preds = %._crit_edge163, %40, %48
  %.0146 = phi double [ %45, %40 ], [ %55, %48 ], [ 0.000000e+00, %._crit_edge163 ]
  %.0141 = phi double [ %47, %40 ], [ %53, %48 ], [ 1.000000e+00, %._crit_edge163 ]
  %57 = getelementptr inbounds double, ptr %2, i64 %26
  store double %.0141, ptr %57, align 8, !tbaa !9
  %58 = or disjoint i64 %26, 1
  %59 = getelementptr inbounds double, ptr %2, i64 %58
  store double %.0146, ptr %59, align 8, !tbaa !9
  %60 = fneg double %.0146
  %61 = fmul double %34, %60
  %62 = tail call double @llvm.fmuladd.f64(double %.0141, double %30, double %61)
  store double %62, ptr %29, align 8, !tbaa !9
  %63 = fcmp oeq double %62, 0.000000e+00
  %64 = trunc nuw nsw i64 %indvars.iv.next175 to i32
  %.1 = select i1 %63, i32 %64, i32 %.0165
  %indvars.iv.next171 = add nsw i32 %indvars.iv170, 1
  %exitcond178.not = icmp eq i64 %indvars.iv.next175, %wide.trip.count177
  br i1 %exitcond178.not, label %.loopexit, label %.preheader, !llvm.loop !13

65:                                               ; preds = %4
  %66 = add i32 %0, -1
  %67 = icmp sgt i32 %0, 1
  br i1 %67, label %.lr.ph, label %.._crit_edge_crit_edge

.._crit_edge_crit_edge:                           ; preds = %65
  %.pre183 = sext i32 %66 to i64
  br label %._crit_edge

.lr.ph:                                           ; preds = %65
  %68 = zext nneg i32 %66 to i64
  %wide.trip.count = zext nneg i32 %66 to i64
  %.pre = load ptr, ptr %1, align 8, !tbaa !5
  %.phi.trans.insert = getelementptr inbounds double, ptr %.pre, i64 %68
  %.pre179 = load double, ptr %.phi.trans.insert, align 8, !tbaa !9
  br label %69

69:                                               ; preds = %.lr.ph, %69
  %70 = phi double [ %.pre179, %.lr.ph ], [ %87, %69 ]
  %71 = phi ptr [ %.pre, %.lr.ph ], [ %75, %69 ]
  %indvars.iv = phi i64 [ 0, %.lr.ph ], [ %indvars.iv.next, %69 ]
  %72 = shl nuw nsw i64 %indvars.iv, 1
  %73 = getelementptr inbounds double, ptr %71, i64 %68
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %74 = getelementptr inbounds ptr, ptr %1, i64 %indvars.iv.next
  %75 = load ptr, ptr %74, align 8, !tbaa !5
  %76 = getelementptr inbounds double, ptr %75, i64 %68
  %77 = load double, ptr %76, align 8, !tbaa !9
  %78 = getelementptr inbounds double, ptr %2, i64 %72
  %79 = load double, ptr %78, align 8, !tbaa !9
  %80 = or disjoint i64 %72, 1
  %81 = getelementptr inbounds double, ptr %2, i64 %80
  %82 = load double, ptr %81, align 8, !tbaa !9
  %83 = fneg double %82
  %84 = fmul double %77, %83
  %85 = tail call double @llvm.fmuladd.f64(double %79, double %70, double %84)
  store double %85, ptr %73, align 8, !tbaa !9
  %86 = fmul double %77, %79
  %87 = tail call double @llvm.fmuladd.f64(double %82, double %70, double %86)
  store double %87, ptr %76, align 8, !tbaa !9
  %exitcond.not = icmp eq i64 %indvars.iv.next, %wide.trip.count
  br i1 %exitcond.not, label %._crit_edge, label %69, !llvm.loop !14

._crit_edge:                                      ; preds = %69, %.._crit_edge_crit_edge
  %.pre-phi = phi i64 [ %.pre183, %.._crit_edge_crit_edge ], [ %68, %69 ]
  %88 = getelementptr inbounds ptr, ptr %1, i64 %.pre-phi
  %89 = load ptr, ptr %88, align 8, !tbaa !5
  %90 = getelementptr inbounds double, ptr %89, i64 %.pre-phi
  %91 = load double, ptr %90, align 8, !tbaa !9
  %92 = sext i32 %0 to i64
  %93 = getelementptr inbounds ptr, ptr %1, i64 %92
  %94 = load ptr, ptr %93, align 8, !tbaa !5
  %95 = getelementptr inbounds double, ptr %94, i64 %.pre-phi
  %96 = load double, ptr %95, align 8, !tbaa !9
  %97 = fcmp oeq double %96, 0.000000e+00
  br i1 %97, label %118, label %98

98:                                               ; preds = %._crit_edge
  %99 = tail call double @llvm.fabs.f64(double %96)
  %100 = tail call double @llvm.fabs.f64(double %91)
  %101 = fcmp ult double %99, %100
  br i1 %101, label %110, label %102

102:                                              ; preds = %98
  %103 = fdiv double %91, %96
  %104 = tail call double @llvm.fmuladd.f64(double %103, double %103, double 1.000000e+00)
  %105 = fcmp ugt double %104, 0.000000e+00
  %sqrt157 = tail call double @llvm.sqrt.f64(double %104)
  %106 = fdiv double -1.000000e+00, %sqrt157
  %107 = select i1 %105, double %106, double 0xFFF0000000000000
  %108 = fneg double %107
  %109 = fmul double %103, %108
  br label %118

110:                                              ; preds = %98
  %111 = fdiv double %96, %91
  %112 = tail call double @llvm.fmuladd.f64(double %111, double %111, double 1.000000e+00)
  %113 = fcmp ugt double %112, 0.000000e+00
  %sqrt158 = tail call double @llvm.sqrt.f64(double %112)
  %114 = fdiv double 1.000000e+00, %sqrt158
  %115 = select i1 %113, double %114, double 0x7FF0000000000000
  %116 = fneg double %115
  %117 = fmul double %111, %116
  br label %118

118:                                              ; preds = %._crit_edge, %102, %110
  %.1147 = phi double [ %107, %102 ], [ %117, %110 ], [ 0.000000e+00, %._crit_edge ]
  %.1142 = phi double [ %109, %102 ], [ %115, %110 ], [ 1.000000e+00, %._crit_edge ]
  %119 = shl nsw i32 %66, 1
  %120 = sext i32 %119 to i64
  %121 = getelementptr inbounds double, ptr %2, i64 %120
  store double %.1142, ptr %121, align 8, !tbaa !9
  %122 = or disjoint i32 %119, 1
  %123 = sext i32 %122 to i64
  %124 = getelementptr inbounds double, ptr %2, i64 %123
  store double %.1147, ptr %124, align 8, !tbaa !9
  %125 = fneg double %.1147
  %126 = fmul double %96, %125
  %127 = tail call double @llvm.fmuladd.f64(double %.1142, double %91, double %126)
  store double %127, ptr %90, align 8, !tbaa !9
  %128 = fcmp oeq double %127, 0.000000e+00
  br i1 %128, label %129, label %.loopexit

129:                                              ; preds = %118
  br label %.loopexit

.loopexit:                                        ; preds = %56, %.preheader159, %118, %129
  %.2 = phi i32 [ %0, %129 ], [ 0, %118 ], [ 0, %.preheader159 ], [ %.1, %56 ]
  ret i32 %.2
}

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

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

attributes #0 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #1 = { nofree norecurse nosync nounwind memory(readwrite, inaccessiblemem: none) 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 = !{!"any pointer", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"double", !7, i64 0}
!11 = distinct !{!11, !12}
!12 = !{!"llvm.loop.mustprogress"}
!13 = distinct !{!13, !12}
!14 = distinct !{!14, !12}


; Effective O3 passes: AggressiveInstCombinePass,BDCEPass,CorrelatedValuePropagationPass,EarlyCSEPass,GVNPass,GlobalOptPass,IndVarSimplifyPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,JumpThreadingPass,LCSSAPass,LICMPass,LoopRotatePass,LoopSimplifyPass,LoopUnrollPass,PostOrderFunctionAttrsPass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass