; LLVM IR from CLASS-SZ/class_sz :: class-sz/tools (hyperspherical)

; 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.HypersphericalInterpolationStructure = type { i32, double, double, i32, i32, ptr, ptr, i32, ptr, ptr, ptr, ptr, ptr }

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

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

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

; Function Attrs: nounwind uwtable
declare dso_local i32 @ClosedModY(i32 noundef, i32 noundef, ptr noundef, ptr noundef, ptr noundef) #2

; Function Attrs: nocallback nofree nounwind willreturn memory(argmem: write)
declare void @llvm.memset.p0.i64(ptr nocapture writeonly, i8, i64, i1 immarg) #3

; Function Attrs: nounwind uwtable
define dso_local i32 @hyperspherical_Hermite3_interpolation_vector_PhidPhid2Phi(ptr noundef %0, i32 noundef %1, i32 noundef %2, ptr noundef %3, ptr noundef %4, ptr noundef %5, ptr noundef %6, ptr noundef %7) #2 {
  %9 = alloca ptr, align 8
  %10 = alloca i32, align 4
  %11 = alloca i32, align 4
  %12 = alloca ptr, align 8
  %13 = alloca ptr, align 8
  %14 = alloca ptr, align 8
  %15 = alloca ptr, align 8
  %16 = alloca ptr, align 8
  %17 = alloca i32, align 4
  %18 = alloca double, align 8
  %19 = alloca double, align 8
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca [2 x double], align 16
  %23 = alloca [3 x double], align 16
  %24 = alloca [3 x double], align 16
  %25 = alloca [3 x double], align 16
  %26 = alloca double, align 8
  %27 = alloca double, align 8
  %28 = alloca double, align 8
  %29 = alloca double, align 8
  %30 = alloca double, align 8
  %31 = alloca double, align 8
  %32 = alloca double, align 8
  %33 = alloca double, align 8
  %34 = alloca double, align 8
  %35 = alloca double, align 8
  %36 = alloca ptr, align 8
  %37 = alloca ptr, align 8
  %38 = alloca i32, align 4
  %39 = alloca double, align 8
  %40 = alloca double, align 8
  %41 = alloca double, align 8
  %42 = alloca ptr, align 8
  %43 = alloca double, align 8
  %44 = alloca double, align 8
  %45 = alloca double, align 8
  %46 = alloca double, align 8
  %47 = alloca double, align 8
  %48 = alloca double, align 8
  %49 = alloca i32, align 4
  %50 = alloca i32, align 4
  %51 = alloca i32, align 4
  %52 = alloca ptr, align 8
  %53 = alloca ptr, align 8
  %54 = alloca i32, align 4
  %55 = alloca i32, align 4
  store ptr %0, ptr %9, align 8, !tbaa !5
  store i32 %1, ptr %10, align 4, !tbaa !9
  store i32 %2, ptr %11, align 4, !tbaa !9
  store ptr %3, ptr %12, align 8, !tbaa !5
  store ptr %4, ptr %13, align 8, !tbaa !5
  store ptr %5, ptr %14, align 8, !tbaa !5
  store ptr %6, ptr %15, align 8, !tbaa !5
  store ptr %7, ptr %16, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %17) #4
  %56 = load ptr, ptr %9, align 8, !tbaa !5
  %57 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %56, i32 0, i32 5
  %58 = load ptr, ptr %57, align 8, !tbaa !11
  %59 = load i32, ptr %11, align 4, !tbaa !9
  %60 = sext i32 %59 to i64
  %61 = getelementptr inbounds i32, ptr %58, i64 %60
  %62 = load i32, ptr %61, align 4, !tbaa !9
  store i32 %62, ptr %17, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #4
  store double 0.000000e+00, ptr %18, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %19) #4
  store double 0.000000e+00, ptr %19, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #4
  store double 0.000000e+00, ptr %20, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #4
  call void @llvm.lifetime.start.p0(i64 16, ptr %22) #4
  call void @llvm.memset.p0.i64(ptr align 16 %22, i8 0, i64 16, i1 false)
  call void @llvm.lifetime.start.p0(i64 24, ptr %23) #4
  call void @llvm.memset.p0.i64(ptr align 16 %23, i8 0, i64 24, i1 false)
  call void @llvm.lifetime.start.p0(i64 24, ptr %24) #4
  call void @llvm.memset.p0.i64(ptr align 16 %24, i8 0, i64 24, i1 false)
  call void @llvm.lifetime.start.p0(i64 24, ptr %25) #4
  call void @llvm.memset.p0.i64(ptr align 16 %25, i8 0, i64 24, i1 false)
  call void @llvm.lifetime.start.p0(i64 8, ptr %26) #4
  store double 0.000000e+00, ptr %26, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %27) #4
  store double 0.000000e+00, ptr %27, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %28) #4
  store double 0.000000e+00, ptr %28, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %29) #4
  store double 0.000000e+00, ptr %29, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %30) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %31) #4
  store double 0.000000e+00, ptr %31, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %32) #4
  store double 0.000000e+00, ptr %32, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %33) #4
  store double 0.000000e+00, ptr %33, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %34) #4
  store double 0.000000e+00, ptr %34, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %35) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %36) #4
  %63 = load ptr, ptr %9, align 8, !tbaa !5
  %64 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %63, i32 0, i32 9
  %65 = load ptr, ptr %64, align 8, !tbaa !15
  store ptr %65, ptr %36, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %37) #4
  %66 = load ptr, ptr %9, align 8, !tbaa !5
  %67 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %66, i32 0, i32 10
  %68 = load ptr, ptr %67, align 8, !tbaa !16
  store ptr %68, ptr %37, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %38) #4
  %69 = load ptr, ptr %9, align 8, !tbaa !5
  %70 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %69, i32 0, i32 0
  %71 = load i32, ptr %70, align 8, !tbaa !17
  store i32 %71, ptr %38, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %39) #4
  %72 = load i32, ptr %17, align 4, !tbaa !9
  %73 = sitofp i32 %72 to double
  %74 = load i32, ptr %17, align 4, !tbaa !9
  %75 = sitofp i32 %74 to double
  %76 = fadd double %75, 1.000000e+00
  %77 = fmul double %73, %76
  store double %77, ptr %39, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %40) #4
  %78 = load ptr, ptr %9, align 8, !tbaa !5
  %79 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %78, i32 0, i32 1
  %80 = load double, ptr %79, align 8, !tbaa !18
  store double %80, ptr %40, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %41) #4
  %81 = load double, ptr %40, align 8, !tbaa !14
  %82 = load double, ptr %40, align 8, !tbaa !14
  %83 = fmul double %81, %82
  store double %83, ptr %41, align 8, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 8, ptr %42) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %43) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %44) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %45) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %46) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %47) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %48) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %49) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %50) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %51) #4
  store i32 0, ptr %51, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %52) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %53) #4
  call void @llvm.lifetime.start.p0(i64 4, ptr %54) #4
  store i32 1, ptr %54, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 4, ptr %55) #4
  store i32 1, ptr %55, align 4, !tbaa !9
  %84 = load ptr, ptr %9, align 8, !tbaa !5
  %85 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %84, i32 0, i32 8
  %86 = load ptr, ptr %85, align 8, !tbaa !19
  store ptr %86, ptr %42, align 8, !tbaa !5
  %87 = load ptr, ptr %9, align 8, !tbaa !5
  %88 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %87, i32 0, i32 2
  %89 = load double, ptr %88, align 8, !tbaa !20
  store double %89, ptr %45, align 8, !tbaa !14
  %90 = load ptr, ptr %9, align 8, !tbaa !5
  %91 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %90, i32 0, i32 7
  %92 = load i32, ptr %91, align 8, !tbaa !21
  store i32 %92, ptr %50, align 4, !tbaa !9
  %93 = load ptr, ptr %9, align 8, !tbaa !5
  %94 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %93, i32 0, i32 11
  %95 = load ptr, ptr %94, align 8, !tbaa !22
  %96 = load i32, ptr %11, align 4, !tbaa !9
  %97 = load i32, ptr %50, align 4, !tbaa !9
  %98 = mul nsw i32 %96, %97
  %99 = sext i32 %98 to i64
  %100 = getelementptr inbounds double, ptr %95, i64 %99
  store ptr %100, ptr %52, align 8, !tbaa !5
  %101 = load ptr, ptr %9, align 8, !tbaa !5
  %102 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %101, i32 0, i32 12
  %103 = load ptr, ptr %102, align 8, !tbaa !23
  %104 = load i32, ptr %11, align 4, !tbaa !9
  %105 = load i32, ptr %50, align 4, !tbaa !9
  %106 = mul nsw i32 %104, %105
  %107 = sext i32 %106 to i64
  %108 = getelementptr inbounds double, ptr %103, i64 %107
  store ptr %108, ptr %53, align 8, !tbaa !5
  %109 = load ptr, ptr %42, align 8, !tbaa !5
  %110 = getelementptr inbounds double, ptr %109, i64 0
  %111 = load double, ptr %110, align 8, !tbaa !14
  store double %111, ptr %43, align 8, !tbaa !14
  %112 = load ptr, ptr %42, align 8, !tbaa !5
  %113 = load i32, ptr %50, align 4, !tbaa !9
  %114 = sub nsw i32 %113, 1
  %115 = sext i32 %114 to i64
  %116 = getelementptr inbounds double, ptr %112, i64 %115
  %117 = load double, ptr %116, align 8, !tbaa !14
  store double %117, ptr %44, align 8, !tbaa !14
  %118 = load double, ptr %44, align 8, !tbaa !14
  store double %118, ptr %46, align 8, !tbaa !14
  %119 = load double, ptr %43, align 8, !tbaa !14
  store double %119, ptr %47, align 8, !tbaa !14
  %120 = load double, ptr %43, align 8, !tbaa !14
  store double %120, ptr %48, align 8, !tbaa !14
  store i32 0, ptr %49, align 4, !tbaa !9
  br label %121

121:                                              ; preds = %470, %8
  %122 = load i32, ptr %49, align 4, !tbaa !9
  %123 = load i32, ptr %10, align 4, !tbaa !9
  %124 = icmp slt i32 %122, %123
  br i1 %124, label %125, label %473

125:                                              ; preds = %121
  %126 = load ptr, ptr %12, align 8, !tbaa !5
  %127 = load i32, ptr %49, align 4, !tbaa !9
  %128 = sext i32 %127 to i64
  %129 = getelementptr inbounds double, ptr %126, i64 %128
  %130 = load double, ptr %129, align 8, !tbaa !14
  store double %130, ptr %21, align 8, !tbaa !14
  %131 = load ptr, ptr %9, align 8, !tbaa !5
  %132 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %131, i32 0, i32 0
  %133 = load i32, ptr %132, align 8, !tbaa !17
  %134 = icmp eq i32 %133, 1
  br i1 %134, label %135, label %143

135:                                              ; preds = %125
  %136 = load i32, ptr %17, align 4, !tbaa !9
  %137 = load ptr, ptr %9, align 8, !tbaa !5
  %138 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %137, i32 0, i32 1
  %139 = load double, ptr %138, align 8, !tbaa !18
  %140 = fadd double %139, 2.000000e-01
  %141 = fptosi double %140 to i32
  %142 = call i32 @ClosedModY(i32 noundef %136, i32 noundef %141, ptr noundef %21, ptr noundef %54, ptr noundef %55)
  br label %143

143:                                              ; preds = %135, %125
  %144 = load double, ptr %21, align 8, !tbaa !14
  %145 = load double, ptr %43, align 8, !tbaa !14
  %146 = fcmp olt double %144, %145
  br i1 %146, label %151, label %147

147:                                              ; preds = %143
  %148 = load double, ptr %21, align 8, !tbaa !14
  %149 = load double, ptr %44, align 8, !tbaa !14
  %150 = fcmp ogt double %148, %149
  br i1 %150, label %151, label %164

151:                                              ; preds = %147, %143
  %152 = load ptr, ptr %13, align 8, !tbaa !5
  %153 = load i32, ptr %49, align 4, !tbaa !9
  %154 = sext i32 %153 to i64
  %155 = getelementptr inbounds double, ptr %152, i64 %154
  store double 0.000000e+00, ptr %155, align 8, !tbaa !14
  %156 = load ptr, ptr %14, align 8, !tbaa !5
  %157 = load i32, ptr %49, align 4, !tbaa !9
  %158 = sext i32 %157 to i64
  %159 = getelementptr inbounds double, ptr %156, i64 %158
  store double 0.000000e+00, ptr %159, align 8, !tbaa !14
  %160 = load ptr, ptr %15, align 8, !tbaa !5
  %161 = load i32, ptr %49, align 4, !tbaa !9
  %162 = sext i32 %161 to i64
  %163 = getelementptr inbounds double, ptr %160, i64 %162
  store double 0.000000e+00, ptr %163, align 8, !tbaa !14
  br label %470

164:                                              ; preds = %147
  %165 = load double, ptr %21, align 8, !tbaa !14
  %166 = load double, ptr %47, align 8, !tbaa !14
  %167 = fcmp ogt double %165, %166
  br i1 %167, label %172, label %168

168:                                              ; preds = %164
  %169 = load double, ptr %21, align 8, !tbaa !14
  %170 = load double, ptr %46, align 8, !tbaa !14
  %171 = fcmp olt double %169, %170
  br i1 %171, label %172, label %403

172:                                              ; preds = %168, %164
  %173 = load double, ptr %21, align 8, !tbaa !14
  %174 = load double, ptr %48, align 8, !tbaa !14
  %175 = fcmp ogt double %173, %174
  br i1 %175, label %180, label %176

176:                                              ; preds = %172
  %177 = load double, ptr %21, align 8, !tbaa !14
  %178 = load double, ptr %46, align 8, !tbaa !14
  %179 = fcmp olt double %177, %178
  br i1 %179, label %180, label %247

180:                                              ; preds = %176, %172
  %181 = load double, ptr %21, align 8, !tbaa !14
  %182 = load double, ptr %43, align 8, !tbaa !14
  %183 = fsub double %181, %182
  %184 = load double, ptr %45, align 8, !tbaa !14
  %185 = fdiv double %183, %184
  %186 = fptosi double %185 to i32
  %187 = add nsw i32 %186, 1
  store i32 %187, ptr %51, align 4, !tbaa !9
  %188 = load i32, ptr %51, align 4, !tbaa !9
  %189 = icmp slt i32 1, %188
  br i1 %189, label %190, label %192

190:                                              ; preds = %180
  %191 = load i32, ptr %51, align 4, !tbaa !9
  br label %193

192:                                              ; preds = %180
  br label %193

193:                                              ; preds = %192, %190
  %194 = phi i32 [ %191, %190 ], [ 1, %192 ]
  store i32 %194, ptr %51, align 4, !tbaa !9
  %195 = load i32, ptr %50, align 4, !tbaa !9
  %196 = sub nsw i32 %195, 1
  %197 = load i32, ptr %51, align 4, !tbaa !9
  %198 = icmp slt i32 %196, %197
  br i1 %198, label %199, label %202

199:                                              ; preds = %193
  %200 = load i32, ptr %50, align 4, !tbaa !9
  %201 = sub nsw i32 %200, 1
  br label %204

202:                                              ; preds = %193
  %203 = load i32, ptr %51, align 4, !tbaa !9
  br label %204

204:                                              ; preds = %202, %199
  %205 = phi i32 [ %201, %199 ], [ %203, %202 ]
  store i32 %205, ptr %51, align 4, !tbaa !9
  %206 = load ptr, ptr %52, align 8, !tbaa !5
  %207 = load i32, ptr %51, align 4, !tbaa !9
  %208 = sub nsw i32 %207, 1
  %209 = sext i32 %208 to i64
  %210 = getelementptr inbounds double, ptr %206, i64 %209
  %211 = load double, ptr %210, align 8, !tbaa !14
  store double %211, ptr %18, align 8, !tbaa !14
  %212 = load ptr, ptr %53, align 8, !tbaa !5
  %213 = load i32, ptr %51, align 4, !tbaa !9
  %214 = sub nsw i32 %213, 1
  %215 = sext i32 %214 to i64
  %216 = getelementptr inbounds double, ptr %212, i64 %215
  %217 = load double, ptr %216, align 8, !tbaa !14
  store double %217, ptr %31, align 8, !tbaa !14
  %218 = load ptr, ptr %37, align 8, !tbaa !5
  %219 = load i32, ptr %51, align 4, !tbaa !9
  %220 = sub nsw i32 %219, 1
  %221 = sext i32 %220 to i64
  %222 = getelementptr inbounds double, ptr %218, i64 %221
  %223 = load double, ptr %222, align 8, !tbaa !14
  store double %223, ptr %28, align 8, !tbaa !14
  %224 = load ptr, ptr %36, align 8, !tbaa !5
  %225 = load i32, ptr %51, align 4, !tbaa !9
  %226 = sub nsw i32 %225, 1
  %227 = sext i32 %226 to i64
  %228 = getelementptr inbounds double, ptr %224, i64 %227
  %229 = load double, ptr %228, align 8, !tbaa !14
  store double %229, ptr %29, align 8, !tbaa !14
  %230 = load double, ptr %29, align 8, !tbaa !14
  %231 = load double, ptr %29, align 8, !tbaa !14
  %232 = fmul double %230, %231
  store double %232, ptr %30, align 8, !tbaa !14
  %233 = load double, ptr %31, align 8, !tbaa !14
  %234 = fmul double -2.000000e+00, %233
  %235 = load double, ptr %28, align 8, !tbaa !14
  %236 = load double, ptr %18, align 8, !tbaa !14
  %237 = load double, ptr %39, align 8, !tbaa !14
  %238 = load double, ptr %30, align 8, !tbaa !14
  %239 = fdiv double %237, %238
  %240 = load double, ptr %41, align 8, !tbaa !14
  %241 = fsub double %239, %240
  %242 = load i32, ptr %38, align 4, !tbaa !9
  %243 = sitofp i32 %242 to double
  %244 = fadd double %241, %243
  %245 = fmul double %236, %244
  %246 = call double @llvm.fmuladd.f64(double %234, double %235, double %245)
  store double %246, ptr %26, align 8, !tbaa !14
  br label %255

247:                                              ; preds = %176
  %248 = load i32, ptr %51, align 4, !tbaa !9
  %249 = add nsw i32 %248, 1
  store i32 %249, ptr %51, align 4, !tbaa !9
  %250 = load double, ptr %19, align 8, !tbaa !14
  store double %250, ptr %18, align 8, !tbaa !14
  %251 = load double, ptr %20, align 8, !tbaa !14
  store double %251, ptr %31, align 8, !tbaa !14
  %252 = load double, ptr %32, align 8, !tbaa !14
  store double %252, ptr %26, align 8, !tbaa !14
  %253 = load double, ptr %34, align 8, !tbaa !14
  store double %253, ptr %29, align 8, !tbaa !14
  %254 = load double, ptr %33, align 8, !tbaa !14
  store double %254, ptr %28, align 8, !tbaa !14
  br label %255

255:                                              ; preds = %247, %204
  %256 = load ptr, ptr %42, align 8, !tbaa !5
  %257 = load i32, ptr %51, align 4, !tbaa !9
  %258 = sub nsw i32 %257, 1
  %259 = icmp slt i32 0, %258
  br i1 %259, label %260, label %263

260:                                              ; preds = %255
  %261 = load i32, ptr %51, align 4, !tbaa !9
  %262 = sub nsw i32 %261, 1
  br label %264

263:                                              ; preds = %255
  br label %264

264:                                              ; preds = %263, %260
  %265 = phi i32 [ %262, %260 ], [ 0, %263 ]
  %266 = sext i32 %265 to i64
  %267 = getelementptr inbounds double, ptr %256, i64 %266
  %268 = load double, ptr %267, align 8, !tbaa !14
  store double %268, ptr %46, align 8, !tbaa !14
  %269 = load ptr, ptr %42, align 8, !tbaa !5
  %270 = load i32, ptr %51, align 4, !tbaa !9
  %271 = sext i32 %270 to i64
  %272 = getelementptr inbounds double, ptr %269, i64 %271
  %273 = load double, ptr %272, align 8, !tbaa !14
  store double %273, ptr %47, align 8, !tbaa !14
  %274 = load ptr, ptr %42, align 8, !tbaa !5
  %275 = load i32, ptr %50, align 4, !tbaa !9
  %276 = sub nsw i32 %275, 1
  %277 = load i32, ptr %51, align 4, !tbaa !9
  %278 = add nsw i32 %277, 1
  %279 = icmp slt i32 %276, %278
  br i1 %279, label %280, label %283

280:                                              ; preds = %264
  %281 = load i32, ptr %50, align 4, !tbaa !9
  %282 = sub nsw i32 %281, 1
  br label %286

283:                                              ; preds = %264
  %284 = load i32, ptr %51, align 4, !tbaa !9
  %285 = add nsw i32 %284, 1
  br label %286

286:                                              ; preds = %283, %280
  %287 = phi i32 [ %282, %280 ], [ %285, %283 ]
  %288 = sext i32 %287 to i64
  %289 = getelementptr inbounds double, ptr %274, i64 %288
  %290 = load double, ptr %289, align 8, !tbaa !14
  store double %290, ptr %48, align 8, !tbaa !14
  %291 = load ptr, ptr %52, align 8, !tbaa !5
  %292 = load i32, ptr %51, align 4, !tbaa !9
  %293 = sext i32 %292 to i64
  %294 = getelementptr inbounds double, ptr %291, i64 %293
  %295 = load double, ptr %294, align 8, !tbaa !14
  store double %295, ptr %19, align 8, !tbaa !14
  %296 = load ptr, ptr %53, align 8, !tbaa !5
  %297 = load i32, ptr %51, align 4, !tbaa !9
  %298 = sext i32 %297 to i64
  %299 = getelementptr inbounds double, ptr %296, i64 %298
  %300 = load double, ptr %299, align 8, !tbaa !14
  store double %300, ptr %20, align 8, !tbaa !14
  %301 = load ptr, ptr %37, align 8, !tbaa !5
  %302 = load i32, ptr %51, align 4, !tbaa !9
  %303 = sext i32 %302 to i64
  %304 = getelementptr inbounds double, ptr %301, i64 %303
  %305 = load double, ptr %304, align 8, !tbaa !14
  store double %305, ptr %33, align 8, !tbaa !14
  %306 = load ptr, ptr %36, align 8, !tbaa !5
  %307 = load i32, ptr %51, align 4, !tbaa !9
  %308 = sext i32 %307 to i64
  %309 = getelementptr inbounds double, ptr %306, i64 %308
  %310 = load double, ptr %309, align 8, !tbaa !14
  store double %310, ptr %34, align 8, !tbaa !14
  %311 = load double, ptr %34, align 8, !tbaa !14
  %312 = load double, ptr %34, align 8, !tbaa !14
  %313 = fmul double %311, %312
  store double %313, ptr %35, align 8, !tbaa !14
  %314 = load double, ptr %20, align 8, !tbaa !14
  %315 = fmul double -2.000000e+00, %314
  %316 = load double, ptr %33, align 8, !tbaa !14
  %317 = load double, ptr %19, align 8, !tbaa !14
  %318 = load double, ptr %39, align 8, !tbaa !14
  %319 = load double, ptr %35, align 8, !tbaa !14
  %320 = fdiv double %318, %319
  %321 = load double, ptr %41, align 8, !tbaa !14
  %322 = fsub double %320, %321
  %323 = load i32, ptr %38, align 4, !tbaa !9
  %324 = sitofp i32 %323 to double
  %325 = fadd double %322, %324
  %326 = fmul double %317, %325
  %327 = call double @llvm.fmuladd.f64(double %315, double %316, double %326)
  store double %327, ptr %32, align 8, !tbaa !14
  %328 = load double, ptr %33, align 8, !tbaa !14
  %329 = fmul double -2.000000e+00, %328
  %330 = load double, ptr %32, align 8, !tbaa !14
  %331 = load double, ptr %19, align 8, !tbaa !14
  %332 = fmul double 2.000000e+00, %331
  %333 = load double, ptr %39, align 8, !tbaa !14
  %334 = fmul double %332, %333
  %335 = load double, ptr %33, align 8, !tbaa !14
  %336 = fmul double %334, %335
  %337 = load double, ptr %35, align 8, !tbaa !14
  %338 = fdiv double %336, %337
  %339 = fneg double %338
  %340 = call double @llvm.fmuladd.f64(double %329, double %330, double %339)
  %341 = load double, ptr %20, align 8, !tbaa !14
  %342 = load i32, ptr %38, align 4, !tbaa !9
  %343 = sitofp i32 %342 to double
  %344 = load double, ptr %41, align 8, !tbaa !14
  %345 = fsub double %343, %344
  %346 = load double, ptr %39, align 8, !tbaa !14
  %347 = fadd double 2.000000e+00, %346
  %348 = load double, ptr %35, align 8, !tbaa !14
  %349 = fdiv double %347, %348
  %350 = fadd double %345, %349
  %351 = call double @llvm.fmuladd.f64(double %341, double %350, double %340)
  store double %351, ptr %27, align 8, !tbaa !14
  %352 = load double, ptr %20, align 8, !tbaa !14
  %353 = fneg double %352
  %354 = load double, ptr %45, align 8, !tbaa !14
  %355 = load double, ptr %18, align 8, !tbaa !14
  %356 = fmul double 2.000000e+00, %355
  %357 = fneg double %356
  %358 = call double @llvm.fmuladd.f64(double %353, double %354, double %357)
  %359 = load double, ptr %19, align 8, !tbaa !14
  %360 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %359, double %358)
  %361 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 0
  store double %360, ptr %361, align 16, !tbaa !14
  %362 = load double, ptr %20, align 8, !tbaa !14
  %363 = load double, ptr %45, align 8, !tbaa !14
  %364 = load double, ptr %18, align 8, !tbaa !14
  %365 = call double @llvm.fmuladd.f64(double %362, double %363, double %364)
  %366 = load double, ptr %19, align 8, !tbaa !14
  %367 = fsub double %365, %366
  %368 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 1
  store double %367, ptr %368, align 8, !tbaa !14
  %369 = load double, ptr %32, align 8, !tbaa !14
  %370 = fneg double %369
  %371 = load double, ptr %45, align 8, !tbaa !14
  %372 = load double, ptr %31, align 8, !tbaa !14
  %373 = fmul double 2.000000e+00, %372
  %374 = fneg double %373
  %375 = call double @llvm.fmuladd.f64(double %370, double %371, double %374)
  %376 = load double, ptr %20, align 8, !tbaa !14
  %377 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %376, double %375)
  %378 = getelementptr inbounds [3 x double], ptr %24, i64 0, i64 0
  store double %377, ptr %378, align 16, !tbaa !14
  %379 = load double, ptr %32, align 8, !tbaa !14
  %380 = load double, ptr %45, align 8, !tbaa !14
  %381 = load double, ptr %31, align 8, !tbaa !14
  %382 = call double @llvm.fmuladd.f64(double %379, double %380, double %381)
  %383 = load double, ptr %20, align 8, !tbaa !14
  %384 = fsub double %382, %383
  %385 = getelementptr inbounds [3 x double], ptr %24, i64 0, i64 1
  store double %384, ptr %385, align 8, !tbaa !14
  %386 = load double, ptr %27, align 8, !tbaa !14
  %387 = fneg double %386
  %388 = load double, ptr %45, align 8, !tbaa !14
  %389 = load double, ptr %26, align 8, !tbaa !14
  %390 = fmul double 2.000000e+00, %389
  %391 = fneg double %390
  %392 = call double @llvm.fmuladd.f64(double %387, double %388, double %391)
  %393 = load double, ptr %32, align 8, !tbaa !14
  %394 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %393, double %392)
  %395 = getelementptr inbounds [3 x double], ptr %25, i64 0, i64 0
  store double %394, ptr %395, align 16, !tbaa !14
  %396 = load double, ptr %27, align 8, !tbaa !14
  %397 = load double, ptr %45, align 8, !tbaa !14
  %398 = load double, ptr %26, align 8, !tbaa !14
  %399 = call double @llvm.fmuladd.f64(double %396, double %397, double %398)
  %400 = load double, ptr %32, align 8, !tbaa !14
  %401 = fsub double %399, %400
  %402 = getelementptr inbounds [3 x double], ptr %25, i64 0, i64 1
  store double %401, ptr %402, align 8, !tbaa !14
  br label %403

403:                                              ; preds = %286, %168
  %404 = load double, ptr %21, align 8, !tbaa !14
  %405 = load double, ptr %46, align 8, !tbaa !14
  %406 = fsub double %404, %405
  %407 = load double, ptr %45, align 8, !tbaa !14
  %408 = fdiv double %406, %407
  %409 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 0
  store double %408, ptr %409, align 16, !tbaa !14
  %410 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 0
  %411 = load double, ptr %410, align 16, !tbaa !14
  %412 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 0
  %413 = load double, ptr %412, align 16, !tbaa !14
  %414 = fmul double %411, %413
  %415 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 1
  store double %414, ptr %415, align 8, !tbaa !14
  %416 = load double, ptr %18, align 8, !tbaa !14
  %417 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 0
  %418 = load double, ptr %417, align 16, !tbaa !14
  %419 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 0
  %420 = load double, ptr %419, align 16, !tbaa !14
  %421 = call double @llvm.fmuladd.f64(double %418, double %420, double %416)
  %422 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 1
  %423 = load double, ptr %422, align 8, !tbaa !14
  %424 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 1
  %425 = load double, ptr %424, align 8, !tbaa !14
  %426 = call double @llvm.fmuladd.f64(double %423, double %425, double %421)
  %427 = load i32, ptr %54, align 4, !tbaa !9
  %428 = sitofp i32 %427 to double
  %429 = fmul double %426, %428
  %430 = load ptr, ptr %13, align 8, !tbaa !5
  %431 = load i32, ptr %49, align 4, !tbaa !9
  %432 = sext i32 %431 to i64
  %433 = getelementptr inbounds double, ptr %430, i64 %432
  store double %429, ptr %433, align 8, !tbaa !14
  %434 = load double, ptr %31, align 8, !tbaa !14
  %435 = getelementptr inbounds [3 x double], ptr %24, i64 0, i64 0
  %436 = load double, ptr %435, align 16, !tbaa !14
  %437 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 0
  %438 = load double, ptr %437, align 16, !tbaa !14
  %439 = call double @llvm.fmuladd.f64(double %436, double %438, double %434)
  %440 = getelementptr inbounds [3 x double], ptr %24, i64 0, i64 1
  %441 = load double, ptr %440, align 8, !tbaa !14
  %442 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 1
  %443 = load double, ptr %442, align 8, !tbaa !14
  %444 = call double @llvm.fmuladd.f64(double %441, double %443, double %439)
  %445 = load i32, ptr %55, align 4, !tbaa !9
  %446 = sitofp i32 %445 to double
  %447 = fmul double %444, %446
  %448 = load ptr, ptr %14, align 8, !tbaa !5
  %449 = load i32, ptr %49, align 4, !tbaa !9
  %450 = sext i32 %449 to i64
  %451 = getelementptr inbounds double, ptr %448, i64 %450
  store double %447, ptr %451, align 8, !tbaa !14
  %452 = load double, ptr %26, align 8, !tbaa !14
  %453 = getelementptr inbounds [3 x double], ptr %25, i64 0, i64 0
  %454 = load double, ptr %453, align 16, !tbaa !14
  %455 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 0
  %456 = load double, ptr %455, align 16, !tbaa !14
  %457 = call double @llvm.fmuladd.f64(double %454, double %456, double %452)
  %458 = getelementptr inbounds [3 x double], ptr %25, i64 0, i64 1
  %459 = load double, ptr %458, align 8, !tbaa !14
  %460 = getelementptr inbounds [2 x double], ptr %22, i64 0, i64 1
  %461 = load double, ptr %460, align 8, !tbaa !14
  %462 = call double @llvm.fmuladd.f64(double %459, double %461, double %457)
  %463 = load i32, ptr %54, align 4, !tbaa !9
  %464 = sitofp i32 %463 to double
  %465 = fmul double %462, %464
  %466 = load ptr, ptr %15, align 8, !tbaa !5
  %467 = load i32, ptr %49, align 4, !tbaa !9
  %468 = sext i32 %467 to i64
  %469 = getelementptr inbounds double, ptr %466, i64 %468
  store double %465, ptr %469, align 8, !tbaa !14
  br label %470

470:                                              ; preds = %403, %151
  %471 = load i32, ptr %49, align 4, !tbaa !9
  %472 = add nsw i32 %471, 1
  store i32 %472, ptr %49, align 4, !tbaa !9
  br label %121, !llvm.loop !24

473:                                              ; preds = %121
  call void @llvm.lifetime.end.p0(i64 4, ptr %55) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %54) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %53) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %52) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %51) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %50) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %49) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %48) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %47) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %46) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %45) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %44) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %43) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %42) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %41) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %40) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %39) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %38) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %37) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %36) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %35) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %34) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %33) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %32) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %31) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %30) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %29) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %26) #4
  call void @llvm.lifetime.end.p0(i64 24, ptr %25) #4
  call void @llvm.lifetime.end.p0(i64 24, ptr %24) #4
  call void @llvm.lifetime.end.p0(i64 24, ptr %23) #4
  call void @llvm.lifetime.end.p0(i64 16, ptr %22) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %19) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %18) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %17) #4
  ret i32 0
}

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #2 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { nocallback nofree nounwind willreturn memory(argmem: write) }
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 = !{!"any pointer", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"int", !7, i64 0}
!11 = !{!12, !6, i64 32}
!12 = !{!"HypersphericalInterpolationStructure", !10, i64 0, !13, i64 8, !13, i64 16, !10, i64 24, !10, i64 28, !6, i64 32, !6, i64 40, !10, i64 48, !6, i64 56, !6, i64 64, !6, i64 72, !6, i64 80, !6, i64 88}
!13 = !{!"double", !7, i64 0}
!14 = !{!13, !13, i64 0}
!15 = !{!12, !6, i64 64}
!16 = !{!12, !6, i64 72}
!17 = !{!12, !10, i64 0}
!18 = !{!12, !13, i64 8}
!19 = !{!12, !6, i64 56}
!20 = !{!12, !13, i64 16}
!21 = !{!12, !10, i64 48}
!22 = !{!12, !6, i64 80}
!23 = !{!12, !6, i64 88}
!24 = distinct !{!24, !25}
!25 = !{!"llvm.loop.mustprogress"}


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

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

; Function Attrs: nounwind uwtable
declare dso_local i32 @ClosedModY(i32 noundef, i32 noundef, ptr noundef, ptr noundef, ptr noundef) local_unnamed_addr #2

; Function Attrs: nounwind uwtable
define dso_local noundef i32 @hyperspherical_Hermite3_interpolation_vector_PhidPhid2Phi(ptr nocapture noundef readonly %0, i32 noundef %1, i32 noundef %2, ptr nocapture noundef readonly %3, ptr nocapture noundef writeonly %4, ptr nocapture noundef writeonly %5, ptr nocapture noundef writeonly %6, ptr nocapture noundef readnone %7) local_unnamed_addr #2 {
  %9 = alloca double, align 8
  %10 = alloca i32, align 4
  %11 = alloca i32, align 4
  %12 = getelementptr inbounds i8, ptr %0, i64 32
  %13 = load ptr, ptr %12, align 8, !tbaa !5
  %14 = sext i32 %2 to i64
  %15 = getelementptr inbounds i32, ptr %13, i64 %14
  %16 = load i32, ptr %15, align 4, !tbaa !12
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %9) #4
  %17 = getelementptr inbounds i8, ptr %0, i64 64
  %18 = load ptr, ptr %17, align 8, !tbaa !13
  %19 = getelementptr inbounds i8, ptr %0, i64 72
  %20 = load ptr, ptr %19, align 8, !tbaa !14
  %21 = load i32, ptr %0, align 8, !tbaa !15
  %22 = sitofp i32 %16 to double
  %23 = fadd double %22, 1.000000e+00
  %24 = fmul double %23, %22
  %25 = getelementptr inbounds i8, ptr %0, i64 8
  %26 = load double, ptr %25, align 8, !tbaa !16
  %27 = fmul double %26, %26
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %10) #4
  store i32 1, ptr %10, align 4, !tbaa !12
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %11) #4
  store i32 1, ptr %11, align 4, !tbaa !12
  %28 = getelementptr inbounds i8, ptr %0, i64 56
  %29 = load ptr, ptr %28, align 8, !tbaa !17
  %30 = getelementptr inbounds i8, ptr %0, i64 16
  %31 = load double, ptr %30, align 8, !tbaa !18
  %32 = getelementptr inbounds i8, ptr %0, i64 48
  %33 = load i32, ptr %32, align 8, !tbaa !19
  %34 = getelementptr inbounds i8, ptr %0, i64 80
  %35 = load ptr, ptr %34, align 8, !tbaa !20
  %36 = mul nsw i32 %33, %2
  %37 = sext i32 %36 to i64
  %38 = getelementptr inbounds double, ptr %35, i64 %37
  %39 = getelementptr inbounds i8, ptr %0, i64 88
  %40 = load ptr, ptr %39, align 8, !tbaa !21
  %41 = getelementptr inbounds double, ptr %40, i64 %37
  %42 = load double, ptr %29, align 8, !tbaa !22
  %43 = add nsw i32 %33, -1
  %44 = sext i32 %43 to i64
  %45 = getelementptr inbounds double, ptr %29, i64 %44
  %46 = load double, ptr %45, align 8, !tbaa !22
  %invariant.gep = getelementptr i8, ptr %29, i64 -8
  %47 = icmp sgt i32 %1, 0
  br i1 %47, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %8
  %48 = sitofp i32 %21 to double
  %49 = fsub double %48, %27
  %50 = fadd double %24, 2.000000e+00
  %wide.trip.count = zext nneg i32 %1 to i64
  br label %51

51:                                               ; preds = %.lr.ph, %171
  %indvars.iv = phi i64 [ 0, %.lr.ph ], [ %indvars.iv.next, %171 ]
  %.0221 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3, %171 ]
  %.0161220 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2163, %171 ]
  %.0164219 = phi i32 [ 0, %.lr.ph ], [ %.3167, %171 ]
  %.0168218 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2170, %171 ]
  %.0172216 = phi double [ %42, %.lr.ph ], [ %.2174, %171 ]
  %.0175215 = phi double [ %42, %.lr.ph ], [ %.2177, %171 ]
  %.0178214 = phi double [ %46, %.lr.ph ], [ %.2180, %171 ]
  %.sroa.0109.0213 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.0109.2, %171 ]
  %.sroa.5110.0212 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.5110.2, %171 ]
  %.sroa.0106.0211 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.0106.2, %171 ]
  %.sroa.5107.0210 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.5107.2, %171 ]
  %.sroa.0.0209 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.0.2, %171 ]
  %.sroa.5.0208 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.5.2, %171 ]
  %.0181207 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3184, %171 ]
  %.0191204 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2193, %171 ]
  %.0194203 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3197, %171 ]
  %52 = getelementptr inbounds double, ptr %3, i64 %indvars.iv
  %53 = load double, ptr %52, align 8, !tbaa !22
  store double %53, ptr %9, align 8, !tbaa !22
  %54 = load i32, ptr %0, align 8, !tbaa !15
  %55 = icmp eq i32 %54, 1
  br i1 %55, label %56, label %61

56:                                               ; preds = %51
  %57 = load double, ptr %25, align 8, !tbaa !16
  %58 = fadd double %57, 2.000000e-01
  %59 = fptosi double %58 to i32
  %60 = call i32 @ClosedModY(i32 noundef %16, i32 noundef %59, ptr noundef nonnull %9, ptr noundef nonnull %10, ptr noundef nonnull %11)
  %.pre = load double, ptr %9, align 8, !tbaa !22
  br label %61

61:                                               ; preds = %56, %51
  %62 = phi double [ %.pre, %56 ], [ %53, %51 ]
  %63 = fcmp olt double %62, %42
  %64 = fcmp ogt double %62, %46
  %or.cond = select i1 %63, i1 true, i1 %64
  br i1 %or.cond, label %65, label %68

65:                                               ; preds = %61
  %66 = getelementptr inbounds double, ptr %4, i64 %indvars.iv
  store double 0.000000e+00, ptr %66, align 8, !tbaa !22
  %67 = getelementptr inbounds double, ptr %5, i64 %indvars.iv
  store double 0.000000e+00, ptr %67, align 8, !tbaa !22
  br label %171

68:                                               ; preds = %61
  %69 = fcmp ogt double %62, %.0175215
  %70 = fcmp olt double %62, %.0178214
  %or.cond201 = select i1 %69, i1 true, i1 %70
  br i1 %or.cond201, label %71, label %152

71:                                               ; preds = %68
  %72 = fcmp ogt double %62, %.0172216
  %or.cond202 = select i1 %72, i1 true, i1 %70
  br i1 %or.cond202, label %73, label %97

73:                                               ; preds = %71
  %74 = fsub double %62, %42
  %75 = fdiv double %74, %31
  %76 = fptosi double %75 to i32
  %77 = call i32 @llvm.smax.i32(i32 %76, i32 0)
  %78 = add nuw nsw i32 %77, 1
  %79 = call i32 @llvm.smin.i32(i32 %43, i32 %78)
  %80 = add nsw i32 %79, -1
  %81 = sext i32 %80 to i64
  %82 = getelementptr inbounds double, ptr %38, i64 %81
  %83 = load double, ptr %82, align 8, !tbaa !22
  %84 = getelementptr inbounds double, ptr %41, i64 %81
  %85 = load double, ptr %84, align 8, !tbaa !22
  %86 = getelementptr inbounds double, ptr %20, i64 %81
  %87 = load double, ptr %86, align 8, !tbaa !22
  %88 = getelementptr inbounds double, ptr %18, i64 %81
  %89 = load double, ptr %88, align 8, !tbaa !22
  %90 = fmul double %89, %89
  %91 = fmul double %85, -2.000000e+00
  %92 = fdiv double %24, %90
  %93 = fsub double %92, %27
  %94 = fadd double %93, %48
  %95 = fmul double %83, %94
  %96 = call double @llvm.fmuladd.f64(double %91, double %87, double %95)
  br label %99

97:                                               ; preds = %71
  %98 = add nsw i32 %.0164219, 1
  %.pre223 = fmul double %.0168218, -2.000000e+00
  br label %99

99:                                               ; preds = %97, %73
  %.pre-phi = phi double [ %.pre223, %97 ], [ %91, %73 ]
  %.1195 = phi double [ %.0168218, %97 ], [ %85, %73 ]
  %.1182 = phi double [ %.0191204, %97 ], [ %96, %73 ]
  %.1165 = phi i32 [ %98, %97 ], [ %79, %73 ]
  %.1 = phi double [ %.0161220, %97 ], [ %83, %73 ]
  %100 = call i32 @llvm.smax.i32(i32 %.1165, i32 1)
  %101 = zext nneg i32 %100 to i64
  %gep = getelementptr double, ptr %invariant.gep, i64 %101
  %102 = load double, ptr %gep, align 8, !tbaa !22
  %103 = sext i32 %.1165 to i64
  %104 = getelementptr inbounds double, ptr %29, i64 %103
  %105 = load double, ptr %104, align 8, !tbaa !22
  %106 = add nsw i32 %.1165, 1
  %107 = call i32 @llvm.smin.i32(i32 %43, i32 %106)
  %108 = sext i32 %107 to i64
  %109 = getelementptr inbounds double, ptr %29, i64 %108
  %110 = load double, ptr %109, align 8, !tbaa !22
  %111 = getelementptr inbounds double, ptr %38, i64 %103
  %112 = load double, ptr %111, align 8, !tbaa !22
  %113 = getelementptr inbounds double, ptr %41, i64 %103
  %114 = load double, ptr %113, align 8, !tbaa !22
  %115 = getelementptr inbounds double, ptr %20, i64 %103
  %116 = load double, ptr %115, align 8, !tbaa !22
  %117 = getelementptr inbounds double, ptr %18, i64 %103
  %118 = load double, ptr %117, align 8, !tbaa !22
  %119 = fmul double %118, %118
  %120 = fmul double %114, -2.000000e+00
  %121 = fdiv double %24, %119
  %122 = fsub double %121, %27
  %123 = fadd double %122, %48
  %124 = fmul double %112, %123
  %125 = call double @llvm.fmuladd.f64(double %120, double %116, double %124)
  %126 = fmul double %116, -2.000000e+00
  %127 = fmul double %112, -2.000000e+00
  %128 = fmul double %24, %127
  %129 = fmul double %128, %116
  %130 = fdiv double %129, %119
  %131 = call double @llvm.fmuladd.f64(double %126, double %125, double %130)
  %132 = fdiv double %50, %119
  %133 = fadd double %49, %132
  %134 = call double @llvm.fmuladd.f64(double %114, double %133, double %131)
  %135 = fneg double %114
  %136 = fmul double %.1, -2.000000e+00
  %137 = call double @llvm.fmuladd.f64(double %135, double %31, double %136)
  %138 = call double @llvm.fmuladd.f64(double %112, double 2.000000e+00, double %137)
  %139 = call double @llvm.fmuladd.f64(double %114, double %31, double %.1)
  %140 = fsub double %139, %112
  %141 = fneg double %125
  %142 = call double @llvm.fmuladd.f64(double %141, double %31, double %.pre-phi)
  %143 = call double @llvm.fmuladd.f64(double %114, double 2.000000e+00, double %142)
  %144 = call double @llvm.fmuladd.f64(double %125, double %31, double %.1195)
  %145 = fsub double %144, %114
  %146 = fneg double %134
  %147 = fmul double %.1182, -2.000000e+00
  %148 = call double @llvm.fmuladd.f64(double %146, double %31, double %147)
  %149 = call double @llvm.fmuladd.f64(double %125, double 2.000000e+00, double %148)
  %150 = call double @llvm.fmuladd.f64(double %134, double %31, double %.1182)
  %151 = fsub double %150, %125
  br label %152

152:                                              ; preds = %68, %99
  %.2196 = phi double [ %.1195, %99 ], [ %.0194203, %68 ]
  %.1192 = phi double [ %125, %99 ], [ %.0191204, %68 ]
  %.2183 = phi double [ %.1182, %99 ], [ %.0181207, %68 ]
  %.sroa.5.1 = phi double [ %151, %99 ], [ %.sroa.5.0208, %68 ]
  %.sroa.0.1 = phi double [ %149, %99 ], [ %.sroa.0.0209, %68 ]
  %.sroa.5107.1 = phi double [ %145, %99 ], [ %.sroa.5107.0210, %68 ]
  %.sroa.0106.1 = phi double [ %143, %99 ], [ %.sroa.0106.0211, %68 ]
  %.sroa.5110.1 = phi double [ %140, %99 ], [ %.sroa.5110.0212, %68 ]
  %.sroa.0109.1 = phi double [ %138, %99 ], [ %.sroa.0109.0213, %68 ]
  %.1179 = phi double [ %102, %99 ], [ %.0178214, %68 ]
  %.1176 = phi double [ %105, %99 ], [ %.0175215, %68 ]
  %.1173 = phi double [ %110, %99 ], [ %.0172216, %68 ]
  %.1169 = phi double [ %114, %99 ], [ %.0168218, %68 ]
  %.2166 = phi i32 [ %.1165, %99 ], [ %.0164219, %68 ]
  %.1162 = phi double [ %112, %99 ], [ %.0161220, %68 ]
  %.2 = phi double [ %.1, %99 ], [ %.0221, %68 ]
  %153 = fsub double %62, %.1179
  %154 = fdiv double %153, %31
  %155 = fmul double %154, %154
  %156 = call double @llvm.fmuladd.f64(double %.sroa.0109.1, double %154, double %.2)
  %157 = call double @llvm.fmuladd.f64(double %.sroa.5110.1, double %155, double %156)
  %158 = load i32, ptr %10, align 4, !tbaa !12
  %159 = sitofp i32 %158 to double
  %160 = fmul double %157, %159
  %161 = getelementptr inbounds double, ptr %4, i64 %indvars.iv
  store double %160, ptr %161, align 8, !tbaa !22
  %162 = call double @llvm.fmuladd.f64(double %.sroa.0106.1, double %154, double %.2196)
  %163 = call double @llvm.fmuladd.f64(double %.sroa.5107.1, double %155, double %162)
  %164 = load i32, ptr %11, align 4, !tbaa !12
  %165 = sitofp i32 %164 to double
  %166 = fmul double %163, %165
  %167 = getelementptr inbounds double, ptr %5, i64 %indvars.iv
  store double %166, ptr %167, align 8, !tbaa !22
  %168 = call double @llvm.fmuladd.f64(double %.sroa.0.1, double %154, double %.2183)
  %169 = call double @llvm.fmuladd.f64(double %.sroa.5.1, double %155, double %168)
  %170 = fmul double %169, %159
  br label %171

171:                                              ; preds = %152, %65
  %.sink = phi double [ %170, %152 ], [ 0.000000e+00, %65 ]
  %.3197 = phi double [ %.2196, %152 ], [ %.0194203, %65 ]
  %.2193 = phi double [ %.1192, %152 ], [ %.0191204, %65 ]
  %.3184 = phi double [ %.2183, %152 ], [ %.0181207, %65 ]
  %.sroa.5.2 = phi double [ %.sroa.5.1, %152 ], [ %.sroa.5.0208, %65 ]
  %.sroa.0.2 = phi double [ %.sroa.0.1, %152 ], [ %.sroa.0.0209, %65 ]
  %.sroa.5107.2 = phi double [ %.sroa.5107.1, %152 ], [ %.sroa.5107.0210, %65 ]
  %.sroa.0106.2 = phi double [ %.sroa.0106.1, %152 ], [ %.sroa.0106.0211, %65 ]
  %.sroa.5110.2 = phi double [ %.sroa.5110.1, %152 ], [ %.sroa.5110.0212, %65 ]
  %.sroa.0109.2 = phi double [ %.sroa.0109.1, %152 ], [ %.sroa.0109.0213, %65 ]
  %.2180 = phi double [ %.1179, %152 ], [ %.0178214, %65 ]
  %.2177 = phi double [ %.1176, %152 ], [ %.0175215, %65 ]
  %.2174 = phi double [ %.1173, %152 ], [ %.0172216, %65 ]
  %.2170 = phi double [ %.1169, %152 ], [ %.0168218, %65 ]
  %.3167 = phi i32 [ %.2166, %152 ], [ %.0164219, %65 ]
  %.2163 = phi double [ %.1162, %152 ], [ %.0161220, %65 ]
  %.3 = phi double [ %.2, %152 ], [ %.0221, %65 ]
  %172 = getelementptr inbounds double, ptr %6, i64 %indvars.iv
  store double %.sink, ptr %172, align 8, !tbaa !22
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, %wide.trip.count
  br i1 %exitcond.not, label %._crit_edge, label %51, !llvm.loop !23

._crit_edge:                                      ; preds = %171, %8
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %11) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %10) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %9) #4
  ret i32 0
}

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

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

attributes #0 = { mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #2 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
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, !11, i64 32}
!6 = !{!"HypersphericalInterpolationStructure", !7, i64 0, !10, i64 8, !10, i64 16, !7, i64 24, !7, i64 28, !11, i64 32, !11, i64 40, !7, i64 48, !11, i64 56, !11, i64 64, !11, i64 72, !11, i64 80, !11, i64 88}
!7 = !{!"int", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C/C++ TBAA"}
!10 = !{!"double", !8, i64 0}
!11 = !{!"any pointer", !8, i64 0}
!12 = !{!7, !7, i64 0}
!13 = !{!6, !11, i64 64}
!14 = !{!6, !11, i64 72}
!15 = !{!6, !7, i64 0}
!16 = !{!6, !10, i64 8}
!17 = !{!6, !11, i64 56}
!18 = !{!6, !10, i64 16}
!19 = !{!6, !7, i64 48}
!20 = !{!6, !11, i64 80}
!21 = !{!6, !11, i64 88}
!22 = !{!10, !10, i64 0}
!23 = distinct !{!23, !24}
!24 = !{!"llvm.loop.mustprogress"}


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