; LLVM IR from QUILL-PIC/Quill :: src/init.ll (ZN14spatial_region8add_beamEdddddddd)

; 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.int_vector3d = type { i32, i32, i32 }
%class.spatial_region = type { i32, i32, i32, double, double, double, double, i32, i32, i32, double, ptr, i32, i32, i32, ptr, i32, i32, i32, i32, i32, i32, i32, %"class.std::unique_ptr", %"class.std::function", double, double, double, i32, i32, i32, ptr, double, double, double, double, ptr, double, double, double, ptr, %"class.std::__cxx11::basic_string", ptr, ptr, ptr, ptr, %class.field3d, %class.field3d.26, %class.field3d.51, %class.field3d.76, %class.field3d.101, %"class.std::vector", %"class.std::vector.129" }
%"class.std::unique_ptr" = type { %"struct.std::__uniq_ptr_data" }
%"struct.std::__uniq_ptr_data" = type { %"class.std::__uniq_ptr_impl" }
%"class.std::__uniq_ptr_impl" = type { %"class.std::tuple" }
%"class.std::tuple" = type { %"struct.std::_Tuple_impl" }
%"struct.std::_Tuple_impl" = type { %"struct.std::_Head_base.1" }
%"struct.std::_Head_base.1" = type { ptr }
%"class.std::function" = type { %"class.std::_Function_base", ptr }
%"class.std::_Function_base" = type { %"union.std::_Any_data", ptr }
%"union.std::_Any_data" = type { %"union.std::_Nocopy_types" }
%"union.std::_Nocopy_types" = type { { i64, i64 } }
%"class.std::__cxx11::basic_string" = type { %"struct.std::__cxx11::basic_string<char>::_Alloc_hider", i64, %union.anon }
%"struct.std::__cxx11::basic_string<char>::_Alloc_hider" = type { ptr }
%union.anon = type { i64, [8 x i8] }
%class.field3d = type { i32, i32, i32, %"class.std::unique_ptr.2", %"class.std::unique_ptr.10", %"class.std::unique_ptr.18" }
%"class.std::unique_ptr.2" = type { %"struct.std::__uniq_ptr_data.3" }
%"struct.std::__uniq_ptr_data.3" = type { %"class.std::__uniq_ptr_impl.4" }
%"class.std::__uniq_ptr_impl.4" = type { %"class.std::tuple.5" }
%"class.std::tuple.5" = type { %"struct.std::_Tuple_impl.6" }
%"struct.std::_Tuple_impl.6" = type { %"struct.std::_Head_base.9" }
%"struct.std::_Head_base.9" = type { ptr }
%"class.std::unique_ptr.10" = type { %"struct.std::__uniq_ptr_data.11" }
%"struct.std::__uniq_ptr_data.11" = type { %"class.std::__uniq_ptr_impl.12" }
%"class.std::__uniq_ptr_impl.12" = type { %"class.std::tuple.13" }
%"class.std::tuple.13" = type { %"struct.std::_Tuple_impl.14" }
%"struct.std::_Tuple_impl.14" = type { %"struct.std::_Head_base.17" }
%"struct.std::_Head_base.17" = type { ptr }
%"class.std::unique_ptr.18" = type { %"struct.std::__uniq_ptr_data.19" }
%"struct.std::__uniq_ptr_data.19" = type { %"class.std::__uniq_ptr_impl.20" }
%"class.std::__uniq_ptr_impl.20" = type { %"class.std::tuple.21" }
%"class.std::tuple.21" = type { %"struct.std::_Tuple_impl.22" }
%"struct.std::_Tuple_impl.22" = type { %"struct.std::_Head_base.25" }
%"struct.std::_Head_base.25" = type { ptr }
%class.field3d.26 = type { i32, i32, i32, %"class.std::unique_ptr.27", %"class.std::unique_ptr.35", %"class.std::unique_ptr.43" }
%"class.std::unique_ptr.27" = type { %"struct.std::__uniq_ptr_data.28" }
%"struct.std::__uniq_ptr_data.28" = type { %"class.std::__uniq_ptr_impl.29" }
%"class.std::__uniq_ptr_impl.29" = type { %"class.std::tuple.30" }
%"class.std::tuple.30" = type { %"struct.std::_Tuple_impl.31" }
%"struct.std::_Tuple_impl.31" = type { %"struct.std::_Head_base.34" }
%"struct.std::_Head_base.34" = type { ptr }
%"class.std::unique_ptr.35" = type { %"struct.std::__uniq_ptr_data.36" }
%"struct.std::__uniq_ptr_data.36" = type { %"class.std::__uniq_ptr_impl.37" }
%"class.std::__uniq_ptr_impl.37" = type { %"class.std::tuple.38" }
%"class.std::tuple.38" = type { %"struct.std::_Tuple_impl.39" }
%"struct.std::_Tuple_impl.39" = type { %"struct.std::_Head_base.42" }
%"struct.std::_Head_base.42" = type { ptr }
%"class.std::unique_ptr.43" = type { %"struct.std::__uniq_ptr_data.44" }
%"struct.std::__uniq_ptr_data.44" = type { %"class.std::__uniq_ptr_impl.45" }
%"class.std::__uniq_ptr_impl.45" = type { %"class.std::tuple.46" }
%"class.std::tuple.46" = type { %"struct.std::_Tuple_impl.47" }
%"struct.std::_Tuple_impl.47" = type { %"struct.std::_Head_base.50" }
%"struct.std::_Head_base.50" = type { ptr }
%class.field3d.51 = type { i32, i32, i32, %"class.std::unique_ptr.52", %"class.std::unique_ptr.60", %"class.std::unique_ptr.68" }
%"class.std::unique_ptr.52" = type { %"struct.std::__uniq_ptr_data.53" }
%"struct.std::__uniq_ptr_data.53" = type { %"class.std::__uniq_ptr_impl.54" }
%"class.std::__uniq_ptr_impl.54" = type { %"class.std::tuple.55" }
%"class.std::tuple.55" = type { %"struct.std::_Tuple_impl.56" }
%"struct.std::_Tuple_impl.56" = type { %"struct.std::_Head_base.59" }
%"struct.std::_Head_base.59" = type { ptr }
%"class.std::unique_ptr.60" = type { %"struct.std::__uniq_ptr_data.61" }
%"struct.std::__uniq_ptr_data.61" = type { %"class.std::__uniq_ptr_impl.62" }
%"class.std::__uniq_ptr_impl.62" = type { %"class.std::tuple.63" }
%"class.std::tuple.63" = type { %"struct.std::_Tuple_impl.64" }
%"struct.std::_Tuple_impl.64" = type { %"struct.std::_Head_base.67" }
%"struct.std::_Head_base.67" = type { ptr }
%"class.std::unique_ptr.68" = type { %"struct.std::__uniq_ptr_data.69" }
%"struct.std::__uniq_ptr_data.69" = type { %"class.std::__uniq_ptr_impl.70" }
%"class.std::__uniq_ptr_impl.70" = type { %"class.std::tuple.71" }
%"class.std::tuple.71" = type { %"struct.std::_Tuple_impl.72" }
%"struct.std::_Tuple_impl.72" = type { %"struct.std::_Head_base.75" }
%"struct.std::_Head_base.75" = type { ptr }
%class.field3d.76 = type { i32, i32, i32, %"class.std::unique_ptr.77", %"class.std::unique_ptr.85", %"class.std::unique_ptr.93" }
%"class.std::unique_ptr.77" = type { %"struct.std::__uniq_ptr_data.78" }
%"struct.std::__uniq_ptr_data.78" = type { %"class.std::__uniq_ptr_impl.79" }
%"class.std::__uniq_ptr_impl.79" = type { %"class.std::tuple.80" }
%"class.std::tuple.80" = type { %"struct.std::_Tuple_impl.81" }
%"struct.std::_Tuple_impl.81" = type { %"struct.std::_Head_base.84" }
%"struct.std::_Head_base.84" = type { ptr }
%"class.std::unique_ptr.85" = type { %"struct.std::__uniq_ptr_data.86" }
%"struct.std::__uniq_ptr_data.86" = type { %"class.std::__uniq_ptr_impl.87" }
%"class.std::__uniq_ptr_impl.87" = type { %"class.std::tuple.88" }
%"class.std::tuple.88" = type { %"struct.std::_Tuple_impl.89" }
%"struct.std::_Tuple_impl.89" = type { %"struct.std::_Head_base.92" }
%"struct.std::_Head_base.92" = type { ptr }
%"class.std::unique_ptr.93" = type { %"struct.std::__uniq_ptr_data.94" }
%"struct.std::__uniq_ptr_data.94" = type { %"class.std::__uniq_ptr_impl.95" }
%"class.std::__uniq_ptr_impl.95" = type { %"class.std::tuple.96" }
%"class.std::tuple.96" = type { %"struct.std::_Tuple_impl.97" }
%"struct.std::_Tuple_impl.97" = type { %"struct.std::_Head_base.100" }
%"struct.std::_Head_base.100" = type { ptr }
%class.field3d.101 = type { i32, i32, i32, %"class.std::unique_ptr.102", %"class.std::unique_ptr.110", %"class.std::unique_ptr.118" }
%"class.std::unique_ptr.102" = type { %"struct.std::__uniq_ptr_data.103" }
%"struct.std::__uniq_ptr_data.103" = type { %"class.std::__uniq_ptr_impl.104" }
%"class.std::__uniq_ptr_impl.104" = type { %"class.std::tuple.105" }
%"class.std::tuple.105" = type { %"struct.std::_Tuple_impl.106" }
%"struct.std::_Tuple_impl.106" = type { %"struct.std::_Head_base.109" }
%"struct.std::_Head_base.109" = type { ptr }
%"class.std::unique_ptr.110" = type { %"struct.std::__uniq_ptr_data.111" }
%"struct.std::__uniq_ptr_data.111" = type { %"class.std::__uniq_ptr_impl.112" }
%"class.std::__uniq_ptr_impl.112" = type { %"class.std::tuple.113" }
%"class.std::tuple.113" = type { %"struct.std::_Tuple_impl.114" }
%"struct.std::_Tuple_impl.114" = type { %"struct.std::_Head_base.117" }
%"struct.std::_Head_base.117" = type { ptr }
%"class.std::unique_ptr.118" = type { %"struct.std::__uniq_ptr_data.119" }
%"struct.std::__uniq_ptr_data.119" = type { %"class.std::__uniq_ptr_impl.120" }
%"class.std::__uniq_ptr_impl.120" = type { %"class.std::tuple.121" }
%"class.std::tuple.121" = type { %"struct.std::_Tuple_impl.122" }
%"struct.std::_Tuple_impl.122" = type { %"struct.std::_Head_base.125" }
%"struct.std::_Head_base.125" = type { ptr }
%"class.std::vector" = type { %"struct.std::_Vector_base" }
%"struct.std::_Vector_base" = type { %"struct.std::_Vector_base<field3d<double>, std::allocator<field3d<double>>>::_Vector_impl" }
%"struct.std::_Vector_base<field3d<double>, std::allocator<field3d<double>>>::_Vector_impl" = type { %"struct.std::_Vector_base<field3d<double>, std::allocator<field3d<double>>>::_Vector_impl_data" }
%"struct.std::_Vector_base<field3d<double>, std::allocator<field3d<double>>>::_Vector_impl_data" = type { ptr, ptr, ptr }
%"class.std::vector.129" = type { %"struct.std::_Vector_base.130" }
%"struct.std::_Vector_base.130" = type { %"struct.std::_Vector_base<spatial_region::deleted_particle, std::allocator<spatial_region::deleted_particle>>::_Vector_impl" }
%"struct.std::_Vector_base<spatial_region::deleted_particle, std::allocator<spatial_region::deleted_particle>>::_Vector_impl" = type { %"struct.std::_Vector_base<spatial_region::deleted_particle, std::allocator<spatial_region::deleted_particle>>::_Vector_impl_data" }
%"struct.std::_Vector_base<spatial_region::deleted_particle, std::allocator<spatial_region::deleted_particle>>::_Vector_impl_data" = type { ptr, ptr, ptr }
%struct.celle = type { double, double, double }
%struct.cellb = type { double, double, double }

; 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: nounwind
declare double @sqrt(double noundef) #1

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

; Function Attrs: nounwind
declare double @cos(double noundef) #1

; Function Attrs: nounwind
declare double @sin(double noundef) #1

; Function Attrs: inlinehint mustprogress uwtable
declare dso_local noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40), i32 noundef) #3 align 2

; Function Attrs: inlinehint mustprogress uwtable
declare dso_local noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40), i32 noundef) #3 align 2

; Function Attrs: mustprogress uwtable
declare dso_local void @_ZN14spatial_region22fill_cell_by_particlesEdR12int_vector3dS1_ddddd(ptr noundef nonnull align 8 dereferenceable(608), double noundef, ptr noundef nonnull align 4 dereferenceable(12), ptr noundef nonnull align 4 dereferenceable(12), double noundef, double noundef, double noundef, double noundef, double noundef) #4 align 2

; Function Attrs: mustprogress uwtable
define dso_local void @_ZN14spatial_region8add_beamEdddddddd(ptr noundef nonnull align 8 dereferenceable(608) %0, double noundef %1, double noundef %2, double noundef %3, double noundef %4, double noundef %5, double noundef %6, double noundef %7, double noundef %8) #4 align 2 {
  %10 = alloca ptr, align 8
  %11 = alloca double, align 8
  %12 = alloca double, align 8
  %13 = alloca double, align 8
  %14 = alloca double, align 8
  %15 = alloca double, align 8
  %16 = alloca double, align 8
  %17 = alloca double, align 8
  %18 = alloca double, align 8
  %19 = alloca %struct.int_vector3d, align 4
  %20 = alloca %struct.int_vector3d, align 4
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %25 = alloca double, align 8
  %26 = alloca double, align 8
  %27 = alloca double, align 8
  %28 = alloca double, align 8
  %29 = alloca double, align 8
  %30 = alloca i32, align 4
  %31 = alloca i32, align 4
  %32 = alloca i32, align 4
  %33 = alloca i32, align 4
  store ptr %0, ptr %10, align 8, !tbaa !5
  store double %1, ptr %11, align 8, !tbaa !9
  store double %2, ptr %12, align 8, !tbaa !9
  store double %3, ptr %13, align 8, !tbaa !9
  store double %4, ptr %14, align 8, !tbaa !9
  store double %5, ptr %15, align 8, !tbaa !9
  store double %6, ptr %16, align 8, !tbaa !9
  store double %7, ptr %17, align 8, !tbaa !9
  store double %8, ptr %18, align 8, !tbaa !9
  %34 = load ptr, ptr %10, align 8
  call void @llvm.lifetime.start.p0(i64 12, ptr %19) #6
  call void @_ZN12int_vector3dC2Ev(ptr noundef nonnull align 4 dereferenceable(12) %19)
  call void @llvm.lifetime.start.p0(i64 12, ptr %20) #6
  call void @_ZN12int_vector3dC2Ev(ptr noundef nonnull align 4 dereferenceable(12) %20)
  %35 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 7
  %36 = load i32, ptr %35, align 8, !tbaa !11
  %37 = getelementptr inbounds %struct.int_vector3d, ptr %20, i32 0, i32 0
  store i32 %36, ptr %37, align 4, !tbaa !128
  %38 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 8
  %39 = load i32, ptr %38, align 4, !tbaa !130
  %40 = getelementptr inbounds %struct.int_vector3d, ptr %20, i32 0, i32 1
  store i32 %39, ptr %40, align 4, !tbaa !131
  %41 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 9
  %42 = load i32, ptr %41, align 8, !tbaa !132
  %43 = getelementptr inbounds %struct.int_vector3d, ptr %20, i32 0, i32 2
  store i32 %42, ptr %43, align 4, !tbaa !133
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %22) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %25) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %26) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %27) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %28) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %29) #6
  %44 = load double, ptr %13, align 8, !tbaa !9
  %45 = fcmp ogt double %44, 0.000000e+00
  %46 = zext i1 %45 to i32
  %47 = load double, ptr %13, align 8, !tbaa !9
  %48 = fcmp olt double %47, 0.000000e+00
  %49 = zext i1 %48 to i32
  %50 = sub nsw i32 %46, %49
  %51 = sitofp i32 %50 to double
  store double %51, ptr %27, align 8, !tbaa !9
  %52 = load double, ptr %11, align 8, !tbaa !9
  %53 = fcmp ogt double %52, 0.000000e+00
  %54 = zext i1 %53 to i32
  %55 = load double, ptr %11, align 8, !tbaa !9
  %56 = fcmp olt double %55, 0.000000e+00
  %57 = zext i1 %56 to i32
  %58 = sub nsw i32 %54, %57
  %59 = sitofp i32 %58 to double
  store double %59, ptr %28, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 4, ptr %30) #6
  store i32 0, ptr %30, align 4, !tbaa !134
  br label %60

60:                                               ; preds = %571, %9
  %61 = load i32, ptr %30, align 4, !tbaa !134
  %62 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 0
  %63 = load i32, ptr %62, align 8, !tbaa !135
  %64 = icmp slt i32 %61, %63
  br i1 %64, label %66, label %65

65:                                               ; preds = %60
  store i32 2, ptr %31, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %30) #6
  br label %574

66:                                               ; preds = %60
  call void @llvm.lifetime.start.p0(i64 4, ptr %32) #6
  store i32 0, ptr %32, align 4, !tbaa !134
  br label %67

67:                                               ; preds = %567, %66
  %68 = load i32, ptr %32, align 4, !tbaa !134
  %69 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 1
  %70 = load i32, ptr %69, align 4, !tbaa !136
  %71 = icmp slt i32 %68, %70
  br i1 %71, label %73, label %72

72:                                               ; preds = %67
  store i32 5, ptr %31, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %32) #6
  br label %570

73:                                               ; preds = %67
  call void @llvm.lifetime.start.p0(i64 4, ptr %33) #6
  store i32 0, ptr %33, align 4, !tbaa !134
  br label %74

74:                                               ; preds = %563, %73
  %75 = load i32, ptr %33, align 4, !tbaa !134
  %76 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 2
  %77 = load i32, ptr %76, align 8, !tbaa !137
  %78 = icmp slt i32 %75, %77
  br i1 %78, label %80, label %79

79:                                               ; preds = %74
  store i32 8, ptr %31, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %33) #6
  br label %566

80:                                               ; preds = %74
  %81 = load i32, ptr %30, align 4, !tbaa !134
  %82 = getelementptr inbounds %struct.int_vector3d, ptr %19, i32 0, i32 0
  store i32 %81, ptr %82, align 4, !tbaa !128
  %83 = load i32, ptr %32, align 4, !tbaa !134
  %84 = getelementptr inbounds %struct.int_vector3d, ptr %19, i32 0, i32 1
  store i32 %83, ptr %84, align 4, !tbaa !131
  %85 = load i32, ptr %33, align 4, !tbaa !134
  %86 = getelementptr inbounds %struct.int_vector3d, ptr %19, i32 0, i32 2
  store i32 %85, ptr %86, align 4, !tbaa !133
  %87 = load i32, ptr %32, align 4, !tbaa !134
  %88 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 1
  %89 = load i32, ptr %88, align 4, !tbaa !136
  %90 = sdiv i32 %89, 2
  %91 = sub nsw i32 %87, %90
  %92 = sitofp i32 %91 to double
  %93 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 4
  %94 = load double, ptr %93, align 8, !tbaa !138
  %95 = load double, ptr %17, align 8, !tbaa !9
  %96 = fneg double %95
  %97 = call double @llvm.fmuladd.f64(double %92, double %94, double %96)
  %98 = load i32, ptr %30, align 4, !tbaa !134
  %99 = sitofp i32 %98 to double
  %100 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 3
  %101 = load double, ptr %100, align 8, !tbaa !139
  %102 = load double, ptr %16, align 8, !tbaa !9
  %103 = fneg double %102
  %104 = call double @llvm.fmuladd.f64(double %99, double %101, double %103)
  %105 = call double @atan2(double noundef %97, double noundef %104) #6
  store double %105, ptr %29, align 8, !tbaa !9
  %106 = load i32, ptr %30, align 4, !tbaa !134
  %107 = sitofp i32 %106 to double
  %108 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 3
  %109 = load double, ptr %108, align 8, !tbaa !139
  %110 = load double, ptr %16, align 8, !tbaa !9
  %111 = fneg double %110
  %112 = call double @llvm.fmuladd.f64(double %107, double %109, double %111)
  %113 = load i32, ptr %30, align 4, !tbaa !134
  %114 = sitofp i32 %113 to double
  %115 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 3
  %116 = load double, ptr %115, align 8, !tbaa !139
  %117 = load double, ptr %16, align 8, !tbaa !9
  %118 = fneg double %117
  %119 = call double @llvm.fmuladd.f64(double %114, double %116, double %118)
  %120 = load i32, ptr %32, align 4, !tbaa !134
  %121 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 1
  %122 = load i32, ptr %121, align 4, !tbaa !136
  %123 = sdiv i32 %122, 2
  %124 = sub nsw i32 %120, %123
  %125 = sitofp i32 %124 to double
  %126 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 4
  %127 = load double, ptr %126, align 8, !tbaa !138
  %128 = load double, ptr %17, align 8, !tbaa !9
  %129 = fneg double %128
  %130 = call double @llvm.fmuladd.f64(double %125, double %127, double %129)
  %131 = load i32, ptr %32, align 4, !tbaa !134
  %132 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 1
  %133 = load i32, ptr %132, align 4, !tbaa !136
  %134 = sdiv i32 %133, 2
  %135 = sub nsw i32 %131, %134
  %136 = sitofp i32 %135 to double
  %137 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 4
  %138 = load double, ptr %137, align 8, !tbaa !138
  %139 = load double, ptr %17, align 8, !tbaa !9
  %140 = fneg double %139
  %141 = call double @llvm.fmuladd.f64(double %136, double %138, double %140)
  %142 = fmul double %130, %141
  %143 = call double @llvm.fmuladd.f64(double %112, double %119, double %142)
  %144 = call double @sqrt(double noundef %143) #6
  %145 = load double, ptr %29, align 8, !tbaa !9
  %146 = load double, ptr %18, align 8, !tbaa !9
  %147 = fsub double %145, %146
  %148 = call double @cos(double noundef %147) #6
  %149 = fmul double %144, %148
  store double %149, ptr %21, align 8, !tbaa !9
  %150 = load i32, ptr %30, align 4, !tbaa !134
  %151 = sitofp i32 %150 to double
  %152 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 3
  %153 = load double, ptr %152, align 8, !tbaa !139
  %154 = load double, ptr %16, align 8, !tbaa !9
  %155 = fneg double %154
  %156 = call double @llvm.fmuladd.f64(double %151, double %153, double %155)
  %157 = load i32, ptr %30, align 4, !tbaa !134
  %158 = sitofp i32 %157 to double
  %159 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 3
  %160 = load double, ptr %159, align 8, !tbaa !139
  %161 = load double, ptr %16, align 8, !tbaa !9
  %162 = fneg double %161
  %163 = call double @llvm.fmuladd.f64(double %158, double %160, double %162)
  %164 = load i32, ptr %32, align 4, !tbaa !134
  %165 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 1
  %166 = load i32, ptr %165, align 4, !tbaa !136
  %167 = sdiv i32 %166, 2
  %168 = sub nsw i32 %164, %167
  %169 = sitofp i32 %168 to double
  %170 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 4
  %171 = load double, ptr %170, align 8, !tbaa !138
  %172 = load double, ptr %17, align 8, !tbaa !9
  %173 = fneg double %172
  %174 = call double @llvm.fmuladd.f64(double %169, double %171, double %173)
  %175 = load i32, ptr %32, align 4, !tbaa !134
  %176 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 1
  %177 = load i32, ptr %176, align 4, !tbaa !136
  %178 = sdiv i32 %177, 2
  %179 = sub nsw i32 %175, %178
  %180 = sitofp i32 %179 to double
  %181 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 4
  %182 = load double, ptr %181, align 8, !tbaa !138
  %183 = load double, ptr %17, align 8, !tbaa !9
  %184 = fneg double %183
  %185 = call double @llvm.fmuladd.f64(double %180, double %182, double %184)
  %186 = fmul double %174, %185
  %187 = call double @llvm.fmuladd.f64(double %156, double %163, double %186)
  %188 = call double @sqrt(double noundef %187) #6
  %189 = load double, ptr %29, align 8, !tbaa !9
  %190 = load double, ptr %18, align 8, !tbaa !9
  %191 = fsub double %189, %190
  %192 = call double @sin(double noundef %191) #6
  %193 = fmul double %188, %192
  store double %193, ptr %22, align 8, !tbaa !9
  %194 = load i32, ptr %33, align 4, !tbaa !134
  %195 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 2
  %196 = load i32, ptr %195, align 8, !tbaa !137
  %197 = sdiv i32 %196, 2
  %198 = sub nsw i32 %194, %197
  %199 = sitofp i32 %198 to double
  %200 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 5
  %201 = load double, ptr %200, align 8, !tbaa !140
  %202 = fmul double %199, %201
  store double %202, ptr %23, align 8, !tbaa !9
  %203 = load double, ptr %22, align 8, !tbaa !9
  %204 = load double, ptr %22, align 8, !tbaa !9
  %205 = load double, ptr %23, align 8, !tbaa !9
  %206 = load double, ptr %23, align 8, !tbaa !9
  %207 = fmul double %205, %206
  %208 = call double @llvm.fmuladd.f64(double %203, double %204, double %207)
  %209 = call double @sqrt(double noundef %208) #6
  store double %209, ptr %24, align 8, !tbaa !9
  %210 = load double, ptr %21, align 8, !tbaa !9
  %211 = load double, ptr %14, align 8, !tbaa !9
  %212 = fneg double %211
  %213 = fcmp ogt double %210, %212
  br i1 %213, label %214, label %401

214:                                              ; preds = %80
  %215 = load double, ptr %21, align 8, !tbaa !9
  %216 = load double, ptr %14, align 8, !tbaa !9
  %217 = fcmp olt double %215, %216
  br i1 %217, label %218, label %401

218:                                              ; preds = %214
  %219 = load double, ptr %24, align 8, !tbaa !9
  %220 = load double, ptr %15, align 8, !tbaa !9
  %221 = fcmp olt double %219, %220
  br i1 %221, label %222, label %401

222:                                              ; preds = %218
  %223 = load double, ptr %12, align 8, !tbaa !9
  %224 = load double, ptr %21, align 8, !tbaa !9
  %225 = load double, ptr %21, align 8, !tbaa !9
  %226 = fmul double %224, %225
  %227 = load double, ptr %14, align 8, !tbaa !9
  %228 = load double, ptr %14, align 8, !tbaa !9
  %229 = fmul double %227, %228
  %230 = fdiv double %226, %229
  %231 = fsub double 1.000000e+00, %230
  %232 = fmul double %223, %231
  %233 = load double, ptr %24, align 8, !tbaa !9
  %234 = load double, ptr %24, align 8, !tbaa !9
  %235 = fmul double %233, %234
  %236 = load double, ptr %15, align 8, !tbaa !9
  %237 = load double, ptr %15, align 8, !tbaa !9
  %238 = fmul double %236, %237
  %239 = fdiv double %235, %238
  %240 = fsub double 1.000000e+00, %239
  %241 = fmul double %232, %240
  store double %241, ptr %25, align 8, !tbaa !9
  %242 = load double, ptr %11, align 8, !tbaa !9
  %243 = load double, ptr %25, align 8, !tbaa !9
  %244 = load double, ptr %13, align 8, !tbaa !9
  %245 = load double, ptr %18, align 8, !tbaa !9
  %246 = call double @cos(double noundef %245) #6
  %247 = fmul double %244, %246
  %248 = load double, ptr %13, align 8, !tbaa !9
  %249 = load double, ptr %18, align 8, !tbaa !9
  %250 = call double @sin(double noundef %249) #6
  %251 = fmul double %248, %250
  call void @_ZN14spatial_region22fill_cell_by_particlesEdR12int_vector3dS1_ddddd(ptr noundef nonnull align 8 dereferenceable(608) %34, double noundef %242, ptr noundef nonnull align 4 dereferenceable(12) %19, ptr noundef nonnull align 4 dereferenceable(12) %20, double noundef %243, double noundef %247, double noundef %251, double noundef 0.000000e+00, double noundef 0.000000e+00)
  %252 = load double, ptr %12, align 8, !tbaa !9
  %253 = fmul double 5.000000e-01, %252
  %254 = load double, ptr %21, align 8, !tbaa !9
  %255 = load double, ptr %21, align 8, !tbaa !9
  %256 = fmul double %254, %255
  %257 = load double, ptr %14, align 8, !tbaa !9
  %258 = load double, ptr %14, align 8, !tbaa !9
  %259 = fmul double %257, %258
  %260 = fdiv double %256, %259
  %261 = fsub double 1.000000e+00, %260
  %262 = fmul double %253, %261
  %263 = load double, ptr %24, align 8, !tbaa !9
  %264 = load double, ptr %24, align 8, !tbaa !9
  %265 = fmul double %263, %264
  %266 = load double, ptr %15, align 8, !tbaa !9
  %267 = fmul double 2.000000e+00, %266
  %268 = load double, ptr %15, align 8, !tbaa !9
  %269 = fmul double %267, %268
  %270 = fdiv double %265, %269
  %271 = fsub double 1.000000e+00, %270
  %272 = fmul double %262, %271
  store double %272, ptr %26, align 8, !tbaa !9
  %273 = load double, ptr %28, align 8, !tbaa !9
  %274 = load double, ptr %26, align 8, !tbaa !9
  %275 = fmul double %273, %274
  %276 = load double, ptr %22, align 8, !tbaa !9
  %277 = fmul double %275, %276
  %278 = load double, ptr %18, align 8, !tbaa !9
  %279 = call double @sin(double noundef %278) #6
  %280 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 46
  %281 = load i32, ptr %30, align 4, !tbaa !134
  %282 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %280, i32 noundef %281)
  %283 = load ptr, ptr %282, align 8, !tbaa !5
  %284 = load i32, ptr %32, align 4, !tbaa !134
  %285 = sext i32 %284 to i64
  %286 = getelementptr inbounds ptr, ptr %283, i64 %285
  %287 = load ptr, ptr %286, align 8, !tbaa !5
  %288 = load i32, ptr %33, align 4, !tbaa !134
  %289 = sext i32 %288 to i64
  %290 = getelementptr inbounds %struct.celle, ptr %287, i64 %289
  %291 = getelementptr inbounds %struct.celle, ptr %290, i32 0, i32 0
  %292 = load double, ptr %291, align 8, !tbaa !141
  %293 = fneg double %277
  %294 = call double @llvm.fmuladd.f64(double %293, double %279, double %292)
  store double %294, ptr %291, align 8, !tbaa !141
  %295 = load double, ptr %28, align 8, !tbaa !9
  %296 = load double, ptr %26, align 8, !tbaa !9
  %297 = fmul double %295, %296
  %298 = load double, ptr %22, align 8, !tbaa !9
  %299 = fmul double %297, %298
  %300 = load double, ptr %18, align 8, !tbaa !9
  %301 = call double @cos(double noundef %300) #6
  %302 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 46
  %303 = load i32, ptr %30, align 4, !tbaa !134
  %304 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %302, i32 noundef %303)
  %305 = load ptr, ptr %304, align 8, !tbaa !5
  %306 = load i32, ptr %32, align 4, !tbaa !134
  %307 = sext i32 %306 to i64
  %308 = getelementptr inbounds ptr, ptr %305, i64 %307
  %309 = load ptr, ptr %308, align 8, !tbaa !5
  %310 = load i32, ptr %33, align 4, !tbaa !134
  %311 = sext i32 %310 to i64
  %312 = getelementptr inbounds %struct.celle, ptr %309, i64 %311
  %313 = getelementptr inbounds %struct.celle, ptr %312, i32 0, i32 1
  %314 = load double, ptr %313, align 8, !tbaa !143
  %315 = call double @llvm.fmuladd.f64(double %299, double %301, double %314)
  store double %315, ptr %313, align 8, !tbaa !143
  %316 = load double, ptr %28, align 8, !tbaa !9
  %317 = load double, ptr %26, align 8, !tbaa !9
  %318 = fmul double %316, %317
  %319 = load double, ptr %23, align 8, !tbaa !9
  %320 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 46
  %321 = load i32, ptr %30, align 4, !tbaa !134
  %322 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %320, i32 noundef %321)
  %323 = load ptr, ptr %322, align 8, !tbaa !5
  %324 = load i32, ptr %32, align 4, !tbaa !134
  %325 = sext i32 %324 to i64
  %326 = getelementptr inbounds ptr, ptr %323, i64 %325
  %327 = load ptr, ptr %326, align 8, !tbaa !5
  %328 = load i32, ptr %33, align 4, !tbaa !134
  %329 = sext i32 %328 to i64
  %330 = getelementptr inbounds %struct.celle, ptr %327, i64 %329
  %331 = getelementptr inbounds %struct.celle, ptr %330, i32 0, i32 2
  %332 = load double, ptr %331, align 8, !tbaa !144
  %333 = call double @llvm.fmuladd.f64(double %318, double %319, double %332)
  store double %333, ptr %331, align 8, !tbaa !144
  %334 = load double, ptr %28, align 8, !tbaa !9
  %335 = load double, ptr %27, align 8, !tbaa !9
  %336 = fmul double %334, %335
  %337 = load double, ptr %26, align 8, !tbaa !9
  %338 = fmul double %336, %337
  %339 = load double, ptr %23, align 8, !tbaa !9
  %340 = fmul double %338, %339
  %341 = load double, ptr %18, align 8, !tbaa !9
  %342 = call double @sin(double noundef %341) #6
  %343 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 47
  %344 = load i32, ptr %30, align 4, !tbaa !134
  %345 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %343, i32 noundef %344)
  %346 = load ptr, ptr %345, align 8, !tbaa !5
  %347 = load i32, ptr %32, align 4, !tbaa !134
  %348 = sext i32 %347 to i64
  %349 = getelementptr inbounds ptr, ptr %346, i64 %348
  %350 = load ptr, ptr %349, align 8, !tbaa !5
  %351 = load i32, ptr %33, align 4, !tbaa !134
  %352 = sext i32 %351 to i64
  %353 = getelementptr inbounds %struct.cellb, ptr %350, i64 %352
  %354 = getelementptr inbounds %struct.cellb, ptr %353, i32 0, i32 0
  %355 = load double, ptr %354, align 8, !tbaa !145
  %356 = call double @llvm.fmuladd.f64(double %340, double %342, double %355)
  store double %356, ptr %354, align 8, !tbaa !145
  %357 = load double, ptr %28, align 8, !tbaa !9
  %358 = load double, ptr %27, align 8, !tbaa !9
  %359 = fmul double %357, %358
  %360 = load double, ptr %26, align 8, !tbaa !9
  %361 = fmul double %359, %360
  %362 = load double, ptr %23, align 8, !tbaa !9
  %363 = fmul double %361, %362
  %364 = load double, ptr %18, align 8, !tbaa !9
  %365 = call double @cos(double noundef %364) #6
  %366 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 47
  %367 = load i32, ptr %30, align 4, !tbaa !134
  %368 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %366, i32 noundef %367)
  %369 = load ptr, ptr %368, align 8, !tbaa !5
  %370 = load i32, ptr %32, align 4, !tbaa !134
  %371 = sext i32 %370 to i64
  %372 = getelementptr inbounds ptr, ptr %369, i64 %371
  %373 = load ptr, ptr %372, align 8, !tbaa !5
  %374 = load i32, ptr %33, align 4, !tbaa !134
  %375 = sext i32 %374 to i64
  %376 = getelementptr inbounds %struct.cellb, ptr %373, i64 %375
  %377 = getelementptr inbounds %struct.cellb, ptr %376, i32 0, i32 1
  %378 = load double, ptr %377, align 8, !tbaa !147
  %379 = fneg double %363
  %380 = call double @llvm.fmuladd.f64(double %379, double %365, double %378)
  store double %380, ptr %377, align 8, !tbaa !147
  %381 = load double, ptr %28, align 8, !tbaa !9
  %382 = load double, ptr %27, align 8, !tbaa !9
  %383 = fmul double %381, %382
  %384 = load double, ptr %26, align 8, !tbaa !9
  %385 = fmul double %383, %384
  %386 = load double, ptr %22, align 8, !tbaa !9
  %387 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 47
  %388 = load i32, ptr %30, align 4, !tbaa !134
  %389 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %387, i32 noundef %388)
  %390 = load ptr, ptr %389, align 8, !tbaa !5
  %391 = load i32, ptr %32, align 4, !tbaa !134
  %392 = sext i32 %391 to i64
  %393 = getelementptr inbounds ptr, ptr %390, i64 %392
  %394 = load ptr, ptr %393, align 8, !tbaa !5
  %395 = load i32, ptr %33, align 4, !tbaa !134
  %396 = sext i32 %395 to i64
  %397 = getelementptr inbounds %struct.cellb, ptr %394, i64 %396
  %398 = getelementptr inbounds %struct.cellb, ptr %397, i32 0, i32 2
  %399 = load double, ptr %398, align 8, !tbaa !148
  %400 = call double @llvm.fmuladd.f64(double %385, double %386, double %399)
  store double %400, ptr %398, align 8, !tbaa !148
  br label %401

401:                                              ; preds = %222, %218, %214, %80
  %402 = load double, ptr %21, align 8, !tbaa !9
  %403 = load double, ptr %14, align 8, !tbaa !9
  %404 = fneg double %403
  %405 = fcmp ogt double %402, %404
  br i1 %405, label %406, label %562

406:                                              ; preds = %401
  %407 = load double, ptr %21, align 8, !tbaa !9
  %408 = load double, ptr %14, align 8, !tbaa !9
  %409 = fcmp olt double %407, %408
  br i1 %409, label %410, label %562

410:                                              ; preds = %406
  %411 = load double, ptr %24, align 8, !tbaa !9
  %412 = load double, ptr %15, align 8, !tbaa !9
  %413 = fcmp oge double %411, %412
  br i1 %413, label %414, label %562

414:                                              ; preds = %410
  %415 = load double, ptr %12, align 8, !tbaa !9
  %416 = fmul double 2.500000e-01, %415
  %417 = load double, ptr %21, align 8, !tbaa !9
  %418 = load double, ptr %21, align 8, !tbaa !9
  %419 = fmul double %417, %418
  %420 = load double, ptr %14, align 8, !tbaa !9
  %421 = load double, ptr %14, align 8, !tbaa !9
  %422 = fmul double %420, %421
  %423 = fdiv double %419, %422
  %424 = fsub double 1.000000e+00, %423
  %425 = fmul double %416, %424
  %426 = load double, ptr %15, align 8, !tbaa !9
  %427 = fmul double %425, %426
  %428 = load double, ptr %15, align 8, !tbaa !9
  %429 = fmul double %427, %428
  %430 = load double, ptr %24, align 8, !tbaa !9
  %431 = load double, ptr %24, align 8, !tbaa !9
  %432 = fmul double %430, %431
  %433 = fdiv double %429, %432
  store double %433, ptr %26, align 8, !tbaa !9
  %434 = load double, ptr %28, align 8, !tbaa !9
  %435 = load double, ptr %26, align 8, !tbaa !9
  %436 = fmul double %434, %435
  %437 = load double, ptr %22, align 8, !tbaa !9
  %438 = fmul double %436, %437
  %439 = load double, ptr %18, align 8, !tbaa !9
  %440 = call double @sin(double noundef %439) #6
  %441 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 46
  %442 = load i32, ptr %30, align 4, !tbaa !134
  %443 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %441, i32 noundef %442)
  %444 = load ptr, ptr %443, align 8, !tbaa !5
  %445 = load i32, ptr %32, align 4, !tbaa !134
  %446 = sext i32 %445 to i64
  %447 = getelementptr inbounds ptr, ptr %444, i64 %446
  %448 = load ptr, ptr %447, align 8, !tbaa !5
  %449 = load i32, ptr %33, align 4, !tbaa !134
  %450 = sext i32 %449 to i64
  %451 = getelementptr inbounds %struct.celle, ptr %448, i64 %450
  %452 = getelementptr inbounds %struct.celle, ptr %451, i32 0, i32 0
  %453 = load double, ptr %452, align 8, !tbaa !141
  %454 = fneg double %438
  %455 = call double @llvm.fmuladd.f64(double %454, double %440, double %453)
  store double %455, ptr %452, align 8, !tbaa !141
  %456 = load double, ptr %28, align 8, !tbaa !9
  %457 = load double, ptr %26, align 8, !tbaa !9
  %458 = fmul double %456, %457
  %459 = load double, ptr %22, align 8, !tbaa !9
  %460 = fmul double %458, %459
  %461 = load double, ptr %18, align 8, !tbaa !9
  %462 = call double @cos(double noundef %461) #6
  %463 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 46
  %464 = load i32, ptr %30, align 4, !tbaa !134
  %465 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %463, i32 noundef %464)
  %466 = load ptr, ptr %465, align 8, !tbaa !5
  %467 = load i32, ptr %32, align 4, !tbaa !134
  %468 = sext i32 %467 to i64
  %469 = getelementptr inbounds ptr, ptr %466, i64 %468
  %470 = load ptr, ptr %469, align 8, !tbaa !5
  %471 = load i32, ptr %33, align 4, !tbaa !134
  %472 = sext i32 %471 to i64
  %473 = getelementptr inbounds %struct.celle, ptr %470, i64 %472
  %474 = getelementptr inbounds %struct.celle, ptr %473, i32 0, i32 1
  %475 = load double, ptr %474, align 8, !tbaa !143
  %476 = call double @llvm.fmuladd.f64(double %460, double %462, double %475)
  store double %476, ptr %474, align 8, !tbaa !143
  %477 = load double, ptr %28, align 8, !tbaa !9
  %478 = load double, ptr %26, align 8, !tbaa !9
  %479 = fmul double %477, %478
  %480 = load double, ptr %23, align 8, !tbaa !9
  %481 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 46
  %482 = load i32, ptr %30, align 4, !tbaa !134
  %483 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %481, i32 noundef %482)
  %484 = load ptr, ptr %483, align 8, !tbaa !5
  %485 = load i32, ptr %32, align 4, !tbaa !134
  %486 = sext i32 %485 to i64
  %487 = getelementptr inbounds ptr, ptr %484, i64 %486
  %488 = load ptr, ptr %487, align 8, !tbaa !5
  %489 = load i32, ptr %33, align 4, !tbaa !134
  %490 = sext i32 %489 to i64
  %491 = getelementptr inbounds %struct.celle, ptr %488, i64 %490
  %492 = getelementptr inbounds %struct.celle, ptr %491, i32 0, i32 2
  %493 = load double, ptr %492, align 8, !tbaa !144
  %494 = call double @llvm.fmuladd.f64(double %479, double %480, double %493)
  store double %494, ptr %492, align 8, !tbaa !144
  %495 = load double, ptr %28, align 8, !tbaa !9
  %496 = load double, ptr %27, align 8, !tbaa !9
  %497 = fmul double %495, %496
  %498 = load double, ptr %26, align 8, !tbaa !9
  %499 = fmul double %497, %498
  %500 = load double, ptr %23, align 8, !tbaa !9
  %501 = fmul double %499, %500
  %502 = load double, ptr %18, align 8, !tbaa !9
  %503 = call double @sin(double noundef %502) #6
  %504 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 47
  %505 = load i32, ptr %30, align 4, !tbaa !134
  %506 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %504, i32 noundef %505)
  %507 = load ptr, ptr %506, align 8, !tbaa !5
  %508 = load i32, ptr %32, align 4, !tbaa !134
  %509 = sext i32 %508 to i64
  %510 = getelementptr inbounds ptr, ptr %507, i64 %509
  %511 = load ptr, ptr %510, align 8, !tbaa !5
  %512 = load i32, ptr %33, align 4, !tbaa !134
  %513 = sext i32 %512 to i64
  %514 = getelementptr inbounds %struct.cellb, ptr %511, i64 %513
  %515 = getelementptr inbounds %struct.cellb, ptr %514, i32 0, i32 0
  %516 = load double, ptr %515, align 8, !tbaa !145
  %517 = call double @llvm.fmuladd.f64(double %501, double %503, double %516)
  store double %517, ptr %515, align 8, !tbaa !145
  %518 = load double, ptr %28, align 8, !tbaa !9
  %519 = load double, ptr %27, align 8, !tbaa !9
  %520 = fmul double %518, %519
  %521 = load double, ptr %26, align 8, !tbaa !9
  %522 = fmul double %520, %521
  %523 = load double, ptr %23, align 8, !tbaa !9
  %524 = fmul double %522, %523
  %525 = load double, ptr %18, align 8, !tbaa !9
  %526 = call double @cos(double noundef %525) #6
  %527 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 47
  %528 = load i32, ptr %30, align 4, !tbaa !134
  %529 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %527, i32 noundef %528)
  %530 = load ptr, ptr %529, align 8, !tbaa !5
  %531 = load i32, ptr %32, align 4, !tbaa !134
  %532 = sext i32 %531 to i64
  %533 = getelementptr inbounds ptr, ptr %530, i64 %532
  %534 = load ptr, ptr %533, align 8, !tbaa !5
  %535 = load i32, ptr %33, align 4, !tbaa !134
  %536 = sext i32 %535 to i64
  %537 = getelementptr inbounds %struct.cellb, ptr %534, i64 %536
  %538 = getelementptr inbounds %struct.cellb, ptr %537, i32 0, i32 1
  %539 = load double, ptr %538, align 8, !tbaa !147
  %540 = fneg double %524
  %541 = call double @llvm.fmuladd.f64(double %540, double %526, double %539)
  store double %541, ptr %538, align 8, !tbaa !147
  %542 = load double, ptr %28, align 8, !tbaa !9
  %543 = load double, ptr %27, align 8, !tbaa !9
  %544 = fmul double %542, %543
  %545 = load double, ptr %26, align 8, !tbaa !9
  %546 = fmul double %544, %545
  %547 = load double, ptr %22, align 8, !tbaa !9
  %548 = getelementptr inbounds %class.spatial_region, ptr %34, i32 0, i32 47
  %549 = load i32, ptr %30, align 4, !tbaa !134
  %550 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %548, i32 noundef %549)
  %551 = load ptr, ptr %550, align 8, !tbaa !5
  %552 = load i32, ptr %32, align 4, !tbaa !134
  %553 = sext i32 %552 to i64
  %554 = getelementptr inbounds ptr, ptr %551, i64 %553
  %555 = load ptr, ptr %554, align 8, !tbaa !5
  %556 = load i32, ptr %33, align 4, !tbaa !134
  %557 = sext i32 %556 to i64
  %558 = getelementptr inbounds %struct.cellb, ptr %555, i64 %557
  %559 = getelementptr inbounds %struct.cellb, ptr %558, i32 0, i32 2
  %560 = load double, ptr %559, align 8, !tbaa !148
  %561 = call double @llvm.fmuladd.f64(double %546, double %547, double %560)
  store double %561, ptr %559, align 8, !tbaa !148
  br label %562

562:                                              ; preds = %414, %410, %406, %401
  br label %563

563:                                              ; preds = %562
  %564 = load i32, ptr %33, align 4, !tbaa !134
  %565 = add nsw i32 %564, 1
  store i32 %565, ptr %33, align 4, !tbaa !134
  br label %74, !llvm.loop !149

566:                                              ; preds = %79
  br label %567

567:                                              ; preds = %566
  %568 = load i32, ptr %32, align 4, !tbaa !134
  %569 = add nsw i32 %568, 1
  store i32 %569, ptr %32, align 4, !tbaa !134
  br label %67, !llvm.loop !151

570:                                              ; preds = %72
  br label %571

571:                                              ; preds = %570
  %572 = load i32, ptr %30, align 4, !tbaa !134
  %573 = add nsw i32 %572, 1
  store i32 %573, ptr %30, align 4, !tbaa !134
  br label %60, !llvm.loop !152

574:                                              ; preds = %65
  call void @llvm.lifetime.end.p0(i64 8, ptr %29) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %26) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %22) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #6
  call void @llvm.lifetime.end.p0(i64 12, ptr %20) #6
  call void @llvm.lifetime.end.p0(i64 12, ptr %19) #6
  ret void
}

; Function Attrs: mustprogress nounwind uwtable
declare dso_local void @_ZN12int_vector3dC2Ev(ptr noundef nonnull align 4 dereferenceable(12)) unnamed_addr #5 align 2

; Function Attrs: nounwind
declare double @atan2(double noundef, double noundef) #1

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nounwind "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #3 = { inlinehint 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 #4 = { 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 #5 = { 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 #6 = { 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 = !{!"double", !7, i64 0}
!11 = !{!12, !13, i64 48}
!12 = !{!"_ZTS14spatial_region", !13, i64 0, !13, i64 4, !13, i64 8, !10, i64 16, !10, i64 24, !10, i64 32, !10, i64 40, !13, i64 48, !13, i64 52, !13, i64 56, !10, i64 64, !6, i64 72, !13, i64 80, !13, i64 84, !13, i64 88, !6, i64 96, !13, i64 104, !13, i64 108, !13, i64 112, !13, i64 116, !13, i64 120, !13, i64 124, !13, i64 128, !14, i64 136, !20, i64 144, !10, i64 176, !10, i64 184, !10, i64 192, !13, i64 200, !13, i64 204, !13, i64 208, !6, i64 216, !10, i64 224, !10, i64 232, !10, i64 240, !10, i64 248, !6, i64 256, !10, i64 264, !10, i64 272, !10, i64 280, !6, i64 288, !22, i64 296, !6, i64 328, !6, i64 336, !6, i64 344, !6, i64 352, !25, i64 360, !44, i64 400, !63, i64 440, !82, i64 480, !101, i64 520, !120, i64 560, !124, i64 584}
!13 = !{!"int", !7, i64 0}
!14 = !{!"_ZTSSt10unique_ptrI14maxwell_solverSt14default_deleteIS0_EE", !15, i64 0}
!15 = !{!"_ZTSSt15__uniq_ptr_dataI14maxwell_solverSt14default_deleteIS0_ELb1ELb1EE", !16, i64 0}
!16 = !{!"_ZTSSt15__uniq_ptr_implI14maxwell_solverSt14default_deleteIS0_EE", !17, i64 0}
!17 = !{!"_ZTSSt5tupleIJP14maxwell_solverSt14default_deleteIS0_EEE", !18, i64 0}
!18 = !{!"_ZTSSt11_Tuple_implILm0EJP14maxwell_solverSt14default_deleteIS0_EEE", !19, i64 0}
!19 = !{!"_ZTSSt10_Head_baseILm0EP14maxwell_solverLb0EE", !6, i64 0}
!20 = !{!"_ZTSSt8functionIFvR8particleR8vector3dS3_RdEE", !21, i64 0, !6, i64 24}
!21 = !{!"_ZTSSt14_Function_base", !7, i64 0, !6, i64 16}
!22 = !{!"_ZTSNSt7__cxx1112basic_stringIcSt11char_traitsIcESaIcEEE", !23, i64 0, !24, i64 8, !7, i64 16}
!23 = !{!"_ZTSNSt7__cxx1112basic_stringIcSt11char_traitsIcESaIcEE12_Alloc_hiderE", !6, i64 0}
!24 = !{!"long", !7, i64 0}
!25 = !{!"_ZTS7field3dI5celleE", !13, i64 0, !13, i64 4, !13, i64 8, !26, i64 16, !32, i64 24, !38, i64 32}
!26 = !{!"_ZTSSt10unique_ptrIA_PP5celleSt14default_deleteIS3_EE", !27, i64 0}
!27 = !{!"_ZTSSt15__uniq_ptr_dataIPP5celleSt14default_deleteIA_S2_ELb1ELb1EE", !28, i64 0}
!28 = !{!"_ZTSSt15__uniq_ptr_implIPP5celleSt14default_deleteIA_S2_EE", !29, i64 0}
!29 = !{!"_ZTSSt5tupleIJPPP5celleSt14default_deleteIA_S2_EEE", !30, i64 0}
!30 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5celleSt14default_deleteIA_S2_EEE", !31, i64 0}
!31 = !{!"_ZTSSt10_Head_baseILm0EPPP5celleLb0EE", !6, i64 0}
!32 = !{!"_ZTSSt10unique_ptrIA_P5celleSt14default_deleteIS2_EE", !33, i64 0}
!33 = !{!"_ZTSSt15__uniq_ptr_dataIP5celleSt14default_deleteIA_S1_ELb1ELb1EE", !34, i64 0}
!34 = !{!"_ZTSSt15__uniq_ptr_implIP5celleSt14default_deleteIA_S1_EE", !35, i64 0}
!35 = !{!"_ZTSSt5tupleIJPP5celleSt14default_deleteIA_S1_EEE", !36, i64 0}
!36 = !{!"_ZTSSt11_Tuple_implILm0EJPP5celleSt14default_deleteIA_S1_EEE", !37, i64 0}
!37 = !{!"_ZTSSt10_Head_baseILm0EPP5celleLb0EE", !6, i64 0}
!38 = !{!"_ZTSSt10unique_ptrIA_5celleSt14default_deleteIS1_EE", !39, i64 0}
!39 = !{!"_ZTSSt15__uniq_ptr_dataI5celleSt14default_deleteIA_S0_ELb1ELb1EE", !40, i64 0}
!40 = !{!"_ZTSSt15__uniq_ptr_implI5celleSt14default_deleteIA_S0_EE", !41, i64 0}
!41 = !{!"_ZTSSt5tupleIJP5celleSt14default_deleteIA_S0_EEE", !42, i64 0}
!42 = !{!"_ZTSSt11_Tuple_implILm0EJP5celleSt14default_deleteIA_S0_EEE", !43, i64 0}
!43 = !{!"_ZTSSt10_Head_baseILm0EP5celleLb0EE", !6, i64 0}
!44 = !{!"_ZTS7field3dI5cellbE", !13, i64 0, !13, i64 4, !13, i64 8, !45, i64 16, !51, i64 24, !57, i64 32}
!45 = !{!"_ZTSSt10unique_ptrIA_PP5cellbSt14default_deleteIS3_EE", !46, i64 0}
!46 = !{!"_ZTSSt15__uniq_ptr_dataIPP5cellbSt14default_deleteIA_S2_ELb1ELb1EE", !47, i64 0}
!47 = !{!"_ZTSSt15__uniq_ptr_implIPP5cellbSt14default_deleteIA_S2_EE", !48, i64 0}
!48 = !{!"_ZTSSt5tupleIJPPP5cellbSt14default_deleteIA_S2_EEE", !49, i64 0}
!49 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5cellbSt14default_deleteIA_S2_EEE", !50, i64 0}
!50 = !{!"_ZTSSt10_Head_baseILm0EPPP5cellbLb0EE", !6, i64 0}
!51 = !{!"_ZTSSt10unique_ptrIA_P5cellbSt14default_deleteIS2_EE", !52, i64 0}
!52 = !{!"_ZTSSt15__uniq_ptr_dataIP5cellbSt14default_deleteIA_S1_ELb1ELb1EE", !53, i64 0}
!53 = !{!"_ZTSSt15__uniq_ptr_implIP5cellbSt14default_deleteIA_S1_EE", !54, i64 0}
!54 = !{!"_ZTSSt5tupleIJPP5cellbSt14default_deleteIA_S1_EEE", !55, i64 0}
!55 = !{!"_ZTSSt11_Tuple_implILm0EJPP5cellbSt14default_deleteIA_S1_EEE", !56, i64 0}
!56 = !{!"_ZTSSt10_Head_baseILm0EPP5cellbLb0EE", !6, i64 0}
!57 = !{!"_ZTSSt10unique_ptrIA_5cellbSt14default_deleteIS1_EE", !58, i64 0}
!58 = !{!"_ZTSSt15__uniq_ptr_dataI5cellbSt14default_deleteIA_S0_ELb1ELb1EE", !59, i64 0}
!59 = !{!"_ZTSSt15__uniq_ptr_implI5cellbSt14default_deleteIA_S0_EE", !60, i64 0}
!60 = !{!"_ZTSSt5tupleIJP5cellbSt14default_deleteIA_S0_EEE", !61, i64 0}
!61 = !{!"_ZTSSt11_Tuple_implILm0EJP5cellbSt14default_deleteIA_S0_EEE", !62, i64 0}
!62 = !{!"_ZTSSt10_Head_baseILm0EP5cellbLb0EE", !6, i64 0}
!63 = !{!"_ZTS7field3dI5celljE", !13, i64 0, !13, i64 4, !13, i64 8, !64, i64 16, !70, i64 24, !76, i64 32}
!64 = !{!"_ZTSSt10unique_ptrIA_PP5celljSt14default_deleteIS3_EE", !65, i64 0}
!65 = !{!"_ZTSSt15__uniq_ptr_dataIPP5celljSt14default_deleteIA_S2_ELb1ELb1EE", !66, i64 0}
!66 = !{!"_ZTSSt15__uniq_ptr_implIPP5celljSt14default_deleteIA_S2_EE", !67, i64 0}
!67 = !{!"_ZTSSt5tupleIJPPP5celljSt14default_deleteIA_S2_EEE", !68, i64 0}
!68 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5celljSt14default_deleteIA_S2_EEE", !69, i64 0}
!69 = !{!"_ZTSSt10_Head_baseILm0EPPP5celljLb0EE", !6, i64 0}
!70 = !{!"_ZTSSt10unique_ptrIA_P5celljSt14default_deleteIS2_EE", !71, i64 0}
!71 = !{!"_ZTSSt15__uniq_ptr_dataIP5celljSt14default_deleteIA_S1_ELb1ELb1EE", !72, i64 0}
!72 = !{!"_ZTSSt15__uniq_ptr_implIP5celljSt14default_deleteIA_S1_EE", !73, i64 0}
!73 = !{!"_ZTSSt5tupleIJPP5celljSt14default_deleteIA_S1_EEE", !74, i64 0}
!74 = !{!"_ZTSSt11_Tuple_implILm0EJPP5celljSt14default_deleteIA_S1_EEE", !75, i64 0}
!75 = !{!"_ZTSSt10_Head_baseILm0EPP5celljLb0EE", !6, i64 0}
!76 = !{!"_ZTSSt10unique_ptrIA_5celljSt14default_deleteIS1_EE", !77, i64 0}
!77 = !{!"_ZTSSt15__uniq_ptr_dataI5celljSt14default_deleteIA_S0_ELb1ELb1EE", !78, i64 0}
!78 = !{!"_ZTSSt15__uniq_ptr_implI5celljSt14default_deleteIA_S0_EE", !79, i64 0}
!79 = !{!"_ZTSSt5tupleIJP5celljSt14default_deleteIA_S0_EEE", !80, i64 0}
!80 = !{!"_ZTSSt11_Tuple_implILm0EJP5celljSt14default_deleteIA_S0_EEE", !81, i64 0}
!81 = !{!"_ZTSSt10_Head_baseILm0EP5celljLb0EE", !6, i64 0}
!82 = !{!"_ZTS7field3dI6cellbeE", !13, i64 0, !13, i64 4, !13, i64 8, !83, i64 16, !89, i64 24, !95, i64 32}
!83 = !{!"_ZTSSt10unique_ptrIA_PP6cellbeSt14default_deleteIS3_EE", !84, i64 0}
!84 = !{!"_ZTSSt15__uniq_ptr_dataIPP6cellbeSt14default_deleteIA_S2_ELb1ELb1EE", !85, i64 0}
!85 = !{!"_ZTSSt15__uniq_ptr_implIPP6cellbeSt14default_deleteIA_S2_EE", !86, i64 0}
!86 = !{!"_ZTSSt5tupleIJPPP6cellbeSt14default_deleteIA_S2_EEE", !87, i64 0}
!87 = !{!"_ZTSSt11_Tuple_implILm0EJPPP6cellbeSt14default_deleteIA_S2_EEE", !88, i64 0}
!88 = !{!"_ZTSSt10_Head_baseILm0EPPP6cellbeLb0EE", !6, i64 0}
!89 = !{!"_ZTSSt10unique_ptrIA_P6cellbeSt14default_deleteIS2_EE", !90, i64 0}
!90 = !{!"_ZTSSt15__uniq_ptr_dataIP6cellbeSt14default_deleteIA_S1_ELb1ELb1EE", !91, i64 0}
!91 = !{!"_ZTSSt15__uniq_ptr_implIP6cellbeSt14default_deleteIA_S1_EE", !92, i64 0}
!92 = !{!"_ZTSSt5tupleIJPP6cellbeSt14default_deleteIA_S1_EEE", !93, i64 0}
!93 = !{!"_ZTSSt11_Tuple_implILm0EJPP6cellbeSt14default_deleteIA_S1_EEE", !94, i64 0}
!94 = !{!"_ZTSSt10_Head_baseILm0EPP6cellbeLb0EE", !6, i64 0}
!95 = !{!"_ZTSSt10unique_ptrIA_6cellbeSt14default_deleteIS1_EE", !96, i64 0}
!96 = !{!"_ZTSSt15__uniq_ptr_dataI6cellbeSt14default_deleteIA_S0_ELb1ELb1EE", !97, i64 0}
!97 = !{!"_ZTSSt15__uniq_ptr_implI6cellbeSt14default_deleteIA_S0_EE", !98, i64 0}
!98 = !{!"_ZTSSt5tupleIJP6cellbeSt14default_deleteIA_S0_EEE", !99, i64 0}
!99 = !{!"_ZTSSt11_Tuple_implILm0EJP6cellbeSt14default_deleteIA_S0_EEE", !100, i64 0}
!100 = !{!"_ZTSSt10_Head_baseILm0EP6cellbeLb0EE", !6, i64 0}
!101 = !{!"_ZTS7field3dI5cellpE", !13, i64 0, !13, i64 4, !13, i64 8, !102, i64 16, !108, i64 24, !114, i64 32}
!102 = !{!"_ZTSSt10unique_ptrIA_PP5cellpSt14default_deleteIS3_EE", !103, i64 0}
!103 = !{!"_ZTSSt15__uniq_ptr_dataIPP5cellpSt14default_deleteIA_S2_ELb1ELb1EE", !104, i64 0}
!104 = !{!"_ZTSSt15__uniq_ptr_implIPP5cellpSt14default_deleteIA_S2_EE", !105, i64 0}
!105 = !{!"_ZTSSt5tupleIJPPP5cellpSt14default_deleteIA_S2_EEE", !106, i64 0}
!106 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5cellpSt14default_deleteIA_S2_EEE", !107, i64 0}
!107 = !{!"_ZTSSt10_Head_baseILm0EPPP5cellpLb0EE", !6, i64 0}
!108 = !{!"_ZTSSt10unique_ptrIA_P5cellpSt14default_deleteIS2_EE", !109, i64 0}
!109 = !{!"_ZTSSt15__uniq_ptr_dataIP5cellpSt14default_deleteIA_S1_ELb1ELb1EE", !110, i64 0}
!110 = !{!"_ZTSSt15__uniq_ptr_implIP5cellpSt14default_deleteIA_S1_EE", !111, i64 0}
!111 = !{!"_ZTSSt5tupleIJPP5cellpSt14default_deleteIA_S1_EEE", !112, i64 0}
!112 = !{!"_ZTSSt11_Tuple_implILm0EJPP5cellpSt14default_deleteIA_S1_EEE", !113, i64 0}
!113 = !{!"_ZTSSt10_Head_baseILm0EPP5cellpLb0EE", !6, i64 0}
!114 = !{!"_ZTSSt10unique_ptrIA_5cellpSt14default_deleteIS1_EE", !115, i64 0}
!115 = !{!"_ZTSSt15__uniq_ptr_dataI5cellpSt14default_deleteIA_S0_ELb1ELb1EE", !116, i64 0}
!116 = !{!"_ZTSSt15__uniq_ptr_implI5cellpSt14default_deleteIA_S0_EE", !117, i64 0}
!117 = !{!"_ZTSSt5tupleIJP5cellpSt14default_deleteIA_S0_EEE", !118, i64 0}
!118 = !{!"_ZTSSt11_Tuple_implILm0EJP5cellpSt14default_deleteIA_S0_EEE", !119, i64 0}
!119 = !{!"_ZTSSt10_Head_baseILm0EP5cellpLb0EE", !6, i64 0}
!120 = !{!"_ZTSSt6vectorI7field3dIdESaIS1_EE", !121, i64 0}
!121 = !{!"_ZTSSt12_Vector_baseI7field3dIdESaIS1_EE", !122, i64 0}
!122 = !{!"_ZTSNSt12_Vector_baseI7field3dIdESaIS1_EE12_Vector_implE", !123, i64 0}
!123 = !{!"_ZTSNSt12_Vector_baseI7field3dIdESaIS1_EE17_Vector_impl_dataE", !6, i64 0, !6, i64 8, !6, i64 16}
!124 = !{!"_ZTSSt6vectorIN14spatial_region16deleted_particleESaIS1_EE", !125, i64 0}
!125 = !{!"_ZTSSt12_Vector_baseIN14spatial_region16deleted_particleESaIS1_EE", !126, i64 0}
!126 = !{!"_ZTSNSt12_Vector_baseIN14spatial_region16deleted_particleESaIS1_EE12_Vector_implE", !127, i64 0}
!127 = !{!"_ZTSNSt12_Vector_baseIN14spatial_region16deleted_particleESaIS1_EE17_Vector_impl_dataE", !6, i64 0, !6, i64 8, !6, i64 16}
!128 = !{!129, !13, i64 0}
!129 = !{!"_ZTS12int_vector3d", !13, i64 0, !13, i64 4, !13, i64 8}
!130 = !{!12, !13, i64 52}
!131 = !{!129, !13, i64 4}
!132 = !{!12, !13, i64 56}
!133 = !{!129, !13, i64 8}
!134 = !{!13, !13, i64 0}
!135 = !{!12, !13, i64 0}
!136 = !{!12, !13, i64 4}
!137 = !{!12, !13, i64 8}
!138 = !{!12, !10, i64 24}
!139 = !{!12, !10, i64 16}
!140 = !{!12, !10, i64 32}
!141 = !{!142, !10, i64 0}
!142 = !{!"_ZTS5celle", !10, i64 0, !10, i64 8, !10, i64 16}
!143 = !{!142, !10, i64 8}
!144 = !{!142, !10, i64 16}
!145 = !{!146, !10, i64 0}
!146 = !{!"_ZTS5cellb", !10, i64 0, !10, i64 8, !10, i64 16}
!147 = !{!146, !10, i64 8}
!148 = !{!146, !10, i64 16}
!149 = distinct !{!149, !150}
!150 = !{!"llvm.loop.mustprogress"}
!151 = distinct !{!151, !150}
!152 = distinct !{!152, !150}


; Optimized IR (opt -passes='default<O3>'):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

%struct.int_vector3d = type { i32, i32, i32 }
%struct.celle = type { double, double, double }
%struct.cellb = type { double, double, double }

; 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: mustprogress nofree nounwind willreturn memory(write)
declare double @cos(double noundef) local_unnamed_addr #2

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

; Function Attrs: inlinehint mustprogress uwtable
declare dso_local noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40), i32 noundef) local_unnamed_addr #3 align 2

; Function Attrs: inlinehint mustprogress uwtable
declare dso_local noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40), i32 noundef) local_unnamed_addr #3 align 2

; Function Attrs: mustprogress uwtable
declare dso_local void @_ZN14spatial_region22fill_cell_by_particlesEdR12int_vector3dS1_ddddd(ptr noundef nonnull align 8 dereferenceable(608), double noundef, ptr noundef nonnull align 4 dereferenceable(12), ptr noundef nonnull align 4 dereferenceable(12), double noundef, double noundef, double noundef, double noundef, double noundef) local_unnamed_addr #4 align 2

; Function Attrs: mustprogress uwtable
define dso_local void @_ZN14spatial_region8add_beamEdddddddd(ptr noundef nonnull align 8 dereferenceable(608) %0, double noundef %1, double noundef %2, double noundef %3, double noundef %4, double noundef %5, double noundef %6, double noundef %7, double noundef %8) local_unnamed_addr #4 align 2 {
  %10 = alloca %struct.int_vector3d, align 4
  %11 = alloca %struct.int_vector3d, align 8
  call void @llvm.lifetime.start.p0(i64 12, ptr nonnull %10) #7
  call void @_ZN12int_vector3dC2Ev(ptr noundef nonnull align 4 dereferenceable(12) %10)
  call void @llvm.lifetime.start.p0(i64 12, ptr nonnull %11) #7
  call void @_ZN12int_vector3dC2Ev(ptr noundef nonnull align 4 dereferenceable(12) %11)
  %12 = getelementptr inbounds i8, ptr %0, i64 48
  %13 = load <2 x i32>, ptr %12, align 8, !tbaa !5
  store <2 x i32> %13, ptr %11, align 8, !tbaa !5
  %14 = getelementptr inbounds i8, ptr %0, i64 56
  %15 = load i32, ptr %14, align 8, !tbaa !9
  %16 = getelementptr inbounds i8, ptr %11, i64 8
  store i32 %15, ptr %16, align 8, !tbaa !127
  %17 = fcmp ogt double %1, 0.000000e+00
  %18 = zext i1 %17 to i32
  %19 = fcmp olt double %1, 0.000000e+00
  %.neg181 = sext i1 %19 to i32
  %20 = add nsw i32 %.neg181, %18
  %21 = sitofp i32 %20 to double
  %22 = load i32, ptr %0, align 8, !tbaa !129
  %23 = icmp sgt i32 %22, 0
  br i1 %23, label %.preheader191.lr.ph, label %._crit_edge197

.preheader191.lr.ph:                              ; preds = %9
  %24 = fcmp olt double %3, 0.000000e+00
  %.neg = sext i1 %24 to i32
  %25 = fcmp ogt double %3, 0.000000e+00
  %26 = zext i1 %25 to i32
  %27 = add nsw i32 %.neg, %26
  %28 = sitofp i32 %27 to double
  %29 = getelementptr inbounds i8, ptr %0, i64 4
  %30 = getelementptr inbounds i8, ptr %0, i64 8
  %31 = getelementptr inbounds i8, ptr %10, i64 4
  %32 = getelementptr inbounds i8, ptr %10, i64 8
  %33 = getelementptr inbounds i8, ptr %0, i64 24
  %34 = fneg double %7
  %35 = getelementptr inbounds i8, ptr %0, i64 16
  %36 = fneg double %6
  %37 = getelementptr inbounds i8, ptr %0, i64 32
  %38 = fneg double %4
  %39 = insertelement <2 x double> poison, double %4, i64 0
  %40 = insertelement <2 x double> %39, double %5, i64 1
  %41 = fmul <2 x double> %40, %40
  %42 = fmul double %2, 5.000000e-01
  %43 = fmul double %5, 2.000000e+00
  %44 = fmul double %43, %5
  %45 = getelementptr inbounds i8, ptr %0, i64 360
  %46 = fmul double %21, %28
  %47 = getelementptr inbounds i8, ptr %0, i64 400
  %48 = fmul double %2, 2.500000e-01
  %49 = load i32, ptr %29, align 4, !tbaa !130
  %50 = icmp sgt i32 %49, 0
  br i1 %50, label %.preheader191.preheader, label %._crit_edge197

.preheader191.preheader:                          ; preds = %.preheader191.lr.ph
  %51 = extractelement <2 x double> %41, i64 0
  br label %.preheader191

.preheader191:                                    ; preds = %.preheader191.preheader, %._crit_edge194
  %52 = phi i32 [ %63, %._crit_edge194 ], [ %22, %.preheader191.preheader ]
  %53 = phi i32 [ %64, %._crit_edge194 ], [ %49, %.preheader191.preheader ]
  %54 = phi i32 [ %65, %._crit_edge194 ], [ %49, %.preheader191.preheader ]
  %.0174195 = phi i32 [ %66, %._crit_edge194 ], [ 0, %.preheader191.preheader ]
  %55 = icmp sgt i32 %54, 0
  br i1 %55, label %.preheader.lr.ph, label %._crit_edge194

.preheader.lr.ph:                                 ; preds = %.preheader191
  %56 = uitofp nneg i32 %.0174195 to double
  %57 = load i32, ptr %30, align 8, !tbaa !131
  %58 = icmp sgt i32 %57, 0
  br i1 %58, label %.preheader, label %._crit_edge194

._crit_edge197:                                   ; preds = %._crit_edge194, %.preheader191.lr.ph, %9
  call void @llvm.lifetime.end.p0(i64 12, ptr nonnull %11) #7
  call void @llvm.lifetime.end.p0(i64 12, ptr nonnull %10) #7
  ret void

.preheader:                                       ; preds = %.preheader.lr.ph, %._crit_edge
  %59 = phi i32 [ %68, %._crit_edge ], [ %53, %.preheader.lr.ph ]
  %60 = phi i32 [ %69, %._crit_edge ], [ %57, %.preheader.lr.ph ]
  %indvars.iv201 = phi i64 [ %indvars.iv.next202, %._crit_edge ], [ 0, %.preheader.lr.ph ]
  %61 = icmp sgt i32 %60, 0
  br i1 %61, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %.preheader
  %62 = trunc nuw nsw i64 %indvars.iv201 to i32
  br label %72

._crit_edge194.loopexit:                          ; preds = %._crit_edge
  %.pre204 = load i32, ptr %0, align 8, !tbaa !129
  br label %._crit_edge194

._crit_edge194:                                   ; preds = %.preheader.lr.ph, %._crit_edge194.loopexit, %.preheader191
  %63 = phi i32 [ %.pre204, %._crit_edge194.loopexit ], [ %52, %.preheader191 ], [ %52, %.preheader.lr.ph ]
  %64 = phi i32 [ %68, %._crit_edge194.loopexit ], [ %53, %.preheader191 ], [ %53, %.preheader.lr.ph ]
  %65 = phi i32 [ %68, %._crit_edge194.loopexit ], [ %54, %.preheader191 ], [ %54, %.preheader.lr.ph ]
  %66 = add nuw nsw i32 %.0174195, 1
  %67 = icmp slt i32 %66, %63
  br i1 %67, label %.preheader191, label %._crit_edge197, !llvm.loop !132

._crit_edge.loopexit:                             ; preds = %249
  %.pre = load i32, ptr %29, align 4, !tbaa !130
  br label %._crit_edge

._crit_edge:                                      ; preds = %._crit_edge.loopexit, %.preheader
  %68 = phi i32 [ %.pre, %._crit_edge.loopexit ], [ %59, %.preheader ]
  %69 = phi i32 [ %250, %._crit_edge.loopexit ], [ %60, %.preheader ]
  %indvars.iv.next202 = add nuw nsw i64 %indvars.iv201, 1
  %70 = sext i32 %68 to i64
  %71 = icmp slt i64 %indvars.iv.next202, %70
  br i1 %71, label %.preheader, label %._crit_edge194.loopexit, !llvm.loop !135

72:                                               ; preds = %.lr.ph, %249
  %indvars.iv = phi i64 [ 0, %.lr.ph ], [ %indvars.iv.next, %249 ]
  store i32 %.0174195, ptr %10, align 4, !tbaa !136
  store i32 %62, ptr %31, align 4, !tbaa !137
  %73 = trunc nuw nsw i64 %indvars.iv to i32
  store i32 %73, ptr %32, align 4, !tbaa !127
  %74 = load i32, ptr %29, align 4, !tbaa !130
  %.neg182 = sdiv i32 %74, -2
  %75 = add i32 %.neg182, %62
  %76 = sitofp i32 %75 to double
  %77 = load double, ptr %33, align 8, !tbaa !138
  %78 = call double @llvm.fmuladd.f64(double %76, double %77, double %34)
  %79 = load double, ptr %35, align 8, !tbaa !139
  %80 = call double @llvm.fmuladd.f64(double %56, double %79, double %36)
  %81 = call double @atan2(double noundef %78, double noundef %80) #7
  %82 = load double, ptr %35, align 8, !tbaa !139
  %83 = call double @llvm.fmuladd.f64(double %56, double %82, double %36)
  %84 = load i32, ptr %29, align 4, !tbaa !130
  %.neg183 = sdiv i32 %84, -2
  %85 = add i32 %.neg183, %62
  %86 = sitofp i32 %85 to double
  %87 = load double, ptr %33, align 8, !tbaa !138
  %88 = call double @llvm.fmuladd.f64(double %86, double %87, double %34)
  %89 = fmul double %88, %88
  %90 = call double @llvm.fmuladd.f64(double %83, double %83, double %89)
  %sqrt190 = call double @llvm.sqrt.f64(double %90)
  %91 = fsub double %81, %8
  %92 = call double @cos(double noundef %91) #7
  %93 = fmul double %92, %sqrt190
  %94 = load double, ptr %35, align 8, !tbaa !139
  %95 = call double @llvm.fmuladd.f64(double %56, double %94, double %36)
  %96 = load i32, ptr %29, align 4, !tbaa !130
  %.neg184 = sdiv i32 %96, -2
  %97 = add i32 %.neg184, %62
  %98 = sitofp i32 %97 to double
  %99 = load double, ptr %33, align 8, !tbaa !138
  %100 = call double @llvm.fmuladd.f64(double %98, double %99, double %34)
  %101 = fmul double %100, %100
  %102 = call double @llvm.fmuladd.f64(double %95, double %95, double %101)
  %sqrt189 = call double @llvm.sqrt.f64(double %102)
  %103 = call double @sin(double noundef %91) #7
  %104 = fmul double %103, %sqrt189
  %105 = load i32, ptr %30, align 8, !tbaa !131
  %.neg185 = sdiv i32 %105, -2
  %106 = add i32 %.neg185, %73
  %107 = sitofp i32 %106 to double
  %108 = load double, ptr %37, align 8, !tbaa !140
  %109 = fmul double %108, %107
  %110 = fmul double %109, %109
  %111 = call double @llvm.fmuladd.f64(double %104, double %104, double %110)
  %sqrt = call double @llvm.sqrt.f64(double %111)
  %112 = fcmp ogt double %93, %38
  %113 = fcmp olt double %93, %4
  %or.cond = and i1 %112, %113
  %114 = fcmp olt double %sqrt, %5
  %or.cond186 = select i1 %or.cond, i1 %114, i1 false
  br i1 %or.cond186, label %115, label %186

115:                                              ; preds = %72
  %116 = insertelement <2 x double> poison, double %93, i64 0
  %117 = insertelement <2 x double> %116, double %sqrt, i64 1
  %118 = fmul <2 x double> %117, %117
  %119 = fdiv <2 x double> %118, %41
  %120 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %119
  %121 = extractelement <2 x double> %120, i64 0
  %122 = fmul double %121, %2
  %123 = extractelement <2 x double> %120, i64 1
  %124 = fmul double %122, %123
  %125 = call double @cos(double noundef %8) #7
  %126 = fmul double %125, %3
  %127 = call double @sin(double noundef %8) #7
  %128 = fmul double %127, %3
  call void @_ZN14spatial_region22fill_cell_by_particlesEdR12int_vector3dS1_ddddd(ptr noundef nonnull align 8 dereferenceable(608) %0, double noundef %1, ptr noundef nonnull align 4 dereferenceable(12) %10, ptr noundef nonnull align 4 dereferenceable(12) %11, double noundef %124, double noundef %126, double noundef %128, double noundef 0.000000e+00, double noundef 0.000000e+00)
  %129 = fmul double %42, %121
  %130 = extractelement <2 x double> %118, i64 1
  %131 = fdiv double %130, %44
  %132 = fsub double 1.000000e+00, %131
  %133 = fmul double %129, %132
  %134 = fmul double %133, %21
  %135 = fmul double %104, %134
  %136 = call double @sin(double noundef %8) #7
  %137 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %45, i32 noundef %.0174195)
  %138 = load ptr, ptr %137, align 8, !tbaa !141
  %139 = getelementptr inbounds ptr, ptr %138, i64 %indvars.iv201
  %140 = load ptr, ptr %139, align 8, !tbaa !141
  %141 = getelementptr inbounds %struct.celle, ptr %140, i64 %indvars.iv
  %142 = load double, ptr %141, align 8, !tbaa !142
  %143 = fneg double %135
  %144 = call double @llvm.fmuladd.f64(double %143, double %136, double %142)
  store double %144, ptr %141, align 8, !tbaa !142
  %145 = call double @cos(double noundef %8) #7
  %146 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %45, i32 noundef %.0174195)
  %147 = load ptr, ptr %146, align 8, !tbaa !141
  %148 = getelementptr inbounds ptr, ptr %147, i64 %indvars.iv201
  %149 = load ptr, ptr %148, align 8, !tbaa !141
  %150 = getelementptr inbounds %struct.celle, ptr %149, i64 %indvars.iv, i32 1
  %151 = load double, ptr %150, align 8, !tbaa !144
  %152 = call double @llvm.fmuladd.f64(double %135, double %145, double %151)
  store double %152, ptr %150, align 8, !tbaa !144
  %153 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %45, i32 noundef %.0174195)
  %154 = load ptr, ptr %153, align 8, !tbaa !141
  %155 = getelementptr inbounds ptr, ptr %154, i64 %indvars.iv201
  %156 = load ptr, ptr %155, align 8, !tbaa !141
  %157 = getelementptr inbounds %struct.celle, ptr %156, i64 %indvars.iv, i32 2
  %158 = load double, ptr %157, align 8, !tbaa !145
  %159 = call double @llvm.fmuladd.f64(double %134, double %109, double %158)
  store double %159, ptr %157, align 8, !tbaa !145
  %160 = fmul double %46, %133
  %161 = fmul double %109, %160
  %162 = call double @sin(double noundef %8) #7
  %163 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %47, i32 noundef %.0174195)
  %164 = load ptr, ptr %163, align 8, !tbaa !141
  %165 = getelementptr inbounds ptr, ptr %164, i64 %indvars.iv201
  %166 = load ptr, ptr %165, align 8, !tbaa !141
  %167 = getelementptr inbounds %struct.cellb, ptr %166, i64 %indvars.iv
  %168 = load double, ptr %167, align 8, !tbaa !146
  %169 = call double @llvm.fmuladd.f64(double %161, double %162, double %168)
  store double %169, ptr %167, align 8, !tbaa !146
  %170 = call double @cos(double noundef %8) #7
  %171 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %47, i32 noundef %.0174195)
  %172 = load ptr, ptr %171, align 8, !tbaa !141
  %173 = getelementptr inbounds ptr, ptr %172, i64 %indvars.iv201
  %174 = load ptr, ptr %173, align 8, !tbaa !141
  %175 = getelementptr inbounds %struct.cellb, ptr %174, i64 %indvars.iv, i32 1
  %176 = load double, ptr %175, align 8, !tbaa !148
  %177 = fneg double %161
  %178 = call double @llvm.fmuladd.f64(double %177, double %170, double %176)
  store double %178, ptr %175, align 8, !tbaa !148
  %179 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %47, i32 noundef %.0174195)
  %180 = load ptr, ptr %179, align 8, !tbaa !141
  %181 = getelementptr inbounds ptr, ptr %180, i64 %indvars.iv201
  %182 = load ptr, ptr %181, align 8, !tbaa !141
  %183 = getelementptr inbounds %struct.cellb, ptr %182, i64 %indvars.iv, i32 2
  %184 = load double, ptr %183, align 8, !tbaa !149
  %185 = call double @llvm.fmuladd.f64(double %160, double %104, double %184)
  store double %185, ptr %183, align 8, !tbaa !149
  br label %186

186:                                              ; preds = %115, %72
  %187 = fcmp oge double %sqrt, %5
  %or.cond188.not = select i1 %or.cond, i1 %187, i1 false
  br i1 %or.cond188.not, label %188, label %249

188:                                              ; preds = %186
  %189 = fmul double %93, %93
  %190 = fdiv double %189, %51
  %191 = fsub double 1.000000e+00, %190
  %192 = fmul double %48, %191
  %193 = fmul double %192, %5
  %194 = fmul double %193, %5
  %195 = fmul double %sqrt, %sqrt
  %196 = fdiv double %194, %195
  %197 = fmul double %196, %21
  %198 = fmul double %104, %197
  %199 = call double @sin(double noundef %8) #7
  %200 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %45, i32 noundef %.0174195)
  %201 = load ptr, ptr %200, align 8, !tbaa !141
  %202 = getelementptr inbounds ptr, ptr %201, i64 %indvars.iv201
  %203 = load ptr, ptr %202, align 8, !tbaa !141
  %204 = getelementptr inbounds %struct.celle, ptr %203, i64 %indvars.iv
  %205 = load double, ptr %204, align 8, !tbaa !142
  %206 = fneg double %198
  %207 = call double @llvm.fmuladd.f64(double %206, double %199, double %205)
  store double %207, ptr %204, align 8, !tbaa !142
  %208 = call double @cos(double noundef %8) #7
  %209 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %45, i32 noundef %.0174195)
  %210 = load ptr, ptr %209, align 8, !tbaa !141
  %211 = getelementptr inbounds ptr, ptr %210, i64 %indvars.iv201
  %212 = load ptr, ptr %211, align 8, !tbaa !141
  %213 = getelementptr inbounds %struct.celle, ptr %212, i64 %indvars.iv, i32 1
  %214 = load double, ptr %213, align 8, !tbaa !144
  %215 = call double @llvm.fmuladd.f64(double %198, double %208, double %214)
  store double %215, ptr %213, align 8, !tbaa !144
  %216 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5celleEixEi(ptr noundef nonnull align 8 dereferenceable(40) %45, i32 noundef %.0174195)
  %217 = load ptr, ptr %216, align 8, !tbaa !141
  %218 = getelementptr inbounds ptr, ptr %217, i64 %indvars.iv201
  %219 = load ptr, ptr %218, align 8, !tbaa !141
  %220 = getelementptr inbounds %struct.celle, ptr %219, i64 %indvars.iv, i32 2
  %221 = load double, ptr %220, align 8, !tbaa !145
  %222 = call double @llvm.fmuladd.f64(double %197, double %109, double %221)
  store double %222, ptr %220, align 8, !tbaa !145
  %223 = fmul double %46, %196
  %224 = fmul double %109, %223
  %225 = call double @sin(double noundef %8) #7
  %226 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %47, i32 noundef %.0174195)
  %227 = load ptr, ptr %226, align 8, !tbaa !141
  %228 = getelementptr inbounds ptr, ptr %227, i64 %indvars.iv201
  %229 = load ptr, ptr %228, align 8, !tbaa !141
  %230 = getelementptr inbounds %struct.cellb, ptr %229, i64 %indvars.iv
  %231 = load double, ptr %230, align 8, !tbaa !146
  %232 = call double @llvm.fmuladd.f64(double %224, double %225, double %231)
  store double %232, ptr %230, align 8, !tbaa !146
  %233 = call double @cos(double noundef %8) #7
  %234 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %47, i32 noundef %.0174195)
  %235 = load ptr, ptr %234, align 8, !tbaa !141
  %236 = getelementptr inbounds ptr, ptr %235, i64 %indvars.iv201
  %237 = load ptr, ptr %236, align 8, !tbaa !141
  %238 = getelementptr inbounds %struct.cellb, ptr %237, i64 %indvars.iv, i32 1
  %239 = load double, ptr %238, align 8, !tbaa !148
  %240 = fneg double %224
  %241 = call double @llvm.fmuladd.f64(double %240, double %233, double %239)
  store double %241, ptr %238, align 8, !tbaa !148
  %242 = call noundef nonnull align 8 dereferenceable(8) ptr @_ZNK7field3dI5cellbEixEi(ptr noundef nonnull align 8 dereferenceable(40) %47, i32 noundef %.0174195)
  %243 = load ptr, ptr %242, align 8, !tbaa !141
  %244 = getelementptr inbounds ptr, ptr %243, i64 %indvars.iv201
  %245 = load ptr, ptr %244, align 8, !tbaa !141
  %246 = getelementptr inbounds %struct.cellb, ptr %245, i64 %indvars.iv, i32 2
  %247 = load double, ptr %246, align 8, !tbaa !149
  %248 = call double @llvm.fmuladd.f64(double %223, double %104, double %247)
  store double %248, ptr %246, align 8, !tbaa !149
  br label %249

249:                                              ; preds = %186, %188
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %250 = load i32, ptr %30, align 8, !tbaa !131
  %251 = sext i32 %250 to i64
  %252 = icmp slt i64 %indvars.iv.next, %251
  br i1 %252, label %72, label %._crit_edge.loopexit, !llvm.loop !150
}

; Function Attrs: mustprogress nounwind uwtable
declare dso_local void @_ZN12int_vector3dC2Ev(ptr noundef nonnull align 4 dereferenceable(12)) unnamed_addr #5 align 2

; Function Attrs: mustprogress nofree nounwind willreturn memory(write)
declare double @atan2(double noundef, double noundef) local_unnamed_addr #2

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

attributes #0 = { mustprogress nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #2 = { mustprogress nofree nounwind willreturn memory(write) "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { inlinehint 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 #4 = { 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 #5 = { 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 #6 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #7 = { 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 = !{!"int", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C++ TBAA"}
!9 = !{!10, !6, i64 56}
!10 = !{!"_ZTS14spatial_region", !6, i64 0, !6, i64 4, !6, i64 8, !11, i64 16, !11, i64 24, !11, i64 32, !11, i64 40, !6, i64 48, !6, i64 52, !6, i64 56, !11, i64 64, !12, i64 72, !6, i64 80, !6, i64 84, !6, i64 88, !12, i64 96, !6, i64 104, !6, i64 108, !6, i64 112, !6, i64 116, !6, i64 120, !6, i64 124, !6, i64 128, !13, i64 136, !19, i64 144, !11, i64 176, !11, i64 184, !11, i64 192, !6, i64 200, !6, i64 204, !6, i64 208, !12, i64 216, !11, i64 224, !11, i64 232, !11, i64 240, !11, i64 248, !12, i64 256, !11, i64 264, !11, i64 272, !11, i64 280, !12, i64 288, !21, i64 296, !12, i64 328, !12, i64 336, !12, i64 344, !12, i64 352, !24, i64 360, !43, i64 400, !62, i64 440, !81, i64 480, !100, i64 520, !119, i64 560, !123, i64 584}
!11 = !{!"double", !7, i64 0}
!12 = !{!"any pointer", !7, i64 0}
!13 = !{!"_ZTSSt10unique_ptrI14maxwell_solverSt14default_deleteIS0_EE", !14, i64 0}
!14 = !{!"_ZTSSt15__uniq_ptr_dataI14maxwell_solverSt14default_deleteIS0_ELb1ELb1EE", !15, i64 0}
!15 = !{!"_ZTSSt15__uniq_ptr_implI14maxwell_solverSt14default_deleteIS0_EE", !16, i64 0}
!16 = !{!"_ZTSSt5tupleIJP14maxwell_solverSt14default_deleteIS0_EEE", !17, i64 0}
!17 = !{!"_ZTSSt11_Tuple_implILm0EJP14maxwell_solverSt14default_deleteIS0_EEE", !18, i64 0}
!18 = !{!"_ZTSSt10_Head_baseILm0EP14maxwell_solverLb0EE", !12, i64 0}
!19 = !{!"_ZTSSt8functionIFvR8particleR8vector3dS3_RdEE", !20, i64 0, !12, i64 24}
!20 = !{!"_ZTSSt14_Function_base", !7, i64 0, !12, i64 16}
!21 = !{!"_ZTSNSt7__cxx1112basic_stringIcSt11char_traitsIcESaIcEEE", !22, i64 0, !23, i64 8, !7, i64 16}
!22 = !{!"_ZTSNSt7__cxx1112basic_stringIcSt11char_traitsIcESaIcEE12_Alloc_hiderE", !12, i64 0}
!23 = !{!"long", !7, i64 0}
!24 = !{!"_ZTS7field3dI5celleE", !6, i64 0, !6, i64 4, !6, i64 8, !25, i64 16, !31, i64 24, !37, i64 32}
!25 = !{!"_ZTSSt10unique_ptrIA_PP5celleSt14default_deleteIS3_EE", !26, i64 0}
!26 = !{!"_ZTSSt15__uniq_ptr_dataIPP5celleSt14default_deleteIA_S2_ELb1ELb1EE", !27, i64 0}
!27 = !{!"_ZTSSt15__uniq_ptr_implIPP5celleSt14default_deleteIA_S2_EE", !28, i64 0}
!28 = !{!"_ZTSSt5tupleIJPPP5celleSt14default_deleteIA_S2_EEE", !29, i64 0}
!29 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5celleSt14default_deleteIA_S2_EEE", !30, i64 0}
!30 = !{!"_ZTSSt10_Head_baseILm0EPPP5celleLb0EE", !12, i64 0}
!31 = !{!"_ZTSSt10unique_ptrIA_P5celleSt14default_deleteIS2_EE", !32, i64 0}
!32 = !{!"_ZTSSt15__uniq_ptr_dataIP5celleSt14default_deleteIA_S1_ELb1ELb1EE", !33, i64 0}
!33 = !{!"_ZTSSt15__uniq_ptr_implIP5celleSt14default_deleteIA_S1_EE", !34, i64 0}
!34 = !{!"_ZTSSt5tupleIJPP5celleSt14default_deleteIA_S1_EEE", !35, i64 0}
!35 = !{!"_ZTSSt11_Tuple_implILm0EJPP5celleSt14default_deleteIA_S1_EEE", !36, i64 0}
!36 = !{!"_ZTSSt10_Head_baseILm0EPP5celleLb0EE", !12, i64 0}
!37 = !{!"_ZTSSt10unique_ptrIA_5celleSt14default_deleteIS1_EE", !38, i64 0}
!38 = !{!"_ZTSSt15__uniq_ptr_dataI5celleSt14default_deleteIA_S0_ELb1ELb1EE", !39, i64 0}
!39 = !{!"_ZTSSt15__uniq_ptr_implI5celleSt14default_deleteIA_S0_EE", !40, i64 0}
!40 = !{!"_ZTSSt5tupleIJP5celleSt14default_deleteIA_S0_EEE", !41, i64 0}
!41 = !{!"_ZTSSt11_Tuple_implILm0EJP5celleSt14default_deleteIA_S0_EEE", !42, i64 0}
!42 = !{!"_ZTSSt10_Head_baseILm0EP5celleLb0EE", !12, i64 0}
!43 = !{!"_ZTS7field3dI5cellbE", !6, i64 0, !6, i64 4, !6, i64 8, !44, i64 16, !50, i64 24, !56, i64 32}
!44 = !{!"_ZTSSt10unique_ptrIA_PP5cellbSt14default_deleteIS3_EE", !45, i64 0}
!45 = !{!"_ZTSSt15__uniq_ptr_dataIPP5cellbSt14default_deleteIA_S2_ELb1ELb1EE", !46, i64 0}
!46 = !{!"_ZTSSt15__uniq_ptr_implIPP5cellbSt14default_deleteIA_S2_EE", !47, i64 0}
!47 = !{!"_ZTSSt5tupleIJPPP5cellbSt14default_deleteIA_S2_EEE", !48, i64 0}
!48 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5cellbSt14default_deleteIA_S2_EEE", !49, i64 0}
!49 = !{!"_ZTSSt10_Head_baseILm0EPPP5cellbLb0EE", !12, i64 0}
!50 = !{!"_ZTSSt10unique_ptrIA_P5cellbSt14default_deleteIS2_EE", !51, i64 0}
!51 = !{!"_ZTSSt15__uniq_ptr_dataIP5cellbSt14default_deleteIA_S1_ELb1ELb1EE", !52, i64 0}
!52 = !{!"_ZTSSt15__uniq_ptr_implIP5cellbSt14default_deleteIA_S1_EE", !53, i64 0}
!53 = !{!"_ZTSSt5tupleIJPP5cellbSt14default_deleteIA_S1_EEE", !54, i64 0}
!54 = !{!"_ZTSSt11_Tuple_implILm0EJPP5cellbSt14default_deleteIA_S1_EEE", !55, i64 0}
!55 = !{!"_ZTSSt10_Head_baseILm0EPP5cellbLb0EE", !12, i64 0}
!56 = !{!"_ZTSSt10unique_ptrIA_5cellbSt14default_deleteIS1_EE", !57, i64 0}
!57 = !{!"_ZTSSt15__uniq_ptr_dataI5cellbSt14default_deleteIA_S0_ELb1ELb1EE", !58, i64 0}
!58 = !{!"_ZTSSt15__uniq_ptr_implI5cellbSt14default_deleteIA_S0_EE", !59, i64 0}
!59 = !{!"_ZTSSt5tupleIJP5cellbSt14default_deleteIA_S0_EEE", !60, i64 0}
!60 = !{!"_ZTSSt11_Tuple_implILm0EJP5cellbSt14default_deleteIA_S0_EEE", !61, i64 0}
!61 = !{!"_ZTSSt10_Head_baseILm0EP5cellbLb0EE", !12, i64 0}
!62 = !{!"_ZTS7field3dI5celljE", !6, i64 0, !6, i64 4, !6, i64 8, !63, i64 16, !69, i64 24, !75, i64 32}
!63 = !{!"_ZTSSt10unique_ptrIA_PP5celljSt14default_deleteIS3_EE", !64, i64 0}
!64 = !{!"_ZTSSt15__uniq_ptr_dataIPP5celljSt14default_deleteIA_S2_ELb1ELb1EE", !65, i64 0}
!65 = !{!"_ZTSSt15__uniq_ptr_implIPP5celljSt14default_deleteIA_S2_EE", !66, i64 0}
!66 = !{!"_ZTSSt5tupleIJPPP5celljSt14default_deleteIA_S2_EEE", !67, i64 0}
!67 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5celljSt14default_deleteIA_S2_EEE", !68, i64 0}
!68 = !{!"_ZTSSt10_Head_baseILm0EPPP5celljLb0EE", !12, i64 0}
!69 = !{!"_ZTSSt10unique_ptrIA_P5celljSt14default_deleteIS2_EE", !70, i64 0}
!70 = !{!"_ZTSSt15__uniq_ptr_dataIP5celljSt14default_deleteIA_S1_ELb1ELb1EE", !71, i64 0}
!71 = !{!"_ZTSSt15__uniq_ptr_implIP5celljSt14default_deleteIA_S1_EE", !72, i64 0}
!72 = !{!"_ZTSSt5tupleIJPP5celljSt14default_deleteIA_S1_EEE", !73, i64 0}
!73 = !{!"_ZTSSt11_Tuple_implILm0EJPP5celljSt14default_deleteIA_S1_EEE", !74, i64 0}
!74 = !{!"_ZTSSt10_Head_baseILm0EPP5celljLb0EE", !12, i64 0}
!75 = !{!"_ZTSSt10unique_ptrIA_5celljSt14default_deleteIS1_EE", !76, i64 0}
!76 = !{!"_ZTSSt15__uniq_ptr_dataI5celljSt14default_deleteIA_S0_ELb1ELb1EE", !77, i64 0}
!77 = !{!"_ZTSSt15__uniq_ptr_implI5celljSt14default_deleteIA_S0_EE", !78, i64 0}
!78 = !{!"_ZTSSt5tupleIJP5celljSt14default_deleteIA_S0_EEE", !79, i64 0}
!79 = !{!"_ZTSSt11_Tuple_implILm0EJP5celljSt14default_deleteIA_S0_EEE", !80, i64 0}
!80 = !{!"_ZTSSt10_Head_baseILm0EP5celljLb0EE", !12, i64 0}
!81 = !{!"_ZTS7field3dI6cellbeE", !6, i64 0, !6, i64 4, !6, i64 8, !82, i64 16, !88, i64 24, !94, i64 32}
!82 = !{!"_ZTSSt10unique_ptrIA_PP6cellbeSt14default_deleteIS3_EE", !83, i64 0}
!83 = !{!"_ZTSSt15__uniq_ptr_dataIPP6cellbeSt14default_deleteIA_S2_ELb1ELb1EE", !84, i64 0}
!84 = !{!"_ZTSSt15__uniq_ptr_implIPP6cellbeSt14default_deleteIA_S2_EE", !85, i64 0}
!85 = !{!"_ZTSSt5tupleIJPPP6cellbeSt14default_deleteIA_S2_EEE", !86, i64 0}
!86 = !{!"_ZTSSt11_Tuple_implILm0EJPPP6cellbeSt14default_deleteIA_S2_EEE", !87, i64 0}
!87 = !{!"_ZTSSt10_Head_baseILm0EPPP6cellbeLb0EE", !12, i64 0}
!88 = !{!"_ZTSSt10unique_ptrIA_P6cellbeSt14default_deleteIS2_EE", !89, i64 0}
!89 = !{!"_ZTSSt15__uniq_ptr_dataIP6cellbeSt14default_deleteIA_S1_ELb1ELb1EE", !90, i64 0}
!90 = !{!"_ZTSSt15__uniq_ptr_implIP6cellbeSt14default_deleteIA_S1_EE", !91, i64 0}
!91 = !{!"_ZTSSt5tupleIJPP6cellbeSt14default_deleteIA_S1_EEE", !92, i64 0}
!92 = !{!"_ZTSSt11_Tuple_implILm0EJPP6cellbeSt14default_deleteIA_S1_EEE", !93, i64 0}
!93 = !{!"_ZTSSt10_Head_baseILm0EPP6cellbeLb0EE", !12, i64 0}
!94 = !{!"_ZTSSt10unique_ptrIA_6cellbeSt14default_deleteIS1_EE", !95, i64 0}
!95 = !{!"_ZTSSt15__uniq_ptr_dataI6cellbeSt14default_deleteIA_S0_ELb1ELb1EE", !96, i64 0}
!96 = !{!"_ZTSSt15__uniq_ptr_implI6cellbeSt14default_deleteIA_S0_EE", !97, i64 0}
!97 = !{!"_ZTSSt5tupleIJP6cellbeSt14default_deleteIA_S0_EEE", !98, i64 0}
!98 = !{!"_ZTSSt11_Tuple_implILm0EJP6cellbeSt14default_deleteIA_S0_EEE", !99, i64 0}
!99 = !{!"_ZTSSt10_Head_baseILm0EP6cellbeLb0EE", !12, i64 0}
!100 = !{!"_ZTS7field3dI5cellpE", !6, i64 0, !6, i64 4, !6, i64 8, !101, i64 16, !107, i64 24, !113, i64 32}
!101 = !{!"_ZTSSt10unique_ptrIA_PP5cellpSt14default_deleteIS3_EE", !102, i64 0}
!102 = !{!"_ZTSSt15__uniq_ptr_dataIPP5cellpSt14default_deleteIA_S2_ELb1ELb1EE", !103, i64 0}
!103 = !{!"_ZTSSt15__uniq_ptr_implIPP5cellpSt14default_deleteIA_S2_EE", !104, i64 0}
!104 = !{!"_ZTSSt5tupleIJPPP5cellpSt14default_deleteIA_S2_EEE", !105, i64 0}
!105 = !{!"_ZTSSt11_Tuple_implILm0EJPPP5cellpSt14default_deleteIA_S2_EEE", !106, i64 0}
!106 = !{!"_ZTSSt10_Head_baseILm0EPPP5cellpLb0EE", !12, i64 0}
!107 = !{!"_ZTSSt10unique_ptrIA_P5cellpSt14default_deleteIS2_EE", !108, i64 0}
!108 = !{!"_ZTSSt15__uniq_ptr_dataIP5cellpSt14default_deleteIA_S1_ELb1ELb1EE", !109, i64 0}
!109 = !{!"_ZTSSt15__uniq_ptr_implIP5cellpSt14default_deleteIA_S1_EE", !110, i64 0}
!110 = !{!"_ZTSSt5tupleIJPP5cellpSt14default_deleteIA_S1_EEE", !111, i64 0}
!111 = !{!"_ZTSSt11_Tuple_implILm0EJPP5cellpSt14default_deleteIA_S1_EEE", !112, i64 0}
!112 = !{!"_ZTSSt10_Head_baseILm0EPP5cellpLb0EE", !12, i64 0}
!113 = !{!"_ZTSSt10unique_ptrIA_5cellpSt14default_deleteIS1_EE", !114, i64 0}
!114 = !{!"_ZTSSt15__uniq_ptr_dataI5cellpSt14default_deleteIA_S0_ELb1ELb1EE", !115, i64 0}
!115 = !{!"_ZTSSt15__uniq_ptr_implI5cellpSt14default_deleteIA_S0_EE", !116, i64 0}
!116 = !{!"_ZTSSt5tupleIJP5cellpSt14default_deleteIA_S0_EEE", !117, i64 0}
!117 = !{!"_ZTSSt11_Tuple_implILm0EJP5cellpSt14default_deleteIA_S0_EEE", !118, i64 0}
!118 = !{!"_ZTSSt10_Head_baseILm0EP5cellpLb0EE", !12, i64 0}
!119 = !{!"_ZTSSt6vectorI7field3dIdESaIS1_EE", !120, i64 0}
!120 = !{!"_ZTSSt12_Vector_baseI7field3dIdESaIS1_EE", !121, i64 0}
!121 = !{!"_ZTSNSt12_Vector_baseI7field3dIdESaIS1_EE12_Vector_implE", !122, i64 0}
!122 = !{!"_ZTSNSt12_Vector_baseI7field3dIdESaIS1_EE17_Vector_impl_dataE", !12, i64 0, !12, i64 8, !12, i64 16}
!123 = !{!"_ZTSSt6vectorIN14spatial_region16deleted_particleESaIS1_EE", !124, i64 0}
!124 = !{!"_ZTSSt12_Vector_baseIN14spatial_region16deleted_particleESaIS1_EE", !125, i64 0}
!125 = !{!"_ZTSNSt12_Vector_baseIN14spatial_region16deleted_particleESaIS1_EE12_Vector_implE", !126, i64 0}
!126 = !{!"_ZTSNSt12_Vector_baseIN14spatial_region16deleted_particleESaIS1_EE17_Vector_impl_dataE", !12, i64 0, !12, i64 8, !12, i64 16}
!127 = !{!128, !6, i64 8}
!128 = !{!"_ZTS12int_vector3d", !6, i64 0, !6, i64 4, !6, i64 8}
!129 = !{!10, !6, i64 0}
!130 = !{!10, !6, i64 4}
!131 = !{!10, !6, i64 8}
!132 = distinct !{!132, !133, !134}
!133 = !{!"llvm.loop.mustprogress"}
!134 = !{!"llvm.loop.unswitch.partial.disable"}
!135 = distinct !{!135, !133, !134}
!136 = !{!128, !6, i64 0}
!137 = !{!128, !6, i64 4}
!138 = !{!10, !11, i64 24}
!139 = !{!10, !11, i64 16}
!140 = !{!10, !11, i64 32}
!141 = !{!12, !12, i64 0}
!142 = !{!143, !11, i64 0}
!143 = !{!"_ZTS5celle", !11, i64 0, !11, i64 8, !11, i64 16}
!144 = !{!143, !11, i64 8}
!145 = !{!143, !11, i64 16}
!146 = !{!147, !11, i64 0}
!147 = !{!"_ZTS5cellb", !11, i64 0, !11, i64 8, !11, i64 16}
!148 = !{!147, !11, i64 8}
!149 = !{!147, !11, i64 16}
!150 = distinct !{!150, !133}


; Effective O3 passes: AggressiveInstCombinePass,EarlyCSEPass,GVNPass,GlobalOptPass,IndVarSimplifyPass,InferAlignmentPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,JumpThreadingPass,LCSSAPass,LICMPass,LoopInstSimplifyPass,LoopRotatePass,LoopSimplifyCFGPass,LoopSimplifyPass,LoopUnrollPass,ReassociatePass,SLPVectorizerPass,SROAPass,SimpleLoopUnswitchPass,SimplifyCFGPass