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

119:                                              ; preds = %492, %7
  %120 = load i32, ptr %47, align 4, !tbaa !9
  %121 = load i32, ptr %9, align 4, !tbaa !9
  %122 = icmp slt i32 %120, %121
  br i1 %122, label %123, label %495

123:                                              ; preds = %119
  %124 = load ptr, ptr %11, align 8, !tbaa !5
  %125 = load i32, ptr %47, align 4, !tbaa !9
  %126 = sext i32 %125 to i64
  %127 = getelementptr inbounds double, ptr %124, i64 %126
  %128 = load double, ptr %127, align 8, !tbaa !14
  store double %128, ptr %20, align 8, !tbaa !14
  %129 = load ptr, ptr %8, align 8, !tbaa !5
  %130 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %129, i32 0, i32 0
  %131 = load i32, ptr %130, align 8, !tbaa !17
  %132 = icmp eq i32 %131, 1
  br i1 %132, label %133, label %141

133:                                              ; preds = %123
  %134 = load i32, ptr %15, align 4, !tbaa !9
  %135 = load ptr, ptr %8, align 8, !tbaa !5
  %136 = getelementptr inbounds %struct.HypersphericalInterpolationStructure, ptr %135, i32 0, i32 1
  %137 = load double, ptr %136, align 8, !tbaa !18
  %138 = fadd double %137, 2.000000e-01
  %139 = fptosi double %138 to i32
  %140 = call i32 @ClosedModY(i32 noundef %134, i32 noundef %139, ptr noundef %20, ptr noundef %52, ptr noundef %53)
  br label %141

141:                                              ; preds = %133, %123
  %142 = load double, ptr %20, align 8, !tbaa !14
  %143 = load double, ptr %41, align 8, !tbaa !14
  %144 = fcmp olt double %142, %143
  br i1 %144, label %149, label %145

145:                                              ; preds = %141
  %146 = load double, ptr %20, align 8, !tbaa !14
  %147 = load double, ptr %42, align 8, !tbaa !14
  %148 = fcmp ogt double %146, %147
  br i1 %148, label %149, label %158

149:                                              ; preds = %145, %141
  %150 = load ptr, ptr %12, align 8, !tbaa !5
  %151 = load i32, ptr %47, align 4, !tbaa !9
  %152 = sext i32 %151 to i64
  %153 = getelementptr inbounds double, ptr %150, i64 %152
  store double 0.000000e+00, ptr %153, align 8, !tbaa !14
  %154 = load ptr, ptr %13, align 8, !tbaa !5
  %155 = load i32, ptr %47, align 4, !tbaa !9
  %156 = sext i32 %155 to i64
  %157 = getelementptr inbounds double, ptr %154, i64 %156
  store double 0.000000e+00, ptr %157, align 8, !tbaa !14
  br label %492

158:                                              ; preds = %145
  %159 = load double, ptr %20, align 8, !tbaa !14
  %160 = load double, ptr %45, align 8, !tbaa !14
  %161 = fcmp ogt double %159, %160
  br i1 %161, label %166, label %162

162:                                              ; preds = %158
  %163 = load double, ptr %20, align 8, !tbaa !14
  %164 = load double, ptr %44, align 8, !tbaa !14
  %165 = fcmp olt double %163, %164
  br i1 %165, label %166, label %427

166:                                              ; preds = %162, %158
  %167 = load double, ptr %20, align 8, !tbaa !14
  %168 = load double, ptr %46, align 8, !tbaa !14
  %169 = fcmp ogt double %167, %168
  br i1 %169, label %174, label %170

170:                                              ; preds = %166
  %171 = load double, ptr %20, align 8, !tbaa !14
  %172 = load double, ptr %44, align 8, !tbaa !14
  %173 = fcmp olt double %171, %172
  br i1 %173, label %174, label %265

174:                                              ; preds = %170, %166
  %175 = load double, ptr %20, align 8, !tbaa !14
  %176 = load double, ptr %41, align 8, !tbaa !14
  %177 = fsub double %175, %176
  %178 = load double, ptr %43, align 8, !tbaa !14
  %179 = fdiv double %177, %178
  %180 = fptosi double %179 to i32
  %181 = add nsw i32 %180, 1
  store i32 %181, ptr %49, align 4, !tbaa !9
  %182 = load i32, ptr %49, align 4, !tbaa !9
  %183 = icmp slt i32 1, %182
  br i1 %183, label %184, label %186

184:                                              ; preds = %174
  %185 = load i32, ptr %49, align 4, !tbaa !9
  br label %187

186:                                              ; preds = %174
  br label %187

187:                                              ; preds = %186, %184
  %188 = phi i32 [ %185, %184 ], [ 1, %186 ]
  store i32 %188, ptr %49, align 4, !tbaa !9
  %189 = load i32, ptr %48, align 4, !tbaa !9
  %190 = sub nsw i32 %189, 1
  %191 = load i32, ptr %49, align 4, !tbaa !9
  %192 = icmp slt i32 %190, %191
  br i1 %192, label %193, label %196

193:                                              ; preds = %187
  %194 = load i32, ptr %48, align 4, !tbaa !9
  %195 = sub nsw i32 %194, 1
  br label %198

196:                                              ; preds = %187
  %197 = load i32, ptr %49, align 4, !tbaa !9
  br label %198

198:                                              ; preds = %196, %193
  %199 = phi i32 [ %195, %193 ], [ %197, %196 ]
  store i32 %199, ptr %49, align 4, !tbaa !9
  %200 = load ptr, ptr %50, align 8, !tbaa !5
  %201 = load i32, ptr %49, align 4, !tbaa !9
  %202 = sub nsw i32 %201, 1
  %203 = sext i32 %202 to i64
  %204 = getelementptr inbounds double, ptr %200, i64 %203
  %205 = load double, ptr %204, align 8, !tbaa !14
  store double %205, ptr %16, align 8, !tbaa !14
  %206 = load ptr, ptr %51, align 8, !tbaa !5
  %207 = load i32, ptr %49, 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 %27, align 8, !tbaa !5
  %213 = load i32, ptr %49, 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 %28, align 8, !tbaa !14
  %218 = load ptr, ptr %26, align 8, !tbaa !5
  %219 = load i32, ptr %49, 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 %30, align 8, !tbaa !14
  %224 = load double, ptr %30, align 8, !tbaa !14
  %225 = load double, ptr %30, align 8, !tbaa !14
  %226 = fmul double %224, %225
  store double %226, ptr %32, align 8, !tbaa !14
  %227 = load double, ptr %18, align 8, !tbaa !14
  %228 = fmul double -2.000000e+00, %227
  %229 = load double, ptr %28, align 8, !tbaa !14
  %230 = load double, ptr %16, align 8, !tbaa !14
  %231 = load double, ptr %37, align 8, !tbaa !14
  %232 = load double, ptr %32, align 8, !tbaa !14
  %233 = fdiv double %231, %232
  %234 = load double, ptr %39, align 8, !tbaa !14
  %235 = fsub double %233, %234
  %236 = load i32, ptr %36, align 4, !tbaa !9
  %237 = sitofp i32 %236 to double
  %238 = fadd double %235, %237
  %239 = fmul double %230, %238
  %240 = call double @llvm.fmuladd.f64(double %228, double %229, double %239)
  store double %240, ptr %34, align 8, !tbaa !14
  %241 = load double, ptr %28, align 8, !tbaa !14
  %242 = fmul double -2.000000e+00, %241
  %243 = load double, ptr %34, align 8, !tbaa !14
  %244 = load double, ptr %16, align 8, !tbaa !14
  %245 = fmul double 2.000000e+00, %244
  %246 = load double, ptr %37, align 8, !tbaa !14
  %247 = fmul double %245, %246
  %248 = load double, ptr %28, align 8, !tbaa !14
  %249 = fmul double %247, %248
  %250 = load double, ptr %32, align 8, !tbaa !14
  %251 = fdiv double %249, %250
  %252 = fneg double %251
  %253 = call double @llvm.fmuladd.f64(double %242, double %243, double %252)
  %254 = load double, ptr %18, align 8, !tbaa !14
  %255 = load i32, ptr %36, align 4, !tbaa !9
  %256 = sitofp i32 %255 to double
  %257 = load double, ptr %39, align 8, !tbaa !14
  %258 = fsub double %256, %257
  %259 = load double, ptr %37, align 8, !tbaa !14
  %260 = fadd double 2.000000e+00, %259
  %261 = load double, ptr %32, align 8, !tbaa !14
  %262 = fdiv double %260, %261
  %263 = fadd double %258, %262
  %264 = call double @llvm.fmuladd.f64(double %254, double %263, double %253)
  store double %264, ptr %24, align 8, !tbaa !14
  br label %274

265:                                              ; preds = %170
  %266 = load i32, ptr %49, align 4, !tbaa !9
  %267 = add nsw i32 %266, 1
  store i32 %267, ptr %49, align 4, !tbaa !9
  %268 = load double, ptr %17, align 8, !tbaa !14
  store double %268, ptr %16, align 8, !tbaa !14
  %269 = load double, ptr %19, align 8, !tbaa !14
  store double %269, ptr %18, align 8, !tbaa !14
  %270 = load double, ptr %35, align 8, !tbaa !14
  store double %270, ptr %34, align 8, !tbaa !14
  %271 = load double, ptr %31, align 8, !tbaa !14
  store double %271, ptr %30, align 8, !tbaa !14
  %272 = load double, ptr %29, align 8, !tbaa !14
  store double %272, ptr %28, align 8, !tbaa !14
  %273 = load double, ptr %25, align 8, !tbaa !14
  store double %273, ptr %24, align 8, !tbaa !14
  br label %274

274:                                              ; preds = %265, %198
  %275 = load ptr, ptr %40, align 8, !tbaa !5
  %276 = load i32, ptr %49, align 4, !tbaa !9
  %277 = sub nsw i32 %276, 1
  %278 = icmp slt i32 0, %277
  br i1 %278, label %279, label %282

279:                                              ; preds = %274
  %280 = load i32, ptr %49, align 4, !tbaa !9
  %281 = sub nsw i32 %280, 1
  br label %283

282:                                              ; preds = %274
  br label %283

283:                                              ; preds = %282, %279
  %284 = phi i32 [ %281, %279 ], [ 0, %282 ]
  %285 = sext i32 %284 to i64
  %286 = getelementptr inbounds double, ptr %275, i64 %285
  %287 = load double, ptr %286, align 8, !tbaa !14
  store double %287, ptr %44, align 8, !tbaa !14
  %288 = load ptr, ptr %40, align 8, !tbaa !5
  %289 = load i32, ptr %49, align 4, !tbaa !9
  %290 = sext i32 %289 to i64
  %291 = getelementptr inbounds double, ptr %288, i64 %290
  %292 = load double, ptr %291, align 8, !tbaa !14
  store double %292, ptr %45, align 8, !tbaa !14
  %293 = load ptr, ptr %40, align 8, !tbaa !5
  %294 = load i32, ptr %48, align 4, !tbaa !9
  %295 = sub nsw i32 %294, 1
  %296 = load i32, ptr %49, align 4, !tbaa !9
  %297 = add nsw i32 %296, 1
  %298 = icmp slt i32 %295, %297
  br i1 %298, label %299, label %302

299:                                              ; preds = %283
  %300 = load i32, ptr %48, align 4, !tbaa !9
  %301 = sub nsw i32 %300, 1
  br label %305

302:                                              ; preds = %283
  %303 = load i32, ptr %49, align 4, !tbaa !9
  %304 = add nsw i32 %303, 1
  br label %305

305:                                              ; preds = %302, %299
  %306 = phi i32 [ %301, %299 ], [ %304, %302 ]
  %307 = sext i32 %306 to i64
  %308 = getelementptr inbounds double, ptr %293, i64 %307
  %309 = load double, ptr %308, align 8, !tbaa !14
  store double %309, ptr %46, align 8, !tbaa !14
  %310 = load ptr, ptr %50, align 8, !tbaa !5
  %311 = load i32, ptr %49, align 4, !tbaa !9
  %312 = sext i32 %311 to i64
  %313 = getelementptr inbounds double, ptr %310, i64 %312
  %314 = load double, ptr %313, align 8, !tbaa !14
  store double %314, ptr %17, align 8, !tbaa !14
  %315 = load ptr, ptr %51, align 8, !tbaa !5
  %316 = load i32, ptr %49, align 4, !tbaa !9
  %317 = sext i32 %316 to i64
  %318 = getelementptr inbounds double, ptr %315, i64 %317
  %319 = load double, ptr %318, align 8, !tbaa !14
  store double %319, ptr %19, align 8, !tbaa !14
  %320 = load ptr, ptr %27, align 8, !tbaa !5
  %321 = load i32, ptr %49, align 4, !tbaa !9
  %322 = sext i32 %321 to i64
  %323 = getelementptr inbounds double, ptr %320, i64 %322
  %324 = load double, ptr %323, align 8, !tbaa !14
  store double %324, ptr %29, align 8, !tbaa !14
  %325 = load ptr, ptr %26, align 8, !tbaa !5
  %326 = load i32, ptr %49, align 4, !tbaa !9
  %327 = sext i32 %326 to i64
  %328 = getelementptr inbounds double, ptr %325, i64 %327
  %329 = load double, ptr %328, align 8, !tbaa !14
  store double %329, ptr %31, align 8, !tbaa !14
  %330 = load double, ptr %31, align 8, !tbaa !14
  %331 = load double, ptr %31, align 8, !tbaa !14
  %332 = fmul double %330, %331
  store double %332, ptr %33, align 8, !tbaa !14
  %333 = load double, ptr %19, align 8, !tbaa !14
  %334 = fmul double -2.000000e+00, %333
  %335 = load double, ptr %29, align 8, !tbaa !14
  %336 = load double, ptr %17, align 8, !tbaa !14
  %337 = load double, ptr %37, align 8, !tbaa !14
  %338 = load double, ptr %33, align 8, !tbaa !14
  %339 = fdiv double %337, %338
  %340 = load double, ptr %39, align 8, !tbaa !14
  %341 = fsub double %339, %340
  %342 = load i32, ptr %36, align 4, !tbaa !9
  %343 = sitofp i32 %342 to double
  %344 = fadd double %341, %343
  %345 = fmul double %336, %344
  %346 = call double @llvm.fmuladd.f64(double %334, double %335, double %345)
  store double %346, ptr %35, align 8, !tbaa !14
  %347 = load double, ptr %29, align 8, !tbaa !14
  %348 = fmul double -2.000000e+00, %347
  %349 = load double, ptr %35, align 8, !tbaa !14
  %350 = load double, ptr %17, align 8, !tbaa !14
  %351 = fmul double 2.000000e+00, %350
  %352 = load double, ptr %37, align 8, !tbaa !14
  %353 = fmul double %351, %352
  %354 = load double, ptr %29, align 8, !tbaa !14
  %355 = fmul double %353, %354
  %356 = load double, ptr %33, align 8, !tbaa !14
  %357 = fdiv double %355, %356
  %358 = fneg double %357
  %359 = call double @llvm.fmuladd.f64(double %348, double %349, double %358)
  %360 = load double, ptr %19, align 8, !tbaa !14
  %361 = load i32, ptr %36, align 4, !tbaa !9
  %362 = sitofp i32 %361 to double
  %363 = load double, ptr %39, align 8, !tbaa !14
  %364 = fsub double %362, %363
  %365 = load double, ptr %37, align 8, !tbaa !14
  %366 = fadd double 2.000000e+00, %365
  %367 = load double, ptr %33, align 8, !tbaa !14
  %368 = fdiv double %366, %367
  %369 = fadd double %364, %368
  %370 = call double @llvm.fmuladd.f64(double %360, double %369, double %359)
  store double %370, ptr %25, align 8, !tbaa !14
  %371 = load double, ptr %18, align 8, !tbaa !14
  %372 = load double, ptr %43, align 8, !tbaa !14
  %373 = fmul double %371, %372
  %374 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 0
  store double %373, ptr %374, align 16, !tbaa !14
  %375 = load double, ptr %18, align 8, !tbaa !14
  %376 = fmul double -2.000000e+00, %375
  %377 = load double, ptr %43, align 8, !tbaa !14
  %378 = load double, ptr %19, align 8, !tbaa !14
  %379 = load double, ptr %43, align 8, !tbaa !14
  %380 = fmul double %378, %379
  %381 = fneg double %380
  %382 = call double @llvm.fmuladd.f64(double %376, double %377, double %381)
  %383 = load double, ptr %16, align 8, !tbaa !14
  %384 = call double @llvm.fmuladd.f64(double -3.000000e+00, double %383, double %382)
  %385 = load double, ptr %17, align 8, !tbaa !14
  %386 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %385, double %384)
  %387 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 1
  store double %386, ptr %387, align 8, !tbaa !14
  %388 = load double, ptr %18, align 8, !tbaa !14
  %389 = load double, ptr %43, align 8, !tbaa !14
  %390 = load double, ptr %19, align 8, !tbaa !14
  %391 = load double, ptr %43, align 8, !tbaa !14
  %392 = fmul double %390, %391
  %393 = call double @llvm.fmuladd.f64(double %388, double %389, double %392)
  %394 = load double, ptr %16, align 8, !tbaa !14
  %395 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %394, double %393)
  %396 = load double, ptr %17, align 8, !tbaa !14
  %397 = call double @llvm.fmuladd.f64(double -2.000000e+00, double %396, double %395)
  %398 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 2
  store double %397, ptr %398, align 16, !tbaa !14
  %399 = load double, ptr %24, align 8, !tbaa !14
  %400 = load double, ptr %43, align 8, !tbaa !14
  %401 = fmul double %399, %400
  %402 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 0
  store double %401, ptr %402, align 16, !tbaa !14
  %403 = load double, ptr %24, align 8, !tbaa !14
  %404 = fmul double -2.000000e+00, %403
  %405 = load double, ptr %43, align 8, !tbaa !14
  %406 = load double, ptr %25, align 8, !tbaa !14
  %407 = load double, ptr %43, align 8, !tbaa !14
  %408 = fmul double %406, %407
  %409 = fneg double %408
  %410 = call double @llvm.fmuladd.f64(double %404, double %405, double %409)
  %411 = load double, ptr %34, align 8, !tbaa !14
  %412 = call double @llvm.fmuladd.f64(double -3.000000e+00, double %411, double %410)
  %413 = load double, ptr %35, align 8, !tbaa !14
  %414 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %413, double %412)
  %415 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 1
  store double %414, ptr %415, align 8, !tbaa !14
  %416 = load double, ptr %24, align 8, !tbaa !14
  %417 = load double, ptr %43, align 8, !tbaa !14
  %418 = load double, ptr %25, align 8, !tbaa !14
  %419 = load double, ptr %43, align 8, !tbaa !14
  %420 = fmul double %418, %419
  %421 = call double @llvm.fmuladd.f64(double %416, double %417, double %420)
  %422 = load double, ptr %34, align 8, !tbaa !14
  %423 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %422, double %421)
  %424 = load double, ptr %35, align 8, !tbaa !14
  %425 = call double @llvm.fmuladd.f64(double -2.000000e+00, double %424, double %423)
  %426 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 2
  store double %425, ptr %426, align 16, !tbaa !14
  br label %427

427:                                              ; preds = %305, %162
  %428 = load double, ptr %20, align 8, !tbaa !14
  %429 = load double, ptr %44, align 8, !tbaa !14
  %430 = fsub double %428, %429
  %431 = load double, ptr %43, align 8, !tbaa !14
  %432 = fdiv double %430, %431
  %433 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  store double %432, ptr %433, align 16, !tbaa !14
  %434 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  %435 = load double, ptr %434, align 16, !tbaa !14
  %436 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  %437 = load double, ptr %436, align 16, !tbaa !14
  %438 = fmul double %435, %437
  %439 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 1
  store double %438, ptr %439, align 8, !tbaa !14
  %440 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 1
  %441 = load double, ptr %440, align 8, !tbaa !14
  %442 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  %443 = load double, ptr %442, align 16, !tbaa !14
  %444 = fmul double %441, %443
  %445 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 2
  store double %444, ptr %445, align 16, !tbaa !14
  %446 = load double, ptr %16, align 8, !tbaa !14
  %447 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 0
  %448 = load double, ptr %447, align 16, !tbaa !14
  %449 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  %450 = load double, ptr %449, align 16, !tbaa !14
  %451 = call double @llvm.fmuladd.f64(double %448, double %450, double %446)
  %452 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 1
  %453 = load double, ptr %452, align 8, !tbaa !14
  %454 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 1
  %455 = load double, ptr %454, align 8, !tbaa !14
  %456 = call double @llvm.fmuladd.f64(double %453, double %455, double %451)
  %457 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 2
  %458 = load double, ptr %457, align 16, !tbaa !14
  %459 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 2
  %460 = load double, ptr %459, align 16, !tbaa !14
  %461 = call double @llvm.fmuladd.f64(double %458, double %460, double %456)
  %462 = load i32, ptr %52, align 4, !tbaa !9
  %463 = sitofp i32 %462 to double
  %464 = fmul double %461, %463
  %465 = load ptr, ptr %12, align 8, !tbaa !5
  %466 = load i32, ptr %47, align 4, !tbaa !9
  %467 = sext i32 %466 to i64
  %468 = getelementptr inbounds double, ptr %465, i64 %467
  store double %464, ptr %468, align 8, !tbaa !14
  %469 = load double, ptr %34, align 8, !tbaa !14
  %470 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 0
  %471 = load double, ptr %470, align 16, !tbaa !14
  %472 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  %473 = load double, ptr %472, align 16, !tbaa !14
  %474 = call double @llvm.fmuladd.f64(double %471, double %473, double %469)
  %475 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 1
  %476 = load double, ptr %475, align 8, !tbaa !14
  %477 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 1
  %478 = load double, ptr %477, align 8, !tbaa !14
  %479 = call double @llvm.fmuladd.f64(double %476, double %478, double %474)
  %480 = getelementptr inbounds [3 x double], ptr %23, i64 0, i64 2
  %481 = load double, ptr %480, align 16, !tbaa !14
  %482 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 2
  %483 = load double, ptr %482, align 16, !tbaa !14
  %484 = call double @llvm.fmuladd.f64(double %481, double %483, double %479)
  %485 = load i32, ptr %52, align 4, !tbaa !9
  %486 = sitofp i32 %485 to double
  %487 = fmul double %484, %486
  %488 = load ptr, ptr %13, align 8, !tbaa !5
  %489 = load i32, ptr %47, align 4, !tbaa !9
  %490 = sext i32 %489 to i64
  %491 = getelementptr inbounds double, ptr %488, i64 %490
  store double %487, ptr %491, align 8, !tbaa !14
  br label %492

492:                                              ; preds = %427, %149
  %493 = load i32, ptr %47, align 4, !tbaa !9
  %494 = add nsw i32 %493, 1
  store i32 %494, ptr %47, align 4, !tbaa !9
  br label %119, !llvm.loop !24

495:                                              ; preds = %119
  call void @llvm.lifetime.end.p0(i64 4, ptr %53) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %52) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %51) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %50) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %49) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %48) #4
  call void @llvm.lifetime.end.p0(i64 4, 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 8, ptr %38) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %37) #4
  call void @llvm.lifetime.end.p0(i64 4, 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 8, ptr %25) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #4
  call void @llvm.lifetime.end.p0(i64 24, ptr %23) #4
  call void @llvm.lifetime.end.p0(i64 24, ptr %22) #4
  call void @llvm.lifetime.end.p0(i64 24, 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 8, ptr %17) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %15) #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_Hermite4_interpolation_vector_Phid2Phi(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 readnone %6) local_unnamed_addr #2 {
  %8 = alloca double, align 8
  %9 = alloca i32, align 4
  %10 = alloca i32, align 4
  %11 = getelementptr inbounds i8, ptr %0, i64 32
  %12 = load ptr, ptr %11, align 8, !tbaa !5
  %13 = sext i32 %2 to i64
  %14 = getelementptr inbounds i32, ptr %12, i64 %13
  %15 = load i32, ptr %14, align 4, !tbaa !12
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %8) #4
  %16 = getelementptr inbounds i8, ptr %0, i64 64
  %17 = load ptr, ptr %16, align 8, !tbaa !13
  %18 = getelementptr inbounds i8, ptr %0, i64 72
  %19 = load ptr, ptr %18, align 8, !tbaa !14
  %20 = load i32, ptr %0, align 8, !tbaa !15
  %21 = sitofp i32 %15 to double
  %22 = fadd double %21, 1.000000e+00
  %23 = fmul double %22, %21
  %24 = getelementptr inbounds i8, ptr %0, i64 8
  %25 = load double, ptr %24, align 8, !tbaa !16
  %26 = fmul double %25, %25
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %9) #4
  store i32 1, ptr %9, align 4, !tbaa !12
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %10) #4
  store i32 1, ptr %10, align 4, !tbaa !12
  %27 = getelementptr inbounds i8, ptr %0, i64 56
  %28 = load ptr, ptr %27, align 8, !tbaa !17
  %29 = getelementptr inbounds i8, ptr %0, i64 16
  %30 = load double, ptr %29, align 8, !tbaa !18
  %31 = getelementptr inbounds i8, ptr %0, i64 48
  %32 = load i32, ptr %31, align 8, !tbaa !19
  %33 = getelementptr inbounds i8, ptr %0, i64 80
  %34 = load ptr, ptr %33, align 8, !tbaa !20
  %35 = mul nsw i32 %32, %2
  %36 = sext i32 %35 to i64
  %37 = getelementptr inbounds double, ptr %34, i64 %36
  %38 = getelementptr inbounds i8, ptr %0, i64 88
  %39 = load ptr, ptr %38, align 8, !tbaa !21
  %40 = getelementptr inbounds double, ptr %39, i64 %36
  %41 = load double, ptr %28, align 8, !tbaa !22
  %42 = add nsw i32 %32, -1
  %43 = sext i32 %42 to i64
  %44 = getelementptr inbounds double, ptr %28, i64 %43
  %45 = load double, ptr %44, align 8, !tbaa !22
  %invariant.gep = getelementptr i8, ptr %28, i64 -8
  %46 = icmp sgt i32 %1, 0
  br i1 %46, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %7
  %47 = sitofp i32 %20 to double
  %48 = fsub double %47, %26
  %49 = fadd double %23, 2.000000e+00
  %wide.trip.count = zext nneg i32 %1 to i64
  br label %50

50:                                               ; preds = %.lr.ph, %177
  %indvars.iv = phi i64 [ 0, %.lr.ph ], [ %indvars.iv.next, %177 ]
  %.0231 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3, %177 ]
  %.0170230 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2172, %177 ]
  %.0173229 = phi i32 [ 0, %.lr.ph ], [ %.3176, %177 ]
  %.0179227 = phi double [ %41, %.lr.ph ], [ %.2181, %177 ]
  %.0182226 = phi double [ %41, %.lr.ph ], [ %.2184, %177 ]
  %.0185225 = phi double [ %45, %.lr.ph ], [ %.2187, %177 ]
  %.0188224 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2190, %177 ]
  %.sroa.0115.0223 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.0115.2, %177 ]
  %.sroa.5116.0222 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.5116.2, %177 ]
  %.sroa.7117.0221 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.7117.2, %177 ]
  %.sroa.0.0220 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.0.2, %177 ]
  %.0191219 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2193, %177 ]
  %.0194218 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3197, %177 ]
  %.sroa.5.0217 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.5.2, %177 ]
  %.sroa.7.0216 = phi double [ 0.000000e+00, %.lr.ph ], [ %.sroa.7.2, %177 ]
  %.0205213 = phi double [ 0.000000e+00, %.lr.ph ], [ %.2207, %177 ]
  %51 = getelementptr inbounds double, ptr %3, i64 %indvars.iv
  %52 = load double, ptr %51, align 8, !tbaa !22
  store double %52, ptr %8, align 8, !tbaa !22
  %53 = load i32, ptr %0, align 8, !tbaa !15
  %54 = icmp eq i32 %53, 1
  br i1 %54, label %55, label %60

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

60:                                               ; preds = %55, %50
  %61 = phi double [ %.pre, %55 ], [ %52, %50 ]
  %62 = fcmp olt double %61, %41
  %63 = fcmp ogt double %61, %45
  %or.cond = select i1 %62, i1 true, i1 %63
  br i1 %or.cond, label %64, label %66

64:                                               ; preds = %60
  %65 = getelementptr inbounds double, ptr %4, i64 %indvars.iv
  store double 0.000000e+00, ptr %65, align 8, !tbaa !22
  br label %177

66:                                               ; preds = %60
  %67 = fcmp ogt double %61, %.0182226
  %68 = fcmp olt double %61, %.0185225
  %or.cond211 = select i1 %67, i1 true, i1 %68
  br i1 %or.cond211, label %69, label %161

69:                                               ; preds = %66
  %70 = fcmp ogt double %61, %.0179227
  %or.cond212 = select i1 %70, i1 true, i1 %68
  br i1 %or.cond212, label %71, label %104

71:                                               ; preds = %69
  %72 = fsub double %61, %41
  %73 = fdiv double %72, %30
  %74 = fptosi double %73 to i32
  %75 = call i32 @llvm.smax.i32(i32 %74, i32 0)
  %76 = add nuw nsw i32 %75, 1
  %77 = call i32 @llvm.smin.i32(i32 %42, i32 %76)
  %78 = add nsw i32 %77, -1
  %79 = sext i32 %78 to i64
  %80 = getelementptr inbounds double, ptr %37, i64 %79
  %81 = load double, ptr %80, align 8, !tbaa !22
  %82 = getelementptr inbounds double, ptr %40, i64 %79
  %83 = load double, ptr %82, align 8, !tbaa !22
  %84 = getelementptr inbounds double, ptr %19, i64 %79
  %85 = load double, ptr %84, align 8, !tbaa !22
  %86 = getelementptr inbounds double, ptr %17, i64 %79
  %87 = load double, ptr %86, align 8, !tbaa !22
  %88 = fmul double %87, %87
  %89 = fmul double %83, -2.000000e+00
  %90 = fdiv double %23, %88
  %91 = fsub double %90, %26
  %92 = fadd double %91, %47
  %93 = fmul double %81, %92
  %94 = call double @llvm.fmuladd.f64(double %89, double %85, double %93)
  %95 = fmul double %85, -2.000000e+00
  %96 = fmul double %81, -2.000000e+00
  %97 = fmul double %23, %96
  %98 = fmul double %97, %85
  %99 = fdiv double %98, %88
  %100 = call double @llvm.fmuladd.f64(double %95, double %94, double %99)
  %101 = fdiv double %49, %88
  %102 = fadd double %48, %101
  %103 = call double @llvm.fmuladd.f64(double %83, double %102, double %100)
  br label %106

104:                                              ; preds = %69
  %105 = add nsw i32 %.0173229, 1
  %.pre233 = fmul double %.0188224, -2.000000e+00
  br label %106

106:                                              ; preds = %104, %71
  %.pre-phi = phi double [ %.pre233, %104 ], [ %89, %71 ]
  %.0204 = phi double [ %.0205213, %104 ], [ %103, %71 ]
  %.1195 = phi double [ %.0191219, %104 ], [ %94, %71 ]
  %.0177 = phi double [ %.0188224, %104 ], [ %83, %71 ]
  %.1174 = phi i32 [ %105, %104 ], [ %77, %71 ]
  %.1 = phi double [ %.0170230, %104 ], [ %81, %71 ]
  %107 = call i32 @llvm.smax.i32(i32 %.1174, i32 1)
  %108 = zext nneg i32 %107 to i64
  %gep = getelementptr double, ptr %invariant.gep, i64 %108
  %109 = load double, ptr %gep, align 8, !tbaa !22
  %110 = sext i32 %.1174 to i64
  %111 = getelementptr inbounds double, ptr %28, i64 %110
  %112 = load double, ptr %111, align 8, !tbaa !22
  %113 = add nsw i32 %.1174, 1
  %114 = call i32 @llvm.smin.i32(i32 %42, i32 %113)
  %115 = sext i32 %114 to i64
  %116 = getelementptr inbounds double, ptr %28, i64 %115
  %117 = load double, ptr %116, align 8, !tbaa !22
  %118 = getelementptr inbounds double, ptr %37, i64 %110
  %119 = load double, ptr %118, align 8, !tbaa !22
  %120 = getelementptr inbounds double, ptr %40, i64 %110
  %121 = load double, ptr %120, align 8, !tbaa !22
  %122 = getelementptr inbounds double, ptr %19, i64 %110
  %123 = load double, ptr %122, align 8, !tbaa !22
  %124 = getelementptr inbounds double, ptr %17, i64 %110
  %125 = load double, ptr %124, align 8, !tbaa !22
  %126 = fmul double %125, %125
  %127 = fmul double %121, -2.000000e+00
  %128 = fdiv double %23, %126
  %129 = fsub double %128, %26
  %130 = fadd double %129, %47
  %131 = fmul double %119, %130
  %132 = call double @llvm.fmuladd.f64(double %127, double %123, double %131)
  %133 = fmul double %123, -2.000000e+00
  %134 = fmul double %119, -2.000000e+00
  %135 = fmul double %23, %134
  %136 = fmul double %135, %123
  %137 = fdiv double %136, %126
  %138 = call double @llvm.fmuladd.f64(double %133, double %132, double %137)
  %139 = fdiv double %49, %126
  %140 = fadd double %48, %139
  %141 = call double @llvm.fmuladd.f64(double %121, double %140, double %138)
  %142 = fmul double %30, %.0177
  %143 = fmul double %30, %121
  %144 = fneg double %143
  %145 = call double @llvm.fmuladd.f64(double %.pre-phi, double %30, double %144)
  %146 = call double @llvm.fmuladd.f64(double %.1, double -3.000000e+00, double %145)
  %147 = call double @llvm.fmuladd.f64(double %119, double 3.000000e+00, double %146)
  %148 = call double @llvm.fmuladd.f64(double %.0177, double %30, double %143)
  %149 = call double @llvm.fmuladd.f64(double %.1, double 2.000000e+00, double %148)
  %150 = call double @llvm.fmuladd.f64(double %119, double -2.000000e+00, double %149)
  %151 = fmul double %30, %.0204
  %152 = fmul double %.0204, -2.000000e+00
  %153 = fmul double %30, %141
  %154 = fneg double %153
  %155 = call double @llvm.fmuladd.f64(double %152, double %30, double %154)
  %156 = call double @llvm.fmuladd.f64(double %.1195, double -3.000000e+00, double %155)
  %157 = call double @llvm.fmuladd.f64(double %132, double 3.000000e+00, double %156)
  %158 = call double @llvm.fmuladd.f64(double %.0204, double %30, double %153)
  %159 = call double @llvm.fmuladd.f64(double %.1195, double 2.000000e+00, double %158)
  %160 = call double @llvm.fmuladd.f64(double %132, double -2.000000e+00, double %159)
  br label %161

161:                                              ; preds = %66, %106
  %.1206 = phi double [ %141, %106 ], [ %.0205213, %66 ]
  %.sroa.7.1 = phi double [ %160, %106 ], [ %.sroa.7.0216, %66 ]
  %.sroa.5.1 = phi double [ %157, %106 ], [ %.sroa.5.0217, %66 ]
  %.2196 = phi double [ %.1195, %106 ], [ %.0194218, %66 ]
  %.1192 = phi double [ %132, %106 ], [ %.0191219, %66 ]
  %.sroa.0.1 = phi double [ %151, %106 ], [ %.sroa.0.0220, %66 ]
  %.sroa.7117.1 = phi double [ %150, %106 ], [ %.sroa.7117.0221, %66 ]
  %.sroa.5116.1 = phi double [ %147, %106 ], [ %.sroa.5116.0222, %66 ]
  %.sroa.0115.1 = phi double [ %142, %106 ], [ %.sroa.0115.0223, %66 ]
  %.1189 = phi double [ %121, %106 ], [ %.0188224, %66 ]
  %.1186 = phi double [ %109, %106 ], [ %.0185225, %66 ]
  %.1183 = phi double [ %112, %106 ], [ %.0182226, %66 ]
  %.1180 = phi double [ %117, %106 ], [ %.0179227, %66 ]
  %.2175 = phi i32 [ %.1174, %106 ], [ %.0173229, %66 ]
  %.1171 = phi double [ %119, %106 ], [ %.0170230, %66 ]
  %.2 = phi double [ %.1, %106 ], [ %.0231, %66 ]
  %162 = fsub double %61, %.1186
  %163 = fdiv double %162, %30
  %164 = fmul double %163, %163
  %165 = fmul double %163, %164
  %166 = call double @llvm.fmuladd.f64(double %.sroa.0115.1, double %163, double %.2)
  %167 = call double @llvm.fmuladd.f64(double %.sroa.5116.1, double %164, double %166)
  %168 = call double @llvm.fmuladd.f64(double %.sroa.7117.1, double %165, double %167)
  %169 = load i32, ptr %9, align 4, !tbaa !12
  %170 = sitofp i32 %169 to double
  %171 = fmul double %168, %170
  %172 = getelementptr inbounds double, ptr %4, i64 %indvars.iv
  store double %171, ptr %172, align 8, !tbaa !22
  %173 = call double @llvm.fmuladd.f64(double %.sroa.0.1, double %163, double %.2196)
  %174 = call double @llvm.fmuladd.f64(double %.sroa.5.1, double %164, double %173)
  %175 = call double @llvm.fmuladd.f64(double %.sroa.7.1, double %165, double %174)
  %176 = fmul double %175, %170
  br label %177

177:                                              ; preds = %161, %64
  %.sink = phi double [ %176, %161 ], [ 0.000000e+00, %64 ]
  %.2207 = phi double [ %.1206, %161 ], [ %.0205213, %64 ]
  %.sroa.7.2 = phi double [ %.sroa.7.1, %161 ], [ %.sroa.7.0216, %64 ]
  %.sroa.5.2 = phi double [ %.sroa.5.1, %161 ], [ %.sroa.5.0217, %64 ]
  %.3197 = phi double [ %.2196, %161 ], [ %.0194218, %64 ]
  %.2193 = phi double [ %.1192, %161 ], [ %.0191219, %64 ]
  %.sroa.0.2 = phi double [ %.sroa.0.1, %161 ], [ %.sroa.0.0220, %64 ]
  %.sroa.7117.2 = phi double [ %.sroa.7117.1, %161 ], [ %.sroa.7117.0221, %64 ]
  %.sroa.5116.2 = phi double [ %.sroa.5116.1, %161 ], [ %.sroa.5116.0222, %64 ]
  %.sroa.0115.2 = phi double [ %.sroa.0115.1, %161 ], [ %.sroa.0115.0223, %64 ]
  %.2190 = phi double [ %.1189, %161 ], [ %.0188224, %64 ]
  %.2187 = phi double [ %.1186, %161 ], [ %.0185225, %64 ]
  %.2184 = phi double [ %.1183, %161 ], [ %.0182226, %64 ]
  %.2181 = phi double [ %.1180, %161 ], [ %.0179227, %64 ]
  %.3176 = phi i32 [ %.2175, %161 ], [ %.0173229, %64 ]
  %.2172 = phi double [ %.1171, %161 ], [ %.0170230, %64 ]
  %.3 = phi double [ %.2, %161 ], [ %.0231, %64 ]
  %178 = getelementptr inbounds double, ptr %5, i64 %indvars.iv
  store double %.sink, ptr %178, 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 %50, !llvm.loop !23

._crit_edge:                                      ; preds = %177, %7
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %10) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %9) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %8) #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