; LLVM IR from QUILL-PIC/Quill :: src/mc.ll (ZN14spatial_region1wERKdRdS1_)

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

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

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

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

; Function Attrs: nounwind
declare double @sqrt(double noundef) #2

; Function Attrs: mustprogress nounwind uwtable
define dso_local noundef double @_ZN14spatial_region1wERKdRdS1_(ptr noundef nonnull align 8 dereferenceable(608) %0, ptr noundef nonnull align 8 dereferenceable(8) %1, ptr noundef nonnull align 8 dereferenceable(8) %2, ptr noundef nonnull align 8 dereferenceable(8) %3) #3 align 2 {
  %5 = alloca double, align 8
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca double, 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
  store ptr %0, ptr %6, align 8, !tbaa !5
  store ptr %1, ptr %7, align 8, !tbaa !5
  store ptr %2, ptr %8, align 8, !tbaa !5
  store ptr %3, ptr %9, align 8, !tbaa !5
  %17 = load ptr, ptr %6, align 8
  %18 = load ptr, ptr %8, align 8, !tbaa !5
  %19 = load double, ptr %18, align 8, !tbaa !9
  %20 = fsub double 1.000000e+00, %19
  %21 = load ptr, ptr %7, align 8, !tbaa !5
  %22 = load double, ptr %21, align 8, !tbaa !9
  %23 = fmul double %20, %22
  %24 = fcmp olt double %23, 1.000000e+00
  br i1 %24, label %33, label %25

25:                                               ; preds = %4
  %26 = load ptr, ptr %8, align 8, !tbaa !5
  %27 = load double, ptr %26, align 8, !tbaa !9
  %28 = fcmp oeq double %27, 0.000000e+00
  br i1 %28, label %33, label %29

29:                                               ; preds = %25
  %30 = load ptr, ptr %9, align 8, !tbaa !5
  %31 = load double, ptr %30, align 8, !tbaa !9
  %32 = fcmp oeq double %31, 0.000000e+00
  br i1 %32, label %33, label %34

33:                                               ; preds = %29, %25, %4
  store double 0.000000e+00, ptr %5, align 8
  br label %194

34:                                               ; preds = %29
  %35 = load ptr, ptr %9, align 8, !tbaa !5
  %36 = load double, ptr %35, align 8, !tbaa !9
  %37 = fcmp ogt double %36, 1.330000e-01
  br i1 %37, label %38, label %109

38:                                               ; preds = %34
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #4
  %39 = load ptr, ptr %8, align 8, !tbaa !5
  %40 = load double, ptr %39, align 8, !tbaa !9
  %41 = load ptr, ptr %8, align 8, !tbaa !5
  %42 = load double, ptr %41, align 8, !tbaa !9
  %43 = fsub double 1.000000e+00, %42
  %44 = load ptr, ptr %9, align 8, !tbaa !5
  %45 = load double, ptr %44, align 8, !tbaa !9
  %46 = fmul double %43, %45
  %47 = fdiv double %40, %46
  store double %47, ptr %10, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %11) #4
  %48 = load double, ptr %10, align 8, !tbaa !9
  %49 = call double @pow(double noundef %48, double noundef 0x3FD5555555555555) #4
  store double %49, ptr %11, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #4
  %50 = load double, ptr %11, align 8, !tbaa !9
  %51 = load double, ptr %11, align 8, !tbaa !9
  %52 = fmul double %50, %51
  store double %52, ptr %12, align 8, !tbaa !9
  %53 = getelementptr inbounds %class.spatial_region, ptr %17, i32 0, i32 10
  %54 = load double, ptr %53, align 8, !tbaa !11
  %55 = fmul double -2.059100e-03, %54
  %56 = load ptr, ptr %7, align 8, !tbaa !5
  %57 = load double, ptr %56, align 8, !tbaa !9
  %58 = load ptr, ptr %7, align 8, !tbaa !5
  %59 = load double, ptr %58, align 8, !tbaa !9
  %60 = fmul double %57, %59
  %61 = fdiv double %55, %60
  %62 = load double, ptr %10, align 8, !tbaa !9
  %63 = fmul double -2.000000e+00, %62
  %64 = fdiv double %63, 3.000000e+00
  %65 = call double @exp(double noundef %64) #4
  %66 = fmul double %61, %65
  %67 = load double, ptr %12, align 8, !tbaa !9
  %68 = fadd double %67, 8.004900e-01
  %69 = load double, ptr %12, align 8, !tbaa !9
  %70 = fadd double %69, 8.004900e-01
  %71 = fmul double %68, %70
  %72 = load double, ptr %12, align 8, !tbaa !9
  %73 = fadd double %72, 8.004900e-01
  %74 = fmul double %71, %73
  %75 = call double @sqrt(double noundef %74) #4
  %76 = call double @sqrt(double noundef %75) #4
  %77 = fdiv double 1.000000e+00, %76
  %78 = load ptr, ptr %8, align 8, !tbaa !5
  %79 = load double, ptr %78, align 8, !tbaa !9
  %80 = fsub double 1.000000e+00, %79
  %81 = load ptr, ptr %8, align 8, !tbaa !5
  %82 = load double, ptr %81, align 8, !tbaa !9
  %83 = load ptr, ptr %8, align 8, !tbaa !5
  %84 = load double, ptr %83, align 8, !tbaa !9
  %85 = fmul double %82, %84
  %86 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %80, double %85)
  %87 = load ptr, ptr %9, align 8, !tbaa !5
  %88 = load double, ptr %87, align 8, !tbaa !9
  %89 = fmul double %86, %88
  %90 = load double, ptr %11, align 8, !tbaa !9
  %91 = fmul double %89, %90
  %92 = load ptr, ptr %8, align 8, !tbaa !5
  %93 = load double, ptr %92, align 8, !tbaa !9
  %94 = fdiv double %91, %93
  %95 = load double, ptr %12, align 8, !tbaa !9
  %96 = load double, ptr %12, align 8, !tbaa !9
  %97 = fmul double 6.500000e-01, %96
  %98 = load double, ptr %12, align 8, !tbaa !9
  %99 = load double, ptr %12, align 8, !tbaa !9
  %100 = call double @llvm.fmuladd.f64(double %98, double %99, double 1.000000e+00)
  %101 = fdiv double %97, %100
  %102 = fsub double 1.000000e+00, %101
  %103 = call double @llvm.fmuladd.f64(double 0x3FE6ACEEE0F3CB3E, double %102, double %95)
  %104 = call double @sqrt(double noundef %103) #4
  %105 = call double @sqrt(double noundef %104) #4
  %106 = fneg double %94
  %107 = call double @llvm.fmuladd.f64(double %106, double %105, double %77)
  %108 = fmul double %66, %107
  store double %108, ptr %5, align 8
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %11) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #4
  br label %194

109:                                              ; preds = %34
  call void @llvm.lifetime.start.p0(i64 8, ptr %13) #4
  %110 = load ptr, ptr %9, align 8, !tbaa !5
  %111 = load double, ptr %110, align 8, !tbaa !9
  %112 = fmul double 1.500000e+01, %111
  %113 = fdiv double 2.000000e+00, %112
  %114 = fadd double 1.000000e+00, %113
  %115 = fdiv double 1.000000e+00, %114
  store double %115, ptr %13, align 8, !tbaa !9
  %116 = load double, ptr %13, align 8, !tbaa !9
  %117 = load ptr, ptr %8, align 8, !tbaa !5
  %118 = load double, ptr %117, align 8, !tbaa !9
  %119 = fmul double %116, %118
  %120 = load ptr, ptr %8, align 8, !tbaa !5
  store double %119, ptr %120, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #4
  %121 = load ptr, ptr %8, align 8, !tbaa !5
  %122 = load double, ptr %121, align 8, !tbaa !9
  %123 = load ptr, ptr %8, align 8, !tbaa !5
  %124 = load double, ptr %123, align 8, !tbaa !9
  %125 = fsub double 1.000000e+00, %124
  %126 = load ptr, ptr %9, align 8, !tbaa !5
  %127 = load double, ptr %126, align 8, !tbaa !9
  %128 = fmul double %125, %127
  %129 = fdiv double %122, %128
  store double %129, ptr %14, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #4
  %130 = load double, ptr %14, align 8, !tbaa !9
  %131 = call double @pow(double noundef %130, double noundef 0x3FD5555555555555) #4
  store double %131, ptr %15, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #4
  %132 = load double, ptr %15, align 8, !tbaa !9
  %133 = load double, ptr %15, align 8, !tbaa !9
  %134 = fmul double %132, %133
  store double %134, ptr %16, align 8, !tbaa !9
  %135 = load double, ptr %13, align 8, !tbaa !9
  %136 = fneg double %135
  %137 = fmul double %136, 2.059100e-03
  %138 = getelementptr inbounds %class.spatial_region, ptr %17, i32 0, i32 10
  %139 = load double, ptr %138, align 8, !tbaa !11
  %140 = fmul double %137, %139
  %141 = load ptr, ptr %7, align 8, !tbaa !5
  %142 = load double, ptr %141, align 8, !tbaa !9
  %143 = load ptr, ptr %7, align 8, !tbaa !5
  %144 = load double, ptr %143, align 8, !tbaa !9
  %145 = fmul double %142, %144
  %146 = fdiv double %140, %145
  %147 = load double, ptr %14, align 8, !tbaa !9
  %148 = fmul double -2.000000e+00, %147
  %149 = fdiv double %148, 3.000000e+00
  %150 = call double @exp(double noundef %149) #4
  %151 = fmul double %146, %150
  %152 = load double, ptr %16, align 8, !tbaa !9
  %153 = fadd double %152, 8.004900e-01
  %154 = load double, ptr %16, align 8, !tbaa !9
  %155 = fadd double %154, 8.004900e-01
  %156 = fmul double %153, %155
  %157 = load double, ptr %16, align 8, !tbaa !9
  %158 = fadd double %157, 8.004900e-01
  %159 = fmul double %156, %158
  %160 = call double @sqrt(double noundef %159) #4
  %161 = call double @sqrt(double noundef %160) #4
  %162 = fdiv double 1.000000e+00, %161
  %163 = load ptr, ptr %8, align 8, !tbaa !5
  %164 = load double, ptr %163, align 8, !tbaa !9
  %165 = fsub double 1.000000e+00, %164
  %166 = load ptr, ptr %8, align 8, !tbaa !5
  %167 = load double, ptr %166, align 8, !tbaa !9
  %168 = load ptr, ptr %8, align 8, !tbaa !5
  %169 = load double, ptr %168, align 8, !tbaa !9
  %170 = fmul double %167, %169
  %171 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %165, double %170)
  %172 = load ptr, ptr %9, align 8, !tbaa !5
  %173 = load double, ptr %172, align 8, !tbaa !9
  %174 = fmul double %171, %173
  %175 = load double, ptr %15, align 8, !tbaa !9
  %176 = fmul double %174, %175
  %177 = load ptr, ptr %8, align 8, !tbaa !5
  %178 = load double, ptr %177, align 8, !tbaa !9
  %179 = fdiv double %176, %178
  %180 = load double, ptr %16, align 8, !tbaa !9
  %181 = load double, ptr %16, align 8, !tbaa !9
  %182 = fmul double 6.500000e-01, %181
  %183 = load double, ptr %16, align 8, !tbaa !9
  %184 = load double, ptr %16, align 8, !tbaa !9
  %185 = call double @llvm.fmuladd.f64(double %183, double %184, double 1.000000e+00)
  %186 = fdiv double %182, %185
  %187 = fsub double 1.000000e+00, %186
  %188 = call double @llvm.fmuladd.f64(double 0x3FE6ACEEE0F3CB3E, double %187, double %180)
  %189 = call double @sqrt(double noundef %188) #4
  %190 = call double @sqrt(double noundef %189) #4
  %191 = fneg double %179
  %192 = call double @llvm.fmuladd.f64(double %191, double %190, double %162)
  %193 = fmul double %151, %192
  store double %193, ptr %5, align 8
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %13) #4
  br label %194

194:                                              ; preds = %109, %38, %33
  %195 = load double, ptr %5, align 8
  ret double %195
}

; Function Attrs: nounwind
declare double @pow(double noundef, double noundef) #2

; Function Attrs: nounwind
declare double @exp(double noundef) #2

attributes #0 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }
attributes #1 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #2 = { 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 #3 = { 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 = { 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, !10, i64 64}
!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}


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

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

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

; Function Attrs: mustprogress nofree nounwind willreturn memory(write, argmem: readwrite) uwtable
define dso_local noundef double @_ZN14spatial_region1wERKdRdS1_(ptr nocapture noundef nonnull readonly align 8 dereferenceable(608) %0, ptr nocapture noundef nonnull readonly align 8 dereferenceable(8) %1, ptr nocapture noundef nonnull align 8 dereferenceable(8) %2, ptr nocapture noundef nonnull readonly align 8 dereferenceable(8) %3) local_unnamed_addr #2 align 2 {
  %5 = load double, ptr %2, align 8, !tbaa !5
  %6 = fsub double 1.000000e+00, %5
  %7 = load double, ptr %1, align 8, !tbaa !5
  %8 = fmul double %6, %7
  %9 = fcmp olt double %8, 1.000000e+00
  %10 = fcmp oeq double %5, 0.000000e+00
  %or.cond = or i1 %10, %9
  br i1 %or.cond, label %88, label %11

11:                                               ; preds = %4
  %12 = load double, ptr %3, align 8, !tbaa !5
  %13 = fcmp oeq double %12, 0.000000e+00
  br i1 %13, label %88, label %14

14:                                               ; preds = %11
  %15 = fcmp ogt double %12, 1.330000e-01
  br i1 %15, label %16, label %43

16:                                               ; preds = %14
  %17 = fmul double %6, %12
  %18 = fdiv double %5, %17
  %19 = tail call double @pow(double noundef %18, double noundef 0x3FD5555555555555) #3
  %20 = fmul double %19, %19
  %21 = getelementptr inbounds i8, ptr %0, i64 64
  %22 = load double, ptr %21, align 8, !tbaa !9
  %23 = fmul double %22, -2.059100e-03
  %24 = load double, ptr %1, align 8, !tbaa !5
  %25 = fmul double %24, %24
  %26 = fdiv double %23, %25
  %27 = fmul double %18, -2.000000e+00
  %28 = fdiv double %27, 3.000000e+00
  %29 = tail call double @exp(double noundef %28) #3
  %30 = fmul double %29, %26
  %31 = fadd double %20, 8.004900e-01
  %32 = fmul double %31, %31
  %33 = fmul double %31, %32
  %34 = tail call double @sqrt(double noundef %33) #3
  %35 = tail call double @sqrt(double noundef %34) #3
  %36 = load double, ptr %2, align 8, !tbaa !5
  %37 = fsub double 1.000000e+00, %36
  %38 = fmul double %36, %36
  %39 = tail call double @llvm.fmuladd.f64(double %37, double 2.000000e+00, double %38)
  %40 = fmul double %20, 6.500000e-01
  %41 = tail call double @llvm.fmuladd.f64(double %20, double %20, double 1.000000e+00)
  %42 = fdiv double %40, %41
  br label %.sink.split

43:                                               ; preds = %14
  %44 = fmul double %12, 1.500000e+01
  %45 = fdiv double 2.000000e+00, %44
  %46 = fadd double %45, 1.000000e+00
  %47 = fdiv double 1.000000e+00, %46
  %48 = fmul double %5, %47
  store double %48, ptr %2, align 8, !tbaa !5
  %49 = fsub double 1.000000e+00, %48
  %50 = load double, ptr %3, align 8, !tbaa !5
  %51 = fmul double %49, %50
  %52 = fdiv double %48, %51
  %53 = tail call double @pow(double noundef %52, double noundef 0x3FD5555555555555) #3
  %54 = fmul double %53, %53
  %55 = fmul double %47, -2.059100e-03
  %56 = getelementptr inbounds i8, ptr %0, i64 64
  %57 = load double, ptr %56, align 8, !tbaa !9
  %58 = fmul double %55, %57
  %59 = load double, ptr %1, align 8, !tbaa !5
  %60 = fmul double %59, %59
  %61 = fdiv double %58, %60
  %62 = fmul double %52, -2.000000e+00
  %63 = fdiv double %62, 3.000000e+00
  %64 = tail call double @exp(double noundef %63) #3
  %65 = fmul double %64, %61
  %66 = fadd double %54, 8.004900e-01
  %67 = fmul double %66, %66
  %68 = fmul double %66, %67
  %69 = tail call double @sqrt(double noundef %68) #3
  %70 = tail call double @sqrt(double noundef %69) #3
  %71 = load double, ptr %2, align 8, !tbaa !5
  %72 = fsub double 1.000000e+00, %71
  %73 = fmul double %71, %71
  %74 = tail call double @llvm.fmuladd.f64(double %72, double 2.000000e+00, double %73)
  %75 = fmul double %54, 6.500000e-01
  %76 = tail call double @llvm.fmuladd.f64(double %54, double %54, double 1.000000e+00)
  %77 = fdiv double %75, %76
  br label %.sink.split

.sink.split:                                      ; preds = %16, %43
  %.sink73 = phi double [ %77, %43 ], [ %42, %16 ]
  %.sink72 = phi double [ %54, %43 ], [ %20, %16 ]
  %.sink68 = phi double [ %74, %43 ], [ %39, %16 ]
  %.sink64 = phi double [ %53, %43 ], [ %19, %16 ]
  %.sink63 = phi double [ %71, %43 ], [ %36, %16 ]
  %.pn = phi double [ %70, %43 ], [ %35, %16 ]
  %.sink = phi double [ %65, %43 ], [ %30, %16 ]
  %.sink61 = fdiv double 1.000000e+00, %.pn
  %.sink66 = load double, ptr %3, align 8, !tbaa !5
  %78 = fsub double 1.000000e+00, %.sink73
  %79 = tail call double @llvm.fmuladd.f64(double %78, double 0x3FE6ACEEE0F3CB3E, double %.sink72)
  %80 = tail call double @sqrt(double noundef %79) #3
  %81 = tail call double @sqrt(double noundef %80) #3
  %82 = fneg double %.sink68
  %83 = fmul double %.sink66, %82
  %84 = fmul double %.sink64, %83
  %85 = fdiv double %84, %.sink63
  %86 = tail call double @llvm.fmuladd.f64(double %85, double %81, double %.sink61)
  %87 = fmul double %.sink, %86
  br label %88

88:                                               ; preds = %.sink.split, %4, %11
  %.0 = phi double [ 0.000000e+00, %11 ], [ 0.000000e+00, %4 ], [ %87, %.sink.split ]
  ret double %.0
}

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

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

attributes #0 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #1 = { 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 #2 = { mustprogress nofree nounwind willreturn memory(write, argmem: readwrite) 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 = { 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 = !{!"double", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C++ TBAA"}
!9 = !{!10, !6, i64 64}
!10 = !{!"_ZTS14spatial_region", !11, i64 0, !11, i64 4, !11, i64 8, !6, i64 16, !6, i64 24, !6, i64 32, !6, i64 40, !11, i64 48, !11, i64 52, !11, i64 56, !6, i64 64, !12, i64 72, !11, i64 80, !11, i64 84, !11, i64 88, !12, i64 96, !11, i64 104, !11, i64 108, !11, i64 112, !11, i64 116, !11, i64 120, !11, i64 124, !11, i64 128, !13, i64 136, !19, i64 144, !6, i64 176, !6, i64 184, !6, i64 192, !11, i64 200, !11, i64 204, !11, i64 208, !12, i64 216, !6, i64 224, !6, i64 232, !6, i64 240, !6, i64 248, !12, i64 256, !6, i64 264, !6, i64 272, !6, 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 = !{!"int", !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", !11, i64 0, !11, i64 4, !11, 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", !11, i64 0, !11, i64 4, !11, 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", !11, i64 0, !11, i64 4, !11, 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", !11, i64 0, !11, i64 4, !11, 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", !11, i64 0, !11, i64 4, !11, 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}


; Effective O3 passes: EarlyCSEPass,GlobalOptPass,InferFunctionAttrsPass,InstCombinePass,PostOrderFunctionAttrsPass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass