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

118:                                              ; preds = %428, %7
  %119 = load i32, ptr %46, align 4, !tbaa !9
  %120 = load i32, ptr %9, align 4, !tbaa !9
  %121 = icmp slt i32 %119, %120
  br i1 %121, label %122, label %431

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

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

140:                                              ; preds = %132, %122
  %141 = load double, ptr %19, align 8, !tbaa !14
  %142 = load double, ptr %40, align 8, !tbaa !14
  %143 = fcmp olt double %141, %142
  br i1 %143, label %148, label %144

144:                                              ; preds = %140
  %145 = load double, ptr %19, align 8, !tbaa !14
  %146 = load double, ptr %41, align 8, !tbaa !14
  %147 = fcmp ogt double %145, %146
  br i1 %147, label %148, label %157

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

157:                                              ; preds = %144
  %158 = load double, ptr %19, align 8, !tbaa !14
  %159 = load double, ptr %44, align 8, !tbaa !14
  %160 = fcmp ogt double %158, %159
  br i1 %160, label %165, label %161

161:                                              ; preds = %157
  %162 = load double, ptr %19, align 8, !tbaa !14
  %163 = load double, ptr %43, align 8, !tbaa !14
  %164 = fcmp olt double %162, %163
  br i1 %164, label %165, label %379

165:                                              ; preds = %161, %157
  %166 = load double, ptr %19, align 8, !tbaa !14
  %167 = load double, ptr %45, align 8, !tbaa !14
  %168 = fcmp ogt double %166, %167
  br i1 %168, label %173, label %169

169:                                              ; preds = %165
  %170 = load double, ptr %19, align 8, !tbaa !14
  %171 = load double, ptr %43, align 8, !tbaa !14
  %172 = fcmp olt double %170, %171
  br i1 %172, label %173, label %240

173:                                              ; preds = %169, %165
  %174 = load double, ptr %19, align 8, !tbaa !14
  %175 = load double, ptr %40, align 8, !tbaa !14
  %176 = fsub double %174, %175
  %177 = load double, ptr %42, align 8, !tbaa !14
  %178 = fdiv double %176, %177
  %179 = fptosi double %178 to i32
  %180 = add nsw i32 %179, 1
  store i32 %180, ptr %48, align 4, !tbaa !9
  %181 = load i32, ptr %48, align 4, !tbaa !9
  %182 = icmp slt i32 1, %181
  br i1 %182, label %183, label %185

183:                                              ; preds = %173
  %184 = load i32, ptr %48, align 4, !tbaa !9
  br label %186

185:                                              ; preds = %173
  br label %186

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

192:                                              ; preds = %186
  %193 = load i32, ptr %47, align 4, !tbaa !9
  %194 = sub nsw i32 %193, 1
  br label %197

195:                                              ; preds = %186
  %196 = load i32, ptr %48, align 4, !tbaa !9
  br label %197

197:                                              ; preds = %195, %192
  %198 = phi i32 [ %194, %192 ], [ %196, %195 ]
  store i32 %198, ptr %48, align 4, !tbaa !9
  %199 = load ptr, ptr %49, align 8, !tbaa !5
  %200 = load i32, ptr %48, align 4, !tbaa !9
  %201 = sub nsw i32 %200, 1
  %202 = sext i32 %201 to i64
  %203 = getelementptr inbounds double, ptr %199, i64 %202
  %204 = load double, ptr %203, align 8, !tbaa !14
  store double %204, ptr %16, align 8, !tbaa !14
  %205 = load ptr, ptr %50, align 8, !tbaa !5
  %206 = load i32, ptr %48, align 4, !tbaa !9
  %207 = sub nsw i32 %206, 1
  %208 = sext i32 %207 to i64
  %209 = getelementptr inbounds double, ptr %205, i64 %208
  %210 = load double, ptr %209, align 8, !tbaa !14
  store double %210, ptr %28, align 8, !tbaa !14
  %211 = load ptr, ptr %34, align 8, !tbaa !5
  %212 = load i32, ptr %48, align 4, !tbaa !9
  %213 = sub nsw i32 %212, 1
  %214 = sext i32 %213 to i64
  %215 = getelementptr inbounds double, ptr %211, i64 %214
  %216 = load double, ptr %215, align 8, !tbaa !14
  store double %216, ptr %25, align 8, !tbaa !14
  %217 = load ptr, ptr %33, align 8, !tbaa !5
  %218 = load i32, ptr %48, align 4, !tbaa !9
  %219 = sub nsw i32 %218, 1
  %220 = sext i32 %219 to i64
  %221 = getelementptr inbounds double, ptr %217, i64 %220
  %222 = load double, ptr %221, align 8, !tbaa !14
  store double %222, ptr %26, align 8, !tbaa !14
  %223 = load double, ptr %26, align 8, !tbaa !14
  %224 = load double, ptr %26, align 8, !tbaa !14
  %225 = fmul double %223, %224
  store double %225, ptr %27, align 8, !tbaa !14
  %226 = load double, ptr %28, align 8, !tbaa !14
  %227 = fmul double -2.000000e+00, %226
  %228 = load double, ptr %25, align 8, !tbaa !14
  %229 = load double, ptr %16, align 8, !tbaa !14
  %230 = load double, ptr %36, align 8, !tbaa !14
  %231 = load double, ptr %27, align 8, !tbaa !14
  %232 = fdiv double %230, %231
  %233 = load double, ptr %38, align 8, !tbaa !14
  %234 = fsub double %232, %233
  %235 = load i32, ptr %35, align 4, !tbaa !9
  %236 = sitofp i32 %235 to double
  %237 = fadd double %234, %236
  %238 = fmul double %229, %237
  %239 = call double @llvm.fmuladd.f64(double %227, double %228, double %238)
  store double %239, ptr %23, align 8, !tbaa !14
  br label %248

240:                                              ; preds = %169
  %241 = load i32, ptr %48, align 4, !tbaa !9
  %242 = add nsw i32 %241, 1
  store i32 %242, ptr %48, align 4, !tbaa !9
  %243 = load double, ptr %17, align 8, !tbaa !14
  store double %243, ptr %16, align 8, !tbaa !14
  %244 = load double, ptr %18, align 8, !tbaa !14
  store double %244, ptr %28, align 8, !tbaa !14
  %245 = load double, ptr %29, align 8, !tbaa !14
  store double %245, ptr %23, align 8, !tbaa !14
  %246 = load double, ptr %31, align 8, !tbaa !14
  store double %246, ptr %26, align 8, !tbaa !14
  %247 = load double, ptr %30, align 8, !tbaa !14
  store double %247, ptr %25, align 8, !tbaa !14
  br label %248

248:                                              ; preds = %240, %197
  %249 = load ptr, ptr %39, align 8, !tbaa !5
  %250 = load i32, ptr %48, align 4, !tbaa !9
  %251 = sub nsw i32 %250, 1
  %252 = icmp slt i32 0, %251
  br i1 %252, label %253, label %256

253:                                              ; preds = %248
  %254 = load i32, ptr %48, align 4, !tbaa !9
  %255 = sub nsw i32 %254, 1
  br label %257

256:                                              ; preds = %248
  br label %257

257:                                              ; preds = %256, %253
  %258 = phi i32 [ %255, %253 ], [ 0, %256 ]
  %259 = sext i32 %258 to i64
  %260 = getelementptr inbounds double, ptr %249, i64 %259
  %261 = load double, ptr %260, align 8, !tbaa !14
  store double %261, ptr %43, align 8, !tbaa !14
  %262 = load ptr, ptr %39, align 8, !tbaa !5
  %263 = load i32, ptr %48, align 4, !tbaa !9
  %264 = sext i32 %263 to i64
  %265 = getelementptr inbounds double, ptr %262, i64 %264
  %266 = load double, ptr %265, align 8, !tbaa !14
  store double %266, ptr %44, align 8, !tbaa !14
  %267 = load ptr, ptr %39, align 8, !tbaa !5
  %268 = load i32, ptr %47, align 4, !tbaa !9
  %269 = sub nsw i32 %268, 1
  %270 = load i32, ptr %48, align 4, !tbaa !9
  %271 = add nsw i32 %270, 1
  %272 = icmp slt i32 %269, %271
  br i1 %272, label %273, label %276

273:                                              ; preds = %257
  %274 = load i32, ptr %47, align 4, !tbaa !9
  %275 = sub nsw i32 %274, 1
  br label %279

276:                                              ; preds = %257
  %277 = load i32, ptr %48, align 4, !tbaa !9
  %278 = add nsw i32 %277, 1
  br label %279

279:                                              ; preds = %276, %273
  %280 = phi i32 [ %275, %273 ], [ %278, %276 ]
  %281 = sext i32 %280 to i64
  %282 = getelementptr inbounds double, ptr %267, i64 %281
  %283 = load double, ptr %282, align 8, !tbaa !14
  store double %283, ptr %45, align 8, !tbaa !14
  %284 = load ptr, ptr %49, align 8, !tbaa !5
  %285 = load i32, ptr %48, align 4, !tbaa !9
  %286 = sext i32 %285 to i64
  %287 = getelementptr inbounds double, ptr %284, i64 %286
  %288 = load double, ptr %287, align 8, !tbaa !14
  store double %288, ptr %17, align 8, !tbaa !14
  %289 = load ptr, ptr %50, align 8, !tbaa !5
  %290 = load i32, ptr %48, align 4, !tbaa !9
  %291 = sext i32 %290 to i64
  %292 = getelementptr inbounds double, ptr %289, i64 %291
  %293 = load double, ptr %292, align 8, !tbaa !14
  store double %293, ptr %18, align 8, !tbaa !14
  %294 = load ptr, ptr %34, align 8, !tbaa !5
  %295 = load i32, ptr %48, align 4, !tbaa !9
  %296 = sext i32 %295 to i64
  %297 = getelementptr inbounds double, ptr %294, i64 %296
  %298 = load double, ptr %297, align 8, !tbaa !14
  store double %298, ptr %30, align 8, !tbaa !14
  %299 = load ptr, ptr %33, align 8, !tbaa !5
  %300 = load i32, ptr %48, align 4, !tbaa !9
  %301 = sext i32 %300 to i64
  %302 = getelementptr inbounds double, ptr %299, i64 %301
  %303 = load double, ptr %302, align 8, !tbaa !14
  store double %303, ptr %31, align 8, !tbaa !14
  %304 = load double, ptr %31, align 8, !tbaa !14
  %305 = load double, ptr %31, align 8, !tbaa !14
  %306 = fmul double %304, %305
  store double %306, ptr %32, align 8, !tbaa !14
  %307 = load double, ptr %18, align 8, !tbaa !14
  %308 = fmul double -2.000000e+00, %307
  %309 = load double, ptr %30, align 8, !tbaa !14
  %310 = load double, ptr %17, align 8, !tbaa !14
  %311 = load double, ptr %36, align 8, !tbaa !14
  %312 = load double, ptr %32, align 8, !tbaa !14
  %313 = fdiv double %311, %312
  %314 = load double, ptr %38, align 8, !tbaa !14
  %315 = fsub double %313, %314
  %316 = load i32, ptr %35, align 4, !tbaa !9
  %317 = sitofp i32 %316 to double
  %318 = fadd double %315, %317
  %319 = fmul double %310, %318
  %320 = call double @llvm.fmuladd.f64(double %308, double %309, double %319)
  store double %320, ptr %29, align 8, !tbaa !14
  %321 = load double, ptr %30, align 8, !tbaa !14
  %322 = fmul double -2.000000e+00, %321
  %323 = load double, ptr %29, align 8, !tbaa !14
  %324 = load double, ptr %17, align 8, !tbaa !14
  %325 = fmul double 2.000000e+00, %324
  %326 = load double, ptr %36, align 8, !tbaa !14
  %327 = fmul double %325, %326
  %328 = load double, ptr %30, align 8, !tbaa !14
  %329 = fmul double %327, %328
  %330 = load double, ptr %32, align 8, !tbaa !14
  %331 = fdiv double %329, %330
  %332 = fneg double %331
  %333 = call double @llvm.fmuladd.f64(double %322, double %323, double %332)
  %334 = load double, ptr %18, align 8, !tbaa !14
  %335 = load i32, ptr %35, align 4, !tbaa !9
  %336 = sitofp i32 %335 to double
  %337 = load double, ptr %38, align 8, !tbaa !14
  %338 = fsub double %336, %337
  %339 = load double, ptr %36, align 8, !tbaa !14
  %340 = fadd double 2.000000e+00, %339
  %341 = load double, ptr %32, align 8, !tbaa !14
  %342 = fdiv double %340, %341
  %343 = fadd double %338, %342
  %344 = call double @llvm.fmuladd.f64(double %334, double %343, double %333)
  store double %344, ptr %24, align 8, !tbaa !14
  %345 = load double, ptr %18, align 8, !tbaa !14
  %346 = fneg double %345
  %347 = load double, ptr %42, align 8, !tbaa !14
  %348 = load double, ptr %16, align 8, !tbaa !14
  %349 = fmul double 2.000000e+00, %348
  %350 = fneg double %349
  %351 = call double @llvm.fmuladd.f64(double %346, double %347, double %350)
  %352 = load double, ptr %17, align 8, !tbaa !14
  %353 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %352, double %351)
  %354 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  store double %353, ptr %354, align 16, !tbaa !14
  %355 = load double, ptr %18, align 8, !tbaa !14
  %356 = load double, ptr %42, align 8, !tbaa !14
  %357 = load double, ptr %16, align 8, !tbaa !14
  %358 = call double @llvm.fmuladd.f64(double %355, double %356, double %357)
  %359 = load double, ptr %17, align 8, !tbaa !14
  %360 = fsub double %358, %359
  %361 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 1
  store double %360, ptr %361, align 8, !tbaa !14
  %362 = load double, ptr %24, align 8, !tbaa !14
  %363 = fneg double %362
  %364 = load double, ptr %42, align 8, !tbaa !14
  %365 = load double, ptr %23, align 8, !tbaa !14
  %366 = fmul double 2.000000e+00, %365
  %367 = fneg double %366
  %368 = call double @llvm.fmuladd.f64(double %363, double %364, double %367)
  %369 = load double, ptr %29, align 8, !tbaa !14
  %370 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %369, double %368)
  %371 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 0
  store double %370, ptr %371, align 16, !tbaa !14
  %372 = load double, ptr %24, align 8, !tbaa !14
  %373 = load double, ptr %42, align 8, !tbaa !14
  %374 = load double, ptr %23, align 8, !tbaa !14
  %375 = call double @llvm.fmuladd.f64(double %372, double %373, double %374)
  %376 = load double, ptr %29, align 8, !tbaa !14
  %377 = fsub double %375, %376
  %378 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 1
  store double %377, ptr %378, align 8, !tbaa !14
  br label %379

379:                                              ; preds = %279, %161
  %380 = load double, ptr %19, align 8, !tbaa !14
  %381 = load double, ptr %43, align 8, !tbaa !14
  %382 = fsub double %380, %381
  %383 = load double, ptr %42, align 8, !tbaa !14
  %384 = fdiv double %382, %383
  %385 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 0
  store double %384, ptr %385, align 16, !tbaa !14
  %386 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 0
  %387 = load double, ptr %386, align 16, !tbaa !14
  %388 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 0
  %389 = load double, ptr %388, align 16, !tbaa !14
  %390 = fmul double %387, %389
  %391 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 1
  store double %390, ptr %391, align 8, !tbaa !14
  %392 = load double, ptr %16, align 8, !tbaa !14
  %393 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 0
  %394 = load double, ptr %393, align 16, !tbaa !14
  %395 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 0
  %396 = load double, ptr %395, align 16, !tbaa !14
  %397 = call double @llvm.fmuladd.f64(double %394, double %396, double %392)
  %398 = getelementptr inbounds [3 x double], ptr %21, i64 0, i64 1
  %399 = load double, ptr %398, align 8, !tbaa !14
  %400 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 1
  %401 = load double, ptr %400, align 8, !tbaa !14
  %402 = call double @llvm.fmuladd.f64(double %399, double %401, double %397)
  %403 = load i32, ptr %51, align 4, !tbaa !9
  %404 = sitofp i32 %403 to double
  %405 = fmul double %402, %404
  %406 = load ptr, ptr %12, align 8, !tbaa !5
  %407 = load i32, ptr %46, align 4, !tbaa !9
  %408 = sext i32 %407 to i64
  %409 = getelementptr inbounds double, ptr %406, i64 %408
  store double %405, ptr %409, align 8, !tbaa !14
  %410 = load double, ptr %23, align 8, !tbaa !14
  %411 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 0
  %412 = load double, ptr %411, align 16, !tbaa !14
  %413 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 0
  %414 = load double, ptr %413, align 16, !tbaa !14
  %415 = call double @llvm.fmuladd.f64(double %412, double %414, double %410)
  %416 = getelementptr inbounds [3 x double], ptr %22, i64 0, i64 1
  %417 = load double, ptr %416, align 8, !tbaa !14
  %418 = getelementptr inbounds [2 x double], ptr %20, i64 0, i64 1
  %419 = load double, ptr %418, align 8, !tbaa !14
  %420 = call double @llvm.fmuladd.f64(double %417, double %419, double %415)
  %421 = load i32, ptr %51, align 4, !tbaa !9
  %422 = sitofp i32 %421 to double
  %423 = fmul double %420, %422
  %424 = load ptr, ptr %13, align 8, !tbaa !5
  %425 = load i32, ptr %46, align 4, !tbaa !9
  %426 = sext i32 %425 to i64
  %427 = getelementptr inbounds double, ptr %424, i64 %426
  store double %423, ptr %427, align 8, !tbaa !14
  br label %428

428:                                              ; preds = %379, %148
  %429 = load i32, ptr %46, align 4, !tbaa !9
  %430 = add nsw i32 %429, 1
  store i32 %430, ptr %46, align 4, !tbaa !9
  br label %118, !llvm.loop !24

431:                                              ; preds = %118
  call void @llvm.lifetime.end.p0(i64 4, ptr %52) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %51) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %50) #4
  call void @llvm.lifetime.end.p0(i64 8, 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 4, 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 8, ptr %36) #4
  call void @llvm.lifetime.end.p0(i64 4, 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 8, 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 16, 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_Hermite3_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
  %50 = insertelement <2 x double> poison, double %30, i64 0
  %51 = shufflevector <2 x double> %50, <2 x double> poison, <2 x i32> zeroinitializer
  br label %52

52:                                               ; preds = %.lr.ph, %173
  %indvars.iv = phi i64 [ 0, %.lr.ph ], [ %indvars.iv.next, %173 ]
  %.0196 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3, %173 ]
  %.0146194 = phi i32 [ 0, %.lr.ph ], [ %.3149, %173 ]
  %.0154191 = phi double [ %41, %.lr.ph ], [ %.2156, %173 ]
  %.0157190 = phi double [ %41, %.lr.ph ], [ %.2159, %173 ]
  %.0160189 = phi double [ %45, %.lr.ph ], [ %.2162, %173 ]
  %.0163184 = phi double [ 0.000000e+00, %.lr.ph ], [ %.3166, %173 ]
  %53 = phi <2 x double> [ zeroinitializer, %.lr.ph ], [ %174, %173 ]
  %54 = phi <2 x double> [ zeroinitializer, %.lr.ph ], [ %175, %173 ]
  %55 = phi <2 x double> [ zeroinitializer, %.lr.ph ], [ %176, %173 ]
  %56 = getelementptr inbounds double, ptr %3, i64 %indvars.iv
  %57 = load double, ptr %56, align 8, !tbaa !22
  store double %57, ptr %8, align 8, !tbaa !22
  %58 = load i32, ptr %0, align 8, !tbaa !15
  %59 = icmp eq i32 %58, 1
  br i1 %59, label %60, label %65

60:                                               ; preds = %52
  %61 = load double, ptr %24, align 8, !tbaa !16
  %62 = fadd double %61, 2.000000e-01
  %63 = fptosi double %62 to i32
  %64 = call i32 @ClosedModY(i32 noundef %15, i32 noundef %63, ptr noundef nonnull %8, ptr noundef nonnull %9, ptr noundef nonnull %10)
  %.pre = load double, ptr %8, align 8, !tbaa !22
  br label %65

65:                                               ; preds = %60, %52
  %66 = phi double [ %.pre, %60 ], [ %57, %52 ]
  %67 = fcmp olt double %66, %41
  %68 = fcmp ogt double %66, %45
  %or.cond = select i1 %67, i1 true, i1 %68
  br i1 %or.cond, label %69, label %71

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

71:                                               ; preds = %65
  %72 = fcmp ogt double %66, %.0157190
  %73 = fcmp olt double %66, %.0160189
  %or.cond179 = select i1 %72, i1 true, i1 %73
  br i1 %or.cond179, label %74, label %153

74:                                               ; preds = %71
  %75 = fcmp ogt double %66, %.0154191
  %or.cond180 = select i1 %75, i1 true, i1 %73
  br i1 %or.cond180, label %76, label %102

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

102:                                              ; preds = %74
  %103 = add nsw i32 %.0146194, 1
  br label %104

104:                                              ; preds = %102, %76
  %.1147 = phi i32 [ %82, %76 ], [ %103, %102 ]
  %105 = phi <2 x double> [ %101, %76 ], [ %54, %102 ]
  %106 = call i32 @llvm.smax.i32(i32 %.1147, i32 1)
  %107 = zext nneg i32 %106 to i64
  %gep = getelementptr double, ptr %invariant.gep, i64 %107
  %108 = load double, ptr %gep, align 8, !tbaa !22
  %109 = sext i32 %.1147 to i64
  %110 = getelementptr inbounds double, ptr %28, i64 %109
  %111 = load double, ptr %110, align 8, !tbaa !22
  %112 = add nsw i32 %.1147, 1
  %113 = call i32 @llvm.smin.i32(i32 %42, i32 %112)
  %114 = sext i32 %113 to i64
  %115 = getelementptr inbounds double, ptr %28, i64 %114
  %116 = load double, ptr %115, align 8, !tbaa !22
  %117 = getelementptr inbounds double, ptr %37, i64 %109
  %118 = load double, ptr %117, align 8, !tbaa !22
  %119 = getelementptr inbounds double, ptr %40, i64 %109
  %120 = load double, ptr %119, align 8, !tbaa !22
  %121 = getelementptr inbounds double, ptr %19, i64 %109
  %122 = load double, ptr %121, align 8, !tbaa !22
  %123 = getelementptr inbounds double, ptr %17, i64 %109
  %124 = load double, ptr %123, align 8, !tbaa !22
  %125 = fmul double %124, %124
  %126 = fmul double %120, -2.000000e+00
  %127 = fdiv double %23, %125
  %128 = fsub double %127, %26
  %129 = fadd double %128, %47
  %130 = fmul double %118, %129
  %131 = call double @llvm.fmuladd.f64(double %126, double %122, double %130)
  %132 = fmul double %122, -2.000000e+00
  %133 = fmul double %118, -2.000000e+00
  %134 = fmul double %23, %133
  %135 = fmul double %134, %122
  %136 = fdiv double %135, %125
  %137 = call double @llvm.fmuladd.f64(double %132, double %131, double %136)
  %138 = fdiv double %49, %125
  %139 = fadd double %48, %138
  %140 = call double @llvm.fmuladd.f64(double %120, double %139, double %137)
  %141 = fmul <2 x double> %105, <double -2.000000e+00, double -2.000000e+00>
  %142 = insertelement <2 x double> poison, double %140, i64 0
  %143 = insertelement <2 x double> %142, double %120, i64 1
  %144 = fneg <2 x double> %143
  %145 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %144, <2 x double> %51, <2 x double> %141)
  %146 = insertelement <2 x double> poison, double %131, i64 0
  %147 = insertelement <2 x double> %146, double %118, i64 1
  %148 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %147, <2 x double> <double 2.000000e+00, double 2.000000e+00>, <2 x double> %145)
  %149 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %143, <2 x double> %51, <2 x double> %105)
  %150 = fsub <2 x double> %149, %147
  %151 = extractelement <2 x double> %105, i64 0
  %152 = extractelement <2 x double> %105, i64 1
  br label %153

153:                                              ; preds = %71, %104
  %.2165 = phi double [ %151, %104 ], [ %.0163184, %71 ]
  %.1161 = phi double [ %108, %104 ], [ %.0160189, %71 ]
  %.1158 = phi double [ %111, %104 ], [ %.0157190, %71 ]
  %.1155 = phi double [ %116, %104 ], [ %.0154191, %71 ]
  %.2148 = phi i32 [ %.1147, %104 ], [ %.0146194, %71 ]
  %.2 = phi double [ %152, %104 ], [ %.0196, %71 ]
  %154 = phi <2 x double> [ %148, %104 ], [ %53, %71 ]
  %155 = phi <2 x double> [ %147, %104 ], [ %54, %71 ]
  %156 = phi <2 x double> [ %150, %104 ], [ %55, %71 ]
  %157 = fsub double %66, %.1161
  %158 = fdiv double %157, %30
  %159 = fmul double %158, %158
  %160 = extractelement <2 x double> %154, i64 1
  %161 = call double @llvm.fmuladd.f64(double %160, double %158, double %.2)
  %162 = extractelement <2 x double> %156, i64 1
  %163 = call double @llvm.fmuladd.f64(double %162, double %159, double %161)
  %164 = load i32, ptr %9, align 4, !tbaa !12
  %165 = sitofp i32 %164 to double
  %166 = fmul double %163, %165
  %167 = getelementptr inbounds double, ptr %4, i64 %indvars.iv
  store double %166, ptr %167, align 8, !tbaa !22
  %168 = extractelement <2 x double> %154, i64 0
  %169 = call double @llvm.fmuladd.f64(double %168, double %158, double %.2165)
  %170 = extractelement <2 x double> %156, i64 0
  %171 = call double @llvm.fmuladd.f64(double %170, double %159, double %169)
  %172 = fmul double %171, %165
  br label %173

173:                                              ; preds = %153, %69
  %.sink = phi double [ %172, %153 ], [ 0.000000e+00, %69 ]
  %.3166 = phi double [ %.2165, %153 ], [ %.0163184, %69 ]
  %.2162 = phi double [ %.1161, %153 ], [ %.0160189, %69 ]
  %.2159 = phi double [ %.1158, %153 ], [ %.0157190, %69 ]
  %.2156 = phi double [ %.1155, %153 ], [ %.0154191, %69 ]
  %.3149 = phi i32 [ %.2148, %153 ], [ %.0146194, %69 ]
  %.3 = phi double [ %.2, %153 ], [ %.0196, %69 ]
  %174 = phi <2 x double> [ %154, %153 ], [ %53, %69 ]
  %175 = phi <2 x double> [ %155, %153 ], [ %54, %69 ]
  %176 = phi <2 x double> [ %156, %153 ], [ %55, %69 ]
  %177 = getelementptr inbounds double, ptr %5, i64 %indvars.iv
  store double %.sink, ptr %177, 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 %52, !llvm.loop !23

._crit_edge:                                      ; preds = %173, %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

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare <2 x double> @llvm.fmuladd.v2f64(<2 x double>, <2 x double>, <2 x double>) #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,SLPVectorizerPass,SROAPass,SimplifyCFGPass