; LLVM IR from CLASS-SZ/class_sz :: class-sz/hyrec (hydrogen)

; 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.TWO_PHOTON_PARAMS = type { [311 x double], [311 x double], [311 x double], [311 x double], [311 x 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 speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #1

; 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 @exp(double noundef) #2

; Function Attrs: nounwind uwtable
declare dso_local void @interpolate_rates(ptr noundef, ptr noundef, ptr noundef, double noundef, double noundef, ptr noundef) #3

; Function Attrs: nounwind uwtable
declare dso_local double @cube(double noundef) #3

; Function Attrs: nounwind uwtable
declare dso_local void @populate_Diffusion(ptr noundef, ptr noundef, ptr noundef, ptr noundef, double noundef, ptr noundef, ptr noundef) #3

; Function Attrs: nounwind uwtable
define dso_local void @populateTS_2photon(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, ptr noundef %4, ptr noundef %5, ptr noundef %6, double noundef %7, double noundef %8, double noundef %9, double noundef %10, double noundef %11, ptr noundef %12, ptr noundef %13, ptr noundef %14, ptr noundef %15, ptr noundef %16, ptr noundef %17, double noundef %18) #3 {
  %20 = alloca ptr, align 8
  %21 = alloca ptr, align 8
  %22 = alloca ptr, align 8
  %23 = alloca ptr, align 8
  %24 = alloca ptr, align 8
  %25 = alloca ptr, align 8
  %26 = alloca ptr, align 8
  %27 = alloca double, align 8
  %28 = alloca double, align 8
  %29 = alloca double, align 8
  %30 = alloca double, align 8
  %31 = alloca double, align 8
  %32 = alloca ptr, align 8
  %33 = alloca ptr, align 8
  %34 = alloca ptr, align 8
  %35 = alloca ptr, align 8
  %36 = alloca ptr, align 8
  %37 = alloca ptr, align 8
  %38 = alloca double, align 8
  %39 = alloca double, align 8
  %40 = alloca i32, align 4
  %41 = alloca double, align 8
  %42 = alloca double, align 8
  %43 = alloca double, align 8
  %44 = alloca double, align 8
  %45 = alloca ptr, align 8
  %46 = alloca ptr, align 8
  %47 = alloca double, align 8
  %48 = alloca double, align 8
  store ptr %0, ptr %20, align 8, !tbaa !5
  store ptr %1, ptr %21, align 8, !tbaa !5
  store ptr %2, ptr %22, align 8, !tbaa !5
  store ptr %3, ptr %23, align 8, !tbaa !5
  store ptr %4, ptr %24, align 8, !tbaa !5
  store ptr %5, ptr %25, align 8, !tbaa !5
  store ptr %6, ptr %26, align 8, !tbaa !5
  store double %7, ptr %27, align 8, !tbaa !9
  store double %8, ptr %28, align 8, !tbaa !9
  store double %9, ptr %29, align 8, !tbaa !9
  store double %10, ptr %30, align 8, !tbaa !9
  store double %11, ptr %31, align 8, !tbaa !9
  store ptr %12, ptr %32, align 8, !tbaa !5
  store ptr %13, ptr %33, align 8, !tbaa !5
  store ptr %14, ptr %34, align 8, !tbaa !5
  store ptr %15, ptr %35, align 8, !tbaa !5
  store ptr %16, ptr %36, align 8, !tbaa !5
  store ptr %17, ptr %37, align 8, !tbaa !5
  store double %18, ptr %38, align 8, !tbaa !9
  call void @llvm.lifetime.start.p0(i64 8, ptr %39) #5
  call void @llvm.lifetime.start.p0(i64 4, ptr %40) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %41) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %42) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %43) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %44) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %45) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %46) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %47) #5
  call void @llvm.lifetime.start.p0(i64 8, ptr %48) #5
  %49 = call ptr @create_1D_array(i32 noundef 311)
  store ptr %49, ptr %45, align 8, !tbaa !5
  %50 = call ptr @create_1D_array(i32 noundef 311)
  store ptr %50, ptr %46, align 8, !tbaa !5
  %51 = load double, ptr %31, align 8, !tbaa !9
  %52 = fmul double 0x433090E1CA8B9039, %51
  %53 = load double, ptr %30, align 8, !tbaa !9
  %54 = fdiv double %52, %53
  %55 = load double, ptr %27, align 8, !tbaa !9
  %56 = fsub double 1.000000e+00, %55
  %57 = fdiv double %54, %56
  store double %57, ptr %41, align 8, !tbaa !9
  %58 = load ptr, ptr %36, align 8, !tbaa !5
  %59 = load ptr, ptr %37, align 8, !tbaa !5
  %60 = load double, ptr %29, align 8, !tbaa !9
  %61 = load double, ptr %28, align 8, !tbaa !9
  %62 = load double, ptr %29, align 8, !tbaa !9
  %63 = fdiv double %61, %62
  %64 = load ptr, ptr %32, align 8, !tbaa !5
  call void @interpolate_rates(ptr noundef %58, ptr noundef %59, ptr noundef %39, double noundef %60, double noundef %63, ptr noundef %64)
  %65 = load ptr, ptr %37, align 8, !tbaa !5
  %66 = getelementptr inbounds double, ptr %65, i64 0
  %67 = load double, ptr %66, align 8, !tbaa !9
  %68 = load double, ptr %39, align 8, !tbaa !9
  %69 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %68, double %67)
  %70 = load double, ptr %41, align 8, !tbaa !9
  %71 = fmul double 3.000000e+00, %70
  %72 = load double, ptr %29, align 8, !tbaa !9
  %73 = fdiv double 0xBFFE37E9F25A4877, %72
  %74 = call double @exp(double noundef %73) #5
  %75 = load double, ptr %29, align 8, !tbaa !9
  %76 = fdiv double 0xC00465BDE9FCF0EC, %75
  %77 = call double @exp(double noundef %76) #5
  %78 = fmul double 0x3FFF400000000000, %77
  %79 = call double @llvm.fmuladd.f64(double 0x3FFAA2F78F1B4CC7, double %74, double %78)
  %80 = call double @llvm.fmuladd.f64(double %71, double %79, double %69)
  %81 = load ptr, ptr %20, align 8, !tbaa !5
  %82 = getelementptr inbounds [2 x double], ptr %81, i64 0
  %83 = getelementptr inbounds [2 x double], ptr %82, i64 0, i64 0
  store double %80, ptr %83, align 8, !tbaa !9
  %84 = load double, ptr %39, align 8, !tbaa !9
  %85 = fneg double %84
  %86 = load ptr, ptr %20, align 8, !tbaa !5
  %87 = getelementptr inbounds [2 x double], ptr %86, i64 0
  %88 = getelementptr inbounds [2 x double], ptr %87, i64 0, i64 1
  store double %85, ptr %88, align 8, !tbaa !9
  %89 = load double, ptr %30, align 8, !tbaa !9
  %90 = load ptr, ptr %36, align 8, !tbaa !5
  %91 = getelementptr inbounds double, ptr %90, i64 0
  %92 = load double, ptr %91, align 8, !tbaa !9
  %93 = fmul double %89, %92
  %94 = load double, ptr %27, align 8, !tbaa !9
  %95 = fmul double %93, %94
  %96 = load double, ptr %27, align 8, !tbaa !9
  %97 = load double, ptr %41, align 8, !tbaa !9
  %98 = fmul double 3.000000e+00, %97
  %99 = load double, ptr %27, align 8, !tbaa !9
  %100 = fsub double 1.000000e+00, %99
  %101 = fmul double %98, %100
  %102 = load ptr, ptr %35, align 8, !tbaa !5
  %103 = getelementptr inbounds double, ptr %102, i64 1
  %104 = load double, ptr %103, align 8, !tbaa !9
  %105 = load double, ptr %29, align 8, !tbaa !9
  %106 = fdiv double 0xC0297F2D647C2D27, %105
  %107 = call double @exp(double noundef %106) #5
  %108 = fmul double 0x3FFF400000000000, %107
  %109 = call double @llvm.fmuladd.f64(double 0x3FFAA2F78F1B4CC7, double %104, double %108)
  %110 = fmul double %101, %109
  %111 = call double @llvm.fmuladd.f64(double %95, double %96, double %110)
  %112 = load ptr, ptr %24, align 8, !tbaa !5
  %113 = getelementptr inbounds double, ptr %112, i64 0
  store double %111, ptr %113, align 8, !tbaa !9
  %114 = load ptr, ptr %37, align 8, !tbaa !5
  %115 = getelementptr inbounds double, ptr %114, i64 1
  %116 = load double, ptr %115, align 8, !tbaa !9
  %117 = load double, ptr %39, align 8, !tbaa !9
  %118 = fadd double %116, %117
  %119 = load double, ptr %41, align 8, !tbaa !9
  %120 = fadd double %118, %119
  %121 = load ptr, ptr %20, align 8, !tbaa !5
  %122 = getelementptr inbounds [2 x double], ptr %121, i64 1
  %123 = getelementptr inbounds [2 x double], ptr %122, i64 0, i64 1
  store double %120, ptr %123, align 8, !tbaa !9
  %124 = load double, ptr %39, align 8, !tbaa !9
  %125 = fmul double -3.000000e+00, %124
  %126 = load ptr, ptr %20, align 8, !tbaa !5
  %127 = getelementptr inbounds [2 x double], ptr %126, i64 1
  %128 = getelementptr inbounds [2 x double], ptr %127, i64 0, i64 0
  store double %125, ptr %128, align 8, !tbaa !9
  %129 = load double, ptr %30, align 8, !tbaa !9
  %130 = load ptr, ptr %36, align 8, !tbaa !5
  %131 = getelementptr inbounds double, ptr %130, i64 1
  %132 = load double, ptr %131, align 8, !tbaa !9
  %133 = fmul double %129, %132
  %134 = load double, ptr %27, align 8, !tbaa !9
  %135 = fmul double %133, %134
  %136 = load double, ptr %27, align 8, !tbaa !9
  %137 = load double, ptr %41, align 8, !tbaa !9
  %138 = fmul double 3.000000e+00, %137
  %139 = load double, ptr %27, align 8, !tbaa !9
  %140 = fsub double 1.000000e+00, %139
  %141 = fmul double %138, %140
  %142 = load ptr, ptr %35, align 8, !tbaa !5
  %143 = getelementptr inbounds double, ptr %142, i64 0
  %144 = load double, ptr %143, align 8, !tbaa !9
  %145 = fmul double %141, %144
  %146 = call double @llvm.fmuladd.f64(double %135, double %136, double %145)
  %147 = load ptr, ptr %24, align 8, !tbaa !5
  %148 = getelementptr inbounds double, ptr %147, i64 1
  store double %146, ptr %148, align 8, !tbaa !9
  store i32 0, ptr %40, align 4, !tbaa !11
  br label %149

149:                                              ; preds = %286, %19
  %150 = load i32, ptr %40, align 4, !tbaa !11
  %151 = icmp ult i32 %150, 311
  br i1 %151, label %152, label %289

152:                                              ; preds = %149
  %153 = load ptr, ptr %33, align 8, !tbaa !5
  %154 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %153, i32 0, i32 0
  %155 = load i32, ptr %40, align 4, !tbaa !11
  %156 = zext i32 %155 to i64
  %157 = getelementptr inbounds [311 x double], ptr %154, i64 0, i64 %156
  %158 = load double, ptr %157, align 8, !tbaa !9
  %159 = fsub double %158, 0x402465BDE9FCF0EC
  %160 = load double, ptr %29, align 8, !tbaa !9
  %161 = fdiv double %159, %160
  %162 = call double @exp(double noundef %161) #5
  store double %162, ptr %44, align 8, !tbaa !9
  %163 = load ptr, ptr %33, align 8, !tbaa !5
  %164 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %163, i32 0, i32 2
  %165 = load i32, ptr %40, align 4, !tbaa !11
  %166 = zext i32 %165 to i64
  %167 = getelementptr inbounds [311 x double], ptr %164, i64 0, i64 %166
  %168 = load double, ptr %167, align 8, !tbaa !9
  %169 = fneg double %168
  %170 = load ptr, ptr %33, align 8, !tbaa !5
  %171 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %170, i32 0, i32 0
  %172 = load i32, ptr %40, align 4, !tbaa !11
  %173 = zext i32 %172 to i64
  %174 = getelementptr inbounds [311 x double], ptr %171, i64 0, i64 %173
  %175 = load double, ptr %174, align 8, !tbaa !9
  %176 = fsub double %175, 0x402465BDE9FCF0EC
  %177 = load double, ptr %29, align 8, !tbaa !9
  %178 = fdiv double %176, %177
  %179 = call double @exp(double noundef %178) #5
  %180 = fsub double %179, 1.000000e+00
  %181 = call double @llvm.fabs.f64(double %180)
  %182 = fdiv double %169, %181
  %183 = load ptr, ptr %22, align 8, !tbaa !5
  %184 = getelementptr inbounds ptr, ptr %183, i64 0
  %185 = load ptr, ptr %184, align 8, !tbaa !5
  %186 = load i32, ptr %40, align 4, !tbaa !11
  %187 = zext i32 %186 to i64
  %188 = getelementptr inbounds double, ptr %185, i64 %187
  store double %182, ptr %188, align 8, !tbaa !9
  %189 = load ptr, ptr %20, align 8, !tbaa !5
  %190 = getelementptr inbounds [2 x double], ptr %189, i64 0
  %191 = getelementptr inbounds [2 x double], ptr %190, i64 0, i64 0
  %192 = load double, ptr %191, align 8, !tbaa !9
  %193 = fsub double %192, %182
  store double %193, ptr %191, align 8, !tbaa !9
  %194 = load ptr, ptr %22, align 8, !tbaa !5
  %195 = getelementptr inbounds ptr, ptr %194, i64 0
  %196 = load ptr, ptr %195, align 8, !tbaa !5
  %197 = load i32, ptr %40, align 4, !tbaa !11
  %198 = zext i32 %197 to i64
  %199 = getelementptr inbounds double, ptr %196, i64 %198
  %200 = load double, ptr %199, align 8, !tbaa !9
  %201 = load double, ptr %44, align 8, !tbaa !9
  %202 = fmul double %200, %201
  %203 = load ptr, ptr %21, align 8, !tbaa !5
  %204 = getelementptr inbounds ptr, ptr %203, i64 0
  %205 = load ptr, ptr %204, align 8, !tbaa !5
  %206 = load i32, ptr %40, align 4, !tbaa !11
  %207 = zext i32 %206 to i64
  %208 = getelementptr inbounds double, ptr %205, i64 %207
  store double %202, ptr %208, align 8, !tbaa !9
  %209 = load double, ptr %29, align 8, !tbaa !9
  %210 = fdiv double 0xBFFE37E9F25A4877, %209
  %211 = call double @exp(double noundef %210) #5
  %212 = fneg double %211
  %213 = fdiv double %212, 3.000000e+00
  %214 = load ptr, ptr %33, align 8, !tbaa !5
  %215 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %214, i32 0, i32 3
  %216 = load i32, ptr %40, align 4, !tbaa !11
  %217 = zext i32 %216 to i64
  %218 = getelementptr inbounds [311 x double], ptr %215, i64 0, i64 %217
  %219 = load double, ptr %218, align 8, !tbaa !9
  %220 = fmul double %213, %219
  %221 = load ptr, ptr %33, align 8, !tbaa !5
  %222 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %221, i32 0, i32 0
  %223 = load i32, ptr %40, align 4, !tbaa !11
  %224 = zext i32 %223 to i64
  %225 = getelementptr inbounds [311 x double], ptr %222, i64 0, i64 %224
  %226 = load double, ptr %225, align 8, !tbaa !9
  %227 = fsub double %226, 0x40282CBB284839FB
  %228 = load double, ptr %29, align 8, !tbaa !9
  %229 = fdiv double %227, %228
  %230 = call double @exp(double noundef %229) #5
  %231 = fsub double %230, 1.000000e+00
  %232 = call double @llvm.fabs.f64(double %231)
  %233 = fdiv double %220, %232
  %234 = load double, ptr %29, align 8, !tbaa !9
  %235 = fdiv double 0xC00465BDE9FCF0EC, %234
  %236 = call double @exp(double noundef %235) #5
  %237 = fdiv double %236, 3.000000e+00
  %238 = load ptr, ptr %33, align 8, !tbaa !5
  %239 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %238, i32 0, i32 4
  %240 = load i32, ptr %40, align 4, !tbaa !11
  %241 = zext i32 %240 to i64
  %242 = getelementptr inbounds [311 x double], ptr %239, i64 0, i64 %241
  %243 = load double, ptr %242, align 8, !tbaa !9
  %244 = fmul double %237, %243
  %245 = load ptr, ptr %33, align 8, !tbaa !5
  %246 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %245, i32 0, i32 0
  %247 = load i32, ptr %40, align 4, !tbaa !11
  %248 = zext i32 %247 to i64
  %249 = getelementptr inbounds [311 x double], ptr %246, i64 0, i64 %248
  %250 = load double, ptr %249, align 8, !tbaa !9
  %251 = fsub double %250, 0x40297F2D647C2D27
  %252 = load double, ptr %29, align 8, !tbaa !9
  %253 = fdiv double %251, %252
  %254 = call double @exp(double noundef %253) #5
  %255 = fsub double %254, 1.000000e+00
  %256 = call double @llvm.fabs.f64(double %255)
  %257 = fdiv double %244, %256
  %258 = fsub double %233, %257
  %259 = load ptr, ptr %22, align 8, !tbaa !5
  %260 = getelementptr inbounds ptr, ptr %259, i64 1
  %261 = load ptr, ptr %260, align 8, !tbaa !5
  %262 = load i32, ptr %40, align 4, !tbaa !11
  %263 = zext i32 %262 to i64
  %264 = getelementptr inbounds double, ptr %261, i64 %263
  store double %258, ptr %264, align 8, !tbaa !9
  %265 = load ptr, ptr %20, align 8, !tbaa !5
  %266 = getelementptr inbounds [2 x double], ptr %265, i64 1
  %267 = getelementptr inbounds [2 x double], ptr %266, i64 0, i64 1
  %268 = load double, ptr %267, align 8, !tbaa !9
  %269 = fsub double %268, %258
  store double %269, ptr %267, align 8, !tbaa !9
  %270 = load ptr, ptr %22, align 8, !tbaa !5
  %271 = getelementptr inbounds ptr, ptr %270, i64 1
  %272 = load ptr, ptr %271, align 8, !tbaa !5
  %273 = load i32, ptr %40, align 4, !tbaa !11
  %274 = zext i32 %273 to i64
  %275 = getelementptr inbounds double, ptr %272, i64 %274
  %276 = load double, ptr %275, align 8, !tbaa !9
  %277 = fmul double %276, 3.000000e+00
  %278 = load double, ptr %44, align 8, !tbaa !9
  %279 = fmul double %277, %278
  %280 = load ptr, ptr %21, align 8, !tbaa !5
  %281 = getelementptr inbounds ptr, ptr %280, i64 1
  %282 = load ptr, ptr %281, align 8, !tbaa !5
  %283 = load i32, ptr %40, align 4, !tbaa !11
  %284 = zext i32 %283 to i64
  %285 = getelementptr inbounds double, ptr %282, i64 %284
  store double %279, ptr %285, align 8, !tbaa !9
  br label %286

286:                                              ; preds = %152
  %287 = load i32, ptr %40, align 4, !tbaa !11
  %288 = add i32 %287, 1
  store i32 %288, ptr %40, align 4, !tbaa !11
  br label %149, !llvm.loop !13

289:                                              ; preds = %149
  %290 = load ptr, ptr %45, align 8, !tbaa !5
  %291 = load ptr, ptr %46, align 8, !tbaa !5
  %292 = load double, ptr %28, align 8, !tbaa !9
  %293 = load ptr, ptr %33, align 8, !tbaa !5
  %294 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %293, i32 0, i32 0
  %295 = getelementptr inbounds [311 x double], ptr %294, i64 0, i64 0
  %296 = load ptr, ptr %33, align 8, !tbaa !5
  %297 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %296, i32 0, i32 1
  %298 = getelementptr inbounds [311 x double], ptr %297, i64 0, i64 0
  call void @populate_Diffusion(ptr noundef %290, ptr noundef %291, ptr noundef %47, ptr noundef %48, double noundef %292, ptr noundef %295, ptr noundef %298)
  %299 = load double, ptr %48, align 8, !tbaa !9
  %300 = load double, ptr %47, align 8, !tbaa !9
  %301 = fadd double %299, %300
  %302 = load ptr, ptr %20, align 8, !tbaa !5
  %303 = getelementptr inbounds [2 x double], ptr %302, i64 1
  %304 = getelementptr inbounds [2 x double], ptr %303, i64 0, i64 1
  %305 = load double, ptr %304, align 8, !tbaa !9
  %306 = fadd double %305, %301
  store double %306, ptr %304, align 8, !tbaa !9
  store i32 0, ptr %40, align 4, !tbaa !11
  br label %307

307:                                              ; preds = %588, %289
  %308 = load i32, ptr %40, align 4, !tbaa !11
  %309 = icmp ult i32 %308, 311
  br i1 %309, label %310, label %591

310:                                              ; preds = %307
  %311 = load ptr, ptr %21, align 8, !tbaa !5
  %312 = getelementptr inbounds ptr, ptr %311, i64 0
  %313 = load ptr, ptr %312, align 8, !tbaa !5
  %314 = load i32, ptr %40, align 4, !tbaa !11
  %315 = zext i32 %314 to i64
  %316 = getelementptr inbounds double, ptr %313, i64 %315
  %317 = load double, ptr %316, align 8, !tbaa !9
  %318 = load ptr, ptr %21, align 8, !tbaa !5
  %319 = getelementptr inbounds ptr, ptr %318, i64 1
  %320 = load ptr, ptr %319, align 8, !tbaa !5
  %321 = load i32, ptr %40, align 4, !tbaa !11
  %322 = zext i32 %321 to i64
  %323 = getelementptr inbounds double, ptr %320, i64 %322
  %324 = load double, ptr %323, align 8, !tbaa !9
  %325 = fadd double %317, %324
  %326 = fneg double %325
  %327 = load ptr, ptr %45, align 8, !tbaa !5
  %328 = load i32, ptr %40, align 4, !tbaa !11
  %329 = zext i32 %328 to i64
  %330 = getelementptr inbounds double, ptr %327, i64 %329
  %331 = load double, ptr %330, align 8, !tbaa !9
  %332 = fadd double %326, %331
  %333 = load ptr, ptr %46, align 8, !tbaa !5
  %334 = load i32, ptr %40, align 4, !tbaa !11
  %335 = zext i32 %334 to i64
  %336 = getelementptr inbounds double, ptr %333, i64 %335
  %337 = load double, ptr %336, align 8, !tbaa !9
  %338 = fadd double %332, %337
  store double %338, ptr %42, align 8, !tbaa !9
  %339 = load i32, ptr %40, align 4, !tbaa !11
  %340 = icmp uge i32 %339, 100
  br i1 %340, label %341, label %344

341:                                              ; preds = %310
  %342 = load i32, ptr %40, align 4, !tbaa !11
  %343 = icmp ult i32 %342, 139
  br i1 %343, label %350, label %344

344:                                              ; preds = %341, %310
  %345 = load i32, ptr %40, align 4, !tbaa !11
  %346 = icmp ugt i32 %345, 140
  br i1 %346, label %347, label %377

347:                                              ; preds = %344
  %348 = load i32, ptr %40, align 4, !tbaa !11
  %349 = icmp ult i32 %348, 180
  br i1 %349, label %350, label %377

350:                                              ; preds = %347, %341
  %351 = load ptr, ptr %45, align 8, !tbaa !5
  %352 = load i32, ptr %40, align 4, !tbaa !11
  %353 = sub i32 %352, 1
  %354 = zext i32 %353 to i64
  %355 = getelementptr inbounds double, ptr %351, i64 %354
  %356 = load double, ptr %355, align 8, !tbaa !9
  %357 = fneg double %356
  %358 = load ptr, ptr %23, align 8, !tbaa !5
  %359 = getelementptr inbounds ptr, ptr %358, i64 1
  %360 = load ptr, ptr %359, align 8, !tbaa !5
  %361 = load i32, ptr %40, align 4, !tbaa !11
  %362 = zext i32 %361 to i64
  %363 = getelementptr inbounds double, ptr %360, i64 %362
  store double %357, ptr %363, align 8, !tbaa !9
  %364 = load ptr, ptr %46, align 8, !tbaa !5
  %365 = load i32, ptr %40, align 4, !tbaa !11
  %366 = add i32 %365, 1
  %367 = zext i32 %366 to i64
  %368 = getelementptr inbounds double, ptr %364, i64 %367
  %369 = load double, ptr %368, align 8, !tbaa !9
  %370 = fneg double %369
  %371 = load ptr, ptr %23, align 8, !tbaa !5
  %372 = getelementptr inbounds ptr, ptr %371, i64 2
  %373 = load ptr, ptr %372, align 8, !tbaa !5
  %374 = load i32, ptr %40, align 4, !tbaa !11
  %375 = zext i32 %374 to i64
  %376 = getelementptr inbounds double, ptr %373, i64 %375
  store double %370, ptr %376, align 8, !tbaa !9
  br label %377

377:                                              ; preds = %350, %347, %344
  %378 = load i32, ptr %40, align 4, !tbaa !11
  %379 = icmp eq i32 %378, 139
  br i1 %379, label %380, label %422

380:                                              ; preds = %377
  %381 = load ptr, ptr %45, align 8, !tbaa !5
  %382 = load i32, ptr %40, align 4, !tbaa !11
  %383 = sub i32 %382, 1
  %384 = zext i32 %383 to i64
  %385 = getelementptr inbounds double, ptr %381, i64 %384
  %386 = load double, ptr %385, align 8, !tbaa !9
  %387 = fneg double %386
  %388 = load ptr, ptr %23, align 8, !tbaa !5
  %389 = getelementptr inbounds ptr, ptr %388, i64 1
  %390 = load ptr, ptr %389, align 8, !tbaa !5
  %391 = load i32, ptr %40, align 4, !tbaa !11
  %392 = zext i32 %391 to i64
  %393 = getelementptr inbounds double, ptr %390, i64 %392
  store double %387, ptr %393, align 8, !tbaa !9
  %394 = load ptr, ptr %23, align 8, !tbaa !5
  %395 = getelementptr inbounds ptr, ptr %394, i64 2
  %396 = load ptr, ptr %395, align 8, !tbaa !5
  %397 = load i32, ptr %40, align 4, !tbaa !11
  %398 = zext i32 %397 to i64
  %399 = getelementptr inbounds double, ptr %396, i64 %398
  store double 0.000000e+00, ptr %399, align 8, !tbaa !9
  %400 = load double, ptr %48, align 8, !tbaa !9
  %401 = load ptr, ptr %22, align 8, !tbaa !5
  %402 = getelementptr inbounds ptr, ptr %401, i64 1
  %403 = load ptr, ptr %402, align 8, !tbaa !5
  %404 = load i32, ptr %40, align 4, !tbaa !11
  %405 = zext i32 %404 to i64
  %406 = getelementptr inbounds double, ptr %403, i64 %405
  %407 = load double, ptr %406, align 8, !tbaa !9
  %408 = fsub double %407, %400
  store double %408, ptr %406, align 8, !tbaa !9
  %409 = load ptr, ptr %45, align 8, !tbaa !5
  %410 = load i32, ptr %40, align 4, !tbaa !11
  %411 = zext i32 %410 to i64
  %412 = getelementptr inbounds double, ptr %409, i64 %411
  %413 = load double, ptr %412, align 8, !tbaa !9
  %414 = load ptr, ptr %21, align 8, !tbaa !5
  %415 = getelementptr inbounds ptr, ptr %414, i64 1
  %416 = load ptr, ptr %415, align 8, !tbaa !5
  %417 = load i32, ptr %40, align 4, !tbaa !11
  %418 = zext i32 %417 to i64
  %419 = getelementptr inbounds double, ptr %416, i64 %418
  %420 = load double, ptr %419, align 8, !tbaa !9
  %421 = fsub double %420, %413
  store double %421, ptr %419, align 8, !tbaa !9
  br label %422

422:                                              ; preds = %380, %377
  %423 = load i32, ptr %40, align 4, !tbaa !11
  %424 = icmp eq i32 %423, 140
  br i1 %424, label %425, label %467

425:                                              ; preds = %422
  %426 = load ptr, ptr %23, align 8, !tbaa !5
  %427 = getelementptr inbounds ptr, ptr %426, i64 1
  %428 = load ptr, ptr %427, align 8, !tbaa !5
  %429 = load i32, ptr %40, align 4, !tbaa !11
  %430 = zext i32 %429 to i64
  %431 = getelementptr inbounds double, ptr %428, i64 %430
  store double 0.000000e+00, ptr %431, align 8, !tbaa !9
  %432 = load ptr, ptr %46, align 8, !tbaa !5
  %433 = load i32, ptr %40, align 4, !tbaa !11
  %434 = add i32 %433, 1
  %435 = zext i32 %434 to i64
  %436 = getelementptr inbounds double, ptr %432, i64 %435
  %437 = load double, ptr %436, align 8, !tbaa !9
  %438 = fneg double %437
  %439 = load ptr, ptr %23, align 8, !tbaa !5
  %440 = getelementptr inbounds ptr, ptr %439, i64 2
  %441 = load ptr, ptr %440, align 8, !tbaa !5
  %442 = load i32, ptr %40, align 4, !tbaa !11
  %443 = zext i32 %442 to i64
  %444 = getelementptr inbounds double, ptr %441, i64 %443
  store double %438, ptr %444, align 8, !tbaa !9
  %445 = load double, ptr %47, align 8, !tbaa !9
  %446 = load ptr, ptr %22, align 8, !tbaa !5
  %447 = getelementptr inbounds ptr, ptr %446, i64 1
  %448 = load ptr, ptr %447, align 8, !tbaa !5
  %449 = load i32, ptr %40, align 4, !tbaa !11
  %450 = zext i32 %449 to i64
  %451 = getelementptr inbounds double, ptr %448, i64 %450
  %452 = load double, ptr %451, align 8, !tbaa !9
  %453 = fsub double %452, %445
  store double %453, ptr %451, align 8, !tbaa !9
  %454 = load ptr, ptr %46, align 8, !tbaa !5
  %455 = load i32, ptr %40, align 4, !tbaa !11
  %456 = zext i32 %455 to i64
  %457 = getelementptr inbounds double, ptr %454, i64 %456
  %458 = load double, ptr %457, align 8, !tbaa !9
  %459 = load ptr, ptr %21, align 8, !tbaa !5
  %460 = getelementptr inbounds ptr, ptr %459, i64 1
  %461 = load ptr, ptr %460, align 8, !tbaa !5
  %462 = load i32, ptr %40, align 4, !tbaa !11
  %463 = zext i32 %462 to i64
  %464 = getelementptr inbounds double, ptr %461, i64 %463
  %465 = load double, ptr %464, align 8, !tbaa !9
  %466 = fsub double %465, %458
  store double %466, ptr %464, align 8, !tbaa !9
  br label %467

467:                                              ; preds = %425, %422
  %468 = load double, ptr %42, align 8, !tbaa !9
  %469 = load double, ptr %27, align 8, !tbaa !9
  %470 = fsub double 1.000000e+00, %469
  %471 = fmul double %468, %470
  %472 = load ptr, ptr %33, align 8, !tbaa !5
  %473 = getelementptr inbounds %struct.TWO_PHOTON_PARAMS, ptr %472, i32 0, i32 0
  %474 = load i32, ptr %40, align 4, !tbaa !11
  %475 = zext i32 %474 to i64
  %476 = getelementptr inbounds [311 x double], ptr %473, i64 0, i64 %475
  %477 = load double, ptr %476, align 8, !tbaa !9
  %478 = fdiv double 0x3F2040380F7C715F, %477
  %479 = call double @cube(double noundef %478)
  %480 = fmul double %471, %479
  %481 = load double, ptr %30, align 8, !tbaa !9
  %482 = fmul double %480, %481
  %483 = fdiv double %482, 8.000000e+00
  %484 = fdiv double %483, 0x400921FB54442D18
  %485 = load double, ptr %31, align 8, !tbaa !9
  %486 = fdiv double %484, %485
  %487 = load ptr, ptr %26, align 8, !tbaa !5
  %488 = load i32, ptr %40, align 4, !tbaa !11
  %489 = zext i32 %488 to i64
  %490 = getelementptr inbounds double, ptr %487, i64 %489
  store double %486, ptr %490, align 8, !tbaa !9
  %491 = load ptr, ptr %26, align 8, !tbaa !5
  %492 = load i32, ptr %40, align 4, !tbaa !11
  %493 = zext i32 %492 to i64
  %494 = getelementptr inbounds double, ptr %491, i64 %493
  %495 = load double, ptr %494, align 8, !tbaa !9
  %496 = fcmp ogt double %495, 1.000000e-30
  br i1 %496, label %497, label %544

497:                                              ; preds = %467
  %498 = load ptr, ptr %26, align 8, !tbaa !5
  %499 = load i32, ptr %40, align 4, !tbaa !11
  %500 = zext i32 %499 to i64
  %501 = getelementptr inbounds double, ptr %498, i64 %500
  %502 = load double, ptr %501, align 8, !tbaa !9
  %503 = fneg double %502
  %504 = call double @exp(double noundef %503) #5
  %505 = fsub double 1.000000e+00, %504
  %506 = load ptr, ptr %26, align 8, !tbaa !5
  %507 = load i32, ptr %40, align 4, !tbaa !11
  %508 = zext i32 %507 to i64
  %509 = getelementptr inbounds double, ptr %506, i64 %508
  %510 = load double, ptr %509, align 8, !tbaa !9
  %511 = fdiv double %505, %510
  store double %511, ptr %43, align 8, !tbaa !9
  %512 = load double, ptr %42, align 8, !tbaa !9
  %513 = load double, ptr %43, align 8, !tbaa !9
  %514 = fsub double 1.000000e+00, %513
  %515 = fdiv double %512, %514
  %516 = load ptr, ptr %23, align 8, !tbaa !5
  %517 = getelementptr inbounds ptr, ptr %516, i64 0
  %518 = load ptr, ptr %517, align 8, !tbaa !5
  %519 = load i32, ptr %40, align 4, !tbaa !11
  %520 = zext i32 %519 to i64
  %521 = getelementptr inbounds double, ptr %518, i64 %520
  store double %515, ptr %521, align 8, !tbaa !9
  %522 = load ptr, ptr %23, align 8, !tbaa !5
  %523 = getelementptr inbounds ptr, ptr %522, i64 0
  %524 = load ptr, ptr %523, align 8, !tbaa !5
  %525 = load i32, ptr %40, align 4, !tbaa !11
  %526 = zext i32 %525 to i64
  %527 = getelementptr inbounds double, ptr %524, i64 %526
  %528 = load double, ptr %527, align 8, !tbaa !9
  %529 = load double, ptr %27, align 8, !tbaa !9
  %530 = fsub double 1.000000e+00, %529
  %531 = fmul double %528, %530
  %532 = load ptr, ptr %34, align 8, !tbaa !5
  %533 = load i32, ptr %40, align 4, !tbaa !11
  %534 = zext i32 %533 to i64
  %535 = getelementptr inbounds double, ptr %532, i64 %534
  %536 = load double, ptr %535, align 8, !tbaa !9
  %537 = fmul double %531, %536
  %538 = load double, ptr %43, align 8, !tbaa !9
  %539 = fmul double %537, %538
  %540 = load ptr, ptr %25, align 8, !tbaa !5
  %541 = load i32, ptr %40, align 4, !tbaa !11
  %542 = zext i32 %541 to i64
  %543 = getelementptr inbounds double, ptr %540, i64 %542
  store double %539, ptr %543, align 8, !tbaa !9
  br label %587

544:                                              ; preds = %467
  %545 = load ptr, ptr %23, align 8, !tbaa !5
  %546 = getelementptr inbounds ptr, ptr %545, i64 0
  %547 = load ptr, ptr %546, align 8, !tbaa !5
  %548 = load i32, ptr %40, align 4, !tbaa !11
  %549 = zext i32 %548 to i64
  %550 = getelementptr inbounds double, ptr %547, i64 %549
  store double 1.000000e+00, ptr %550, align 8, !tbaa !9
  %551 = load ptr, ptr %22, align 8, !tbaa !5
  %552 = getelementptr inbounds ptr, ptr %551, i64 1
  %553 = load ptr, ptr %552, align 8, !tbaa !5
  %554 = load i32, ptr %40, align 4, !tbaa !11
  %555 = zext i32 %554 to i64
  %556 = getelementptr inbounds double, ptr %553, i64 %555
  store double 0.000000e+00, ptr %556, align 8, !tbaa !9
  %557 = load ptr, ptr %22, align 8, !tbaa !5
  %558 = getelementptr inbounds ptr, ptr %557, i64 0
  %559 = load ptr, ptr %558, align 8, !tbaa !5
  %560 = load i32, ptr %40, align 4, !tbaa !11
  %561 = zext i32 %560 to i64
  %562 = getelementptr inbounds double, ptr %559, i64 %561
  store double 0.000000e+00, ptr %562, align 8, !tbaa !9
  %563 = load ptr, ptr %21, align 8, !tbaa !5
  %564 = getelementptr inbounds ptr, ptr %563, i64 1
  %565 = load ptr, ptr %564, align 8, !tbaa !5
  %566 = load i32, ptr %40, align 4, !tbaa !11
  %567 = zext i32 %566 to i64
  %568 = getelementptr inbounds double, ptr %565, i64 %567
  store double 0.000000e+00, ptr %568, align 8, !tbaa !9
  %569 = load ptr, ptr %21, align 8, !tbaa !5
  %570 = getelementptr inbounds ptr, ptr %569, i64 0
  %571 = load ptr, ptr %570, align 8, !tbaa !5
  %572 = load i32, ptr %40, align 4, !tbaa !11
  %573 = zext i32 %572 to i64
  %574 = getelementptr inbounds double, ptr %571, i64 %573
  store double 0.000000e+00, ptr %574, align 8, !tbaa !9
  %575 = load double, ptr %27, align 8, !tbaa !9
  %576 = fsub double 1.000000e+00, %575
  %577 = load ptr, ptr %34, align 8, !tbaa !5
  %578 = load i32, ptr %40, align 4, !tbaa !11
  %579 = zext i32 %578 to i64
  %580 = getelementptr inbounds double, ptr %577, i64 %579
  %581 = load double, ptr %580, align 8, !tbaa !9
  %582 = fmul double %576, %581
  %583 = load ptr, ptr %25, align 8, !tbaa !5
  %584 = load i32, ptr %40, align 4, !tbaa !11
  %585 = zext i32 %584 to i64
  %586 = getelementptr inbounds double, ptr %583, i64 %585
  store double %582, ptr %586, align 8, !tbaa !9
  br label %587

587:                                              ; preds = %544, %497
  br label %588

588:                                              ; preds = %587
  %589 = load i32, ptr %40, align 4, !tbaa !11
  %590 = add i32 %589, 1
  store i32 %590, ptr %40, align 4, !tbaa !11
  br label %307, !llvm.loop !15

591:                                              ; preds = %307
  %592 = load ptr, ptr %45, align 8, !tbaa !5
  call void @free(ptr noundef %592) #5
  %593 = load ptr, ptr %46, align 8, !tbaa !5
  call void @free(ptr noundef %593) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %48) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %47) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %46) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %45) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %44) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %43) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %42) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %41) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %40) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %39) #5
  ret void
}

declare ptr @create_1D_array(i32 noundef) #4

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

; Function Attrs: nounwind
declare void @free(ptr 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 = { 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 = { "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 = { nounwind }

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"any pointer", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"double", !7, i64 0}
!11 = !{!12, !12, i64 0}
!12 = !{!"int", !7, i64 0}
!13 = distinct !{!13, !14}
!14 = !{!"llvm.loop.mustprogress"}
!15 = distinct !{!15, !14}


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

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

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

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

; Function Attrs: nounwind uwtable
declare dso_local void @interpolate_rates(ptr noundef, ptr noundef, ptr noundef, double noundef, double noundef, ptr noundef) local_unnamed_addr #3

; Function Attrs: nounwind uwtable
declare dso_local double @cube(double noundef) local_unnamed_addr #3

; Function Attrs: nounwind uwtable
declare dso_local void @populate_Diffusion(ptr noundef, ptr noundef, ptr noundef, ptr noundef, double noundef, ptr noundef, ptr noundef) local_unnamed_addr #3

; Function Attrs: nounwind uwtable
define dso_local void @populateTS_2photon(ptr nocapture noundef %0, ptr nocapture noundef readonly %1, ptr nocapture noundef readonly %2, ptr nocapture noundef readonly %3, ptr nocapture noundef writeonly %4, ptr nocapture noundef writeonly %5, ptr nocapture noundef %6, double noundef %7, double noundef %8, double noundef %9, double noundef %10, double noundef %11, ptr noundef %12, ptr noundef %13, ptr nocapture noundef readonly %14, ptr nocapture noundef readonly %15, ptr noundef %16, ptr noundef %17, double noundef %18) local_unnamed_addr #3 {
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %20) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %21) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %22) #6
  %23 = tail call ptr @create_1D_array(i32 noundef 311) #6
  %24 = tail call ptr @create_1D_array(i32 noundef 311) #6
  %25 = fmul double %11, 0x433090E1CA8B9039
  %26 = fdiv double %25, %10
  %27 = fsub double 1.000000e+00, %7
  %28 = fdiv double %26, %27
  %29 = fdiv double %8, %9
  call void @interpolate_rates(ptr noundef %16, ptr noundef %17, ptr noundef nonnull %20, double noundef %9, double noundef %29, ptr noundef %12)
  %30 = load double, ptr %17, align 8, !tbaa !5
  %31 = load double, ptr %20, align 8, !tbaa !5
  %32 = call double @llvm.fmuladd.f64(double %31, double 3.000000e+00, double %30)
  %33 = fmul double %28, 3.000000e+00
  %34 = fdiv double 0xBFFE37E9F25A4877, %9
  %35 = call double @exp(double noundef %34) #6
  %36 = fdiv double 0xC00465BDE9FCF0EC, %9
  %37 = call double @exp(double noundef %36) #6
  %38 = fmul double %37, 0x3FFF400000000000
  %39 = call double @llvm.fmuladd.f64(double %35, double 0x3FFAA2F78F1B4CC7, double %38)
  %40 = call double @llvm.fmuladd.f64(double %33, double %39, double %32)
  store double %40, ptr %0, align 8, !tbaa !5
  %41 = load double, ptr %20, align 8, !tbaa !5
  %42 = fneg double %41
  %43 = getelementptr inbounds i8, ptr %0, i64 8
  store double %42, ptr %43, align 8, !tbaa !5
  %44 = load double, ptr %16, align 8, !tbaa !5
  %45 = fmul double %44, %10
  %46 = fmul double %45, %7
  %47 = fmul double %27, %33
  %48 = getelementptr inbounds i8, ptr %15, i64 8
  %49 = load double, ptr %48, align 8, !tbaa !5
  %50 = fdiv double 0xC0297F2D647C2D27, %9
  %51 = call double @exp(double noundef %50) #6
  %52 = fmul double %51, 0x3FFF400000000000
  %53 = call double @llvm.fmuladd.f64(double %49, double 0x3FFAA2F78F1B4CC7, double %52)
  %54 = fmul double %47, %53
  %55 = call double @llvm.fmuladd.f64(double %46, double %7, double %54)
  store double %55, ptr %4, align 8, !tbaa !5
  %56 = getelementptr inbounds i8, ptr %17, i64 8
  %57 = load double, ptr %56, align 8, !tbaa !5
  %58 = load double, ptr %20, align 8, !tbaa !5
  %59 = fadd double %57, %58
  %60 = fadd double %28, %59
  %61 = getelementptr inbounds i8, ptr %0, i64 16
  %62 = getelementptr inbounds i8, ptr %0, i64 24
  store double %60, ptr %62, align 8, !tbaa !5
  %63 = fmul double %58, -3.000000e+00
  store double %63, ptr %61, align 8, !tbaa !5
  %64 = getelementptr inbounds i8, ptr %16, i64 8
  %65 = load double, ptr %64, align 8, !tbaa !5
  %66 = fmul double %65, %10
  %67 = fmul double %66, %7
  %68 = load double, ptr %15, align 8, !tbaa !5
  %69 = fmul double %47, %68
  %70 = call double @llvm.fmuladd.f64(double %67, double %7, double %69)
  %71 = getelementptr inbounds i8, ptr %4, i64 8
  store double %70, ptr %71, align 8, !tbaa !5
  %72 = getelementptr inbounds i8, ptr %13, i64 4976
  %73 = getelementptr inbounds i8, ptr %13, i64 7464
  %74 = getelementptr inbounds i8, ptr %13, i64 9952
  %75 = getelementptr inbounds i8, ptr %2, i64 8
  %76 = getelementptr inbounds i8, ptr %1, i64 8
  br label %77

77:                                               ; preds = %19, %77
  %indvars.iv = phi i64 [ 0, %19 ], [ %indvars.iv.next, %77 ]
  %78 = getelementptr inbounds [311 x double], ptr %13, i64 0, i64 %indvars.iv
  %79 = load double, ptr %78, align 8, !tbaa !5
  %80 = fadd double %79, 0xC02465BDE9FCF0EC
  %81 = fdiv double %80, %9
  %82 = call double @exp(double noundef %81) #6
  %83 = getelementptr inbounds [311 x double], ptr %72, i64 0, i64 %indvars.iv
  %84 = load double, ptr %83, align 8, !tbaa !5
  %85 = fneg double %84
  %86 = load double, ptr %78, align 8, !tbaa !5
  %87 = fadd double %86, 0xC02465BDE9FCF0EC
  %88 = fdiv double %87, %9
  %89 = call double @exp(double noundef %88) #6
  %90 = fadd double %89, -1.000000e+00
  %91 = call double @llvm.fabs.f64(double %90)
  %92 = fdiv double %85, %91
  %93 = load ptr, ptr %2, align 8, !tbaa !9
  %94 = getelementptr inbounds double, ptr %93, i64 %indvars.iv
  store double %92, ptr %94, align 8, !tbaa !5
  %95 = load double, ptr %0, align 8, !tbaa !5
  %96 = fsub double %95, %92
  store double %96, ptr %0, align 8, !tbaa !5
  %97 = load double, ptr %94, align 8, !tbaa !5
  %98 = fmul double %82, %97
  %99 = load ptr, ptr %1, align 8, !tbaa !9
  %100 = getelementptr inbounds double, ptr %99, i64 %indvars.iv
  store double %98, ptr %100, align 8, !tbaa !5
  %101 = call double @exp(double noundef %34) #6
  %102 = fdiv double %101, -3.000000e+00
  %103 = getelementptr inbounds [311 x double], ptr %73, i64 0, i64 %indvars.iv
  %104 = load double, ptr %103, align 8, !tbaa !5
  %105 = fmul double %102, %104
  %106 = load double, ptr %78, align 8, !tbaa !5
  %107 = fadd double %106, 0xC0282CBB284839FB
  %108 = fdiv double %107, %9
  %109 = call double @exp(double noundef %108) #6
  %110 = fadd double %109, -1.000000e+00
  %111 = call double @llvm.fabs.f64(double %110)
  %112 = fdiv double %105, %111
  %113 = call double @exp(double noundef %36) #6
  %114 = fdiv double %113, 3.000000e+00
  %115 = getelementptr inbounds [311 x double], ptr %74, i64 0, i64 %indvars.iv
  %116 = load double, ptr %115, align 8, !tbaa !5
  %117 = fmul double %114, %116
  %118 = load double, ptr %78, align 8, !tbaa !5
  %119 = fadd double %118, 0xC0297F2D647C2D27
  %120 = fdiv double %119, %9
  %121 = call double @exp(double noundef %120) #6
  %122 = fadd double %121, -1.000000e+00
  %123 = call double @llvm.fabs.f64(double %122)
  %124 = fdiv double %117, %123
  %125 = fsub double %112, %124
  %126 = load ptr, ptr %75, align 8, !tbaa !9
  %127 = getelementptr inbounds double, ptr %126, i64 %indvars.iv
  store double %125, ptr %127, align 8, !tbaa !5
  %128 = load double, ptr %62, align 8, !tbaa !5
  %129 = fsub double %128, %125
  store double %129, ptr %62, align 8, !tbaa !5
  %130 = load double, ptr %127, align 8, !tbaa !5
  %131 = fmul double %130, 3.000000e+00
  %132 = fmul double %82, %131
  %133 = load ptr, ptr %76, align 8, !tbaa !9
  %134 = getelementptr inbounds double, ptr %133, i64 %indvars.iv
  store double %132, ptr %134, align 8, !tbaa !5
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, 311
  br i1 %exitcond.not, label %135, label %77, !llvm.loop !11

135:                                              ; preds = %77
  %136 = getelementptr inbounds i8, ptr %13, i64 2488
  call void @populate_Diffusion(ptr noundef %23, ptr noundef %24, ptr noundef nonnull %21, ptr noundef nonnull %22, double noundef %8, ptr noundef nonnull %13, ptr noundef nonnull %136)
  %137 = load double, ptr %22, align 8, !tbaa !5
  %138 = load double, ptr %21, align 8, !tbaa !5
  %139 = fadd double %137, %138
  %140 = load double, ptr %62, align 8, !tbaa !5
  %141 = fadd double %140, %139
  store double %141, ptr %62, align 8, !tbaa !5
  %invariant.gep = getelementptr i8, ptr %23, i64 -8
  %142 = getelementptr inbounds i8, ptr %3, i64 8
  %143 = getelementptr inbounds i8, ptr %3, i64 16
  %144 = getelementptr inbounds i8, ptr %24, i64 1128
  %145 = getelementptr inbounds i8, ptr %23, i64 1104
  br label %146

146:                                              ; preds = %135, %237
  %indvars.iv178 = phi i64 [ 0, %135 ], [ %indvars.iv.next179, %237 ]
  %147 = load ptr, ptr %1, align 8, !tbaa !9
  %148 = getelementptr inbounds double, ptr %147, i64 %indvars.iv178
  %149 = load double, ptr %148, align 8, !tbaa !5
  %150 = load ptr, ptr %76, align 8, !tbaa !9
  %151 = getelementptr inbounds double, ptr %150, i64 %indvars.iv178
  %152 = load double, ptr %151, align 8, !tbaa !5
  %153 = fadd double %149, %152
  %154 = getelementptr inbounds double, ptr %23, i64 %indvars.iv178
  %155 = load double, ptr %154, align 8, !tbaa !5
  %156 = fsub double %155, %153
  %157 = getelementptr inbounds double, ptr %24, i64 %indvars.iv178
  %158 = load double, ptr %157, align 8, !tbaa !5
  %159 = fadd double %156, %158
  %160 = trunc i64 %indvars.iv178 to i32
  %161 = add i32 %160, -100
  %or.cond = icmp ult i32 %161, 39
  %162 = trunc i64 %indvars.iv178 to i32
  %163 = add i32 %162, -141
  %or.cond3 = icmp ult i32 %163, 39
  %or.cond174 = or i1 %or.cond, %or.cond3
  br i1 %or.cond174, label %164, label %174

164:                                              ; preds = %146
  %gep = getelementptr double, ptr %invariant.gep, i64 %indvars.iv178
  %165 = load double, ptr %gep, align 8, !tbaa !5
  %166 = fneg double %165
  %167 = load ptr, ptr %142, align 8, !tbaa !9
  %168 = getelementptr inbounds double, ptr %167, i64 %indvars.iv178
  store double %166, ptr %168, align 8, !tbaa !5
  %169 = getelementptr inbounds i8, ptr %157, i64 8
  %170 = load double, ptr %169, align 8, !tbaa !5
  %171 = fneg double %170
  %172 = load ptr, ptr %143, align 8, !tbaa !9
  %173 = getelementptr inbounds double, ptr %172, i64 %indvars.iv178
  store double %171, ptr %173, align 8, !tbaa !5
  br label %174

174:                                              ; preds = %146, %164
  %175 = trunc nuw nsw i64 %indvars.iv178 to i32
  switch i32 %175, label %195 [
    i32 139, label %176
    i32 140, label %181
  ]

176:                                              ; preds = %174
  %177 = load double, ptr %145, align 8, !tbaa !5
  %178 = fneg double %177
  %179 = load ptr, ptr %142, align 8, !tbaa !9
  %180 = getelementptr inbounds double, ptr %179, i64 %indvars.iv178
  store double %178, ptr %180, align 8, !tbaa !5
  br label %.sink.split

181:                                              ; preds = %174
  %182 = load ptr, ptr %142, align 8, !tbaa !9
  %183 = getelementptr inbounds double, ptr %182, i64 %indvars.iv178
  store double 0.000000e+00, ptr %183, align 8, !tbaa !5
  %184 = load double, ptr %144, align 8, !tbaa !5
  %185 = fneg double %184
  br label %.sink.split

.sink.split:                                      ; preds = %181, %176
  %.sink192 = phi double [ %185, %181 ], [ 0.000000e+00, %176 ]
  %.sink189.in = phi ptr [ %21, %181 ], [ %22, %176 ]
  %.sink185 = phi ptr [ %157, %181 ], [ %154, %176 ]
  %186 = load ptr, ptr %143, align 8, !tbaa !9
  %187 = getelementptr inbounds double, ptr %186, i64 %indvars.iv178
  store double %.sink192, ptr %187, align 8, !tbaa !5
  %.sink189 = load double, ptr %.sink189.in, align 8, !tbaa !5
  %188 = load ptr, ptr %75, align 8, !tbaa !9
  %189 = getelementptr inbounds double, ptr %188, i64 %indvars.iv178
  %190 = load double, ptr %189, align 8, !tbaa !5
  %191 = fsub double %190, %.sink189
  store double %191, ptr %189, align 8, !tbaa !5
  %192 = load double, ptr %.sink185, align 8, !tbaa !5
  %193 = load double, ptr %151, align 8, !tbaa !5
  %194 = fsub double %193, %192
  store double %194, ptr %151, align 8, !tbaa !5
  br label %195

195:                                              ; preds = %.sink.split, %174
  %196 = fmul double %27, %159
  %197 = getelementptr inbounds [311 x double], ptr %13, i64 0, i64 %indvars.iv178
  %198 = load double, ptr %197, align 8, !tbaa !5
  %199 = fdiv double 0x3F2040380F7C715F, %198
  %200 = call double @cube(double noundef %199)
  %201 = fmul double %196, %200
  %202 = fmul double %201, %10
  %203 = fmul double %202, 1.250000e-01
  %204 = fdiv double %203, 0x400921FB54442D18
  %205 = fdiv double %204, %11
  %206 = getelementptr inbounds double, ptr %6, i64 %indvars.iv178
  store double %205, ptr %206, align 8, !tbaa !5
  %207 = fcmp ogt double %205, 1.000000e-30
  br i1 %207, label %208, label %223

208:                                              ; preds = %195
  %209 = fneg double %205
  %210 = call double @exp(double noundef %209) #6
  %211 = fsub double 1.000000e+00, %210
  %212 = load double, ptr %206, align 8, !tbaa !5
  %213 = fdiv double %211, %212
  %214 = fsub double 1.000000e+00, %213
  %215 = fdiv double %159, %214
  %216 = load ptr, ptr %3, align 8, !tbaa !9
  %217 = getelementptr inbounds double, ptr %216, i64 %indvars.iv178
  store double %215, ptr %217, align 8, !tbaa !5
  %218 = fmul double %27, %215
  %219 = getelementptr inbounds double, ptr %14, i64 %indvars.iv178
  %220 = load double, ptr %219, align 8, !tbaa !5
  %221 = fmul double %220, %218
  %222 = fmul double %213, %221
  br label %237

223:                                              ; preds = %195
  %224 = load ptr, ptr %3, align 8, !tbaa !9
  %225 = getelementptr inbounds double, ptr %224, i64 %indvars.iv178
  store double 1.000000e+00, ptr %225, align 8, !tbaa !5
  %226 = load ptr, ptr %75, align 8, !tbaa !9
  %227 = getelementptr inbounds double, ptr %226, i64 %indvars.iv178
  store double 0.000000e+00, ptr %227, align 8, !tbaa !5
  %228 = load ptr, ptr %2, align 8, !tbaa !9
  %229 = getelementptr inbounds double, ptr %228, i64 %indvars.iv178
  store double 0.000000e+00, ptr %229, align 8, !tbaa !5
  %230 = load ptr, ptr %76, align 8, !tbaa !9
  %231 = getelementptr inbounds double, ptr %230, i64 %indvars.iv178
  store double 0.000000e+00, ptr %231, align 8, !tbaa !5
  %232 = load ptr, ptr %1, align 8, !tbaa !9
  %233 = getelementptr inbounds double, ptr %232, i64 %indvars.iv178
  store double 0.000000e+00, ptr %233, align 8, !tbaa !5
  %234 = getelementptr inbounds double, ptr %14, i64 %indvars.iv178
  %235 = load double, ptr %234, align 8, !tbaa !5
  %236 = fmul double %27, %235
  br label %237

237:                                              ; preds = %208, %223
  %.sink = phi double [ %222, %208 ], [ %236, %223 ]
  %238 = getelementptr inbounds double, ptr %5, i64 %indvars.iv178
  store double %.sink, ptr %238, align 8
  %indvars.iv.next179 = add nuw nsw i64 %indvars.iv178, 1
  %exitcond181.not = icmp eq i64 %indvars.iv.next179, 311
  br i1 %exitcond181.not, label %239, label %146, !llvm.loop !13

239:                                              ; preds = %237
  call void @free(ptr noundef nonnull %23) #6
  call void @free(ptr noundef nonnull %24) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %22) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %21) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %20) #6
  ret void
}

declare ptr @create_1D_array(i32 noundef) local_unnamed_addr #4

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

; Function Attrs: mustprogress nounwind willreturn allockind("free") memory(argmem: readwrite, inaccessiblemem: readwrite)
declare void @free(ptr allocptr nocapture noundef) local_unnamed_addr #5

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 = { 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 = { "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 willreturn allockind("free") memory(argmem: readwrite, inaccessiblemem: readwrite) "alloc-family"="malloc" "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 = !{!"double", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"any pointer", !7, i64 0}
!11 = distinct !{!11, !12}
!12 = !{!"llvm.loop.mustprogress"}
!13 = distinct !{!13, !12}


; Effective O3 passes: CorrelatedValuePropagationPass,EarlyCSEPass,GVNPass,GlobalOptPass,IPSCCPPass,IndVarSimplifyPass,InferFunctionAttrsPass,InstCombinePass,LICMPass,LoopRotatePass,MergedLoadStoreMotionPass,PostOrderFunctionAttrsPass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass