; LLVM IR from VeloxChem/VeloxChem :: src/dft_func (DensityGridCubic)

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

%class.CDensityGridCubic = type { i32, i32, %class.CDensityGridData2D }
%class.CDensityGridData2D = type <{ %"class.std::vector", %"class.std::vector.0", %"class.std::vector.0", %"class.std::vector.0", i32, [4 x i8] }>
%"class.std::vector" = type { %"struct.std::_Vector_base" }
%"struct.std::_Vector_base" = type { %"struct.std::_Vector_base<double, std::allocator<double>>::_Vector_impl" }
%"struct.std::_Vector_base<double, std::allocator<double>>::_Vector_impl" = type { %"struct.std::_Vector_base<double, std::allocator<double>>::_Vector_impl_data" }
%"struct.std::_Vector_base<double, std::allocator<double>>::_Vector_impl_data" = type { ptr, ptr, ptr }
%"class.std::vector.0" = type { %"struct.std::_Vector_base.1" }
%"struct.std::_Vector_base.1" = type { %"struct.std::_Vector_base<int, std::allocator<int>>::_Vector_impl" }
%"struct.std::_Vector_base<int, std::allocator<int>>::_Vector_impl" = type { %"struct.std::_Vector_base<int, std::allocator<int>>::_Vector_impl_data" }
%"struct.std::_Vector_base<int, std::allocator<int>>::_Vector_impl_data" = type { 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: mustprogress uwtable
declare dso_local noundef i32 @_ZNK17CDensityGridCubic21getNumberOfGridPointsEv(ptr noundef nonnull align 8 dereferenceable(112)) #1 align 2

; Function Attrs: mustprogress uwtable
declare dso_local noundef ptr @_ZN17CDensityGridCubic2piEi(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) #1 align 2

; Function Attrs: mustprogress uwtable
declare dso_local noundef ptr @_ZN17CDensityGridCubic3gamEi(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) #1 align 2

; Function Attrs: mustprogress uwtable
declare dso_local noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) #1 align 2

; Function Attrs: mustprogress uwtable
define dso_local void @_ZN17CDensityGridCubic21DensityProdForLDA_CRFERK12CDensityGridS2_i(ptr noundef nonnull align 8 dereferenceable(112) %0, ptr noundef nonnull align 8 dereferenceable(112) %1, ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %3) #1 align 2 {
  %5 = alloca ptr, align 8
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca i32, align 4
  %9 = alloca i32, align 4
  %10 = alloca i32, align 4
  %11 = alloca i32, align 4
  %12 = alloca ptr, align 8
  %13 = alloca ptr, align 8
  %14 = alloca ptr, align 8
  %15 = alloca ptr, align 8
  %16 = alloca ptr, align 8
  %17 = alloca ptr, align 8
  %18 = alloca ptr, align 8
  %19 = alloca ptr, align 8
  %20 = alloca ptr, align 8
  %21 = alloca ptr, align 8
  %22 = alloca ptr, align 8
  %23 = alloca ptr, align 8
  %24 = alloca ptr, align 8
  %25 = alloca ptr, align 8
  %26 = alloca ptr, align 8
  %27 = alloca ptr, align 8
  %28 = alloca ptr, align 8
  %29 = alloca ptr, align 8
  %30 = alloca ptr, align 8
  %31 = alloca ptr, align 8
  %32 = alloca ptr, align 8
  %33 = alloca ptr, align 8
  %34 = alloca i32, align 4
  store ptr %0, ptr %5, align 8, !tbaa !5
  store ptr %1, ptr %6, align 8, !tbaa !5
  store ptr %2, ptr %7, align 8, !tbaa !5
  store i32 %3, ptr %8, align 4, !tbaa !9
  %35 = load ptr, ptr %5, align 8
  %36 = getelementptr inbounds %class.CDensityGridCubic, ptr %35, i32 0, i32 0
  %37 = load i32, ptr %36, align 8, !tbaa !11
  %38 = icmp ne i32 %37, 0
  br i1 %38, label %39, label %40

39:                                               ; preds = %4
  br label %612

40:                                               ; preds = %4
  call void @llvm.lifetime.start.p0(i64 4, ptr %9) #5
  %41 = call noundef i32 @_ZNK17CDensityGridCubic21getNumberOfGridPointsEv(ptr noundef nonnull align 8 dereferenceable(112) %35)
  store i32 %41, ptr %9, align 4, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 4, ptr %10) #5
  store i32 0, ptr %10, align 4, !tbaa !9
  br label %42

42:                                               ; preds = %608, %40
  %43 = load i32, ptr %10, align 4, !tbaa !9
  %44 = load i32, ptr %8, align 4, !tbaa !9
  %45 = sdiv i32 %44, 8
  %46 = icmp slt i32 %43, %45
  br i1 %46, label %48, label %47

47:                                               ; preds = %42
  store i32 2, ptr %11, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %10) #5
  br label %611

48:                                               ; preds = %42
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #5
  %49 = load i32, ptr %10, align 4, !tbaa !9
  %50 = mul nsw i32 6, %49
  %51 = call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %50)
  store ptr %51, ptr %12, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %13) #5
  %52 = load i32, ptr %10, align 4, !tbaa !9
  %53 = mul nsw i32 6, %52
  %54 = add nsw i32 %53, 1
  %55 = call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %54)
  store ptr %55, ptr %13, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #5
  %56 = load i32, ptr %10, align 4, !tbaa !9
  %57 = mul nsw i32 6, %56
  %58 = add nsw i32 %57, 2
  %59 = call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %58)
  store ptr %59, ptr %14, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #5
  %60 = load i32, ptr %10, align 4, !tbaa !9
  %61 = mul nsw i32 6, %60
  %62 = add nsw i32 %61, 3
  %63 = call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %62)
  store ptr %63, ptr %15, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #5
  %64 = load i32, ptr %10, align 4, !tbaa !9
  %65 = mul nsw i32 6, %64
  %66 = add nsw i32 %65, 4
  %67 = call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %66)
  store ptr %67, ptr %16, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #5
  %68 = load i32, ptr %10, align 4, !tbaa !9
  %69 = mul nsw i32 6, %68
  %70 = add nsw i32 %69, 5
  %71 = call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %70)
  store ptr %71, ptr %17, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #5
  %72 = load i32, ptr %10, align 4, !tbaa !9
  %73 = mul nsw i32 2, %72
  %74 = call noundef ptr @_ZN17CDensityGridCubic3gamEi(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %73)
  store ptr %74, ptr %18, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %19) #5
  %75 = load i32, ptr %10, align 4, !tbaa !9
  %76 = mul nsw i32 2, %75
  %77 = add nsw i32 %76, 1
  %78 = call noundef ptr @_ZN17CDensityGridCubic3gamEi(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %77)
  store ptr %78, ptr %19, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #5
  %79 = load i32, ptr %10, align 4, !tbaa !9
  %80 = mul nsw i32 2, %79
  %81 = call noundef ptr @_ZN17CDensityGridCubic2piEi(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %80)
  store ptr %81, ptr %20, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #5
  %82 = load i32, ptr %10, align 4, !tbaa !9
  %83 = mul nsw i32 2, %82
  %84 = add nsw i32 %83, 1
  %85 = call noundef ptr @_ZN17CDensityGridCubic2piEi(ptr noundef nonnull align 8 dereferenceable(112) %35, i32 noundef %84)
  store ptr %85, ptr %21, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %22) #5
  %86 = load ptr, ptr %6, align 8, !tbaa !5
  %87 = load i32, ptr %10, align 4, !tbaa !9
  %88 = mul nsw i32 6, %87
  %89 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %86, i32 noundef %88)
  store ptr %89, ptr %22, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #5
  %90 = load ptr, ptr %6, align 8, !tbaa !5
  %91 = load i32, ptr %10, align 4, !tbaa !9
  %92 = mul nsw i32 6, %91
  %93 = add nsw i32 %92, 1
  %94 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %90, i32 noundef %93)
  store ptr %94, ptr %23, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #5
  %95 = load ptr, ptr %6, align 8, !tbaa !5
  %96 = load i32, ptr %10, align 4, !tbaa !9
  %97 = mul nsw i32 6, %96
  %98 = add nsw i32 %97, 2
  %99 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %95, i32 noundef %98)
  store ptr %99, ptr %24, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %25) #5
  %100 = load ptr, ptr %6, align 8, !tbaa !5
  %101 = load i32, ptr %10, align 4, !tbaa !9
  %102 = mul nsw i32 6, %101
  %103 = add nsw i32 %102, 3
  %104 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %100, i32 noundef %103)
  store ptr %104, ptr %25, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %26) #5
  %105 = load ptr, ptr %6, align 8, !tbaa !5
  %106 = load i32, ptr %10, align 4, !tbaa !9
  %107 = mul nsw i32 6, %106
  %108 = add nsw i32 %107, 4
  %109 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %105, i32 noundef %108)
  store ptr %109, ptr %26, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %27) #5
  %110 = load ptr, ptr %6, align 8, !tbaa !5
  %111 = load i32, ptr %10, align 4, !tbaa !9
  %112 = mul nsw i32 6, %111
  %113 = add nsw i32 %112, 5
  %114 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %110, i32 noundef %113)
  store ptr %114, ptr %27, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %28) #5
  %115 = load ptr, ptr %7, align 8, !tbaa !5
  %116 = load i32, ptr %10, align 4, !tbaa !9
  %117 = mul nsw i32 6, %116
  %118 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %115, i32 noundef %117)
  store ptr %118, ptr %28, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %29) #5
  %119 = load ptr, ptr %7, align 8, !tbaa !5
  %120 = load i32, ptr %10, align 4, !tbaa !9
  %121 = mul nsw i32 6, %120
  %122 = add nsw i32 %121, 1
  %123 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %119, i32 noundef %122)
  store ptr %123, ptr %29, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %30) #5
  %124 = load ptr, ptr %7, align 8, !tbaa !5
  %125 = load i32, ptr %10, align 4, !tbaa !9
  %126 = mul nsw i32 6, %125
  %127 = add nsw i32 %126, 2
  %128 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %124, i32 noundef %127)
  store ptr %128, ptr %30, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %31) #5
  %129 = load ptr, ptr %7, align 8, !tbaa !5
  %130 = load i32, ptr %10, align 4, !tbaa !9
  %131 = mul nsw i32 6, %130
  %132 = add nsw i32 %131, 3
  %133 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %129, i32 noundef %132)
  store ptr %133, ptr %31, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %32) #5
  %134 = load ptr, ptr %7, align 8, !tbaa !5
  %135 = load i32, ptr %10, align 4, !tbaa !9
  %136 = mul nsw i32 6, %135
  %137 = add nsw i32 %136, 4
  %138 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %134, i32 noundef %137)
  store ptr %138, ptr %32, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %33) #5
  %139 = load ptr, ptr %7, align 8, !tbaa !5
  %140 = load i32, ptr %10, align 4, !tbaa !9
  %141 = mul nsw i32 6, %140
  %142 = add nsw i32 %141, 5
  %143 = call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %139, i32 noundef %142)
  store ptr %143, ptr %33, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %34) #5
  store i32 0, ptr %34, align 4, !tbaa !9
  br label %144

144:                                              ; preds = %604, %48
  %145 = load i32, ptr %34, align 4, !tbaa !9
  %146 = load i32, ptr %9, align 4, !tbaa !9
  %147 = icmp slt i32 %145, %146
  br i1 %147, label %149, label %148

148:                                              ; preds = %144
  store i32 5, ptr %11, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %34) #5
  br label %607

149:                                              ; preds = %144
  %150 = load ptr, ptr %22, align 8, !tbaa !5
  %151 = load i32, ptr %34, align 4, !tbaa !9
  %152 = sext i32 %151 to i64
  %153 = getelementptr inbounds double, ptr %150, i64 %152
  %154 = load double, ptr %153, align 8, !tbaa !23
  %155 = load ptr, ptr %23, align 8, !tbaa !5
  %156 = load i32, ptr %34, align 4, !tbaa !9
  %157 = sext i32 %156 to i64
  %158 = getelementptr inbounds double, ptr %155, i64 %157
  %159 = load double, ptr %158, align 8, !tbaa !23
  %160 = load ptr, ptr %24, align 8, !tbaa !5
  %161 = load i32, ptr %34, align 4, !tbaa !9
  %162 = sext i32 %161 to i64
  %163 = getelementptr inbounds double, ptr %160, i64 %162
  %164 = load double, ptr %163, align 8, !tbaa !23
  %165 = load ptr, ptr %25, align 8, !tbaa !5
  %166 = load i32, ptr %34, align 4, !tbaa !9
  %167 = sext i32 %166 to i64
  %168 = getelementptr inbounds double, ptr %165, i64 %167
  %169 = load double, ptr %168, align 8, !tbaa !23
  %170 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %154, double noundef %159, double noundef %164, double noundef %169)
  %171 = fmul double 2.000000e+00, %170
  %172 = load ptr, ptr %12, align 8, !tbaa !5
  %173 = load i32, ptr %34, align 4, !tbaa !9
  %174 = sext i32 %173 to i64
  %175 = getelementptr inbounds double, ptr %172, i64 %174
  store double %171, ptr %175, align 8, !tbaa !23
  %176 = load ptr, ptr %22, align 8, !tbaa !5
  %177 = load i32, ptr %34, align 4, !tbaa !9
  %178 = sext i32 %177 to i64
  %179 = getelementptr inbounds double, ptr %176, i64 %178
  %180 = load double, ptr %179, align 8, !tbaa !23
  %181 = load ptr, ptr %23, align 8, !tbaa !5
  %182 = load i32, ptr %34, align 4, !tbaa !9
  %183 = sext i32 %182 to i64
  %184 = getelementptr inbounds double, ptr %181, i64 %183
  %185 = load double, ptr %184, align 8, !tbaa !23
  %186 = load ptr, ptr %24, align 8, !tbaa !5
  %187 = load i32, ptr %34, align 4, !tbaa !9
  %188 = sext i32 %187 to i64
  %189 = getelementptr inbounds double, ptr %186, i64 %188
  %190 = load double, ptr %189, align 8, !tbaa !23
  %191 = load ptr, ptr %25, align 8, !tbaa !5
  %192 = load i32, ptr %34, align 4, !tbaa !9
  %193 = sext i32 %192 to i64
  %194 = getelementptr inbounds double, ptr %191, i64 %193
  %195 = load double, ptr %194, align 8, !tbaa !23
  %196 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %180, double noundef %185, double noundef %190, double noundef %195)
  %197 = fmul double 2.000000e+00, %196
  %198 = load ptr, ptr %13, align 8, !tbaa !5
  %199 = load i32, ptr %34, align 4, !tbaa !9
  %200 = sext i32 %199 to i64
  %201 = getelementptr inbounds double, ptr %198, i64 %200
  store double %197, ptr %201, align 8, !tbaa !23
  %202 = load ptr, ptr %22, align 8, !tbaa !5
  %203 = load i32, ptr %34, align 4, !tbaa !9
  %204 = sext i32 %203 to i64
  %205 = getelementptr inbounds double, ptr %202, i64 %204
  %206 = load double, ptr %205, align 8, !tbaa !23
  %207 = load ptr, ptr %23, align 8, !tbaa !5
  %208 = load i32, ptr %34, align 4, !tbaa !9
  %209 = sext i32 %208 to i64
  %210 = getelementptr inbounds double, ptr %207, i64 %209
  %211 = load double, ptr %210, align 8, !tbaa !23
  %212 = load ptr, ptr %26, align 8, !tbaa !5
  %213 = load i32, ptr %34, align 4, !tbaa !9
  %214 = sext i32 %213 to i64
  %215 = getelementptr inbounds double, ptr %212, i64 %214
  %216 = load double, ptr %215, align 8, !tbaa !23
  %217 = load ptr, ptr %27, align 8, !tbaa !5
  %218 = load i32, ptr %34, align 4, !tbaa !9
  %219 = sext i32 %218 to i64
  %220 = getelementptr inbounds double, ptr %217, i64 %219
  %221 = load double, ptr %220, align 8, !tbaa !23
  %222 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %206, double noundef %211, double noundef %216, double noundef %221)
  %223 = fmul double 2.000000e+00, %222
  %224 = load ptr, ptr %14, align 8, !tbaa !5
  %225 = load i32, ptr %34, align 4, !tbaa !9
  %226 = sext i32 %225 to i64
  %227 = getelementptr inbounds double, ptr %224, i64 %226
  store double %223, ptr %227, align 8, !tbaa !23
  %228 = load ptr, ptr %22, align 8, !tbaa !5
  %229 = load i32, ptr %34, align 4, !tbaa !9
  %230 = sext i32 %229 to i64
  %231 = getelementptr inbounds double, ptr %228, i64 %230
  %232 = load double, ptr %231, align 8, !tbaa !23
  %233 = load ptr, ptr %23, align 8, !tbaa !5
  %234 = load i32, ptr %34, align 4, !tbaa !9
  %235 = sext i32 %234 to i64
  %236 = getelementptr inbounds double, ptr %233, i64 %235
  %237 = load double, ptr %236, align 8, !tbaa !23
  %238 = load ptr, ptr %26, align 8, !tbaa !5
  %239 = load i32, ptr %34, align 4, !tbaa !9
  %240 = sext i32 %239 to i64
  %241 = getelementptr inbounds double, ptr %238, i64 %240
  %242 = load double, ptr %241, align 8, !tbaa !23
  %243 = load ptr, ptr %27, align 8, !tbaa !5
  %244 = load i32, ptr %34, align 4, !tbaa !9
  %245 = sext i32 %244 to i64
  %246 = getelementptr inbounds double, ptr %243, i64 %245
  %247 = load double, ptr %246, align 8, !tbaa !23
  %248 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %232, double noundef %237, double noundef %242, double noundef %247)
  %249 = fmul double 2.000000e+00, %248
  %250 = load ptr, ptr %15, align 8, !tbaa !5
  %251 = load i32, ptr %34, align 4, !tbaa !9
  %252 = sext i32 %251 to i64
  %253 = getelementptr inbounds double, ptr %250, i64 %252
  store double %249, ptr %253, align 8, !tbaa !23
  %254 = load ptr, ptr %24, align 8, !tbaa !5
  %255 = load i32, ptr %34, align 4, !tbaa !9
  %256 = sext i32 %255 to i64
  %257 = getelementptr inbounds double, ptr %254, i64 %256
  %258 = load double, ptr %257, align 8, !tbaa !23
  %259 = load ptr, ptr %25, align 8, !tbaa !5
  %260 = load i32, ptr %34, align 4, !tbaa !9
  %261 = sext i32 %260 to i64
  %262 = getelementptr inbounds double, ptr %259, i64 %261
  %263 = load double, ptr %262, align 8, !tbaa !23
  %264 = load ptr, ptr %26, align 8, !tbaa !5
  %265 = load i32, ptr %34, align 4, !tbaa !9
  %266 = sext i32 %265 to i64
  %267 = getelementptr inbounds double, ptr %264, i64 %266
  %268 = load double, ptr %267, align 8, !tbaa !23
  %269 = load ptr, ptr %27, align 8, !tbaa !5
  %270 = load i32, ptr %34, align 4, !tbaa !9
  %271 = sext i32 %270 to i64
  %272 = getelementptr inbounds double, ptr %269, i64 %271
  %273 = load double, ptr %272, align 8, !tbaa !23
  %274 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %258, double noundef %263, double noundef %268, double noundef %273)
  %275 = fmul double 2.000000e+00, %274
  %276 = load ptr, ptr %16, align 8, !tbaa !5
  %277 = load i32, ptr %34, align 4, !tbaa !9
  %278 = sext i32 %277 to i64
  %279 = getelementptr inbounds double, ptr %276, i64 %278
  store double %275, ptr %279, align 8, !tbaa !23
  %280 = load ptr, ptr %24, align 8, !tbaa !5
  %281 = load i32, ptr %34, align 4, !tbaa !9
  %282 = sext i32 %281 to i64
  %283 = getelementptr inbounds double, ptr %280, i64 %282
  %284 = load double, ptr %283, align 8, !tbaa !23
  %285 = load ptr, ptr %25, align 8, !tbaa !5
  %286 = load i32, ptr %34, align 4, !tbaa !9
  %287 = sext i32 %286 to i64
  %288 = getelementptr inbounds double, ptr %285, i64 %287
  %289 = load double, ptr %288, align 8, !tbaa !23
  %290 = load ptr, ptr %26, align 8, !tbaa !5
  %291 = load i32, ptr %34, align 4, !tbaa !9
  %292 = sext i32 %291 to i64
  %293 = getelementptr inbounds double, ptr %290, i64 %292
  %294 = load double, ptr %293, align 8, !tbaa !23
  %295 = load ptr, ptr %27, align 8, !tbaa !5
  %296 = load i32, ptr %34, align 4, !tbaa !9
  %297 = sext i32 %296 to i64
  %298 = getelementptr inbounds double, ptr %295, i64 %297
  %299 = load double, ptr %298, align 8, !tbaa !23
  %300 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %284, double noundef %289, double noundef %294, double noundef %299)
  %301 = fmul double 2.000000e+00, %300
  %302 = load ptr, ptr %17, align 8, !tbaa !5
  %303 = load i32, ptr %34, align 4, !tbaa !9
  %304 = sext i32 %303 to i64
  %305 = getelementptr inbounds double, ptr %302, i64 %304
  store double %301, ptr %305, align 8, !tbaa !23
  %306 = load ptr, ptr %26, align 8, !tbaa !5
  %307 = load i32, ptr %34, align 4, !tbaa !9
  %308 = sext i32 %307 to i64
  %309 = getelementptr inbounds double, ptr %306, i64 %308
  %310 = load double, ptr %309, align 8, !tbaa !23
  %311 = load ptr, ptr %27, align 8, !tbaa !5
  %312 = load i32, ptr %34, align 4, !tbaa !9
  %313 = sext i32 %312 to i64
  %314 = getelementptr inbounds double, ptr %311, i64 %313
  %315 = load double, ptr %314, align 8, !tbaa !23
  %316 = load ptr, ptr %28, align 8, !tbaa !5
  %317 = load i32, ptr %34, align 4, !tbaa !9
  %318 = sext i32 %317 to i64
  %319 = getelementptr inbounds double, ptr %316, i64 %318
  %320 = load double, ptr %319, align 8, !tbaa !23
  %321 = load ptr, ptr %29, align 8, !tbaa !5
  %322 = load i32, ptr %34, align 4, !tbaa !9
  %323 = sext i32 %322 to i64
  %324 = getelementptr inbounds double, ptr %321, i64 %323
  %325 = load double, ptr %324, align 8, !tbaa !23
  %326 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %310, double noundef %315, double noundef %320, double noundef %325)
  %327 = fmul double 3.000000e+00, %326
  %328 = load ptr, ptr %18, align 8, !tbaa !5
  %329 = load i32, ptr %34, align 4, !tbaa !9
  %330 = sext i32 %329 to i64
  %331 = getelementptr inbounds double, ptr %328, i64 %330
  store double %327, ptr %331, align 8, !tbaa !23
  %332 = load ptr, ptr %26, align 8, !tbaa !5
  %333 = load i32, ptr %34, align 4, !tbaa !9
  %334 = sext i32 %333 to i64
  %335 = getelementptr inbounds double, ptr %332, i64 %334
  %336 = load double, ptr %335, align 8, !tbaa !23
  %337 = load ptr, ptr %27, align 8, !tbaa !5
  %338 = load i32, ptr %34, align 4, !tbaa !9
  %339 = sext i32 %338 to i64
  %340 = getelementptr inbounds double, ptr %337, i64 %339
  %341 = load double, ptr %340, align 8, !tbaa !23
  %342 = load ptr, ptr %28, align 8, !tbaa !5
  %343 = load i32, ptr %34, align 4, !tbaa !9
  %344 = sext i32 %343 to i64
  %345 = getelementptr inbounds double, ptr %342, i64 %344
  %346 = load double, ptr %345, align 8, !tbaa !23
  %347 = load ptr, ptr %29, align 8, !tbaa !5
  %348 = load i32, ptr %34, align 4, !tbaa !9
  %349 = sext i32 %348 to i64
  %350 = getelementptr inbounds double, ptr %347, i64 %349
  %351 = load double, ptr %350, align 8, !tbaa !23
  %352 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %336, double noundef %341, double noundef %346, double noundef %351)
  %353 = fmul double 3.000000e+00, %352
  %354 = load ptr, ptr %19, align 8, !tbaa !5
  %355 = load i32, ptr %34, align 4, !tbaa !9
  %356 = sext i32 %355 to i64
  %357 = getelementptr inbounds double, ptr %354, i64 %356
  store double %353, ptr %357, align 8, !tbaa !23
  %358 = load ptr, ptr %22, align 8, !tbaa !5
  %359 = load i32, ptr %34, align 4, !tbaa !9
  %360 = sext i32 %359 to i64
  %361 = getelementptr inbounds double, ptr %358, i64 %360
  %362 = load double, ptr %361, align 8, !tbaa !23
  %363 = load ptr, ptr %23, align 8, !tbaa !5
  %364 = load i32, ptr %34, align 4, !tbaa !9
  %365 = sext i32 %364 to i64
  %366 = getelementptr inbounds double, ptr %363, i64 %365
  %367 = load double, ptr %366, align 8, !tbaa !23
  %368 = load ptr, ptr %32, align 8, !tbaa !5
  %369 = load i32, ptr %34, align 4, !tbaa !9
  %370 = sext i32 %369 to i64
  %371 = getelementptr inbounds double, ptr %368, i64 %370
  %372 = load double, ptr %371, align 8, !tbaa !23
  %373 = load ptr, ptr %33, align 8, !tbaa !5
  %374 = load i32, ptr %34, align 4, !tbaa !9
  %375 = sext i32 %374 to i64
  %376 = getelementptr inbounds double, ptr %373, i64 %375
  %377 = load double, ptr %376, align 8, !tbaa !23
  %378 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %362, double noundef %367, double noundef %372, double noundef %377)
  %379 = load ptr, ptr %18, align 8, !tbaa !5
  %380 = load i32, ptr %34, align 4, !tbaa !9
  %381 = sext i32 %380 to i64
  %382 = getelementptr inbounds double, ptr %379, i64 %381
  %383 = load double, ptr %382, align 8, !tbaa !23
  %384 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %378, double %383)
  store double %384, ptr %382, align 8, !tbaa !23
  %385 = load ptr, ptr %22, align 8, !tbaa !5
  %386 = load i32, ptr %34, align 4, !tbaa !9
  %387 = sext i32 %386 to i64
  %388 = getelementptr inbounds double, ptr %385, i64 %387
  %389 = load double, ptr %388, align 8, !tbaa !23
  %390 = load ptr, ptr %23, align 8, !tbaa !5
  %391 = load i32, ptr %34, align 4, !tbaa !9
  %392 = sext i32 %391 to i64
  %393 = getelementptr inbounds double, ptr %390, i64 %392
  %394 = load double, ptr %393, align 8, !tbaa !23
  %395 = load ptr, ptr %32, align 8, !tbaa !5
  %396 = load i32, ptr %34, align 4, !tbaa !9
  %397 = sext i32 %396 to i64
  %398 = getelementptr inbounds double, ptr %395, i64 %397
  %399 = load double, ptr %398, align 8, !tbaa !23
  %400 = load ptr, ptr %33, align 8, !tbaa !5
  %401 = load i32, ptr %34, align 4, !tbaa !9
  %402 = sext i32 %401 to i64
  %403 = getelementptr inbounds double, ptr %400, i64 %402
  %404 = load double, ptr %403, align 8, !tbaa !23
  %405 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %389, double noundef %394, double noundef %399, double noundef %404)
  %406 = load ptr, ptr %19, align 8, !tbaa !5
  %407 = load i32, ptr %34, align 4, !tbaa !9
  %408 = sext i32 %407 to i64
  %409 = getelementptr inbounds double, ptr %406, i64 %408
  %410 = load double, ptr %409, align 8, !tbaa !23
  %411 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %405, double %410)
  store double %411, ptr %409, align 8, !tbaa !23
  %412 = load ptr, ptr %24, align 8, !tbaa !5
  %413 = load i32, ptr %34, align 4, !tbaa !9
  %414 = sext i32 %413 to i64
  %415 = getelementptr inbounds double, ptr %412, i64 %414
  %416 = load double, ptr %415, align 8, !tbaa !23
  %417 = load ptr, ptr %25, align 8, !tbaa !5
  %418 = load i32, ptr %34, align 4, !tbaa !9
  %419 = sext i32 %418 to i64
  %420 = getelementptr inbounds double, ptr %417, i64 %419
  %421 = load double, ptr %420, align 8, !tbaa !23
  %422 = load ptr, ptr %30, align 8, !tbaa !5
  %423 = load i32, ptr %34, align 4, !tbaa !9
  %424 = sext i32 %423 to i64
  %425 = getelementptr inbounds double, ptr %422, i64 %424
  %426 = load double, ptr %425, align 8, !tbaa !23
  %427 = load ptr, ptr %31, align 8, !tbaa !5
  %428 = load i32, ptr %34, align 4, !tbaa !9
  %429 = sext i32 %428 to i64
  %430 = getelementptr inbounds double, ptr %427, i64 %429
  %431 = load double, ptr %430, align 8, !tbaa !23
  %432 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %416, double noundef %421, double noundef %426, double noundef %431)
  %433 = load ptr, ptr %18, align 8, !tbaa !5
  %434 = load i32, ptr %34, align 4, !tbaa !9
  %435 = sext i32 %434 to i64
  %436 = getelementptr inbounds double, ptr %433, i64 %435
  %437 = load double, ptr %436, align 8, !tbaa !23
  %438 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %432, double %437)
  store double %438, ptr %436, align 8, !tbaa !23
  %439 = load ptr, ptr %24, align 8, !tbaa !5
  %440 = load i32, ptr %34, align 4, !tbaa !9
  %441 = sext i32 %440 to i64
  %442 = getelementptr inbounds double, ptr %439, i64 %441
  %443 = load double, ptr %442, align 8, !tbaa !23
  %444 = load ptr, ptr %25, align 8, !tbaa !5
  %445 = load i32, ptr %34, align 4, !tbaa !9
  %446 = sext i32 %445 to i64
  %447 = getelementptr inbounds double, ptr %444, i64 %446
  %448 = load double, ptr %447, align 8, !tbaa !23
  %449 = load ptr, ptr %30, align 8, !tbaa !5
  %450 = load i32, ptr %34, align 4, !tbaa !9
  %451 = sext i32 %450 to i64
  %452 = getelementptr inbounds double, ptr %449, i64 %451
  %453 = load double, ptr %452, align 8, !tbaa !23
  %454 = load ptr, ptr %31, align 8, !tbaa !5
  %455 = load i32, ptr %34, align 4, !tbaa !9
  %456 = sext i32 %455 to i64
  %457 = getelementptr inbounds double, ptr %454, i64 %456
  %458 = load double, ptr %457, align 8, !tbaa !23
  %459 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %443, double noundef %448, double noundef %453, double noundef %458)
  %460 = load ptr, ptr %19, align 8, !tbaa !5
  %461 = load i32, ptr %34, align 4, !tbaa !9
  %462 = sext i32 %461 to i64
  %463 = getelementptr inbounds double, ptr %460, i64 %462
  %464 = load double, ptr %463, align 8, !tbaa !23
  %465 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %459, double %464)
  store double %465, ptr %463, align 8, !tbaa !23
  %466 = load ptr, ptr %22, align 8, !tbaa !5
  %467 = load i32, ptr %34, align 4, !tbaa !9
  %468 = sext i32 %467 to i64
  %469 = getelementptr inbounds double, ptr %466, i64 %468
  %470 = load double, ptr %469, align 8, !tbaa !23
  %471 = load ptr, ptr %23, align 8, !tbaa !5
  %472 = load i32, ptr %34, align 4, !tbaa !9
  %473 = sext i32 %472 to i64
  %474 = getelementptr inbounds double, ptr %471, i64 %473
  %475 = load double, ptr %474, align 8, !tbaa !23
  %476 = load ptr, ptr %16, align 8, !tbaa !5
  %477 = load i32, ptr %34, align 4, !tbaa !9
  %478 = sext i32 %477 to i64
  %479 = getelementptr inbounds double, ptr %476, i64 %478
  %480 = load double, ptr %479, align 8, !tbaa !23
  %481 = load ptr, ptr %17, align 8, !tbaa !5
  %482 = load i32, ptr %34, align 4, !tbaa !9
  %483 = sext i32 %482 to i64
  %484 = getelementptr inbounds double, ptr %481, i64 %483
  %485 = load double, ptr %484, align 8, !tbaa !23
  %486 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %470, double noundef %475, double noundef %480, double noundef %485)
  %487 = load ptr, ptr %24, align 8, !tbaa !5
  %488 = load i32, ptr %34, align 4, !tbaa !9
  %489 = sext i32 %488 to i64
  %490 = getelementptr inbounds double, ptr %487, i64 %489
  %491 = load double, ptr %490, align 8, !tbaa !23
  %492 = load ptr, ptr %25, align 8, !tbaa !5
  %493 = load i32, ptr %34, align 4, !tbaa !9
  %494 = sext i32 %493 to i64
  %495 = getelementptr inbounds double, ptr %492, i64 %494
  %496 = load double, ptr %495, align 8, !tbaa !23
  %497 = load ptr, ptr %14, align 8, !tbaa !5
  %498 = load i32, ptr %34, align 4, !tbaa !9
  %499 = sext i32 %498 to i64
  %500 = getelementptr inbounds double, ptr %497, i64 %499
  %501 = load double, ptr %500, align 8, !tbaa !23
  %502 = load ptr, ptr %15, align 8, !tbaa !5
  %503 = load i32, ptr %34, align 4, !tbaa !9
  %504 = sext i32 %503 to i64
  %505 = getelementptr inbounds double, ptr %502, i64 %504
  %506 = load double, ptr %505, align 8, !tbaa !23
  %507 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %491, double noundef %496, double noundef %501, double noundef %506)
  %508 = fadd double %486, %507
  %509 = load ptr, ptr %26, align 8, !tbaa !5
  %510 = load i32, ptr %34, align 4, !tbaa !9
  %511 = sext i32 %510 to i64
  %512 = getelementptr inbounds double, ptr %509, i64 %511
  %513 = load double, ptr %512, align 8, !tbaa !23
  %514 = load ptr, ptr %27, align 8, !tbaa !5
  %515 = load i32, ptr %34, align 4, !tbaa !9
  %516 = sext i32 %515 to i64
  %517 = getelementptr inbounds double, ptr %514, i64 %516
  %518 = load double, ptr %517, align 8, !tbaa !23
  %519 = load ptr, ptr %12, align 8, !tbaa !5
  %520 = load i32, ptr %34, align 4, !tbaa !9
  %521 = sext i32 %520 to i64
  %522 = getelementptr inbounds double, ptr %519, i64 %521
  %523 = load double, ptr %522, align 8, !tbaa !23
  %524 = load ptr, ptr %13, align 8, !tbaa !5
  %525 = load i32, ptr %34, align 4, !tbaa !9
  %526 = sext i32 %525 to i64
  %527 = getelementptr inbounds double, ptr %524, i64 %526
  %528 = load double, ptr %527, align 8, !tbaa !23
  %529 = call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %513, double noundef %518, double noundef %523, double noundef %528)
  %530 = fadd double %508, %529
  %531 = load ptr, ptr %20, align 8, !tbaa !5
  %532 = load i32, ptr %34, align 4, !tbaa !9
  %533 = sext i32 %532 to i64
  %534 = getelementptr inbounds double, ptr %531, i64 %533
  store double %530, ptr %534, align 8, !tbaa !23
  %535 = load ptr, ptr %22, align 8, !tbaa !5
  %536 = load i32, ptr %34, align 4, !tbaa !9
  %537 = sext i32 %536 to i64
  %538 = getelementptr inbounds double, ptr %535, i64 %537
  %539 = load double, ptr %538, align 8, !tbaa !23
  %540 = load ptr, ptr %23, align 8, !tbaa !5
  %541 = load i32, ptr %34, align 4, !tbaa !9
  %542 = sext i32 %541 to i64
  %543 = getelementptr inbounds double, ptr %540, i64 %542
  %544 = load double, ptr %543, align 8, !tbaa !23
  %545 = load ptr, ptr %16, align 8, !tbaa !5
  %546 = load i32, ptr %34, align 4, !tbaa !9
  %547 = sext i32 %546 to i64
  %548 = getelementptr inbounds double, ptr %545, i64 %547
  %549 = load double, ptr %548, align 8, !tbaa !23
  %550 = load ptr, ptr %17, align 8, !tbaa !5
  %551 = load i32, ptr %34, align 4, !tbaa !9
  %552 = sext i32 %551 to i64
  %553 = getelementptr inbounds double, ptr %550, i64 %552
  %554 = load double, ptr %553, align 8, !tbaa !23
  %555 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %539, double noundef %544, double noundef %549, double noundef %554)
  %556 = load ptr, ptr %24, align 8, !tbaa !5
  %557 = load i32, ptr %34, align 4, !tbaa !9
  %558 = sext i32 %557 to i64
  %559 = getelementptr inbounds double, ptr %556, i64 %558
  %560 = load double, ptr %559, align 8, !tbaa !23
  %561 = load ptr, ptr %25, align 8, !tbaa !5
  %562 = load i32, ptr %34, align 4, !tbaa !9
  %563 = sext i32 %562 to i64
  %564 = getelementptr inbounds double, ptr %561, i64 %563
  %565 = load double, ptr %564, align 8, !tbaa !23
  %566 = load ptr, ptr %14, align 8, !tbaa !5
  %567 = load i32, ptr %34, align 4, !tbaa !9
  %568 = sext i32 %567 to i64
  %569 = getelementptr inbounds double, ptr %566, i64 %568
  %570 = load double, ptr %569, align 8, !tbaa !23
  %571 = load ptr, ptr %15, align 8, !tbaa !5
  %572 = load i32, ptr %34, align 4, !tbaa !9
  %573 = sext i32 %572 to i64
  %574 = getelementptr inbounds double, ptr %571, i64 %573
  %575 = load double, ptr %574, align 8, !tbaa !23
  %576 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %560, double noundef %565, double noundef %570, double noundef %575)
  %577 = fadd double %555, %576
  %578 = load ptr, ptr %26, align 8, !tbaa !5
  %579 = load i32, ptr %34, align 4, !tbaa !9
  %580 = sext i32 %579 to i64
  %581 = getelementptr inbounds double, ptr %578, i64 %580
  %582 = load double, ptr %581, align 8, !tbaa !23
  %583 = load ptr, ptr %27, align 8, !tbaa !5
  %584 = load i32, ptr %34, align 4, !tbaa !9
  %585 = sext i32 %584 to i64
  %586 = getelementptr inbounds double, ptr %583, i64 %585
  %587 = load double, ptr %586, align 8, !tbaa !23
  %588 = load ptr, ptr %12, align 8, !tbaa !5
  %589 = load i32, ptr %34, align 4, !tbaa !9
  %590 = sext i32 %589 to i64
  %591 = getelementptr inbounds double, ptr %588, i64 %590
  %592 = load double, ptr %591, align 8, !tbaa !23
  %593 = load ptr, ptr %13, align 8, !tbaa !5
  %594 = load i32, ptr %34, align 4, !tbaa !9
  %595 = sext i32 %594 to i64
  %596 = getelementptr inbounds double, ptr %593, i64 %595
  %597 = load double, ptr %596, align 8, !tbaa !23
  %598 = call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %35, double noundef %582, double noundef %587, double noundef %592, double noundef %597)
  %599 = fadd double %577, %598
  %600 = load ptr, ptr %21, align 8, !tbaa !5
  %601 = load i32, ptr %34, align 4, !tbaa !9
  %602 = sext i32 %601 to i64
  %603 = getelementptr inbounds double, ptr %600, i64 %602
  store double %599, ptr %603, align 8, !tbaa !23
  br label %604

604:                                              ; preds = %149
  %605 = load i32, ptr %34, align 4, !tbaa !9
  %606 = add nsw i32 %605, 1
  store i32 %606, ptr %34, align 4, !tbaa !9
  br label %144, !llvm.loop !25

607:                                              ; preds = %148
  call void @llvm.lifetime.end.p0(i64 8, ptr %33) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %32) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %31) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %30) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %29) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %26) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %22) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %19) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %18) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %17) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %13) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #5
  br label %608

608:                                              ; preds = %607
  %609 = load i32, ptr %10, align 4, !tbaa !9
  %610 = add nsw i32 %609, 1
  store i32 %610, ptr %10, align 4, !tbaa !9
  br label %42, !llvm.loop !27

611:                                              ; preds = %47
  call void @llvm.lifetime.end.p0(i64 4, ptr %9) #5
  br label %612

612:                                              ; preds = %611, %39
  ret void
}

declare noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) #2

; Function Attrs: inlinehint mustprogress nounwind uwtable
declare dso_local noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112), double noundef, double noundef, double noundef, double noundef) #3 align 2

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

; Function Attrs: inlinehint mustprogress nounwind uwtable
declare dso_local noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112), double noundef, double noundef, double noundef, double noundef) #3 align 2

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { mustprogress uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { "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 = { inlinehint mustprogress nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #4 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #5 = { nounwind }

!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++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"int", !7, i64 0}
!11 = !{!12, !13, i64 0}
!12 = !{!"_ZTS17CDensityGridCubic", !13, i64 0, !10, i64 4, !14, i64 8}
!13 = !{!"_ZTS7dengrid", !7, i64 0}
!14 = !{!"_ZTS18CDensityGridData2D", !15, i64 0, !19, i64 24, !19, i64 48, !19, i64 72, !10, i64 96}
!15 = !{!"_ZTSSt6vectorIdSaIdEE", !16, i64 0}
!16 = !{!"_ZTSSt12_Vector_baseIdSaIdEE", !17, i64 0}
!17 = !{!"_ZTSNSt12_Vector_baseIdSaIdEE12_Vector_implE", !18, i64 0}
!18 = !{!"_ZTSNSt12_Vector_baseIdSaIdEE17_Vector_impl_dataE", !6, i64 0, !6, i64 8, !6, i64 16}
!19 = !{!"_ZTSSt6vectorIiSaIiEE", !20, i64 0}
!20 = !{!"_ZTSSt12_Vector_baseIiSaIiEE", !21, i64 0}
!21 = !{!"_ZTSNSt12_Vector_baseIiSaIiEE12_Vector_implE", !22, i64 0}
!22 = !{!"_ZTSNSt12_Vector_baseIiSaIiEE17_Vector_impl_dataE", !6, i64 0, !6, i64 8, !6, i64 16}
!23 = !{!24, !24, i64 0}
!24 = !{!"double", !7, i64 0}
!25 = distinct !{!25, !26}
!26 = !{!"llvm.loop.mustprogress"}
!27 = distinct !{!27, !26}


; 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 uwtable
declare dso_local noundef i32 @_ZNK17CDensityGridCubic21getNumberOfGridPointsEv(ptr noundef nonnull align 8 dereferenceable(112)) local_unnamed_addr #0 align 2

; Function Attrs: mustprogress uwtable
declare dso_local noundef ptr @_ZN17CDensityGridCubic2piEi(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) local_unnamed_addr #0 align 2

; Function Attrs: mustprogress uwtable
declare dso_local noundef ptr @_ZN17CDensityGridCubic3gamEi(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) local_unnamed_addr #0 align 2

; Function Attrs: mustprogress uwtable
declare dso_local noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) local_unnamed_addr #0 align 2

; Function Attrs: mustprogress uwtable
define dso_local void @_ZN17CDensityGridCubic21DensityProdForLDA_CRFERK12CDensityGridS2_i(ptr noundef nonnull align 8 dereferenceable(112) %0, ptr noundef nonnull align 8 dereferenceable(112) %1, ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %3) local_unnamed_addr #0 align 2 {
  %5 = load i32, ptr %0, align 8, !tbaa !5
  %.not = icmp eq i32 %5, 0
  br i1 %.not, label %6, label %.loopexit

6:                                                ; preds = %4
  %7 = tail call noundef i32 @_ZNK17CDensityGridCubic21getNumberOfGridPointsEv(ptr noundef nonnull align 8 dereferenceable(112) %0)
  %8 = sdiv i32 %3, 8
  %9 = icmp sgt i32 %3, 7
  br i1 %9, label %.lr.ph216, label %.loopexit

.lr.ph216:                                        ; preds = %6
  %10 = icmp sgt i32 %7, 0
  %wide.trip.count = zext nneg i32 %7 to i64
  br label %11

11:                                               ; preds = %.lr.ph216, %._crit_edge
  %.0212214 = phi i32 [ 0, %.lr.ph216 ], [ %42, %._crit_edge ]
  %12 = mul nuw nsw i32 %.0212214, 6
  %13 = tail call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %12)
  %14 = or disjoint i32 %12, 1
  %15 = tail call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %14)
  %16 = add nuw nsw i32 %12, 2
  %17 = tail call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %16)
  %18 = add nuw nsw i32 %12, 3
  %19 = tail call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %18)
  %20 = add nuw nsw i32 %12, 4
  %21 = tail call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %20)
  %22 = add nuw nsw i32 %12, 5
  %23 = tail call noundef ptr @_ZN17CDensityGridCubic4gam2Ei(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %22)
  %24 = shl nuw nsw i32 %.0212214, 1
  %25 = tail call noundef ptr @_ZN17CDensityGridCubic3gamEi(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %24)
  %26 = or disjoint i32 %24, 1
  %27 = tail call noundef ptr @_ZN17CDensityGridCubic3gamEi(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %26)
  %28 = tail call noundef ptr @_ZN17CDensityGridCubic2piEi(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %24)
  %29 = tail call noundef ptr @_ZN17CDensityGridCubic2piEi(ptr noundef nonnull align 8 dereferenceable(112) %0, i32 noundef %26)
  %30 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %1, i32 noundef %12)
  %31 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %1, i32 noundef %14)
  %32 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %1, i32 noundef %16)
  %33 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %1, i32 noundef %18)
  %34 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %1, i32 noundef %20)
  %35 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %1, i32 noundef %22)
  %36 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %12)
  %37 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %14)
  %38 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %16)
  %39 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %18)
  %40 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %20)
  %41 = tail call noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112) %2, i32 noundef %22)
  br i1 %10, label %.lr.ph, label %._crit_edge

._crit_edge:                                      ; preds = %.lr.ph, %11
  %42 = add nuw nsw i32 %.0212214, 1
  %exitcond218.not = icmp eq i32 %42, %8
  br i1 %exitcond218.not, label %.loopexit, label %11, !llvm.loop !21

.lr.ph:                                           ; preds = %11, %.lr.ph
  %indvars.iv = phi i64 [ %indvars.iv.next, %.lr.ph ], [ 0, %11 ]
  %43 = getelementptr inbounds double, ptr %30, i64 %indvars.iv
  %44 = load double, ptr %43, align 8, !tbaa !23
  %45 = getelementptr inbounds double, ptr %31, i64 %indvars.iv
  %46 = load double, ptr %45, align 8, !tbaa !23
  %47 = getelementptr inbounds double, ptr %32, i64 %indvars.iv
  %48 = load double, ptr %47, align 8, !tbaa !23
  %49 = getelementptr inbounds double, ptr %33, i64 %indvars.iv
  %50 = load double, ptr %49, align 8, !tbaa !23
  %51 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %44, double noundef %46, double noundef %48, double noundef %50)
  %52 = fmul double %51, 2.000000e+00
  %53 = getelementptr inbounds double, ptr %13, i64 %indvars.iv
  store double %52, ptr %53, align 8, !tbaa !23
  %54 = load double, ptr %43, align 8, !tbaa !23
  %55 = load double, ptr %45, align 8, !tbaa !23
  %56 = load double, ptr %47, align 8, !tbaa !23
  %57 = load double, ptr %49, align 8, !tbaa !23
  %58 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %54, double noundef %55, double noundef %56, double noundef %57)
  %59 = fmul double %58, 2.000000e+00
  %60 = getelementptr inbounds double, ptr %15, i64 %indvars.iv
  store double %59, ptr %60, align 8, !tbaa !23
  %61 = load double, ptr %43, align 8, !tbaa !23
  %62 = load double, ptr %45, align 8, !tbaa !23
  %63 = getelementptr inbounds double, ptr %34, i64 %indvars.iv
  %64 = load double, ptr %63, align 8, !tbaa !23
  %65 = getelementptr inbounds double, ptr %35, i64 %indvars.iv
  %66 = load double, ptr %65, align 8, !tbaa !23
  %67 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %61, double noundef %62, double noundef %64, double noundef %66)
  %68 = fmul double %67, 2.000000e+00
  %69 = getelementptr inbounds double, ptr %17, i64 %indvars.iv
  store double %68, ptr %69, align 8, !tbaa !23
  %70 = load double, ptr %43, align 8, !tbaa !23
  %71 = load double, ptr %45, align 8, !tbaa !23
  %72 = load double, ptr %63, align 8, !tbaa !23
  %73 = load double, ptr %65, align 8, !tbaa !23
  %74 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %70, double noundef %71, double noundef %72, double noundef %73)
  %75 = fmul double %74, 2.000000e+00
  %76 = getelementptr inbounds double, ptr %19, i64 %indvars.iv
  store double %75, ptr %76, align 8, !tbaa !23
  %77 = load double, ptr %47, align 8, !tbaa !23
  %78 = load double, ptr %49, align 8, !tbaa !23
  %79 = load double, ptr %63, align 8, !tbaa !23
  %80 = load double, ptr %65, align 8, !tbaa !23
  %81 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %77, double noundef %78, double noundef %79, double noundef %80)
  %82 = fmul double %81, 2.000000e+00
  %83 = getelementptr inbounds double, ptr %21, i64 %indvars.iv
  store double %82, ptr %83, align 8, !tbaa !23
  %84 = load double, ptr %47, align 8, !tbaa !23
  %85 = load double, ptr %49, align 8, !tbaa !23
  %86 = load double, ptr %63, align 8, !tbaa !23
  %87 = load double, ptr %65, align 8, !tbaa !23
  %88 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %84, double noundef %85, double noundef %86, double noundef %87)
  %89 = fmul double %88, 2.000000e+00
  %90 = getelementptr inbounds double, ptr %23, i64 %indvars.iv
  store double %89, ptr %90, align 8, !tbaa !23
  %91 = load double, ptr %63, align 8, !tbaa !23
  %92 = load double, ptr %65, align 8, !tbaa !23
  %93 = getelementptr inbounds double, ptr %36, i64 %indvars.iv
  %94 = load double, ptr %93, align 8, !tbaa !23
  %95 = getelementptr inbounds double, ptr %37, i64 %indvars.iv
  %96 = load double, ptr %95, align 8, !tbaa !23
  %97 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %91, double noundef %92, double noundef %94, double noundef %96)
  %98 = fmul double %97, 3.000000e+00
  %99 = getelementptr inbounds double, ptr %25, i64 %indvars.iv
  store double %98, ptr %99, align 8, !tbaa !23
  %100 = load double, ptr %63, align 8, !tbaa !23
  %101 = load double, ptr %65, align 8, !tbaa !23
  %102 = load double, ptr %93, align 8, !tbaa !23
  %103 = load double, ptr %95, align 8, !tbaa !23
  %104 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %100, double noundef %101, double noundef %102, double noundef %103)
  %105 = fmul double %104, 3.000000e+00
  %106 = getelementptr inbounds double, ptr %27, i64 %indvars.iv
  store double %105, ptr %106, align 8, !tbaa !23
  %107 = load double, ptr %43, align 8, !tbaa !23
  %108 = load double, ptr %45, align 8, !tbaa !23
  %109 = getelementptr inbounds double, ptr %40, i64 %indvars.iv
  %110 = load double, ptr %109, align 8, !tbaa !23
  %111 = getelementptr inbounds double, ptr %41, i64 %indvars.iv
  %112 = load double, ptr %111, align 8, !tbaa !23
  %113 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %107, double noundef %108, double noundef %110, double noundef %112)
  %114 = load double, ptr %99, align 8, !tbaa !23
  %115 = tail call double @llvm.fmuladd.f64(double %113, double 3.000000e+00, double %114)
  store double %115, ptr %99, align 8, !tbaa !23
  %116 = load double, ptr %43, align 8, !tbaa !23
  %117 = load double, ptr %45, align 8, !tbaa !23
  %118 = load double, ptr %109, align 8, !tbaa !23
  %119 = load double, ptr %111, align 8, !tbaa !23
  %120 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %116, double noundef %117, double noundef %118, double noundef %119)
  %121 = load double, ptr %106, align 8, !tbaa !23
  %122 = tail call double @llvm.fmuladd.f64(double %120, double 3.000000e+00, double %121)
  store double %122, ptr %106, align 8, !tbaa !23
  %123 = load double, ptr %47, align 8, !tbaa !23
  %124 = load double, ptr %49, align 8, !tbaa !23
  %125 = getelementptr inbounds double, ptr %38, i64 %indvars.iv
  %126 = load double, ptr %125, align 8, !tbaa !23
  %127 = getelementptr inbounds double, ptr %39, i64 %indvars.iv
  %128 = load double, ptr %127, align 8, !tbaa !23
  %129 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %123, double noundef %124, double noundef %126, double noundef %128)
  %130 = load double, ptr %99, align 8, !tbaa !23
  %131 = tail call double @llvm.fmuladd.f64(double %129, double 3.000000e+00, double %130)
  store double %131, ptr %99, align 8, !tbaa !23
  %132 = load double, ptr %47, align 8, !tbaa !23
  %133 = load double, ptr %49, align 8, !tbaa !23
  %134 = load double, ptr %125, align 8, !tbaa !23
  %135 = load double, ptr %127, align 8, !tbaa !23
  %136 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %132, double noundef %133, double noundef %134, double noundef %135)
  %137 = load double, ptr %106, align 8, !tbaa !23
  %138 = tail call double @llvm.fmuladd.f64(double %136, double 3.000000e+00, double %137)
  store double %138, ptr %106, align 8, !tbaa !23
  %139 = load double, ptr %43, align 8, !tbaa !23
  %140 = load double, ptr %45, align 8, !tbaa !23
  %141 = load double, ptr %83, align 8, !tbaa !23
  %142 = load double, ptr %90, align 8, !tbaa !23
  %143 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %139, double noundef %140, double noundef %141, double noundef %142)
  %144 = load double, ptr %47, align 8, !tbaa !23
  %145 = load double, ptr %49, align 8, !tbaa !23
  %146 = load double, ptr %69, align 8, !tbaa !23
  %147 = load double, ptr %76, align 8, !tbaa !23
  %148 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %144, double noundef %145, double noundef %146, double noundef %147)
  %149 = fadd double %143, %148
  %150 = load double, ptr %63, align 8, !tbaa !23
  %151 = load double, ptr %65, align 8, !tbaa !23
  %152 = load double, ptr %53, align 8, !tbaa !23
  %153 = load double, ptr %60, align 8, !tbaa !23
  %154 = tail call noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %150, double noundef %151, double noundef %152, double noundef %153)
  %155 = fadd double %149, %154
  %156 = getelementptr inbounds double, ptr %28, i64 %indvars.iv
  store double %155, ptr %156, align 8, !tbaa !23
  %157 = load double, ptr %43, align 8, !tbaa !23
  %158 = load double, ptr %45, align 8, !tbaa !23
  %159 = load double, ptr %83, align 8, !tbaa !23
  %160 = load double, ptr %90, align 8, !tbaa !23
  %161 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %157, double noundef %158, double noundef %159, double noundef %160)
  %162 = load double, ptr %47, align 8, !tbaa !23
  %163 = load double, ptr %49, align 8, !tbaa !23
  %164 = load double, ptr %69, align 8, !tbaa !23
  %165 = load double, ptr %76, align 8, !tbaa !23
  %166 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %162, double noundef %163, double noundef %164, double noundef %165)
  %167 = fadd double %161, %166
  %168 = load double, ptr %63, align 8, !tbaa !23
  %169 = load double, ptr %65, align 8, !tbaa !23
  %170 = load double, ptr %53, align 8, !tbaa !23
  %171 = load double, ptr %60, align 8, !tbaa !23
  %172 = tail call noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112) %0, double noundef %168, double noundef %169, double noundef %170, double noundef %171)
  %173 = fadd double %167, %172
  %174 = getelementptr inbounds double, ptr %29, i64 %indvars.iv
  store double %173, ptr %174, align 8, !tbaa !23
  %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 %.lr.ph, !llvm.loop !25

.loopexit:                                        ; preds = %._crit_edge, %6, %4
  ret void
}

declare noundef ptr @_ZNK12CDensityGrid12alphaDensityEi(ptr noundef nonnull align 8 dereferenceable(112), i32 noundef) local_unnamed_addr #1

; Function Attrs: inlinehint mustprogress nounwind uwtable
declare dso_local noundef double @_ZN17CDensityGridCubic7prod2_rEdddd(ptr noundef nonnull align 8 dereferenceable(112), double noundef, double noundef, double noundef, double noundef) local_unnamed_addr #2 align 2

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

; Function Attrs: inlinehint mustprogress nounwind uwtable
declare dso_local noundef double @_ZN17CDensityGridCubic7prod2_iEdddd(ptr noundef nonnull align 8 dereferenceable(112), double noundef, double noundef, double noundef, double noundef) local_unnamed_addr #2 align 2

attributes #0 = { mustprogress 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 #1 = { "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { inlinehint mustprogress 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 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }

!llvm.module.flags = !{!0, !1, !2, !3}
!llvm.ident = !{!4}

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !7, i64 0}
!6 = !{!"_ZTS17CDensityGridCubic", !7, i64 0, !10, i64 4, !11, i64 8}
!7 = !{!"_ZTS7dengrid", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C++ TBAA"}
!10 = !{!"int", !8, i64 0}
!11 = !{!"_ZTS18CDensityGridData2D", !12, i64 0, !17, i64 24, !17, i64 48, !17, i64 72, !10, i64 96}
!12 = !{!"_ZTSSt6vectorIdSaIdEE", !13, i64 0}
!13 = !{!"_ZTSSt12_Vector_baseIdSaIdEE", !14, i64 0}
!14 = !{!"_ZTSNSt12_Vector_baseIdSaIdEE12_Vector_implE", !15, i64 0}
!15 = !{!"_ZTSNSt12_Vector_baseIdSaIdEE17_Vector_impl_dataE", !16, i64 0, !16, i64 8, !16, i64 16}
!16 = !{!"any pointer", !8, i64 0}
!17 = !{!"_ZTSSt6vectorIiSaIiEE", !18, i64 0}
!18 = !{!"_ZTSSt12_Vector_baseIiSaIiEE", !19, i64 0}
!19 = !{!"_ZTSNSt12_Vector_baseIiSaIiEE12_Vector_implE", !20, i64 0}
!20 = !{!"_ZTSNSt12_Vector_baseIiSaIiEE17_Vector_impl_dataE", !16, i64 0, !16, i64 8, !16, i64 16}
!21 = distinct !{!21, !22}
!22 = !{!"llvm.loop.mustprogress"}
!23 = !{!24, !24, i64 0}
!24 = !{!"double", !8, i64 0}
!25 = distinct !{!25, !22}


; Effective O3 passes: EarlyCSEPass,GlobalOptPass,IndVarSimplifyPass,InferFunctionAttrsPass,InstCombinePass,JumpThreadingPass,LICMPass,LoopRotatePass,LoopSimplifyPass,LoopUnrollPass,SROAPass,SimplifyCFGPass,TailCallElimPass