; LLVM IR from chrr/libECP :: src (type1)

; 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.doubleArray = type { ptr, ptr }
%struct.QIntegrandParameters = type { double, ptr, ptr, ptr, ptr, double }
%struct.Integrand = type { ptr, ptr }
%struct.Type1 = type { ptr, ptr, ptr, ptr, ptr, i32, ptr, ptr, i32, ptr, i32, ptr, i32, i32, ptr, i32, ptr, double, double, double }
%struct.GCIntegrationTable = type { i32, i32, ptr, ptr, double, double, i32, i32 }
%struct._ECP = type { ptr, i32, i32 }

; 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

declare void @GCIntegrationTable_free(ptr noundef) #1

; Function Attrs: nounwind
declare void @free(ptr noundef) #2

; Function Attrs: nounwind uwtable
declare dso_local double @QIntegrand(double noundef, i32 noundef, ptr noundef) #3

; Function Attrs: nounwind
declare double @exp(double noundef) #2

; Function Attrs: nounwind uwtable
define dso_local { ptr, ptr } @calcQ(ptr noundef %0, ptr noundef %1, ptr noundef %2, double noundef %3, i32 noundef %4, double noundef %5, double noundef %6, double noundef %7, double noundef %8, double noundef %9, double noundef %10, double noundef %11, double noundef %12) #3 {
  %14 = alloca %struct.doubleArray, align 8
  %15 = alloca ptr, align 8
  %16 = alloca ptr, align 8
  %17 = alloca ptr, align 8
  %18 = alloca double, align 8
  %19 = alloca i32, align 4
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %25 = alloca double, align 8
  %26 = alloca double, align 8
  %27 = alloca double, align 8
  %28 = alloca ptr, align 8
  %29 = alloca %struct.QIntegrandParameters, align 8
  %30 = alloca %struct.Integrand, align 8
  %31 = alloca i32, align 4
  %32 = alloca i32, align 4
  %33 = alloca i32, align 4
  %34 = alloca double, align 8
  %35 = alloca double, align 8
  %36 = alloca double, align 8
  %37 = alloca double, align 8
  %38 = alloca ptr, align 8
  %39 = alloca ptr, align 8
  %40 = alloca ptr, align 8
  %41 = alloca i32, align 4
  %42 = alloca double, align 8
  %43 = alloca i32, align 4
  %44 = alloca ptr, align 8
  store ptr %0, ptr %15, align 8, !tbaa !5
  store ptr %1, ptr %16, align 8, !tbaa !5
  store ptr %2, ptr %17, align 8, !tbaa !5
  store double %3, ptr %18, align 8, !tbaa !9
  store i32 %4, ptr %19, align 4, !tbaa !11
  store double %5, ptr %20, align 8, !tbaa !9
  store double %6, ptr %21, align 8, !tbaa !9
  store double %7, ptr %22, align 8, !tbaa !9
  store double %8, ptr %23, align 8, !tbaa !9
  store double %9, ptr %24, align 8, !tbaa !9
  store double %10, ptr %25, align 8, !tbaa !9
  store double %11, ptr %26, align 8, !tbaa !9
  store double %12, ptr %27, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %28) #6
  call void @llvm.lifetime.start.p0(i64 48, ptr %29) #6
  call void @llvm.lifetime.start.p0(i64 16, ptr %30) #6
  %45 = getelementptr inbounds %struct.Integrand, ptr %30, i32 0, i32 0
  store ptr @QIntegrand, ptr %45, align 8, !tbaa !13
  %46 = getelementptr inbounds %struct.Integrand, ptr %30, i32 0, i32 1
  store ptr %29, ptr %46, align 8, !tbaa !15
  call void @llvm.lifetime.start.p0(i64 4, ptr %31) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %32) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %33) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %34) #6
  %47 = load double, ptr %22, align 8, !tbaa !9
  %48 = fneg double %47
  %49 = load double, ptr %20, align 8, !tbaa !9
  %50 = fmul double %48, %49
  %51 = load double, ptr %20, align 8, !tbaa !9
  %52 = load double, ptr %26, align 8, !tbaa !9
  %53 = load double, ptr %24, align 8, !tbaa !9
  %54 = fmul double %52, %53
  %55 = load double, ptr %24, align 8, !tbaa !9
  %56 = fmul double %54, %55
  %57 = fneg double %56
  %58 = call double @llvm.fmuladd.f64(double %50, double %51, double %57)
  store double %58, ptr %34, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %35) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %36) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %37) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %38) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %39) #6
  %59 = load i32, ptr %19, align 4, !tbaa !11
  %60 = add nsw i32 %59, 1
  %61 = load i32, ptr %19, align 4, !tbaa !11
  %62 = add nsw i32 %61, 1
  %63 = call { ptr, ptr } (i64, ...) @doubleArray_new(i64 noundef 2, i32 noundef %60, i32 noundef %62)
  %64 = getelementptr inbounds { ptr, ptr }, ptr %14, i32 0, i32 0
  %65 = extractvalue { ptr, ptr } %63, 0
  store ptr %65, ptr %64, align 8
  %66 = getelementptr inbounds { ptr, ptr }, ptr %14, i32 0, i32 1
  %67 = extractvalue { ptr, ptr } %63, 1
  store ptr %67, ptr %66, align 8
  call void @llvm.lifetime.start.p0(i64 8, ptr %40) #6
  %68 = getelementptr inbounds %struct.doubleArray, ptr %14, i32 0, i32 1
  %69 = load ptr, ptr %68, align 8, !tbaa !15
  store ptr %69, ptr %40, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %41) #6
  %70 = getelementptr inbounds %struct.doubleArray, ptr %14, i32 0, i32 0
  %71 = load ptr, ptr %70, align 8, !tbaa !13
  %72 = getelementptr inbounds i64, ptr %71, i64 2
  %73 = load i64, ptr %72, align 8, !tbaa !16
  %74 = trunc i64 %73 to i32
  store i32 %74, ptr %41, align 4, !tbaa !11
  call void @llvm.lifetime.start.p0(i64 8, ptr %42) #6
  %75 = load double, ptr %22, align 8, !tbaa !9
  %76 = fneg double %75
  %77 = load double, ptr %26, align 8, !tbaa !9
  %78 = fsub double %76, %77
  store double %78, ptr %42, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 4, ptr %43) #6
  store i32 0, ptr %43, align 4, !tbaa !11
  call void @llvm.lifetime.start.p0(i64 8, ptr %44) #6
  %79 = load i32, ptr %19, align 4, !tbaa !11
  %80 = add nsw i32 %79, 1
  %81 = mul nsw i32 2, %80
  %82 = load i32, ptr %19, align 4, !tbaa !11
  %83 = add nsw i32 %82, 1
  %84 = mul nsw i32 %81, %83
  %85 = sext i32 %84 to i64
  %86 = call noalias ptr @calloc(i64 noundef %85, i64 noundef 4) #7
  store ptr %86, ptr %44, align 8, !tbaa !5
  %87 = load double, ptr %21, align 8, !tbaa !9
  %88 = load double, ptr %25, align 8, !tbaa !9
  %89 = fmul double %87, %88
  %90 = load double, ptr %34, align 8, !tbaa !9
  %91 = call double @exp(double noundef %90) #6
  %92 = fmul double %89, %91
  %93 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 0
  store double %92, ptr %93, align 8, !tbaa !18
  %94 = load ptr, ptr %15, align 8, !tbaa !5
  %95 = getelementptr inbounds %struct.Type1, ptr %94, i32 0, i32 14
  %96 = load ptr, ptr %95, align 8, !tbaa !20
  %97 = getelementptr inbounds %struct.GCIntegrationTable, ptr %96, i32 0, i32 1
  %98 = load i32, ptr %97, align 4, !tbaa !22
  %99 = sext i32 %98 to i64
  %100 = call noalias ptr @calloc(i64 noundef %99, i64 noundef 8) #7
  %101 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 1
  store ptr %100, ptr %101, align 8, !tbaa !24
  %102 = load ptr, ptr %15, align 8, !tbaa !5
  %103 = getelementptr inbounds %struct.Type1, ptr %102, i32 0, i32 18
  %104 = load double, ptr %103, align 8, !tbaa !25
  %105 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 5
  store double %104, ptr %105, align 8, !tbaa !26
  %106 = load ptr, ptr %15, align 8, !tbaa !5
  %107 = getelementptr inbounds %struct.Type1, ptr %106, i32 0, i32 14
  %108 = load ptr, ptr %107, align 8, !tbaa !20
  %109 = getelementptr inbounds %struct.GCIntegrationTable, ptr %108, i32 0, i32 1
  %110 = load i32, ptr %109, align 4, !tbaa !22
  %111 = sext i32 %110 to i64
  %112 = call noalias ptr @calloc(i64 noundef %111, i64 noundef 8) #7
  %113 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 3
  store ptr %112, ptr %113, align 8, !tbaa !27
  %114 = load i32, ptr %19, align 4, !tbaa !11
  %115 = add nsw i32 %114, 1
  %116 = load ptr, ptr %15, align 8, !tbaa !5
  %117 = getelementptr inbounds %struct.Type1, ptr %116, i32 0, i32 14
  %118 = load ptr, ptr %117, align 8, !tbaa !20
  %119 = getelementptr inbounds %struct.GCIntegrationTable, ptr %118, i32 0, i32 1
  %120 = load i32, ptr %119, align 4, !tbaa !22
  %121 = mul nsw i32 %115, %120
  %122 = sext i32 %121 to i64
  %123 = call noalias ptr @calloc(i64 noundef %122, i64 noundef 8) #7
  store ptr %123, ptr %38, align 8, !tbaa !5
  %124 = load i32, ptr %19, align 4, !tbaa !11
  %125 = add nsw i32 %124, 1
  %126 = load ptr, ptr %15, align 8, !tbaa !5
  %127 = getelementptr inbounds %struct.Type1, ptr %126, i32 0, i32 14
  %128 = load ptr, ptr %127, align 8, !tbaa !20
  %129 = getelementptr inbounds %struct.GCIntegrationTable, ptr %128, i32 0, i32 1
  %130 = load i32, ptr %129, align 4, !tbaa !22
  %131 = mul nsw i32 %125, %130
  %132 = sext i32 %131 to i64
  %133 = call noalias ptr @calloc(i64 noundef %132, i64 noundef 8) #7
  store ptr %133, ptr %39, align 8, !tbaa !5
  %134 = load ptr, ptr %15, align 8, !tbaa !5
  %135 = getelementptr inbounds %struct.Type1, ptr %134, i32 0, i32 16
  %136 = load ptr, ptr %135, align 8, !tbaa !28
  %137 = getelementptr inbounds %struct.GCIntegrationTable, ptr %136, i32 0, i32 6
  %138 = load i32, ptr %137, align 8, !tbaa !29
  store i32 %138, ptr %33, align 4, !tbaa !11
  br label %139

139:                                              ; preds = %235, %13
  %140 = load i32, ptr %33, align 4, !tbaa !11
  %141 = load ptr, ptr %15, align 8, !tbaa !5
  %142 = getelementptr inbounds %struct.Type1, ptr %141, i32 0, i32 16
  %143 = load ptr, ptr %142, align 8, !tbaa !28
  %144 = getelementptr inbounds %struct.GCIntegrationTable, ptr %143, i32 0, i32 7
  %145 = load i32, ptr %144, align 4, !tbaa !30
  %146 = icmp slt i32 %140, %145
  br i1 %146, label %147, label %238

147:                                              ; preds = %139
  %148 = load ptr, ptr %15, align 8, !tbaa !5
  %149 = getelementptr inbounds %struct.Type1, ptr %148, i32 0, i32 16
  %150 = load ptr, ptr %149, align 8, !tbaa !28
  %151 = getelementptr inbounds %struct.GCIntegrationTable, ptr %150, i32 0, i32 2
  %152 = load ptr, ptr %151, align 8, !tbaa !31
  %153 = load i32, ptr %33, align 4, !tbaa !11
  %154 = sext i32 %153 to i64
  %155 = getelementptr inbounds double, ptr %152, i64 %154
  %156 = load double, ptr %155, align 8, !tbaa !9
  store double %156, ptr %37, align 8, !tbaa !9
  %157 = load double, ptr %42, align 8, !tbaa !9
  %158 = load double, ptr %37, align 8, !tbaa !9
  %159 = load double, ptr %18, align 8, !tbaa !9
  %160 = call double @llvm.fmuladd.f64(double %157, double %158, double %159)
  %161 = load double, ptr %37, align 8, !tbaa !9
  %162 = fmul double %160, %161
  %163 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 1
  %164 = load ptr, ptr %163, align 8, !tbaa !24
  %165 = load i32, ptr %33, align 4, !tbaa !11
  %166 = sext i32 %165 to i64
  %167 = getelementptr inbounds double, ptr %164, i64 %166
  store double %162, ptr %167, align 8, !tbaa !9
  %168 = load ptr, ptr %16, align 8, !tbaa !5
  %169 = load i32, ptr %33, align 4, !tbaa !11
  %170 = sext i32 %169 to i64
  %171 = getelementptr inbounds double, ptr %168, i64 %170
  %172 = load double, ptr %171, align 8, !tbaa !9
  %173 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 3
  %174 = load ptr, ptr %173, align 8, !tbaa !27
  %175 = load i32, ptr %33, align 4, !tbaa !11
  %176 = sext i32 %175 to i64
  %177 = getelementptr inbounds double, ptr %174, i64 %176
  store double %172, ptr %177, align 8, !tbaa !9
  %178 = load ptr, ptr %15, align 8, !tbaa !5
  %179 = getelementptr inbounds %struct.Type1, ptr %178, i32 0, i32 11
  %180 = load ptr, ptr %179, align 8, !tbaa !32
  %181 = load ptr, ptr %15, align 8, !tbaa !5
  %182 = getelementptr inbounds %struct.Type1, ptr %181, i32 0, i32 16
  %183 = load ptr, ptr %182, align 8, !tbaa !28
  %184 = getelementptr inbounds %struct.GCIntegrationTable, ptr %183, i32 0, i32 1
  %185 = load i32, ptr %184, align 4, !tbaa !22
  %186 = load i32, ptr %19, align 4, !tbaa !11
  %187 = load double, ptr %18, align 8, !tbaa !9
  %188 = load double, ptr %37, align 8, !tbaa !9
  %189 = fmul double %187, %188
  %190 = load ptr, ptr %38, align 8, !tbaa !5
  %191 = load i32, ptr %33, align 4, !tbaa !11
  %192 = sext i32 %191 to i64
  %193 = getelementptr inbounds double, ptr %190, i64 %192
  call void @weightedBesselFunction(ptr noundef %180, i32 noundef %185, i32 noundef %186, double noundef %189, ptr noundef %193)
  %194 = load ptr, ptr %39, align 8, !tbaa !5
  %195 = load i32, ptr %33, align 4, !tbaa !11
  %196 = sext i32 %195 to i64
  %197 = getelementptr inbounds double, ptr %194, i64 %196
  store double 1.000000e+00, ptr %197, align 8, !tbaa !9
  store i32 1, ptr %31, align 4, !tbaa !11
  br label %198

198:                                              ; preds = %231, %147
  %199 = load i32, ptr %31, align 4, !tbaa !11
  %200 = load i32, ptr %19, align 4, !tbaa !11
  %201 = icmp sle i32 %199, %200
  br i1 %201, label %202, label %234

202:                                              ; preds = %198
  %203 = load double, ptr %37, align 8, !tbaa !9
  %204 = load ptr, ptr %39, align 8, !tbaa !5
  %205 = load i32, ptr %31, align 4, !tbaa !11
  %206 = sub nsw i32 %205, 1
  %207 = load ptr, ptr %15, align 8, !tbaa !5
  %208 = getelementptr inbounds %struct.Type1, ptr %207, i32 0, i32 16
  %209 = load ptr, ptr %208, align 8, !tbaa !28
  %210 = getelementptr inbounds %struct.GCIntegrationTable, ptr %209, i32 0, i32 1
  %211 = load i32, ptr %210, align 4, !tbaa !22
  %212 = mul nsw i32 %206, %211
  %213 = load i32, ptr %33, align 4, !tbaa !11
  %214 = add nsw i32 %212, %213
  %215 = sext i32 %214 to i64
  %216 = getelementptr inbounds double, ptr %204, i64 %215
  %217 = load double, ptr %216, align 8, !tbaa !9
  %218 = fmul double %203, %217
  %219 = load ptr, ptr %39, align 8, !tbaa !5
  %220 = load i32, ptr %31, align 4, !tbaa !11
  %221 = load ptr, ptr %15, align 8, !tbaa !5
  %222 = getelementptr inbounds %struct.Type1, ptr %221, i32 0, i32 16
  %223 = load ptr, ptr %222, align 8, !tbaa !28
  %224 = getelementptr inbounds %struct.GCIntegrationTable, ptr %223, i32 0, i32 1
  %225 = load i32, ptr %224, align 4, !tbaa !22
  %226 = mul nsw i32 %220, %225
  %227 = load i32, ptr %33, align 4, !tbaa !11
  %228 = add nsw i32 %226, %227
  %229 = sext i32 %228 to i64
  %230 = getelementptr inbounds double, ptr %219, i64 %229
  store double %218, ptr %230, align 8, !tbaa !9
  br label %231

231:                                              ; preds = %202
  %232 = load i32, ptr %31, align 4, !tbaa !11
  %233 = add nsw i32 %232, 1
  store i32 %233, ptr %31, align 4, !tbaa !11
  br label %198, !llvm.loop !33

234:                                              ; preds = %198
  br label %235

235:                                              ; preds = %234
  %236 = load i32, ptr %33, align 4, !tbaa !11
  %237 = add nsw i32 %236, 1
  store i32 %237, ptr %33, align 4, !tbaa !11
  br label %139, !llvm.loop !35

238:                                              ; preds = %139
  store i32 0, ptr %31, align 4, !tbaa !11
  br label %239

239:                                              ; preds = %313, %238
  %240 = load i32, ptr %31, align 4, !tbaa !11
  %241 = load i32, ptr %19, align 4, !tbaa !11
  %242 = icmp sle i32 %240, %241
  br i1 %242, label %243, label %316

243:                                              ; preds = %239
  %244 = load ptr, ptr %39, align 8, !tbaa !5
  %245 = load i32, ptr %31, align 4, !tbaa !11
  %246 = load ptr, ptr %15, align 8, !tbaa !5
  %247 = getelementptr inbounds %struct.Type1, ptr %246, i32 0, i32 16
  %248 = load ptr, ptr %247, align 8, !tbaa !28
  %249 = getelementptr inbounds %struct.GCIntegrationTable, ptr %248, i32 0, i32 1
  %250 = load i32, ptr %249, align 4, !tbaa !22
  %251 = mul nsw i32 %245, %250
  %252 = sext i32 %251 to i64
  %253 = getelementptr inbounds double, ptr %244, i64 %252
  %254 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 4
  store ptr %253, ptr %254, align 8, !tbaa !36
  %255 = load i32, ptr %31, align 4, !tbaa !11
  store i32 %255, ptr %32, align 4, !tbaa !11
  br label %256

256:                                              ; preds = %309, %243
  %257 = load i32, ptr %32, align 4, !tbaa !11
  %258 = icmp sge i32 %257, 0
  br i1 %258, label %259, label %312

259:                                              ; preds = %256
  %260 = load ptr, ptr %38, align 8, !tbaa !5
  %261 = load i32, ptr %32, align 4, !tbaa !11
  %262 = load ptr, ptr %15, align 8, !tbaa !5
  %263 = getelementptr inbounds %struct.Type1, ptr %262, i32 0, i32 16
  %264 = load ptr, ptr %263, align 8, !tbaa !28
  %265 = getelementptr inbounds %struct.GCIntegrationTable, ptr %264, i32 0, i32 1
  %266 = load i32, ptr %265, align 4, !tbaa !22
  %267 = mul nsw i32 %261, %266
  %268 = sext i32 %267 to i64
  %269 = getelementptr inbounds double, ptr %260, i64 %268
  %270 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 2
  store ptr %269, ptr %270, align 8, !tbaa !37
  %271 = load ptr, ptr %15, align 8, !tbaa !5
  %272 = getelementptr inbounds %struct.Type1, ptr %271, i32 0, i32 16
  %273 = load ptr, ptr %272, align 8, !tbaa !28
  %274 = call i32 @integrateGC_PS93(ptr noundef %30, ptr noundef %273)
  %275 = icmp ne i32 0, %274
  br i1 %275, label %276, label %292

276:                                              ; preds = %259
  %277 = load i32, ptr %31, align 4, !tbaa !11
  %278 = load ptr, ptr %44, align 8, !tbaa !5
  %279 = load i32, ptr %43, align 4, !tbaa !11
  %280 = mul nsw i32 %279, 2
  %281 = sext i32 %280 to i64
  %282 = getelementptr inbounds i32, ptr %278, i64 %281
  store i32 %277, ptr %282, align 4, !tbaa !11
  %283 = load i32, ptr %32, align 4, !tbaa !11
  %284 = load ptr, ptr %44, align 8, !tbaa !5
  %285 = load i32, ptr %43, align 4, !tbaa !11
  %286 = mul nsw i32 %285, 2
  %287 = add nsw i32 %286, 1
  %288 = sext i32 %287 to i64
  %289 = getelementptr inbounds i32, ptr %284, i64 %288
  store i32 %283, ptr %289, align 4, !tbaa !11
  %290 = load i32, ptr %43, align 4, !tbaa !11
  %291 = add nsw i32 %290, 1
  store i32 %291, ptr %43, align 4, !tbaa !11
  br label %308

292:                                              ; preds = %259
  %293 = load ptr, ptr %15, align 8, !tbaa !5
  %294 = getelementptr inbounds %struct.Type1, ptr %293, i32 0, i32 16
  %295 = load ptr, ptr %294, align 8, !tbaa !28
  %296 = getelementptr inbounds %struct.GCIntegrationTable, ptr %295, i32 0, i32 5
  %297 = load double, ptr %296, align 8, !tbaa !38
  %298 = load ptr, ptr %40, align 8, !tbaa !5
  %299 = load i32, ptr %31, align 4, !tbaa !11
  %300 = load i32, ptr %41, align 4, !tbaa !11
  %301 = mul nsw i32 %299, %300
  %302 = load i32, ptr %32, align 4, !tbaa !11
  %303 = add nsw i32 %301, %302
  %304 = sext i32 %303 to i64
  %305 = getelementptr inbounds double, ptr %298, i64 %304
  %306 = load double, ptr %305, align 8, !tbaa !9
  %307 = fadd double %306, %297
  store double %307, ptr %305, align 8, !tbaa !9
  br label %308

308:                                              ; preds = %292, %276
  br label %309

309:                                              ; preds = %308
  %310 = load i32, ptr %32, align 4, !tbaa !11
  %311 = sub nsw i32 %310, 2
  store i32 %311, ptr %32, align 4, !tbaa !11
  br label %256, !llvm.loop !39

312:                                              ; preds = %256
  br label %313

313:                                              ; preds = %312
  %314 = load i32, ptr %31, align 4, !tbaa !11
  %315 = add nsw i32 %314, 1
  store i32 %315, ptr %31, align 4, !tbaa !11
  br label %239, !llvm.loop !40

316:                                              ; preds = %239
  %317 = load i32, ptr %43, align 4, !tbaa !11
  %318 = icmp slt i32 0, %317
  br i1 %318, label %319, label %540

319:                                              ; preds = %316
  %320 = load double, ptr %21, align 8, !tbaa !9
  %321 = load double, ptr %25, align 8, !tbaa !9
  %322 = fmul double %320, %321
  %323 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 0
  store double %322, ptr %323, align 8, !tbaa !18
  %324 = load ptr, ptr %15, align 8, !tbaa !5
  %325 = getelementptr inbounds %struct.Type1, ptr %324, i32 0, i32 14
  %326 = load ptr, ptr %325, align 8, !tbaa !20
  %327 = call ptr @GCIntegrationTable_copy(ptr noundef %326)
  store ptr %327, ptr %28, align 8, !tbaa !5
  %328 = load double, ptr %22, align 8, !tbaa !9
  %329 = load double, ptr %26, align 8, !tbaa !9
  %330 = fadd double %328, %329
  store double %330, ptr %35, align 8, !tbaa !9
  %331 = load double, ptr %22, align 8, !tbaa !9
  %332 = load double, ptr %20, align 8, !tbaa !9
  %333 = load double, ptr %26, align 8, !tbaa !9
  %334 = load double, ptr %24, align 8, !tbaa !9
  %335 = fmul double %333, %334
  %336 = call double @llvm.fmuladd.f64(double %331, double %332, double %335)
  %337 = load double, ptr %35, align 8, !tbaa !9
  %338 = fdiv double %336, %337
  store double %338, ptr %36, align 8, !tbaa !9
  %339 = load ptr, ptr %28, align 8, !tbaa !5
  %340 = getelementptr inbounds %struct.GCIntegrationTable, ptr %339, i32 0, i32 1
  %341 = load i32, ptr %340, align 4, !tbaa !22
  %342 = load ptr, ptr %28, align 8, !tbaa !5
  %343 = getelementptr inbounds %struct.GCIntegrationTable, ptr %342, i32 0, i32 2
  %344 = load ptr, ptr %343, align 8, !tbaa !31
  %345 = load ptr, ptr %28, align 8, !tbaa !5
  %346 = getelementptr inbounds %struct.GCIntegrationTable, ptr %345, i32 0, i32 3
  %347 = load ptr, ptr %346, align 8, !tbaa !41
  %348 = load double, ptr %35, align 8, !tbaa !9
  %349 = load double, ptr %36, align 8, !tbaa !9
  call void @transformGrid_FM06(i32 noundef %341, ptr noundef %344, ptr noundef %347, double noundef %348, double noundef %349)
  store i32 0, ptr %33, align 4, !tbaa !11
  br label %350

350:                                              ; preds = %473, %319
  %351 = load i32, ptr %33, align 4, !tbaa !11
  %352 = load ptr, ptr %28, align 8, !tbaa !5
  %353 = getelementptr inbounds %struct.GCIntegrationTable, ptr %352, i32 0, i32 1
  %354 = load i32, ptr %353, align 4, !tbaa !22
  %355 = icmp slt i32 %351, %354
  br i1 %355, label %356, label %476

356:                                              ; preds = %350
  %357 = load ptr, ptr %28, align 8, !tbaa !5
  %358 = getelementptr inbounds %struct.GCIntegrationTable, ptr %357, i32 0, i32 2
  %359 = load ptr, ptr %358, align 8, !tbaa !31
  %360 = load i32, ptr %33, align 4, !tbaa !11
  %361 = sext i32 %360 to i64
  %362 = getelementptr inbounds double, ptr %359, i64 %361
  %363 = load double, ptr %362, align 8, !tbaa !9
  store double %363, ptr %37, align 8, !tbaa !9
  %364 = load double, ptr %42, align 8, !tbaa !9
  %365 = load double, ptr %37, align 8, !tbaa !9
  %366 = load double, ptr %18, align 8, !tbaa !9
  %367 = call double @llvm.fmuladd.f64(double %364, double %365, double %366)
  %368 = load double, ptr %37, align 8, !tbaa !9
  %369 = load double, ptr %34, align 8, !tbaa !9
  %370 = call double @llvm.fmuladd.f64(double %367, double %368, double %369)
  %371 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 1
  %372 = load ptr, ptr %371, align 8, !tbaa !24
  %373 = load i32, ptr %33, align 4, !tbaa !11
  %374 = sext i32 %373 to i64
  %375 = getelementptr inbounds double, ptr %372, i64 %374
  store double %370, ptr %375, align 8, !tbaa !9
  %376 = load double, ptr %37, align 8, !tbaa !9
  %377 = load double, ptr %20, align 8, !tbaa !9
  %378 = fcmp ogt double %376, %377
  br i1 %378, label %379, label %398

379:                                              ; preds = %356
  %380 = load double, ptr %37, align 8, !tbaa !9
  %381 = load double, ptr %24, align 8, !tbaa !9
  %382 = fcmp ogt double %380, %381
  br i1 %382, label %383, label %398

383:                                              ; preds = %379
  %384 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 1
  %385 = load ptr, ptr %384, align 8, !tbaa !24
  %386 = load i32, ptr %33, align 4, !tbaa !11
  %387 = sext i32 %386 to i64
  %388 = getelementptr inbounds double, ptr %385, i64 %387
  %389 = load double, ptr %388, align 8, !tbaa !9
  %390 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 5
  %391 = load double, ptr %390, align 8, !tbaa !26
  %392 = fcmp olt double %389, %391
  br i1 %392, label %393, label %398

393:                                              ; preds = %383
  %394 = load i32, ptr %33, align 4, !tbaa !11
  %395 = sub nsw i32 %394, 1
  %396 = load ptr, ptr %28, align 8, !tbaa !5
  %397 = getelementptr inbounds %struct.GCIntegrationTable, ptr %396, i32 0, i32 7
  store i32 %395, ptr %397, align 4, !tbaa !30
  br label %476

398:                                              ; preds = %383, %379, %356
  %399 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 1
  %400 = load ptr, ptr %399, align 8, !tbaa !24
  %401 = load i32, ptr %33, align 4, !tbaa !11
  %402 = sext i32 %401 to i64
  %403 = getelementptr inbounds double, ptr %400, i64 %402
  %404 = load double, ptr %403, align 8, !tbaa !9
  %405 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 5
  %406 = load double, ptr %405, align 8, !tbaa !26
  %407 = fcmp oge double %404, %406
  br i1 %407, label %408, label %471

408:                                              ; preds = %398
  %409 = load ptr, ptr %15, align 8, !tbaa !5
  %410 = getelementptr inbounds %struct.Type1, ptr %409, i32 0, i32 11
  %411 = load ptr, ptr %410, align 8, !tbaa !32
  %412 = load ptr, ptr %28, align 8, !tbaa !5
  %413 = getelementptr inbounds %struct.GCIntegrationTable, ptr %412, i32 0, i32 1
  %414 = load i32, ptr %413, align 4, !tbaa !22
  %415 = load i32, ptr %19, align 4, !tbaa !11
  %416 = load double, ptr %18, align 8, !tbaa !9
  %417 = load double, ptr %37, align 8, !tbaa !9
  %418 = fmul double %416, %417
  %419 = load ptr, ptr %38, align 8, !tbaa !5
  %420 = load i32, ptr %33, align 4, !tbaa !11
  %421 = sext i32 %420 to i64
  %422 = getelementptr inbounds double, ptr %419, i64 %421
  call void @weightedBesselFunction(ptr noundef %411, i32 noundef %414, i32 noundef %415, double noundef %418, ptr noundef %422)
  %423 = load ptr, ptr %17, align 8, !tbaa !5
  %424 = load ptr, ptr %17, align 8, !tbaa !5
  %425 = getelementptr inbounds %struct._ECP, ptr %424, i32 0, i32 2
  %426 = load i32, ptr %425, align 4, !tbaa !42
  %427 = load double, ptr %37, align 8, !tbaa !9
  %428 = call double @evalECP(ptr noundef %423, i32 noundef %426, double noundef %427)
  %429 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 3
  %430 = load ptr, ptr %429, align 8, !tbaa !27
  %431 = load i32, ptr %33, align 4, !tbaa !11
  %432 = sext i32 %431 to i64
  %433 = getelementptr inbounds double, ptr %430, i64 %432
  store double %428, ptr %433, align 8, !tbaa !9
  %434 = load ptr, ptr %39, align 8, !tbaa !5
  %435 = load i32, ptr %33, align 4, !tbaa !11
  %436 = sext i32 %435 to i64
  %437 = getelementptr inbounds double, ptr %434, i64 %436
  store double 1.000000e+00, ptr %437, align 8, !tbaa !9
  store i32 1, ptr %31, align 4, !tbaa !11
  br label %438

438:                                              ; preds = %467, %408
  %439 = load i32, ptr %31, align 4, !tbaa !11
  %440 = load i32, ptr %19, align 4, !tbaa !11
  %441 = icmp sle i32 %439, %440
  br i1 %441, label %442, label %470

442:                                              ; preds = %438
  %443 = load double, ptr %37, align 8, !tbaa !9
  %444 = load ptr, ptr %39, align 8, !tbaa !5
  %445 = load i32, ptr %31, align 4, !tbaa !11
  %446 = sub nsw i32 %445, 1
  %447 = load ptr, ptr %28, align 8, !tbaa !5
  %448 = getelementptr inbounds %struct.GCIntegrationTable, ptr %447, i32 0, i32 1
  %449 = load i32, ptr %448, align 4, !tbaa !22
  %450 = mul nsw i32 %446, %449
  %451 = load i32, ptr %33, align 4, !tbaa !11
  %452 = add nsw i32 %450, %451
  %453 = sext i32 %452 to i64
  %454 = getelementptr inbounds double, ptr %444, i64 %453
  %455 = load double, ptr %454, align 8, !tbaa !9
  %456 = fmul double %443, %455
  %457 = load ptr, ptr %39, align 8, !tbaa !5
  %458 = load i32, ptr %31, align 4, !tbaa !11
  %459 = load ptr, ptr %28, align 8, !tbaa !5
  %460 = getelementptr inbounds %struct.GCIntegrationTable, ptr %459, i32 0, i32 1
  %461 = load i32, ptr %460, align 4, !tbaa !22
  %462 = mul nsw i32 %458, %461
  %463 = load i32, ptr %33, align 4, !tbaa !11
  %464 = add nsw i32 %462, %463
  %465 = sext i32 %464 to i64
  %466 = getelementptr inbounds double, ptr %457, i64 %465
  store double %456, ptr %466, align 8, !tbaa !9
  br label %467

467:                                              ; preds = %442
  %468 = load i32, ptr %31, align 4, !tbaa !11
  %469 = add nsw i32 %468, 1
  store i32 %469, ptr %31, align 4, !tbaa !11
  br label %438, !llvm.loop !44

470:                                              ; preds = %438
  br label %471

471:                                              ; preds = %470, %398
  br label %472

472:                                              ; preds = %471
  br label %473

473:                                              ; preds = %472
  %474 = load i32, ptr %33, align 4, !tbaa !11
  %475 = add nsw i32 %474, 1
  store i32 %475, ptr %33, align 4, !tbaa !11
  br label %350, !llvm.loop !45

476:                                              ; preds = %393, %350
  store i32 0, ptr %33, align 4, !tbaa !11
  br label %477

477:                                              ; preds = %536, %476
  %478 = load i32, ptr %33, align 4, !tbaa !11
  %479 = load i32, ptr %43, align 4, !tbaa !11
  %480 = icmp slt i32 %478, %479
  br i1 %480, label %481, label %539

481:                                              ; preds = %477
  %482 = load ptr, ptr %44, align 8, !tbaa !5
  %483 = load i32, ptr %33, align 4, !tbaa !11
  %484 = mul nsw i32 %483, 2
  %485 = sext i32 %484 to i64
  %486 = getelementptr inbounds i32, ptr %482, i64 %485
  %487 = load i32, ptr %486, align 4, !tbaa !11
  store i32 %487, ptr %31, align 4, !tbaa !11
  %488 = load ptr, ptr %44, align 8, !tbaa !5
  %489 = load i32, ptr %33, align 4, !tbaa !11
  %490 = mul nsw i32 %489, 2
  %491 = add nsw i32 %490, 1
  %492 = sext i32 %491 to i64
  %493 = getelementptr inbounds i32, ptr %488, i64 %492
  %494 = load i32, ptr %493, align 4, !tbaa !11
  store i32 %494, ptr %32, align 4, !tbaa !11
  %495 = load ptr, ptr %39, align 8, !tbaa !5
  %496 = load i32, ptr %31, align 4, !tbaa !11
  %497 = load ptr, ptr %28, align 8, !tbaa !5
  %498 = getelementptr inbounds %struct.GCIntegrationTable, ptr %497, i32 0, i32 1
  %499 = load i32, ptr %498, align 4, !tbaa !22
  %500 = mul nsw i32 %496, %499
  %501 = sext i32 %500 to i64
  %502 = getelementptr inbounds double, ptr %495, i64 %501
  %503 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 4
  store ptr %502, ptr %503, align 8, !tbaa !36
  %504 = load ptr, ptr %38, align 8, !tbaa !5
  %505 = load i32, ptr %32, align 4, !tbaa !11
  %506 = load ptr, ptr %28, align 8, !tbaa !5
  %507 = getelementptr inbounds %struct.GCIntegrationTable, ptr %506, i32 0, i32 1
  %508 = load i32, ptr %507, align 4, !tbaa !22
  %509 = mul nsw i32 %505, %508
  %510 = sext i32 %509 to i64
  %511 = getelementptr inbounds double, ptr %504, i64 %510
  %512 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 2
  store ptr %511, ptr %512, align 8, !tbaa !37
  %513 = load ptr, ptr %28, align 8, !tbaa !5
  %514 = call i32 @integrateGC_PSM92(ptr noundef %30, ptr noundef %513)
  %515 = icmp ne i32 0, %514
  br i1 %515, label %516, label %521

516:                                              ; preds = %481
  %517 = getelementptr inbounds { ptr, ptr }, ptr %14, i32 0, i32 0
  %518 = load ptr, ptr %517, align 8
  %519 = getelementptr inbounds { ptr, ptr }, ptr %14, i32 0, i32 1
  %520 = load ptr, ptr %519, align 8
  call void @doubleArray_free(ptr %518, ptr %520)
  br label %541

521:                                              ; preds = %481
  %522 = load ptr, ptr %28, align 8, !tbaa !5
  %523 = getelementptr inbounds %struct.GCIntegrationTable, ptr %522, i32 0, i32 5
  %524 = load double, ptr %523, align 8, !tbaa !38
  %525 = load ptr, ptr %40, align 8, !tbaa !5
  %526 = load i32, ptr %31, align 4, !tbaa !11
  %527 = load i32, ptr %41, align 4, !tbaa !11
  %528 = mul nsw i32 %526, %527
  %529 = load i32, ptr %32, align 4, !tbaa !11
  %530 = add nsw i32 %528, %529
  %531 = sext i32 %530 to i64
  %532 = getelementptr inbounds double, ptr %525, i64 %531
  %533 = load double, ptr %532, align 8, !tbaa !9
  %534 = fadd double %533, %524
  store double %534, ptr %532, align 8, !tbaa !9
  br label %535

535:                                              ; preds = %521
  br label %536

536:                                              ; preds = %535
  %537 = load i32, ptr %33, align 4, !tbaa !11
  %538 = add nsw i32 %537, 1
  store i32 %538, ptr %33, align 4, !tbaa !11
  br label %477, !llvm.loop !46

539:                                              ; preds = %477
  br label %540

540:                                              ; preds = %539, %316
  br label %541

541:                                              ; preds = %540, %516
  %542 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 1
  %543 = load ptr, ptr %542, align 8, !tbaa !24
  call void @free(ptr noundef %543) #6
  %544 = getelementptr inbounds %struct.QIntegrandParameters, ptr %29, i32 0, i32 3
  %545 = load ptr, ptr %544, align 8, !tbaa !27
  call void @free(ptr noundef %545) #6
  %546 = load ptr, ptr %38, align 8, !tbaa !5
  call void @free(ptr noundef %546) #6
  %547 = load ptr, ptr %39, align 8, !tbaa !5
  call void @free(ptr noundef %547) #6
  %548 = load ptr, ptr %44, align 8, !tbaa !5
  call void @free(ptr noundef %548) #6
  %549 = load i32, ptr %43, align 4, !tbaa !11
  %550 = icmp slt i32 0, %549
  br i1 %550, label %551, label %553

551:                                              ; preds = %541
  %552 = load ptr, ptr %28, align 8, !tbaa !5
  call void @GCIntegrationTable_free(ptr noundef %552)
  br label %553

553:                                              ; preds = %551, %541
  call void @llvm.lifetime.end.p0(i64 8, ptr %44) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %43) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %42) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %41) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %40) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %39) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %38) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %37) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %36) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %35) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %34) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %33) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %32) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %31) #6
  call void @llvm.lifetime.end.p0(i64 16, ptr %30) #6
  call void @llvm.lifetime.end.p0(i64 48, ptr %29) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #6
  %554 = load { ptr, ptr }, ptr %14, align 8
  ret { ptr, ptr } %554
}

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

declare { ptr, ptr } @doubleArray_new(i64 noundef, ...) #1

; Function Attrs: nounwind allocsize(0,1)
declare noalias ptr @calloc(i64 noundef, i64 noundef) #5

declare void @weightedBesselFunction(ptr noundef, i32 noundef, i32 noundef, double noundef, ptr noundef) #1

declare i32 @integrateGC_PS93(ptr noundef, ptr noundef) #1

declare ptr @GCIntegrationTable_copy(ptr noundef) #1

declare void @transformGrid_FM06(i32 noundef, ptr noundef, ptr noundef, double noundef, double noundef) #1

declare double @evalECP(ptr noundef, i32 noundef, double noundef) #1

declare i32 @integrateGC_PSM92(ptr noundef, ptr noundef) #1

declare void @doubleArray_free(ptr, ptr) #1

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { "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 = { 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 #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 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #5 = { nounwind allocsize(0,1) "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 #6 = { nounwind }
attributes #7 = { nounwind allocsize(0,1) }

!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 = !{!12, !12, i64 0}
!12 = !{!"int", !7, i64 0}
!13 = !{!14, !6, i64 0}
!14 = !{!"", !6, i64 0, !6, i64 8}
!15 = !{!14, !6, i64 8}
!16 = !{!17, !17, i64 0}
!17 = !{!"long", !7, i64 0}
!18 = !{!19, !10, i64 0}
!19 = !{!"", !10, i64 0, !6, i64 8, !6, i64 16, !6, i64 24, !6, i64 32, !10, i64 40}
!20 = !{!21, !6, i64 104}
!21 = !{!"", !6, i64 0, !6, i64 8, !6, i64 16, !6, i64 24, !6, i64 32, !12, i64 40, !6, i64 48, !6, i64 56, !12, i64 64, !6, i64 72, !12, i64 80, !6, i64 88, !12, i64 96, !12, i64 100, !6, i64 104, !12, i64 112, !6, i64 120, !10, i64 128, !10, i64 136, !10, i64 144}
!22 = !{!23, !12, i64 4}
!23 = !{!"", !12, i64 0, !12, i64 4, !6, i64 8, !6, i64 16, !10, i64 24, !10, i64 32, !12, i64 40, !12, i64 44}
!24 = !{!19, !6, i64 8}
!25 = !{!21, !10, i64 136}
!26 = !{!19, !10, i64 40}
!27 = !{!19, !6, i64 24}
!28 = !{!21, !6, i64 120}
!29 = !{!23, !12, i64 40}
!30 = !{!23, !12, i64 44}
!31 = !{!23, !6, i64 8}
!32 = !{!21, !6, i64 88}
!33 = distinct !{!33, !34}
!34 = !{!"llvm.loop.mustprogress"}
!35 = distinct !{!35, !34}
!36 = !{!19, !6, i64 32}
!37 = !{!19, !6, i64 16}
!38 = !{!23, !10, i64 32}
!39 = distinct !{!39, !34}
!40 = distinct !{!40, !34}
!41 = !{!23, !6, i64 16}
!42 = !{!43, !12, i64 12}
!43 = !{!"_ECP", !6, i64 0, !12, i64 8, !12, i64 12}
!44 = distinct !{!44, !34}
!45 = distinct !{!45, !34}
!46 = distinct !{!46, !34}


; 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.QIntegrandParameters = type { double, ptr, ptr, ptr, ptr, double }
%struct.Integrand = type { ptr, ptr }

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

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

declare void @GCIntegrationTable_free(ptr noundef) local_unnamed_addr #1

; Function Attrs: mustprogress nounwind willreturn allockind("free") memory(argmem: readwrite, inaccessiblemem: readwrite)
declare void @free(ptr allocptr nocapture noundef) local_unnamed_addr #2

; Function Attrs: nounwind uwtable
declare dso_local double @QIntegrand(double noundef, i32 noundef, ptr noundef) #3

; Function Attrs: mustprogress nofree nounwind willreturn memory(write)
declare double @exp(double noundef) local_unnamed_addr #4

; Function Attrs: nounwind uwtable
define dso_local { ptr, ptr } @calcQ(ptr nocapture noundef readonly %0, ptr nocapture noundef readonly %1, ptr noundef %2, double noundef %3, i32 noundef %4, double noundef %5, double noundef %6, double noundef %7, double noundef %8, double noundef %9, double noundef %10, double noundef %11, double noundef %12) local_unnamed_addr #3 {
  %14 = alloca %struct.QIntegrandParameters, align 8
  %15 = alloca %struct.Integrand, align 8
  call void @llvm.lifetime.start.p0(i64 48, ptr nonnull %14) #8
  call void @llvm.lifetime.start.p0(i64 16, ptr nonnull %15) #8
  store ptr @QIntegrand, ptr %15, align 8, !tbaa !5
  %16 = getelementptr inbounds i8, ptr %15, i64 8
  store ptr %14, ptr %16, align 8, !tbaa !10
  %17 = fneg double %7
  %18 = fmul double %17, %5
  %19 = fmul double %9, %11
  %20 = fneg double %19
  %21 = fmul double %20, %9
  %22 = call double @llvm.fmuladd.f64(double %18, double %5, double %21)
  %23 = add i32 %4, 1
  %24 = call { ptr, ptr } (i64, ...) @doubleArray_new(i64 noundef 2, i32 noundef %23, i32 noundef %23) #8
  %25 = extractvalue { ptr, ptr } %24, 0
  %26 = extractvalue { ptr, ptr } %24, 1
  %27 = getelementptr inbounds i8, ptr %25, i64 16
  %28 = load i64, ptr %27, align 8, !tbaa !11
  %29 = trunc i64 %28 to i32
  %30 = fsub double %17, %11
  %31 = shl nsw i32 %23, 1
  %32 = mul nsw i32 %31, %23
  %33 = sext i32 %32 to i64
  %34 = call noalias ptr @calloc(i64 noundef %33, i64 noundef 4) #9
  %35 = fmul double %6, %10
  %36 = call double @exp(double noundef %22) #8
  %37 = fmul double %35, %36
  store double %37, ptr %14, align 8, !tbaa !13
  %38 = getelementptr inbounds i8, ptr %0, i64 104
  %39 = load ptr, ptr %38, align 8, !tbaa !16
  %40 = getelementptr inbounds i8, ptr %39, i64 4
  %41 = load i32, ptr %40, align 4, !tbaa !19
  %42 = sext i32 %41 to i64
  %43 = call noalias ptr @calloc(i64 noundef %42, i64 noundef 8) #9
  %44 = getelementptr inbounds i8, ptr %14, i64 8
  store ptr %43, ptr %44, align 8, !tbaa !21
  %45 = getelementptr inbounds i8, ptr %0, i64 136
  %46 = load double, ptr %45, align 8, !tbaa !22
  %47 = getelementptr inbounds i8, ptr %14, i64 40
  store double %46, ptr %47, align 8, !tbaa !23
  %48 = call noalias ptr @calloc(i64 noundef %42, i64 noundef 8) #9
  %49 = getelementptr inbounds i8, ptr %14, i64 24
  store ptr %48, ptr %49, align 8, !tbaa !24
  %50 = mul nsw i32 %41, %23
  %51 = sext i32 %50 to i64
  %52 = call noalias ptr @calloc(i64 noundef %51, i64 noundef 8) #9
  %53 = call noalias ptr @calloc(i64 noundef %51, i64 noundef 8) #9
  %54 = getelementptr inbounds i8, ptr %0, i64 120
  %55 = load ptr, ptr %54, align 8, !tbaa !25
  %56 = getelementptr inbounds i8, ptr %55, i64 40
  %57 = load i32, ptr %56, align 8, !tbaa !26
  %58 = getelementptr inbounds i8, ptr %55, i64 44
  %59 = load i32, ptr %58, align 4, !tbaa !27
  %60 = icmp slt i32 %57, %59
  br i1 %60, label %.lr.ph190, label %.preheader

.lr.ph190:                                        ; preds = %13
  %61 = getelementptr inbounds i8, ptr %0, i64 88
  %.not180186 = icmp slt i32 %4, 1
  %62 = sext i32 %57 to i64
  %wide.trip.count = zext i32 %23 to i64
  %63 = add nsw i64 %wide.trip.count, -1
  %64 = add nsw i64 %wide.trip.count, -2
  %xtraiter = and i64 %63, 1
  %65 = icmp eq i64 %64, 0
  %unroll_iter = and i64 %63, -2
  %lcmp.mod.not = icmp eq i64 %xtraiter, 0
  %xtraiter262 = and i64 %63, 3
  %66 = icmp ult i64 %64, 3
  %unroll_iter264 = and i64 %63, -4
  %lcmp.mod263.not = icmp eq i64 %xtraiter262, 0
  br label %70

.preheader:                                       ; preds = %._crit_edge, %13
  %.not196 = icmp slt i32 %4, 0
  br i1 %.not196, label %.loopexit, label %.lr.ph199

.lr.ph199:                                        ; preds = %.preheader
  %67 = getelementptr inbounds i8, ptr %14, i64 32
  %68 = getelementptr inbounds i8, ptr %14, i64 16
  %sext = shl i64 %28, 32
  %69 = ashr exact i64 %sext, 32
  %wide.trip.count224 = zext i32 %23 to i64
  br label %.lr.ph194

70:                                               ; preds = %.lr.ph190, %._crit_edge
  %indvar = phi i64 [ 0, %.lr.ph190 ], [ %indvar.next, %._crit_edge ]
  %indvars.iv215 = phi i64 [ %62, %.lr.ph190 ], [ %indvars.iv.next216, %._crit_edge ]
  %71 = phi ptr [ %55, %.lr.ph190 ], [ %.pre, %._crit_edge ]
  %72 = add i64 %indvar, %62
  %73 = shl i64 %72, 3
  %scevgep256 = getelementptr i8, ptr %53, i64 %73
  %74 = getelementptr inbounds i8, ptr %71, i64 8
  %75 = load ptr, ptr %74, align 8, !tbaa !28
  %76 = getelementptr inbounds double, ptr %75, i64 %indvars.iv215
  %77 = load double, ptr %76, align 8, !tbaa !29
  %78 = call double @llvm.fmuladd.f64(double %30, double %77, double %3)
  %79 = fmul double %77, %78
  %80 = load ptr, ptr %44, align 8, !tbaa !21
  %81 = getelementptr inbounds double, ptr %80, i64 %indvars.iv215
  store double %79, ptr %81, align 8, !tbaa !29
  %82 = getelementptr inbounds double, ptr %1, i64 %indvars.iv215
  %83 = load double, ptr %82, align 8, !tbaa !29
  %84 = load ptr, ptr %49, align 8, !tbaa !24
  %85 = getelementptr inbounds double, ptr %84, i64 %indvars.iv215
  store double %83, ptr %85, align 8, !tbaa !29
  %86 = load ptr, ptr %61, align 8, !tbaa !30
  %87 = getelementptr inbounds i8, ptr %71, i64 4
  %88 = load i32, ptr %87, align 4, !tbaa !19
  %89 = fmul double %77, %3
  %90 = getelementptr inbounds double, ptr %52, i64 %indvars.iv215
  call void @weightedBesselFunction(ptr noundef %86, i32 noundef %88, i32 noundef %4, double noundef %89, ptr noundef %90) #8
  %91 = getelementptr inbounds double, ptr %53, i64 %indvars.iv215
  store double 1.000000e+00, ptr %91, align 8, !tbaa !29
  %.pre = load ptr, ptr %54, align 8, !tbaa !25
  br i1 %.not180186, label %._crit_edge, label %.lver.check253

.lver.check253:                                   ; preds = %70
  %92 = getelementptr inbounds i8, ptr %.pre, i64 4
  %93 = load i32, ptr %92, align 4, !tbaa !19
  %94 = sext i32 %93 to i64
  %invariant.gep = getelementptr double, ptr %53, i64 %indvars.iv215
  %invariant.gep242 = getelementptr double, ptr %53, i64 %indvars.iv215
  %ident.check252.not = icmp eq i32 %93, 1
  br i1 %ident.check252.not, label %.ph254, label %.ph254.lver.orig.preheader

.ph254.lver.orig.preheader:                       ; preds = %.lver.check253
  br i1 %65, label %._crit_edge.loopexit261.unr-lcssa, label %.ph254.lver.orig

.ph254.lver.orig:                                 ; preds = %.ph254.lver.orig.preheader, %.ph254.lver.orig
  %indvars.iv.lver.orig = phi i64 [ %indvars.iv.next.lver.orig.1, %.ph254.lver.orig ], [ 1, %.ph254.lver.orig.preheader ]
  %niter = phi i64 [ %niter.next.1, %.ph254.lver.orig ], [ 0, %.ph254.lver.orig.preheader ]
  %95 = add nsw i64 %indvars.iv.lver.orig, -1
  %96 = mul nsw i64 %95, %94
  %gep.lver.orig = getelementptr double, ptr %invariant.gep, i64 %96
  %97 = load double, ptr %gep.lver.orig, align 8, !tbaa !29
  %98 = fmul double %77, %97
  %99 = mul nsw i64 %indvars.iv.lver.orig, %94
  %gep243.lver.orig = getelementptr double, ptr %invariant.gep242, i64 %99
  store double %98, ptr %gep243.lver.orig, align 8, !tbaa !29
  %indvars.iv.next.lver.orig = add nuw nsw i64 %indvars.iv.lver.orig, 1
  %100 = mul nsw i64 %indvars.iv.lver.orig, %94
  %gep.lver.orig.1 = getelementptr double, ptr %invariant.gep, i64 %100
  %101 = load double, ptr %gep.lver.orig.1, align 8, !tbaa !29
  %102 = fmul double %77, %101
  %103 = mul nsw i64 %indvars.iv.next.lver.orig, %94
  %gep243.lver.orig.1 = getelementptr double, ptr %invariant.gep242, i64 %103
  store double %102, ptr %gep243.lver.orig.1, align 8, !tbaa !29
  %indvars.iv.next.lver.orig.1 = add nuw nsw i64 %indvars.iv.lver.orig, 2
  %niter.next.1 = add nuw i64 %niter, 2
  %niter.ncmp.1 = icmp eq i64 %niter.next.1, %unroll_iter
  br i1 %niter.ncmp.1, label %._crit_edge.loopexit261.unr-lcssa, label %.ph254.lver.orig, !llvm.loop !31

.ph254:                                           ; preds = %.lver.check253
  %load_initial257 = load double, ptr %scevgep256, align 8
  br i1 %66, label %._crit_edge.loopexit.unr-lcssa, label %.ph254.new

.ph254.new:                                       ; preds = %.ph254, %.ph254.new
  %store_forwarded258 = phi double [ %110, %.ph254.new ], [ %load_initial257, %.ph254 ]
  %indvars.iv = phi i64 [ %indvars.iv.next.3, %.ph254.new ], [ 1, %.ph254 ]
  %niter265 = phi i64 [ %niter265.next.3, %.ph254.new ], [ 0, %.ph254 ]
  %104 = fmul double %77, %store_forwarded258
  %105 = mul nuw nsw i64 %indvars.iv, %94
  %gep243 = getelementptr double, ptr %invariant.gep242, i64 %105
  store double %104, ptr %gep243, align 8, !tbaa !29
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %106 = fmul double %77, %104
  %107 = mul nuw nsw i64 %indvars.iv.next, %94
  %gep243.1 = getelementptr double, ptr %invariant.gep242, i64 %107
  store double %106, ptr %gep243.1, align 8, !tbaa !29
  %indvars.iv.next.1 = add nuw nsw i64 %indvars.iv, 2
  %108 = fmul double %77, %106
  %109 = mul nuw nsw i64 %indvars.iv.next.1, %94
  %gep243.2 = getelementptr double, ptr %invariant.gep242, i64 %109
  store double %108, ptr %gep243.2, align 8, !tbaa !29
  %indvars.iv.next.2 = add nuw nsw i64 %indvars.iv, 3
  %110 = fmul double %77, %108
  %111 = mul nuw nsw i64 %indvars.iv.next.2, %94
  %gep243.3 = getelementptr double, ptr %invariant.gep242, i64 %111
  store double %110, ptr %gep243.3, align 8, !tbaa !29
  %indvars.iv.next.3 = add nuw nsw i64 %indvars.iv, 4
  %niter265.next.3 = add i64 %niter265, 4
  %niter265.ncmp.3 = icmp eq i64 %niter265.next.3, %unroll_iter264
  br i1 %niter265.ncmp.3, label %._crit_edge.loopexit.unr-lcssa, label %.ph254.new, !llvm.loop !31

._crit_edge.loopexit.unr-lcssa:                   ; preds = %.ph254.new, %.ph254
  %store_forwarded258.unr = phi double [ %load_initial257, %.ph254 ], [ %110, %.ph254.new ]
  %indvars.iv.unr = phi i64 [ 1, %.ph254 ], [ %indvars.iv.next.3, %.ph254.new ]
  br i1 %lcmp.mod263.not, label %._crit_edge, label %.epil.preheader

.epil.preheader:                                  ; preds = %._crit_edge.loopexit.unr-lcssa, %.epil.preheader
  %store_forwarded258.epil = phi double [ %112, %.epil.preheader ], [ %store_forwarded258.unr, %._crit_edge.loopexit.unr-lcssa ]
  %indvars.iv.epil = phi i64 [ %indvars.iv.next.epil, %.epil.preheader ], [ %indvars.iv.unr, %._crit_edge.loopexit.unr-lcssa ]
  %epil.iter = phi i64 [ %epil.iter.next, %.epil.preheader ], [ 0, %._crit_edge.loopexit.unr-lcssa ]
  %112 = fmul double %77, %store_forwarded258.epil
  %113 = mul nuw nsw i64 %indvars.iv.epil, %94
  %gep243.epil = getelementptr double, ptr %invariant.gep242, i64 %113
  store double %112, ptr %gep243.epil, align 8, !tbaa !29
  %indvars.iv.next.epil = add nuw nsw i64 %indvars.iv.epil, 1
  %epil.iter.next = add i64 %epil.iter, 1
  %epil.iter.cmp.not = icmp eq i64 %epil.iter.next, %xtraiter262
  br i1 %epil.iter.cmp.not, label %._crit_edge, label %.epil.preheader, !llvm.loop !33

._crit_edge.loopexit261.unr-lcssa:                ; preds = %.ph254.lver.orig, %.ph254.lver.orig.preheader
  %indvars.iv.lver.orig.unr = phi i64 [ 1, %.ph254.lver.orig.preheader ], [ %indvars.iv.next.lver.orig.1, %.ph254.lver.orig ]
  br i1 %lcmp.mod.not, label %._crit_edge, label %.ph254.lver.orig.epil

.ph254.lver.orig.epil:                            ; preds = %._crit_edge.loopexit261.unr-lcssa
  %114 = add nsw i64 %indvars.iv.lver.orig.unr, -1
  %115 = mul nsw i64 %114, %94
  %gep.lver.orig.epil = getelementptr double, ptr %invariant.gep, i64 %115
  %116 = load double, ptr %gep.lver.orig.epil, align 8, !tbaa !29
  %117 = fmul double %77, %116
  %118 = mul nsw i64 %indvars.iv.lver.orig.unr, %94
  %gep243.lver.orig.epil = getelementptr double, ptr %invariant.gep242, i64 %118
  store double %117, ptr %gep243.lver.orig.epil, align 8, !tbaa !29
  br label %._crit_edge

._crit_edge:                                      ; preds = %.ph254.lver.orig.epil, %._crit_edge.loopexit261.unr-lcssa, %._crit_edge.loopexit.unr-lcssa, %.epil.preheader, %70
  %indvars.iv.next216 = add nsw i64 %indvars.iv215, 1
  %119 = getelementptr inbounds i8, ptr %.pre, i64 44
  %120 = load i32, ptr %119, align 4, !tbaa !27
  %121 = sext i32 %120 to i64
  %122 = icmp slt i64 %indvars.iv.next216, %121
  %indvar.next = add i64 %indvar, 1
  br i1 %122, label %70, label %.preheader, !llvm.loop !35

.lr.ph194:                                        ; preds = %._crit_edge195, %.lr.ph199
  %indvars.iv218 = phi i64 [ 0, %.lr.ph199 ], [ %indvars.iv.next219, %._crit_edge195 ]
  %.0167198 = phi i32 [ 0, %.lr.ph199 ], [ %.2, %._crit_edge195 ]
  %123 = load ptr, ptr %54, align 8, !tbaa !25
  %124 = getelementptr inbounds i8, ptr %123, i64 4
  %125 = load i32, ptr %124, align 4, !tbaa !19
  %126 = sext i32 %125 to i64
  %127 = mul nsw i64 %indvars.iv218, %126
  %128 = getelementptr inbounds double, ptr %53, i64 %127
  store ptr %128, ptr %67, align 8, !tbaa !36
  %129 = mul nsw i64 %indvars.iv218, %69
  %130 = trunc nuw nsw i64 %indvars.iv218 to i32
  %invariant.gep244 = getelementptr double, ptr %26, i64 %129
  br label %131

131:                                              ; preds = %.lr.ph194, %154
  %indvars.iv220 = phi i64 [ %indvars.iv218, %.lr.ph194 ], [ %indvars.iv.next221, %154 ]
  %.1192 = phi i32 [ %.0167198, %.lr.ph194 ], [ %.2, %154 ]
  %132 = load ptr, ptr %54, align 8, !tbaa !25
  %133 = getelementptr inbounds i8, ptr %132, i64 4
  %134 = load i32, ptr %133, align 4, !tbaa !19
  %135 = sext i32 %134 to i64
  %136 = mul nsw i64 %indvars.iv220, %135
  %137 = getelementptr inbounds double, ptr %52, i64 %136
  store ptr %137, ptr %68, align 8, !tbaa !37
  %138 = call i32 @integrateGC_PS93(ptr noundef nonnull %15, ptr noundef %132) #8
  %.not179 = icmp eq i32 %138, 0
  br i1 %.not179, label %148, label %139

139:                                              ; preds = %131
  %140 = shl nsw i32 %.1192, 1
  %141 = sext i32 %140 to i64
  %142 = getelementptr inbounds i32, ptr %34, i64 %141
  store i32 %130, ptr %142, align 4, !tbaa !38
  %143 = or disjoint i32 %140, 1
  %144 = sext i32 %143 to i64
  %145 = getelementptr inbounds i32, ptr %34, i64 %144
  %146 = trunc nsw i64 %indvars.iv220 to i32
  store i32 %146, ptr %145, align 4, !tbaa !38
  %147 = add nsw i32 %.1192, 1
  br label %154

148:                                              ; preds = %131
  %149 = load ptr, ptr %54, align 8, !tbaa !25
  %150 = getelementptr inbounds i8, ptr %149, i64 32
  %151 = load double, ptr %150, align 8, !tbaa !39
  %gep245 = getelementptr double, ptr %invariant.gep244, i64 %indvars.iv220
  %152 = load double, ptr %gep245, align 8, !tbaa !29
  %153 = fadd double %151, %152
  store double %153, ptr %gep245, align 8, !tbaa !29
  br label %154

154:                                              ; preds = %139, %148
  %.2 = phi i32 [ %147, %139 ], [ %.1192, %148 ]
  %indvars.iv.next221 = add nsw i64 %indvars.iv220, -2
  %155 = icmp sgt i64 %indvars.iv220, 1
  br i1 %155, label %131, label %._crit_edge195, !llvm.loop !40

._crit_edge195:                                   ; preds = %154
  %indvars.iv.next219 = add nuw nsw i64 %indvars.iv218, 1
  %exitcond225.not = icmp eq i64 %indvars.iv.next219, %wide.trip.count224
  br i1 %exitcond225.not, label %._crit_edge200, label %.lr.ph194, !llvm.loop !41

._crit_edge200:                                   ; preds = %._crit_edge195
  %156 = icmp sgt i32 %.2, 0
  br i1 %156, label %157, label %.loopexit

157:                                              ; preds = %._crit_edge200
  store double %35, ptr %14, align 8, !tbaa !13
  %158 = load ptr, ptr %38, align 8, !tbaa !16
  %159 = call ptr @GCIntegrationTable_copy(ptr noundef %158) #8
  %160 = fadd double %7, %11
  %161 = call double @llvm.fmuladd.f64(double %7, double %5, double %19)
  %162 = fdiv double %161, %160
  %163 = getelementptr inbounds i8, ptr %159, i64 4
  %164 = load i32, ptr %163, align 4, !tbaa !19
  %165 = getelementptr inbounds i8, ptr %159, i64 8
  %166 = load ptr, ptr %165, align 8, !tbaa !28
  %167 = getelementptr inbounds i8, ptr %159, i64 16
  %168 = load ptr, ptr %167, align 8, !tbaa !42
  call void @transformGrid_FM06(i32 noundef %164, ptr noundef %166, ptr noundef %168, double noundef %160, double noundef %162) #8
  %169 = load i32, ptr %163, align 4, !tbaa !19
  %170 = icmp sgt i32 %169, 0
  br i1 %170, label %.lr.ph209, label %.lr.ph212

.lr.ph209:                                        ; preds = %157
  %171 = getelementptr inbounds i8, ptr %0, i64 88
  %172 = getelementptr inbounds i8, ptr %2, i64 12
  %.not177202 = icmp eq i32 %4, 0
  %173 = add nsw i64 %wide.trip.count224, -1
  %174 = add nsw i64 %wide.trip.count224, -2
  %xtraiter266 = and i64 %173, 1
  %175 = icmp eq i64 %174, 0
  %unroll_iter269 = and i64 %173, -2
  %lcmp.mod268.not = icmp eq i64 %xtraiter266, 0
  %xtraiter272 = and i64 %173, 3
  %176 = icmp ult i64 %174, 3
  %unroll_iter275 = and i64 %173, -4
  %lcmp.mod274.not = icmp eq i64 %xtraiter272, 0
  br label %177

177:                                              ; preds = %.lr.ph209, %.loopexit183
  %178 = phi i32 [ %169, %.lr.ph209 ], [ %231, %.loopexit183 ]
  %indvars.iv231 = phi i64 [ 0, %.lr.ph209 ], [ %indvars.iv.next232, %.loopexit183 ]
  %179 = shl nuw nsw i64 %indvars.iv231, 3
  %scevgep = getelementptr i8, ptr %53, i64 %179
  %180 = load ptr, ptr %165, align 8, !tbaa !28
  %181 = getelementptr inbounds double, ptr %180, i64 %indvars.iv231
  %182 = load double, ptr %181, align 8, !tbaa !29
  %183 = call double @llvm.fmuladd.f64(double %30, double %182, double %3)
  %184 = call double @llvm.fmuladd.f64(double %183, double %182, double %22)
  %185 = load ptr, ptr %44, align 8, !tbaa !21
  %186 = getelementptr inbounds double, ptr %185, i64 %indvars.iv231
  store double %184, ptr %186, align 8, !tbaa !29
  %187 = fcmp ogt double %182, %5
  %188 = fcmp ogt double %182, %9
  %or.cond = and i1 %187, %188
  %189 = load double, ptr %47, align 8
  %190 = fcmp olt double %184, %189
  %or.cond182 = select i1 %or.cond, i1 %190, i1 false
  br i1 %or.cond182, label %191, label %195

191:                                              ; preds = %177
  %192 = trunc nuw nsw i64 %indvars.iv231 to i32
  %193 = add nsw i32 %192, -1
  %194 = getelementptr inbounds i8, ptr %159, i64 44
  store i32 %193, ptr %194, align 4, !tbaa !27
  br label %.lr.ph212

195:                                              ; preds = %177
  %196 = fcmp ult double %184, %189
  br i1 %196, label %.loopexit183, label %197

197:                                              ; preds = %195
  %198 = load ptr, ptr %171, align 8, !tbaa !30
  %199 = fmul double %182, %3
  %200 = getelementptr inbounds double, ptr %52, i64 %indvars.iv231
  call void @weightedBesselFunction(ptr noundef %198, i32 noundef %178, i32 noundef %4, double noundef %199, ptr noundef %200) #8
  %201 = load i32, ptr %172, align 4, !tbaa !43
  %202 = call double @evalECP(ptr noundef %2, i32 noundef %201, double noundef %182) #8
  %203 = load ptr, ptr %49, align 8, !tbaa !24
  %204 = getelementptr inbounds double, ptr %203, i64 %indvars.iv231
  store double %202, ptr %204, align 8, !tbaa !29
  %205 = getelementptr inbounds double, ptr %53, i64 %indvars.iv231
  store double 1.000000e+00, ptr %205, align 8, !tbaa !29
  %.pre239 = load i32, ptr %163, align 4, !tbaa !19
  br i1 %.not177202, label %.loopexit183, label %.lver.check

.lver.check:                                      ; preds = %197
  %206 = sext i32 %.pre239 to i64
  %invariant.gep246 = getelementptr double, ptr %53, i64 %indvars.iv231
  %invariant.gep248 = getelementptr double, ptr %53, i64 %indvars.iv231
  %ident.check.not = icmp eq i32 %.pre239, 1
  br i1 %ident.check.not, label %.ph, label %.ph.lver.orig.preheader

.ph.lver.orig.preheader:                          ; preds = %.lver.check
  br i1 %175, label %.loopexit183.loopexit259.unr-lcssa, label %.ph.lver.orig

.ph.lver.orig:                                    ; preds = %.ph.lver.orig.preheader, %.ph.lver.orig
  %indvars.iv226.lver.orig = phi i64 [ %indvars.iv.next227.lver.orig.1, %.ph.lver.orig ], [ 1, %.ph.lver.orig.preheader ]
  %niter270 = phi i64 [ %niter270.next.1, %.ph.lver.orig ], [ 0, %.ph.lver.orig.preheader ]
  %207 = add nsw i64 %indvars.iv226.lver.orig, -1
  %208 = mul nsw i64 %207, %206
  %gep247.lver.orig = getelementptr double, ptr %invariant.gep246, i64 %208
  %209 = load double, ptr %gep247.lver.orig, align 8, !tbaa !29
  %210 = fmul double %182, %209
  %211 = mul nsw i64 %indvars.iv226.lver.orig, %206
  %gep249.lver.orig = getelementptr double, ptr %invariant.gep248, i64 %211
  store double %210, ptr %gep249.lver.orig, align 8, !tbaa !29
  %indvars.iv.next227.lver.orig = add nuw nsw i64 %indvars.iv226.lver.orig, 1
  %212 = mul nsw i64 %indvars.iv226.lver.orig, %206
  %gep247.lver.orig.1 = getelementptr double, ptr %invariant.gep246, i64 %212
  %213 = load double, ptr %gep247.lver.orig.1, align 8, !tbaa !29
  %214 = fmul double %182, %213
  %215 = mul nsw i64 %indvars.iv.next227.lver.orig, %206
  %gep249.lver.orig.1 = getelementptr double, ptr %invariant.gep248, i64 %215
  store double %214, ptr %gep249.lver.orig.1, align 8, !tbaa !29
  %indvars.iv.next227.lver.orig.1 = add nuw nsw i64 %indvars.iv226.lver.orig, 2
  %niter270.next.1 = add nuw i64 %niter270, 2
  %niter270.ncmp.1 = icmp eq i64 %niter270.next.1, %unroll_iter269
  br i1 %niter270.ncmp.1, label %.loopexit183.loopexit259.unr-lcssa, label %.ph.lver.orig, !llvm.loop !45

.ph:                                              ; preds = %.lver.check
  %load_initial = load double, ptr %scevgep, align 8
  br i1 %176, label %.loopexit183.loopexit.unr-lcssa, label %.ph.new

.ph.new:                                          ; preds = %.ph, %.ph.new
  %store_forwarded = phi double [ %222, %.ph.new ], [ %load_initial, %.ph ]
  %indvars.iv226 = phi i64 [ %indvars.iv.next227.3, %.ph.new ], [ 1, %.ph ]
  %niter276 = phi i64 [ %niter276.next.3, %.ph.new ], [ 0, %.ph ]
  %216 = fmul double %182, %store_forwarded
  %217 = mul nuw nsw i64 %indvars.iv226, %206
  %gep249 = getelementptr double, ptr %invariant.gep248, i64 %217
  store double %216, ptr %gep249, align 8, !tbaa !29
  %indvars.iv.next227 = add nuw nsw i64 %indvars.iv226, 1
  %218 = fmul double %182, %216
  %219 = mul nuw nsw i64 %indvars.iv.next227, %206
  %gep249.1 = getelementptr double, ptr %invariant.gep248, i64 %219
  store double %218, ptr %gep249.1, align 8, !tbaa !29
  %indvars.iv.next227.1 = add nuw nsw i64 %indvars.iv226, 2
  %220 = fmul double %182, %218
  %221 = mul nuw nsw i64 %indvars.iv.next227.1, %206
  %gep249.2 = getelementptr double, ptr %invariant.gep248, i64 %221
  store double %220, ptr %gep249.2, align 8, !tbaa !29
  %indvars.iv.next227.2 = add nuw nsw i64 %indvars.iv226, 3
  %222 = fmul double %182, %220
  %223 = mul nuw nsw i64 %indvars.iv.next227.2, %206
  %gep249.3 = getelementptr double, ptr %invariant.gep248, i64 %223
  store double %222, ptr %gep249.3, align 8, !tbaa !29
  %indvars.iv.next227.3 = add nuw nsw i64 %indvars.iv226, 4
  %niter276.next.3 = add i64 %niter276, 4
  %niter276.ncmp.3 = icmp eq i64 %niter276.next.3, %unroll_iter275
  br i1 %niter276.ncmp.3, label %.loopexit183.loopexit.unr-lcssa, label %.ph.new, !llvm.loop !45

.loopexit183.loopexit.unr-lcssa:                  ; preds = %.ph.new, %.ph
  %store_forwarded.unr = phi double [ %load_initial, %.ph ], [ %222, %.ph.new ]
  %indvars.iv226.unr = phi i64 [ 1, %.ph ], [ %indvars.iv.next227.3, %.ph.new ]
  br i1 %lcmp.mod274.not, label %.loopexit183, label %.epil.preheader271

.epil.preheader271:                               ; preds = %.loopexit183.loopexit.unr-lcssa, %.epil.preheader271
  %store_forwarded.epil = phi double [ %224, %.epil.preheader271 ], [ %store_forwarded.unr, %.loopexit183.loopexit.unr-lcssa ]
  %indvars.iv226.epil = phi i64 [ %indvars.iv.next227.epil, %.epil.preheader271 ], [ %indvars.iv226.unr, %.loopexit183.loopexit.unr-lcssa ]
  %epil.iter273 = phi i64 [ %epil.iter273.next, %.epil.preheader271 ], [ 0, %.loopexit183.loopexit.unr-lcssa ]
  %224 = fmul double %182, %store_forwarded.epil
  %225 = mul nuw nsw i64 %indvars.iv226.epil, %206
  %gep249.epil = getelementptr double, ptr %invariant.gep248, i64 %225
  store double %224, ptr %gep249.epil, align 8, !tbaa !29
  %indvars.iv.next227.epil = add nuw nsw i64 %indvars.iv226.epil, 1
  %epil.iter273.next = add i64 %epil.iter273, 1
  %epil.iter273.cmp.not = icmp eq i64 %epil.iter273.next, %xtraiter272
  br i1 %epil.iter273.cmp.not, label %.loopexit183, label %.epil.preheader271, !llvm.loop !46

.loopexit183.loopexit259.unr-lcssa:               ; preds = %.ph.lver.orig, %.ph.lver.orig.preheader
  %indvars.iv226.lver.orig.unr = phi i64 [ 1, %.ph.lver.orig.preheader ], [ %indvars.iv.next227.lver.orig.1, %.ph.lver.orig ]
  br i1 %lcmp.mod268.not, label %.loopexit183, label %.ph.lver.orig.epil

.ph.lver.orig.epil:                               ; preds = %.loopexit183.loopexit259.unr-lcssa
  %226 = add nsw i64 %indvars.iv226.lver.orig.unr, -1
  %227 = mul nsw i64 %226, %206
  %gep247.lver.orig.epil = getelementptr double, ptr %invariant.gep246, i64 %227
  %228 = load double, ptr %gep247.lver.orig.epil, align 8, !tbaa !29
  %229 = fmul double %182, %228
  %230 = mul nsw i64 %indvars.iv226.lver.orig.unr, %206
  %gep249.lver.orig.epil = getelementptr double, ptr %invariant.gep248, i64 %230
  store double %229, ptr %gep249.lver.orig.epil, align 8, !tbaa !29
  br label %.loopexit183

.loopexit183:                                     ; preds = %.ph.lver.orig.epil, %.loopexit183.loopexit259.unr-lcssa, %.loopexit183.loopexit.unr-lcssa, %.epil.preheader271, %197, %195
  %231 = phi i32 [ %.pre239, %197 ], [ %178, %195 ], [ %.pre239, %.epil.preheader271 ], [ %.pre239, %.loopexit183.loopexit.unr-lcssa ], [ %.pre239, %.loopexit183.loopexit259.unr-lcssa ], [ %.pre239, %.ph.lver.orig.epil ]
  %indvars.iv.next232 = add nuw nsw i64 %indvars.iv231, 1
  %232 = sext i32 %231 to i64
  %233 = icmp slt i64 %indvars.iv.next232, %232
  br i1 %233, label %177, label %.lr.ph212, !llvm.loop !47

.lr.ph212:                                        ; preds = %.loopexit183, %191, %157
  %234 = getelementptr inbounds i8, ptr %14, i64 32
  %235 = getelementptr inbounds i8, ptr %14, i64 16
  %236 = getelementptr inbounds i8, ptr %159, i64 32
  %smax = call i32 @llvm.smax.i32(i32 %.2, i32 1)
  %wide.trip.count237 = zext nneg i32 %smax to i64
  br label %237

237:                                              ; preds = %.lr.ph212, %253
  %indvars.iv234 = phi i64 [ 0, %.lr.ph212 ], [ %indvars.iv.next235, %253 ]
  %238 = shl nuw nsw i64 %indvars.iv234, 1
  %239 = getelementptr inbounds i32, ptr %34, i64 %238
  %240 = load i32, ptr %239, align 4, !tbaa !38
  %241 = or disjoint i64 %238, 1
  %242 = getelementptr inbounds i32, ptr %34, i64 %241
  %243 = load i32, ptr %242, align 4, !tbaa !38
  %244 = load i32, ptr %163, align 4, !tbaa !19
  %245 = mul nsw i32 %244, %240
  %246 = sext i32 %245 to i64
  %247 = getelementptr inbounds double, ptr %53, i64 %246
  store ptr %247, ptr %234, align 8, !tbaa !36
  %248 = mul nsw i32 %244, %243
  %249 = sext i32 %248 to i64
  %250 = getelementptr inbounds double, ptr %52, i64 %249
  store ptr %250, ptr %235, align 8, !tbaa !37
  %251 = call i32 @integrateGC_PSM92(ptr noundef nonnull %15, ptr noundef nonnull %159) #8
  %.not178 = icmp eq i32 %251, 0
  br i1 %.not178, label %253, label %252

252:                                              ; preds = %237
  call void @doubleArray_free(ptr %25, ptr %26) #8
  br label %.loopexit

253:                                              ; preds = %237
  %254 = load double, ptr %236, align 8, !tbaa !39
  %255 = mul nsw i32 %240, %29
  %256 = add nsw i32 %255, %243
  %257 = sext i32 %256 to i64
  %258 = getelementptr inbounds double, ptr %26, i64 %257
  %259 = load double, ptr %258, align 8, !tbaa !29
  %260 = fadd double %254, %259
  store double %260, ptr %258, align 8, !tbaa !29
  %indvars.iv.next235 = add nuw nsw i64 %indvars.iv234, 1
  %exitcond238.not = icmp eq i64 %indvars.iv.next235, %wide.trip.count237
  br i1 %exitcond238.not, label %.loopexit, label %237, !llvm.loop !48

.loopexit:                                        ; preds = %253, %.preheader, %._crit_edge200, %252
  %261 = phi i1 [ %156, %252 ], [ false, %._crit_edge200 ], [ false, %.preheader ], [ %156, %253 ]
  %.0 = phi ptr [ %159, %252 ], [ undef, %._crit_edge200 ], [ undef, %.preheader ], [ %159, %253 ]
  %262 = load ptr, ptr %44, align 8, !tbaa !21
  call void @free(ptr noundef %262) #8
  %263 = load ptr, ptr %49, align 8, !tbaa !24
  call void @free(ptr noundef %263) #8
  call void @free(ptr noundef %52) #8
  call void @free(ptr noundef %53) #8
  call void @free(ptr noundef %34) #8
  br i1 %261, label %264, label %265

264:                                              ; preds = %.loopexit
  call void @GCIntegrationTable_free(ptr noundef %.0) #8
  br label %265

265:                                              ; preds = %264, %.loopexit
  call void @llvm.lifetime.end.p0(i64 16, ptr nonnull %15) #8
  call void @llvm.lifetime.end.p0(i64 48, ptr nonnull %14) #8
  ret { ptr, ptr } %24
}

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

declare { ptr, ptr } @doubleArray_new(i64 noundef, ...) local_unnamed_addr #1

; Function Attrs: mustprogress nofree nounwind willreturn allockind("alloc,zeroed") allocsize(0,1) memory(inaccessiblemem: readwrite)
declare noalias noundef ptr @calloc(i64 noundef, i64 noundef) local_unnamed_addr #6

declare void @weightedBesselFunction(ptr noundef, i32 noundef, i32 noundef, double noundef, ptr noundef) local_unnamed_addr #1

declare i32 @integrateGC_PS93(ptr noundef, ptr noundef) local_unnamed_addr #1

declare ptr @GCIntegrationTable_copy(ptr noundef) local_unnamed_addr #1

declare void @transformGrid_FM06(i32 noundef, ptr noundef, ptr noundef, double noundef, double noundef) local_unnamed_addr #1

declare double @evalECP(ptr noundef, i32 noundef, double noundef) local_unnamed_addr #1

declare i32 @integrateGC_PSM92(ptr noundef, ptr noundef) local_unnamed_addr #1

declare void @doubleArray_free(ptr, ptr) local_unnamed_addr #1

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare i32 @llvm.smax.i32(i32, i32) #7

attributes #0 = { mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { "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 = { mustprogress nounwind willreturn allockind("free") memory(argmem: readwrite, inaccessiblemem: readwrite) "alloc-family"="malloc" "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 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 = { mustprogress nofree nounwind willreturn memory(write) "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 #5 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #6 = { mustprogress nofree nounwind willreturn allockind("alloc,zeroed") allocsize(0,1) memory(inaccessiblemem: readwrite) "alloc-family"="malloc" "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 #7 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #8 = { nounwind }
attributes #9 = { nounwind allocsize(0,1) }

!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, !7, i64 0}
!6 = !{!"", !7, i64 0, !7, i64 8}
!7 = !{!"any pointer", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C/C++ TBAA"}
!10 = !{!6, !7, i64 8}
!11 = !{!12, !12, i64 0}
!12 = !{!"long", !8, i64 0}
!13 = !{!14, !15, i64 0}
!14 = !{!"", !15, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !15, i64 40}
!15 = !{!"double", !8, i64 0}
!16 = !{!17, !7, i64 104}
!17 = !{!"", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !18, i64 40, !7, i64 48, !7, i64 56, !18, i64 64, !7, i64 72, !18, i64 80, !7, i64 88, !18, i64 96, !18, i64 100, !7, i64 104, !18, i64 112, !7, i64 120, !15, i64 128, !15, i64 136, !15, i64 144}
!18 = !{!"int", !8, i64 0}
!19 = !{!20, !18, i64 4}
!20 = !{!"", !18, i64 0, !18, i64 4, !7, i64 8, !7, i64 16, !15, i64 24, !15, i64 32, !18, i64 40, !18, i64 44}
!21 = !{!14, !7, i64 8}
!22 = !{!17, !15, i64 136}
!23 = !{!14, !15, i64 40}
!24 = !{!14, !7, i64 24}
!25 = !{!17, !7, i64 120}
!26 = !{!20, !18, i64 40}
!27 = !{!20, !18, i64 44}
!28 = !{!20, !7, i64 8}
!29 = !{!15, !15, i64 0}
!30 = !{!17, !7, i64 88}
!31 = distinct !{!31, !32}
!32 = !{!"llvm.loop.mustprogress"}
!33 = distinct !{!33, !34}
!34 = !{!"llvm.loop.unroll.disable"}
!35 = distinct !{!35, !32}
!36 = !{!14, !7, i64 32}
!37 = !{!14, !7, i64 16}
!38 = !{!18, !18, i64 0}
!39 = !{!20, !15, i64 32}
!40 = distinct !{!40, !32}
!41 = distinct !{!41, !32}
!42 = !{!20, !7, i64 16}
!43 = !{!44, !18, i64 12}
!44 = !{!"_ECP", !7, i64 0, !18, i64 8, !18, i64 12}
!45 = distinct !{!45, !32}
!46 = distinct !{!46, !34}
!47 = distinct !{!47, !32}
!48 = distinct !{!48, !32}


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