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

; 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 double @square(double noundef) #3

; Function Attrs: nounwind uwtable
define dso_local void @populate_Diffusion(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3, double noundef %4, ptr noundef %5, ptr noundef %6) #3 {
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca ptr, align 8
  %11 = alloca ptr, align 8
  %12 = alloca double, align 8
  %13 = alloca ptr, align 8
  %14 = alloca ptr, align 8
  %15 = alloca i32, align 4
  %16 = alloca double, align 8
  store ptr %0, ptr %8, align 8, !tbaa !5
  store ptr %1, ptr %9, align 8, !tbaa !5
  store ptr %2, ptr %10, align 8, !tbaa !5
  store ptr %3, ptr %11, align 8, !tbaa !5
  store double %4, ptr %12, align 8, !tbaa !9
  store ptr %5, ptr %13, align 8, !tbaa !5
  store ptr %6, ptr %14, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %15) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #4
  %17 = load double, ptr %12, align 8, !tbaa !9
  %18 = fmul double 0x406A00E1BF70A905, %17
  %19 = fdiv double %18, 0x41CBFA588B27AE14
  store double %19, ptr %16, align 8, !tbaa !9
  store i32 100, ptr %15, align 4, !tbaa !11
  %20 = load double, ptr %16, align 8, !tbaa !9
  %21 = load ptr, ptr %13, align 8, !tbaa !5
  %22 = load i32, ptr %15, align 4, !tbaa !11
  %23 = add i32 %22, 1
  %24 = zext i32 %23 to i64
  %25 = getelementptr inbounds double, ptr %21, i64 %24
  %26 = load double, ptr %25, align 8, !tbaa !9
  %27 = load ptr, ptr %13, align 8, !tbaa !5
  %28 = load i32, ptr %15, align 4, !tbaa !11
  %29 = zext i32 %28 to i64
  %30 = getelementptr inbounds double, ptr %27, i64 %29
  %31 = load double, ptr %30, align 8, !tbaa !9
  %32 = fsub double %26, %31
  %33 = call double @square(double noundef %32)
  %34 = fdiv double %20, %33
  %35 = load ptr, ptr %14, align 8, !tbaa !5
  %36 = load i32, ptr %15, align 4, !tbaa !11
  %37 = zext i32 %36 to i64
  %38 = getelementptr inbounds double, ptr %35, i64 %37
  %39 = load double, ptr %38, align 8, !tbaa !9
  %40 = fmul double %34, %39
  %41 = load ptr, ptr %8, align 8, !tbaa !5
  %42 = load i32, ptr %15, align 4, !tbaa !11
  %43 = zext i32 %42 to i64
  %44 = getelementptr inbounds double, ptr %41, i64 %43
  store double %40, ptr %44, align 8, !tbaa !9
  store i32 101, ptr %15, align 4, !tbaa !11
  br label %45

45:                                               ; preds = %120, %7
  %46 = load i32, ptr %15, align 4, !tbaa !11
  %47 = icmp ult i32 %46, 139
  br i1 %47, label %48, label %123

48:                                               ; preds = %45
  %49 = load ptr, ptr %13, align 8, !tbaa !5
  %50 = load i32, ptr %15, align 4, !tbaa !11
  %51 = zext i32 %50 to i64
  %52 = getelementptr inbounds double, ptr %49, i64 %51
  %53 = load double, ptr %52, align 8, !tbaa !9
  %54 = load ptr, ptr %13, align 8, !tbaa !5
  %55 = load i32, ptr %15, align 4, !tbaa !11
  %56 = sub i32 %55, 1
  %57 = zext i32 %56 to i64
  %58 = getelementptr inbounds double, ptr %54, i64 %57
  %59 = load double, ptr %58, align 8, !tbaa !9
  %60 = fsub double %53, %59
  %61 = load double, ptr %12, align 8, !tbaa !9
  %62 = fdiv double %60, %61
  %63 = call double @exp(double noundef %62) #4
  %64 = load ptr, ptr %8, align 8, !tbaa !5
  %65 = load i32, ptr %15, align 4, !tbaa !11
  %66 = sub i32 %65, 1
  %67 = zext i32 %66 to i64
  %68 = getelementptr inbounds double, ptr %64, i64 %67
  %69 = load double, ptr %68, align 8, !tbaa !9
  %70 = fmul double %63, %69
  %71 = load ptr, ptr %9, align 8, !tbaa !5
  %72 = load i32, ptr %15, align 4, !tbaa !11
  %73 = zext i32 %72 to i64
  %74 = getelementptr inbounds double, ptr %71, i64 %73
  store double %70, ptr %74, align 8, !tbaa !9
  %75 = load double, ptr %16, align 8, !tbaa !9
  %76 = load ptr, ptr %14, align 8, !tbaa !5
  %77 = load i32, ptr %15, align 4, !tbaa !11
  %78 = zext i32 %77 to i64
  %79 = getelementptr inbounds double, ptr %76, i64 %78
  %80 = load double, ptr %79, align 8, !tbaa !9
  %81 = load ptr, ptr %13, align 8, !tbaa !5
  %82 = load i32, ptr %15, align 4, !tbaa !11
  %83 = zext i32 %82 to i64
  %84 = getelementptr inbounds double, ptr %81, i64 %83
  %85 = load double, ptr %84, align 8, !tbaa !9
  %86 = load ptr, ptr %13, align 8, !tbaa !5
  %87 = load i32, ptr %15, align 4, !tbaa !11
  %88 = sub i32 %87, 1
  %89 = zext i32 %88 to i64
  %90 = getelementptr inbounds double, ptr %86, i64 %89
  %91 = load double, ptr %90, align 8, !tbaa !9
  %92 = fsub double %85, %91
  %93 = call double @square(double noundef %92)
  %94 = load ptr, ptr %9, align 8, !tbaa !5
  %95 = load i32, ptr %15, align 4, !tbaa !11
  %96 = zext i32 %95 to i64
  %97 = getelementptr inbounds double, ptr %94, i64 %96
  %98 = load double, ptr %97, align 8, !tbaa !9
  %99 = fmul double %93, %98
  %100 = fneg double %99
  %101 = call double @llvm.fmuladd.f64(double %75, double %80, double %100)
  %102 = load ptr, ptr %13, align 8, !tbaa !5
  %103 = load i32, ptr %15, align 4, !tbaa !11
  %104 = add i32 %103, 1
  %105 = zext i32 %104 to i64
  %106 = getelementptr inbounds double, ptr %102, i64 %105
  %107 = load double, ptr %106, align 8, !tbaa !9
  %108 = load ptr, ptr %13, align 8, !tbaa !5
  %109 = load i32, ptr %15, align 4, !tbaa !11
  %110 = zext i32 %109 to i64
  %111 = getelementptr inbounds double, ptr %108, i64 %110
  %112 = load double, ptr %111, align 8, !tbaa !9
  %113 = fsub double %107, %112
  %114 = call double @square(double noundef %113)
  %115 = fdiv double %101, %114
  %116 = load ptr, ptr %8, align 8, !tbaa !5
  %117 = load i32, ptr %15, align 4, !tbaa !11
  %118 = zext i32 %117 to i64
  %119 = getelementptr inbounds double, ptr %116, i64 %118
  store double %115, ptr %119, align 8, !tbaa !9
  br label %120

120:                                              ; preds = %48
  %121 = load i32, ptr %15, align 4, !tbaa !11
  %122 = add i32 %121, 1
  store i32 %122, ptr %15, align 4, !tbaa !11
  br label %45, !llvm.loop !13

123:                                              ; preds = %45
  store i32 139, ptr %15, align 4, !tbaa !11
  %124 = load ptr, ptr %13, align 8, !tbaa !5
  %125 = load i32, ptr %15, align 4, !tbaa !11
  %126 = zext i32 %125 to i64
  %127 = getelementptr inbounds double, ptr %124, i64 %126
  %128 = load double, ptr %127, align 8, !tbaa !9
  %129 = load ptr, ptr %13, align 8, !tbaa !5
  %130 = load i32, ptr %15, align 4, !tbaa !11
  %131 = sub i32 %130, 1
  %132 = zext i32 %131 to i64
  %133 = getelementptr inbounds double, ptr %129, i64 %132
  %134 = load double, ptr %133, align 8, !tbaa !9
  %135 = fsub double %128, %134
  %136 = load double, ptr %12, align 8, !tbaa !9
  %137 = fdiv double %135, %136
  %138 = call double @exp(double noundef %137) #4
  %139 = load ptr, ptr %8, align 8, !tbaa !5
  %140 = load i32, ptr %15, align 4, !tbaa !11
  %141 = sub i32 %140, 1
  %142 = zext i32 %141 to i64
  %143 = getelementptr inbounds double, ptr %139, i64 %142
  %144 = load double, ptr %143, align 8, !tbaa !9
  %145 = fmul double %138, %144
  %146 = load ptr, ptr %9, align 8, !tbaa !5
  %147 = load i32, ptr %15, align 4, !tbaa !11
  %148 = zext i32 %147 to i64
  %149 = getelementptr inbounds double, ptr %146, i64 %148
  store double %145, ptr %149, align 8, !tbaa !9
  %150 = load double, ptr %16, align 8, !tbaa !9
  %151 = load ptr, ptr %14, align 8, !tbaa !5
  %152 = load i32, ptr %15, align 4, !tbaa !11
  %153 = zext i32 %152 to i64
  %154 = getelementptr inbounds double, ptr %151, i64 %153
  %155 = load double, ptr %154, align 8, !tbaa !9
  %156 = load ptr, ptr %13, align 8, !tbaa !5
  %157 = load i32, ptr %15, align 4, !tbaa !11
  %158 = zext i32 %157 to i64
  %159 = getelementptr inbounds double, ptr %156, i64 %158
  %160 = load double, ptr %159, align 8, !tbaa !9
  %161 = load ptr, ptr %13, align 8, !tbaa !5
  %162 = load i32, ptr %15, align 4, !tbaa !11
  %163 = sub i32 %162, 1
  %164 = zext i32 %163 to i64
  %165 = getelementptr inbounds double, ptr %161, i64 %164
  %166 = load double, ptr %165, align 8, !tbaa !9
  %167 = fsub double %160, %166
  %168 = call double @square(double noundef %167)
  %169 = load ptr, ptr %9, align 8, !tbaa !5
  %170 = load i32, ptr %15, align 4, !tbaa !11
  %171 = zext i32 %170 to i64
  %172 = getelementptr inbounds double, ptr %169, i64 %171
  %173 = load double, ptr %172, align 8, !tbaa !9
  %174 = fmul double %168, %173
  %175 = fneg double %174
  %176 = call double @llvm.fmuladd.f64(double %150, double %155, double %175)
  %177 = load ptr, ptr %13, align 8, !tbaa !5
  %178 = load i32, ptr %15, align 4, !tbaa !11
  %179 = zext i32 %178 to i64
  %180 = getelementptr inbounds double, ptr %177, i64 %179
  %181 = load double, ptr %180, align 8, !tbaa !9
  %182 = fsub double 0x402465BDE9FCF0EC, %181
  %183 = call double @square(double noundef %182)
  %184 = fdiv double %176, %183
  %185 = load ptr, ptr %8, align 8, !tbaa !5
  %186 = load i32, ptr %15, align 4, !tbaa !11
  %187 = zext i32 %186 to i64
  %188 = getelementptr inbounds double, ptr %185, i64 %187
  store double %184, ptr %188, align 8, !tbaa !9
  %189 = load ptr, ptr %13, align 8, !tbaa !5
  %190 = load i32, ptr %15, align 4, !tbaa !11
  %191 = zext i32 %190 to i64
  %192 = getelementptr inbounds double, ptr %189, i64 %191
  %193 = load double, ptr %192, align 8, !tbaa !9
  %194 = fsub double 0x402465BDE9FCF0EC, %193
  %195 = load double, ptr %12, align 8, !tbaa !9
  %196 = fdiv double %194, %195
  %197 = call double @exp(double noundef %196) #4
  %198 = fdiv double %197, 3.000000e+00
  %199 = load ptr, ptr %8, align 8, !tbaa !5
  %200 = load i32, ptr %15, align 4, !tbaa !11
  %201 = zext i32 %200 to i64
  %202 = getelementptr inbounds double, ptr %199, i64 %201
  %203 = load double, ptr %202, align 8, !tbaa !9
  %204 = fmul double %198, %203
  %205 = load ptr, ptr %11, align 8, !tbaa !5
  store double %204, ptr %205, align 8, !tbaa !9
  store i32 179, ptr %15, align 4, !tbaa !11
  %206 = load double, ptr %16, align 8, !tbaa !9
  %207 = load ptr, ptr %13, align 8, !tbaa !5
  %208 = load i32, ptr %15, align 4, !tbaa !11
  %209 = zext i32 %208 to i64
  %210 = getelementptr inbounds double, ptr %207, i64 %209
  %211 = load double, ptr %210, align 8, !tbaa !9
  %212 = load ptr, ptr %13, align 8, !tbaa !5
  %213 = load i32, ptr %15, align 4, !tbaa !11
  %214 = sub i32 %213, 1
  %215 = zext i32 %214 to i64
  %216 = getelementptr inbounds double, ptr %212, i64 %215
  %217 = load double, ptr %216, align 8, !tbaa !9
  %218 = fsub double %211, %217
  %219 = call double @square(double noundef %218)
  %220 = fdiv double %206, %219
  %221 = load ptr, ptr %14, align 8, !tbaa !5
  %222 = load i32, ptr %15, align 4, !tbaa !11
  %223 = zext i32 %222 to i64
  %224 = getelementptr inbounds double, ptr %221, i64 %223
  %225 = load double, ptr %224, align 8, !tbaa !9
  %226 = fmul double %220, %225
  %227 = load ptr, ptr %9, align 8, !tbaa !5
  %228 = load i32, ptr %15, align 4, !tbaa !11
  %229 = zext i32 %228 to i64
  %230 = getelementptr inbounds double, ptr %227, i64 %229
  store double %226, ptr %230, align 8, !tbaa !9
  store i32 178, ptr %15, align 4, !tbaa !11
  br label %231

231:                                              ; preds = %306, %123
  %232 = load i32, ptr %15, align 4, !tbaa !11
  %233 = icmp ugt i32 %232, 140
  br i1 %233, label %234, label %309

234:                                              ; preds = %231
  %235 = load ptr, ptr %13, align 8, !tbaa !5
  %236 = load i32, ptr %15, align 4, !tbaa !11
  %237 = zext i32 %236 to i64
  %238 = getelementptr inbounds double, ptr %235, i64 %237
  %239 = load double, ptr %238, align 8, !tbaa !9
  %240 = load ptr, ptr %13, align 8, !tbaa !5
  %241 = load i32, ptr %15, align 4, !tbaa !11
  %242 = add i32 %241, 1
  %243 = zext i32 %242 to i64
  %244 = getelementptr inbounds double, ptr %240, i64 %243
  %245 = load double, ptr %244, align 8, !tbaa !9
  %246 = fsub double %239, %245
  %247 = load double, ptr %12, align 8, !tbaa !9
  %248 = fdiv double %246, %247
  %249 = call double @exp(double noundef %248) #4
  %250 = load ptr, ptr %9, align 8, !tbaa !5
  %251 = load i32, ptr %15, align 4, !tbaa !11
  %252 = add i32 %251, 1
  %253 = zext i32 %252 to i64
  %254 = getelementptr inbounds double, ptr %250, i64 %253
  %255 = load double, ptr %254, align 8, !tbaa !9
  %256 = fmul double %249, %255
  %257 = load ptr, ptr %8, align 8, !tbaa !5
  %258 = load i32, ptr %15, align 4, !tbaa !11
  %259 = zext i32 %258 to i64
  %260 = getelementptr inbounds double, ptr %257, i64 %259
  store double %256, ptr %260, align 8, !tbaa !9
  %261 = load double, ptr %16, align 8, !tbaa !9
  %262 = load ptr, ptr %14, align 8, !tbaa !5
  %263 = load i32, ptr %15, align 4, !tbaa !11
  %264 = zext i32 %263 to i64
  %265 = getelementptr inbounds double, ptr %262, i64 %264
  %266 = load double, ptr %265, align 8, !tbaa !9
  %267 = load ptr, ptr %13, align 8, !tbaa !5
  %268 = load i32, ptr %15, align 4, !tbaa !11
  %269 = add i32 %268, 1
  %270 = zext i32 %269 to i64
  %271 = getelementptr inbounds double, ptr %267, i64 %270
  %272 = load double, ptr %271, align 8, !tbaa !9
  %273 = load ptr, ptr %13, align 8, !tbaa !5
  %274 = load i32, ptr %15, align 4, !tbaa !11
  %275 = zext i32 %274 to i64
  %276 = getelementptr inbounds double, ptr %273, i64 %275
  %277 = load double, ptr %276, align 8, !tbaa !9
  %278 = fsub double %272, %277
  %279 = call double @square(double noundef %278)
  %280 = load ptr, ptr %8, align 8, !tbaa !5
  %281 = load i32, ptr %15, align 4, !tbaa !11
  %282 = zext i32 %281 to i64
  %283 = getelementptr inbounds double, ptr %280, i64 %282
  %284 = load double, ptr %283, align 8, !tbaa !9
  %285 = fmul double %279, %284
  %286 = fneg double %285
  %287 = call double @llvm.fmuladd.f64(double %261, double %266, double %286)
  %288 = load ptr, ptr %13, align 8, !tbaa !5
  %289 = load i32, ptr %15, align 4, !tbaa !11
  %290 = zext i32 %289 to i64
  %291 = getelementptr inbounds double, ptr %288, i64 %290
  %292 = load double, ptr %291, align 8, !tbaa !9
  %293 = load ptr, ptr %13, align 8, !tbaa !5
  %294 = load i32, ptr %15, align 4, !tbaa !11
  %295 = sub i32 %294, 1
  %296 = zext i32 %295 to i64
  %297 = getelementptr inbounds double, ptr %293, i64 %296
  %298 = load double, ptr %297, align 8, !tbaa !9
  %299 = fsub double %292, %298
  %300 = call double @square(double noundef %299)
  %301 = fdiv double %287, %300
  %302 = load ptr, ptr %9, align 8, !tbaa !5
  %303 = load i32, ptr %15, align 4, !tbaa !11
  %304 = zext i32 %303 to i64
  %305 = getelementptr inbounds double, ptr %302, i64 %304
  store double %301, ptr %305, align 8, !tbaa !9
  br label %306

306:                                              ; preds = %234
  %307 = load i32, ptr %15, align 4, !tbaa !11
  %308 = add i32 %307, -1
  store i32 %308, ptr %15, align 4, !tbaa !11
  br label %231, !llvm.loop !15

309:                                              ; preds = %231
  store i32 140, ptr %15, align 4, !tbaa !11
  %310 = load ptr, ptr %13, align 8, !tbaa !5
  %311 = load i32, ptr %15, align 4, !tbaa !11
  %312 = zext i32 %311 to i64
  %313 = getelementptr inbounds double, ptr %310, i64 %312
  %314 = load double, ptr %313, align 8, !tbaa !9
  %315 = load ptr, ptr %13, align 8, !tbaa !5
  %316 = load i32, ptr %15, align 4, !tbaa !11
  %317 = add i32 %316, 1
  %318 = zext i32 %317 to i64
  %319 = getelementptr inbounds double, ptr %315, i64 %318
  %320 = load double, ptr %319, align 8, !tbaa !9
  %321 = fsub double %314, %320
  %322 = load double, ptr %12, align 8, !tbaa !9
  %323 = fdiv double %321, %322
  %324 = call double @exp(double noundef %323) #4
  %325 = load ptr, ptr %9, align 8, !tbaa !5
  %326 = load i32, ptr %15, align 4, !tbaa !11
  %327 = add i32 %326, 1
  %328 = zext i32 %327 to i64
  %329 = getelementptr inbounds double, ptr %325, i64 %328
  %330 = load double, ptr %329, align 8, !tbaa !9
  %331 = fmul double %324, %330
  %332 = load ptr, ptr %8, align 8, !tbaa !5
  %333 = load i32, ptr %15, align 4, !tbaa !11
  %334 = zext i32 %333 to i64
  %335 = getelementptr inbounds double, ptr %332, i64 %334
  store double %331, ptr %335, align 8, !tbaa !9
  %336 = load double, ptr %16, align 8, !tbaa !9
  %337 = load ptr, ptr %14, align 8, !tbaa !5
  %338 = load i32, ptr %15, align 4, !tbaa !11
  %339 = zext i32 %338 to i64
  %340 = getelementptr inbounds double, ptr %337, i64 %339
  %341 = load double, ptr %340, align 8, !tbaa !9
  %342 = load ptr, ptr %13, align 8, !tbaa !5
  %343 = load i32, ptr %15, align 4, !tbaa !11
  %344 = add i32 %343, 1
  %345 = zext i32 %344 to i64
  %346 = getelementptr inbounds double, ptr %342, i64 %345
  %347 = load double, ptr %346, align 8, !tbaa !9
  %348 = load ptr, ptr %13, align 8, !tbaa !5
  %349 = load i32, ptr %15, align 4, !tbaa !11
  %350 = zext i32 %349 to i64
  %351 = getelementptr inbounds double, ptr %348, i64 %350
  %352 = load double, ptr %351, align 8, !tbaa !9
  %353 = fsub double %347, %352
  %354 = call double @square(double noundef %353)
  %355 = load ptr, ptr %8, align 8, !tbaa !5
  %356 = load i32, ptr %15, align 4, !tbaa !11
  %357 = zext i32 %356 to i64
  %358 = getelementptr inbounds double, ptr %355, i64 %357
  %359 = load double, ptr %358, align 8, !tbaa !9
  %360 = fmul double %354, %359
  %361 = fneg double %360
  %362 = call double @llvm.fmuladd.f64(double %336, double %341, double %361)
  %363 = load ptr, ptr %13, align 8, !tbaa !5
  %364 = load i32, ptr %15, align 4, !tbaa !11
  %365 = zext i32 %364 to i64
  %366 = getelementptr inbounds double, ptr %363, i64 %365
  %367 = load double, ptr %366, align 8, !tbaa !9
  %368 = fsub double %367, 0x402465BDE9FCF0EC
  %369 = call double @square(double noundef %368)
  %370 = fdiv double %362, %369
  %371 = load ptr, ptr %9, align 8, !tbaa !5
  %372 = load i32, ptr %15, align 4, !tbaa !11
  %373 = zext i32 %372 to i64
  %374 = getelementptr inbounds double, ptr %371, i64 %373
  store double %370, ptr %374, align 8, !tbaa !9
  %375 = load ptr, ptr %13, align 8, !tbaa !5
  %376 = load i32, ptr %15, align 4, !tbaa !11
  %377 = zext i32 %376 to i64
  %378 = getelementptr inbounds double, ptr %375, i64 %377
  %379 = load double, ptr %378, align 8, !tbaa !9
  %380 = fsub double 0x402465BDE9FCF0EC, %379
  %381 = load double, ptr %12, align 8, !tbaa !9
  %382 = fdiv double %380, %381
  %383 = call double @exp(double noundef %382) #4
  %384 = fdiv double %383, 3.000000e+00
  %385 = load ptr, ptr %9, align 8, !tbaa !5
  %386 = load i32, ptr %15, align 4, !tbaa !11
  %387 = zext i32 %386 to i64
  %388 = getelementptr inbounds double, ptr %385, i64 %387
  %389 = load double, ptr %388, align 8, !tbaa !9
  %390 = fmul double %384, %389
  %391 = load ptr, ptr %10, align 8, !tbaa !5
  store double %390, ptr %391, align 8, !tbaa !9
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %15) #4
  ret void
}

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 = { 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 speculatable willreturn memory(none)
declare double @llvm.fmuladd.f64(double, double, double) #0

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

; Function Attrs: nounwind uwtable
declare dso_local double @square(double noundef) local_unnamed_addr #2

; Function Attrs: nounwind uwtable
define dso_local void @populate_Diffusion(ptr nocapture noundef %0, ptr nocapture noundef %1, ptr nocapture noundef writeonly %2, ptr nocapture noundef writeonly %3, double noundef %4, ptr nocapture noundef readonly %5, ptr nocapture noundef readonly %6) local_unnamed_addr #2 {
  %8 = fmul double %4, 0x406A00E1BF70A905
  %9 = fdiv double %8, 0x41CBFA588B27AE14
  %10 = getelementptr inbounds i8, ptr %5, i64 808
  %11 = load double, ptr %10, align 8, !tbaa !5
  %12 = getelementptr inbounds i8, ptr %5, i64 800
  %13 = load double, ptr %12, align 8, !tbaa !5
  %14 = fsub double %11, %13
  %15 = tail call double @square(double noundef %14)
  %16 = fdiv double %9, %15
  %17 = getelementptr inbounds i8, ptr %6, i64 800
  %18 = load double, ptr %17, align 8, !tbaa !5
  %19 = fmul double %16, %18
  %20 = getelementptr inbounds i8, ptr %0, i64 800
  store double %19, ptr %20, align 8, !tbaa !5
  br label %21

21:                                               ; preds = %7, %21
  %indvars.iv = phi i64 [ 101, %7 ], [ %indvars.iv.next, %21 ]
  %22 = getelementptr inbounds double, ptr %5, i64 %indvars.iv
  %23 = load double, ptr %22, align 8, !tbaa !5
  %24 = add nsw i64 %indvars.iv, -1
  %25 = getelementptr inbounds double, ptr %5, i64 %24
  %26 = load double, ptr %25, align 8, !tbaa !5
  %27 = fsub double %23, %26
  %28 = fdiv double %27, %4
  %29 = tail call double @exp(double noundef %28) #3
  %30 = getelementptr inbounds double, ptr %0, i64 %24
  %31 = load double, ptr %30, align 8, !tbaa !5
  %32 = fmul double %29, %31
  %33 = getelementptr inbounds double, ptr %1, i64 %indvars.iv
  store double %32, ptr %33, align 8, !tbaa !5
  %34 = getelementptr inbounds double, ptr %6, i64 %indvars.iv
  %35 = load double, ptr %34, align 8, !tbaa !5
  %36 = load double, ptr %22, align 8, !tbaa !5
  %37 = load double, ptr %25, align 8, !tbaa !5
  %38 = fsub double %36, %37
  %39 = tail call double @square(double noundef %38)
  %40 = load double, ptr %33, align 8, !tbaa !5
  %41 = fneg double %39
  %42 = fmul double %40, %41
  %43 = tail call double @llvm.fmuladd.f64(double %9, double %35, double %42)
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %44 = getelementptr inbounds double, ptr %5, i64 %indvars.iv.next
  %45 = load double, ptr %44, align 8, !tbaa !5
  %46 = load double, ptr %22, align 8, !tbaa !5
  %47 = fsub double %45, %46
  %48 = tail call double @square(double noundef %47)
  %49 = fdiv double %43, %48
  %50 = getelementptr inbounds double, ptr %0, i64 %indvars.iv
  store double %49, ptr %50, align 8, !tbaa !5
  %exitcond.not = icmp eq i64 %indvars.iv.next, 139
  br i1 %exitcond.not, label %51, label %21, !llvm.loop !9

51:                                               ; preds = %21
  %52 = getelementptr inbounds i8, ptr %5, i64 1112
  %53 = load double, ptr %52, align 8, !tbaa !5
  %54 = getelementptr inbounds i8, ptr %5, i64 1104
  %55 = load double, ptr %54, align 8, !tbaa !5
  %56 = fsub double %53, %55
  %57 = fdiv double %56, %4
  %58 = tail call double @exp(double noundef %57) #3
  %59 = getelementptr inbounds i8, ptr %0, i64 1104
  %60 = load double, ptr %59, align 8, !tbaa !5
  %61 = fmul double %58, %60
  %62 = getelementptr inbounds i8, ptr %1, i64 1112
  store double %61, ptr %62, align 8, !tbaa !5
  %63 = getelementptr inbounds i8, ptr %6, i64 1112
  %64 = load double, ptr %63, align 8, !tbaa !5
  %65 = load double, ptr %52, align 8, !tbaa !5
  %66 = load double, ptr %54, align 8, !tbaa !5
  %67 = fsub double %65, %66
  %68 = tail call double @square(double noundef %67)
  %69 = load double, ptr %62, align 8, !tbaa !5
  %70 = fneg double %68
  %71 = fmul double %69, %70
  %72 = tail call double @llvm.fmuladd.f64(double %9, double %64, double %71)
  %73 = load double, ptr %52, align 8, !tbaa !5
  %74 = fsub double 0x402465BDE9FCF0EC, %73
  %75 = tail call double @square(double noundef %74)
  %76 = fdiv double %72, %75
  %77 = getelementptr inbounds i8, ptr %0, i64 1112
  store double %76, ptr %77, align 8, !tbaa !5
  %78 = load double, ptr %52, align 8, !tbaa !5
  %79 = fsub double 0x402465BDE9FCF0EC, %78
  %80 = fdiv double %79, %4
  %81 = tail call double @exp(double noundef %80) #3
  %82 = fdiv double %81, 3.000000e+00
  %83 = load double, ptr %77, align 8, !tbaa !5
  %84 = fmul double %82, %83
  store double %84, ptr %3, align 8, !tbaa !5
  %85 = getelementptr inbounds i8, ptr %5, i64 1432
  %86 = load double, ptr %85, align 8, !tbaa !5
  %87 = getelementptr inbounds i8, ptr %5, i64 1424
  %88 = load double, ptr %87, align 8, !tbaa !5
  %89 = fsub double %86, %88
  %90 = tail call double @square(double noundef %89)
  %91 = fdiv double %9, %90
  %92 = getelementptr inbounds i8, ptr %6, i64 1432
  %93 = load double, ptr %92, align 8, !tbaa !5
  %94 = fmul double %91, %93
  %95 = getelementptr inbounds i8, ptr %1, i64 1432
  store double %94, ptr %95, align 8, !tbaa !5
  br label %96

96:                                               ; preds = %51, %96
  %indvars.iv129 = phi i64 [ 178, %51 ], [ %indvars.iv.next130, %96 ]
  %97 = getelementptr inbounds double, ptr %5, i64 %indvars.iv129
  %98 = load double, ptr %97, align 8, !tbaa !5
  %99 = add nuw nsw i64 %indvars.iv129, 1
  %100 = getelementptr inbounds double, ptr %5, i64 %99
  %101 = load double, ptr %100, align 8, !tbaa !5
  %102 = fsub double %98, %101
  %103 = fdiv double %102, %4
  %104 = tail call double @exp(double noundef %103) #3
  %105 = getelementptr inbounds double, ptr %1, i64 %99
  %106 = load double, ptr %105, align 8, !tbaa !5
  %107 = fmul double %104, %106
  %108 = getelementptr inbounds double, ptr %0, i64 %indvars.iv129
  store double %107, ptr %108, align 8, !tbaa !5
  %109 = getelementptr inbounds double, ptr %6, i64 %indvars.iv129
  %110 = load double, ptr %109, align 8, !tbaa !5
  %111 = load double, ptr %100, align 8, !tbaa !5
  %112 = load double, ptr %97, align 8, !tbaa !5
  %113 = fsub double %111, %112
  %114 = tail call double @square(double noundef %113)
  %115 = load double, ptr %108, align 8, !tbaa !5
  %116 = fneg double %114
  %117 = fmul double %115, %116
  %118 = tail call double @llvm.fmuladd.f64(double %9, double %110, double %117)
  %119 = load double, ptr %97, align 8, !tbaa !5
  %indvars.iv.next130 = add nsw i64 %indvars.iv129, -1
  %120 = getelementptr inbounds double, ptr %5, i64 %indvars.iv.next130
  %121 = load double, ptr %120, align 8, !tbaa !5
  %122 = fsub double %119, %121
  %123 = tail call double @square(double noundef %122)
  %124 = fdiv double %118, %123
  %125 = getelementptr inbounds double, ptr %1, i64 %indvars.iv129
  store double %124, ptr %125, align 8, !tbaa !5
  %126 = icmp ugt i64 %indvars.iv.next130, 140
  br i1 %126, label %96, label %127, !llvm.loop !11

127:                                              ; preds = %96
  %128 = getelementptr inbounds i8, ptr %5, i64 1120
  %129 = load double, ptr %128, align 8, !tbaa !5
  %130 = getelementptr inbounds i8, ptr %5, i64 1128
  %131 = load double, ptr %130, align 8, !tbaa !5
  %132 = fsub double %129, %131
  %133 = fdiv double %132, %4
  %134 = tail call double @exp(double noundef %133) #3
  %135 = getelementptr inbounds i8, ptr %1, i64 1128
  %136 = load double, ptr %135, align 8, !tbaa !5
  %137 = fmul double %134, %136
  %138 = getelementptr inbounds i8, ptr %0, i64 1120
  store double %137, ptr %138, align 8, !tbaa !5
  %139 = getelementptr inbounds i8, ptr %6, i64 1120
  %140 = load double, ptr %139, align 8, !tbaa !5
  %141 = load double, ptr %130, align 8, !tbaa !5
  %142 = load double, ptr %128, align 8, !tbaa !5
  %143 = fsub double %141, %142
  %144 = tail call double @square(double noundef %143)
  %145 = load double, ptr %138, align 8, !tbaa !5
  %146 = fneg double %144
  %147 = fmul double %145, %146
  %148 = tail call double @llvm.fmuladd.f64(double %9, double %140, double %147)
  %149 = load double, ptr %128, align 8, !tbaa !5
  %150 = fadd double %149, 0xC02465BDE9FCF0EC
  %151 = tail call double @square(double noundef %150)
  %152 = fdiv double %148, %151
  %153 = getelementptr inbounds i8, ptr %1, i64 1120
  store double %152, ptr %153, align 8, !tbaa !5
  %154 = load double, ptr %128, align 8, !tbaa !5
  %155 = fsub double 0x402465BDE9FCF0EC, %154
  %156 = fdiv double %155, %4
  %157 = tail call double @exp(double noundef %156) #3
  %158 = fdiv double %157, 3.000000e+00
  %159 = load double, ptr %153, align 8, !tbaa !5
  %160 = fmul double %158, %159
  store double %160, ptr %2, align 8, !tbaa !5
  ret void
}

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 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { 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 = distinct !{!9, !10}
!10 = !{!"llvm.loop.mustprogress"}
!11 = distinct !{!11, !10}


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