; LLVM IR from CLASS-SZ/class_sz :: class-sz/source (perturbations)

; 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.background = type { double, double, double, double, double, double, double, i32, double, double, double, double, i16, double, double, double, double, double, double, double, double, double, i16, double, double, ptr, i32, i32, double, i32, ptr, ptr, double, ptr, double, ptr, double, ptr, double, ptr, ptr, ptr, double, double, double, double, i32, ptr, double, double, double, double, double, double, double, double, double, double, double, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, ptr, ptr, ptr, ptr, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, ptr, i16, i16, i16, i16, i16, i16, [2048 x i8], i16, [2048 x i8] }
%struct.perturbs = type { i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i32, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i32, i32, i32, i32, double, i32, i32, [100 x double], [100 x double], i32, i32, i32, i32, i32, double, i32, i32, [30 x double], double, double, double, ptr, ptr, i32, i16, i16, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i16, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, ptr, ptr, ptr, ptr, double, double, ptr, i32, double, double, double, ptr, ptr, ptr, ptr, ptr, ptr, i32, ptr, ptr, ptr, [8000 x i8], [8000 x i8], [8000 x i8], i32, i32, i32, [30 x ptr], [30 x ptr], [30 x ptr], [30 x i32], [30 x i32], [30 x i32], i16, [2048 x i8] }
%struct.thermo = type { double, i32, i32, i32, double, double, i16, i16, double, double, double, double, i32, ptr, ptr, double, i32, ptr, ptr, double, i32, ptr, ptr, double, i16, double, double, double, double, double, double, double, double, double, double, double, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, i32, ptr, ptr, ptr, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i16, i16, i16, [2048 x i8] }
%struct.precision = type { double, double, double, double, double, double, double, double, double, double, double, double, double, [256 x i8], double, i32, double, i32, i32, double, i32, double, i32, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i32, [256 x i8], [256 x i8], [256 x i8], double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i32, double, double, i32, i32, i32, i32, i32, i32, i32, i32, double, double, double, double, double, double, double, double, i32, double, double, i32, double, double, i32, double, i32, double, double, i32, double, double, double, i32, double, double, double, double, double, i32, double, double, double, double, double, double, double, i32, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, double, i32, double, double, double, double, i32, double, double, i32, i32, i32, double, double, [2048 x i8] }

@.str = external hidden unnamed_addr constant [18 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.4 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.7 = external hidden unnamed_addr constant [27 x i8], align 1
@.str.128 = external hidden unnamed_addr constant [9 x i8], align 1
@__func__.perturb_timesampling_for_sources = external hidden unnamed_addr constant [33 x i8], align 1
@.str.129 = external hidden unnamed_addr constant [11 x i8], align 1
@.str.130 = external hidden unnamed_addr constant [99 x i8], align 1
@.str.131 = external hidden unnamed_addr constant [125 x i8], align 1
@.str.132 = external hidden unnamed_addr constant [274 x i8], align 1
@.str.133 = external hidden unnamed_addr constant [128 x i8], align 1
@.str.134 = external hidden unnamed_addr constant [99 x i8], align 1
@.str.135 = external hidden unnamed_addr constant [171 x i8], align 1
@.str.136 = external hidden unnamed_addr constant [128 x i8], align 1
@.str.137 = external hidden unnamed_addr constant [97 x i8], align 1
@.str.138 = external hidden unnamed_addr constant [52 x i8], align 1
@.str.139 = external hidden unnamed_addr constant [97 x i8], align 1
@.str.140 = external hidden unnamed_addr constant [93 x i8], align 1
@.str.141 = external hidden unnamed_addr constant [125 x i8], align 1
@.str.142 = external hidden unnamed_addr constant [46 x i8], align 1
@.str.143 = external hidden unnamed_addr constant [23 x i8], align 1
@.str.144 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.145 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.146 = external hidden unnamed_addr constant [18 x i8], align 1
@.str.147 = external hidden unnamed_addr constant [33 x i8], align 1
@.str.148 = external hidden unnamed_addr constant [18 x i8], align 1
@.str.149 = external hidden unnamed_addr constant [50 x i8], align 1
@.str.150 = external hidden unnamed_addr constant [155 x i8], align 1
@.str.151 = external hidden unnamed_addr constant [44 x i8], align 1
@.str.152 = external hidden unnamed_addr constant [85 x i8], align 1
@.str.153 = external hidden unnamed_addr constant [12 x i8], align 1
@.str.154 = external hidden unnamed_addr constant [12 x i8], align 1
@.str.155 = external hidden unnamed_addr constant [65 x i8], align 1
@.str.156 = external hidden unnamed_addr constant [72 x i8], align 1

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

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

declare void @class_protect_sprintf(ptr noundef, ptr noundef, ...) #2

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

; Function Attrs: nounwind allocsize(0)
declare noalias ptr @malloc(i64 noundef) #3

declare i32 @background_tau_of_z(ptr noundef, double noundef, ptr noundef) #2

; Function Attrs: nounwind
declare void @free(ptr noundef) #1

; Function Attrs: nounwind uwtable
define dso_local i32 @perturb_timesampling_for_sources(ptr noundef %0, ptr noundef %1, ptr noundef %2, ptr noundef %3) #4 {
  %5 = alloca i32, align 4
  %6 = alloca ptr, align 8
  %7 = alloca ptr, align 8
  %8 = alloca ptr, align 8
  %9 = alloca ptr, align 8
  %10 = alloca i32, align 4
  %11 = alloca i32, align 4
  %12 = alloca i32, align 4
  %13 = alloca i32, align 4
  %14 = alloca i32, align 4
  %15 = alloca i32, align 4
  %16 = alloca i32, align 4
  %17 = alloca i32, align 4
  %18 = alloca i32, align 4
  %19 = alloca double, align 8
  %20 = alloca double, align 8
  %21 = alloca double, align 8
  %22 = alloca double, align 8
  %23 = alloca double, align 8
  %24 = alloca double, align 8
  %25 = alloca double, align 8
  %26 = alloca double, align 8
  %27 = alloca double, align 8
  %28 = alloca double, align 8
  %29 = alloca ptr, align 8
  %30 = alloca ptr, align 8
  %31 = alloca i32, align 4
  %32 = alloca [2048 x i8], align 16
  %33 = alloca i32, align 4
  %34 = alloca i32, align 4
  %35 = alloca [2048 x i8], align 16
  %36 = alloca [2048 x i8], align 16
  %37 = alloca [2048 x i8], align 16
  %38 = alloca [2048 x i8], align 16
  %39 = alloca [2048 x i8], align 16
  %40 = alloca [2048 x i8], align 16
  %41 = alloca [2048 x i8], align 16
  %42 = alloca [2048 x i8], align 16
  %43 = alloca [2048 x i8], align 16
  %44 = alloca [2048 x i8], align 16
  %45 = alloca [2048 x i8], align 16
  %46 = alloca [2048 x i8], align 16
  %47 = alloca [2048 x i8], align 16
  %48 = alloca [2048 x i8], align 16
  %49 = alloca [2048 x i8], align 16
  %50 = alloca [2048 x i8], align 16
  %51 = alloca [2048 x i8], align 16
  %52 = alloca [2048 x i8], align 16
  %53 = alloca [2048 x i8], align 16
  %54 = alloca [2048 x i8], align 16
  %55 = alloca i32, align 4
  %56 = alloca [2048 x i8], align 16
  %57 = alloca [2048 x i8], align 16
  %58 = alloca [2048 x i8], align 16
  %59 = alloca [2048 x i8], align 16
  %60 = alloca [2048 x i8], align 16
  %61 = alloca [2048 x i8], align 16
  %62 = alloca [2048 x i8], align 16
  %63 = alloca [2048 x i8], align 16
  %64 = alloca [2048 x i8], align 16
  %65 = alloca [2048 x i8], align 16
  %66 = alloca [2048 x i8], align 16
  %67 = alloca [2048 x i8], align 16
  %68 = alloca [2048 x i8], align 16
  %69 = alloca [2048 x i8], align 16
  %70 = alloca i32, align 4
  %71 = alloca [2048 x i8], align 16
  %72 = alloca i32, align 4
  %73 = alloca [2048 x i8], align 16
  %74 = alloca i32, align 4
  %75 = alloca [2048 x i8], align 16
  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
  call void @llvm.lifetime.start.p0(i64 4, ptr %10) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %11) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %12) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %13) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %14) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %15) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %16) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %17) #6
  call void @llvm.lifetime.start.p0(i64 4, ptr %18) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %19) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %22) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %25) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %26) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %27) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %28) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %29) #6
  call void @llvm.lifetime.start.p0(i64 8, ptr %30) #6
  %76 = load ptr, ptr %7, align 8, !tbaa !5
  %77 = getelementptr inbounds %struct.background, ptr %76, i32 0, i32 97
  %78 = load i32, ptr %77, align 8, !tbaa !9
  %79 = sext i32 %78 to i64
  %80 = mul i64 %79, 8
  %81 = call noalias ptr @malloc(i64 noundef %80) #7
  store ptr %81, ptr %29, align 8, !tbaa !5
  %82 = load ptr, ptr %29, align 8, !tbaa !5
  %83 = icmp eq ptr %82, null
  br i1 %83, label %84, label %97

84:                                               ; preds = %4
  call void @llvm.lifetime.start.p0(i64 4, ptr %31) #6
  %85 = load ptr, ptr %7, align 8, !tbaa !5
  %86 = getelementptr inbounds %struct.background, ptr %85, i32 0, i32 97
  %87 = load i32, ptr %86, align 8, !tbaa !9
  %88 = sext i32 %87 to i64
  %89 = mul i64 %88, 8
  %90 = trunc i64 %89 to i32
  store i32 %90, ptr %31, align 4, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 2048, ptr %32) #6
  %91 = getelementptr inbounds [2048 x i8], ptr %32, i64 0, i64 0
  %92 = load i32, ptr %31, align 4, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %91, ptr noundef @.str.4, ptr noundef @.str.128, i32 noundef %92)
  %93 = load ptr, ptr %9, align 8, !tbaa !5
  %94 = getelementptr inbounds %struct.perturbs, ptr %93, i32 0, i32 186
  %95 = getelementptr inbounds [2048 x i8], ptr %94, i64 0, i64 0
  %96 = getelementptr inbounds [2048 x i8], ptr %32, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %95, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1431, ptr noundef %96)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %32) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %31) #6
  br label %1345

97:                                               ; preds = %4
  %98 = load ptr, ptr %8, align 8, !tbaa !5
  %99 = getelementptr inbounds %struct.thermo, ptr %98, i32 0, i32 61
  %100 = load i32, ptr %99, align 4, !tbaa !15
  %101 = sext i32 %100 to i64
  %102 = mul i64 %101, 8
  %103 = call noalias ptr @malloc(i64 noundef %102) #7
  store ptr %103, ptr %30, align 8, !tbaa !5
  %104 = load ptr, ptr %30, align 8, !tbaa !5
  %105 = icmp eq ptr %104, null
  br i1 %105, label %106, label %119

106:                                              ; preds = %97
  call void @llvm.lifetime.start.p0(i64 4, ptr %34) #6
  %107 = load ptr, ptr %8, align 8, !tbaa !5
  %108 = getelementptr inbounds %struct.thermo, ptr %107, i32 0, i32 61
  %109 = load i32, ptr %108, align 4, !tbaa !15
  %110 = sext i32 %109 to i64
  %111 = mul i64 %110, 8
  %112 = trunc i64 %111 to i32
  store i32 %112, ptr %34, align 4, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 2048, ptr %35) #6
  %113 = getelementptr inbounds [2048 x i8], ptr %35, i64 0, i64 0
  %114 = load i32, ptr %34, align 4, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %113, ptr noundef @.str.4, ptr noundef @.str.129, i32 noundef %114)
  %115 = load ptr, ptr %9, align 8, !tbaa !5
  %116 = getelementptr inbounds %struct.perturbs, ptr %115, i32 0, i32 186
  %117 = getelementptr inbounds [2048 x i8], ptr %116, i64 0, i64 0
  %118 = getelementptr inbounds [2048 x i8], ptr %35, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %117, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1432, ptr noundef %118)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %35) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %34) #6
  br label %1345

119:                                              ; preds = %97
  %120 = load ptr, ptr %9, align 8, !tbaa !5
  %121 = getelementptr inbounds %struct.perturbs, ptr %120, i32 0, i32 53
  %122 = load i16, ptr %121, align 4, !tbaa !17
  %123 = sext i16 %122 to i32
  %124 = icmp eq i32 %123, 1
  br i1 %124, label %131, label %125

125:                                              ; preds = %119
  %126 = load ptr, ptr %9, align 8, !tbaa !5
  %127 = getelementptr inbounds %struct.perturbs, ptr %126, i32 0, i32 10
  %128 = load i16, ptr %127, align 4, !tbaa !19
  %129 = sext i16 %128 to i32
  %130 = icmp eq i32 %129, 1
  br i1 %130, label %131, label %409

131:                                              ; preds = %125, %119
  %132 = load ptr, ptr %8, align 8, !tbaa !5
  %133 = getelementptr inbounds %struct.thermo, ptr %132, i32 0, i32 88
  %134 = load double, ptr %133, align 8, !tbaa !20
  store double %134, ptr %21, align 8, !tbaa !21
  %135 = load ptr, ptr %7, align 8, !tbaa !5
  %136 = load double, ptr %21, align 8, !tbaa !21
  %137 = load ptr, ptr %7, align 8, !tbaa !5
  %138 = getelementptr inbounds %struct.background, ptr %137, i32 0, i32 141
  %139 = load i16, ptr %138, align 8, !tbaa !22
  %140 = load ptr, ptr %7, align 8, !tbaa !5
  %141 = getelementptr inbounds %struct.background, ptr %140, i32 0, i32 144
  %142 = load i16, ptr %141, align 2, !tbaa !23
  %143 = load ptr, ptr %29, align 8, !tbaa !5
  %144 = call i32 @background_at_tau(ptr noundef %135, double noundef %136, i16 noundef signext %139, i16 noundef signext %142, ptr noundef %17, ptr noundef %143)
  %145 = icmp eq i32 %144, 1
  br i1 %145, label %146, label %155

146:                                              ; preds = %131
  call void @llvm.lifetime.start.p0(i64 2048, ptr %36) #6
  %147 = getelementptr inbounds [2048 x i8], ptr %36, i64 0, i64 0
  %148 = load ptr, ptr %7, align 8, !tbaa !5
  %149 = getelementptr inbounds %struct.background, ptr %148, i32 0, i32 149
  %150 = getelementptr inbounds [2048 x i8], ptr %149, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %147, ptr noundef @.str, ptr noundef @.str.130, ptr noundef %150)
  %151 = load ptr, ptr %9, align 8, !tbaa !5
  %152 = getelementptr inbounds %struct.perturbs, ptr %151, i32 0, i32 186
  %153 = getelementptr inbounds [2048 x i8], ptr %152, i64 0, i64 0
  %154 = getelementptr inbounds [2048 x i8], ptr %36, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %153, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1459, ptr noundef %154)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %36) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

155:                                              ; preds = %131
  %156 = load ptr, ptr %7, align 8, !tbaa !5
  %157 = load ptr, ptr %8, align 8, !tbaa !5
  %158 = load ptr, ptr %29, align 8, !tbaa !5
  %159 = load ptr, ptr %7, align 8, !tbaa !5
  %160 = getelementptr inbounds %struct.background, ptr %159, i32 0, i32 59
  %161 = load i32, ptr %160, align 8, !tbaa !24
  %162 = sext i32 %161 to i64
  %163 = getelementptr inbounds double, ptr %158, i64 %162
  %164 = load double, ptr %163, align 8, !tbaa !21
  %165 = fdiv double 1.000000e+00, %164
  %166 = fsub double %165, 1.000000e+00
  %167 = load ptr, ptr %8, align 8, !tbaa !5
  %168 = getelementptr inbounds %struct.thermo, ptr %167, i32 0, i32 90
  %169 = load i16, ptr %168, align 8, !tbaa !25
  %170 = load ptr, ptr %29, align 8, !tbaa !5
  %171 = load ptr, ptr %30, align 8, !tbaa !5
  %172 = call i32 @thermodynamics_at_z(ptr noundef %156, ptr noundef %157, double noundef %166, i16 noundef signext %169, ptr noundef %18, ptr noundef %170, ptr noundef %171)
  %173 = icmp eq i32 %172, 1
  br i1 %173, label %174, label %183

174:                                              ; preds = %155
  call void @llvm.lifetime.start.p0(i64 2048, ptr %37) #6
  %175 = getelementptr inbounds [2048 x i8], ptr %37, i64 0, i64 0
  %176 = load ptr, ptr %8, align 8, !tbaa !5
  %177 = getelementptr inbounds %struct.thermo, ptr %176, i32 0, i32 93
  %178 = getelementptr inbounds [2048 x i8], ptr %177, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %175, ptr noundef @.str, ptr noundef @.str.131, ptr noundef %178)
  %179 = load ptr, ptr %9, align 8, !tbaa !5
  %180 = getelementptr inbounds %struct.perturbs, ptr %179, i32 0, i32 186
  %181 = getelementptr inbounds [2048 x i8], ptr %180, i64 0, i64 0
  %182 = getelementptr inbounds [2048 x i8], ptr %37, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %181, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1469, ptr noundef %182)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %37) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

183:                                              ; preds = %155
  %184 = load ptr, ptr %29, align 8, !tbaa !5
  %185 = load ptr, ptr %7, align 8, !tbaa !5
  %186 = getelementptr inbounds %struct.background, ptr %185, i32 0, i32 59
  %187 = load i32, ptr %186, align 8, !tbaa !24
  %188 = sext i32 %187 to i64
  %189 = getelementptr inbounds double, ptr %184, i64 %188
  %190 = load double, ptr %189, align 8, !tbaa !21
  %191 = load ptr, ptr %29, align 8, !tbaa !5
  %192 = load ptr, ptr %7, align 8, !tbaa !5
  %193 = getelementptr inbounds %struct.background, ptr %192, i32 0, i32 60
  %194 = load i32, ptr %193, align 4, !tbaa !26
  %195 = sext i32 %194 to i64
  %196 = getelementptr inbounds double, ptr %191, i64 %195
  %197 = load double, ptr %196, align 8, !tbaa !21
  %198 = fmul double %190, %197
  %199 = load ptr, ptr %30, align 8, !tbaa !5
  %200 = load ptr, ptr %8, align 8, !tbaa !5
  %201 = getelementptr inbounds %struct.thermo, ptr %200, i32 0, i32 37
  %202 = load i32, ptr %201, align 4, !tbaa !27
  %203 = sext i32 %202 to i64
  %204 = getelementptr inbounds double, ptr %199, i64 %203
  %205 = load double, ptr %204, align 8, !tbaa !21
  %206 = fdiv double %198, %205
  %207 = load ptr, ptr %6, align 8, !tbaa !5
  %208 = getelementptr inbounds %struct.precision, ptr %207, i32 0, i32 67
  %209 = load double, ptr %208, align 8, !tbaa !28
  %210 = fcmp ogt double %206, %209
  br i1 %210, label %211, label %220

211:                                              ; preds = %183
  call void @llvm.lifetime.start.p0(i64 2048, ptr %38) #6
  %212 = getelementptr inbounds [2048 x i8], ptr %38, i64 0, i64 0
  %213 = load double, ptr %21, align 8, !tbaa !21
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %212, ptr noundef @.str.132, double noundef %213)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %39) #6
  %214 = getelementptr inbounds [2048 x i8], ptr %39, i64 0, i64 0
  %215 = getelementptr inbounds [2048 x i8], ptr %38, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %214, ptr noundef @.str.7, ptr noundef @.str.133, ptr noundef %215)
  %216 = load ptr, ptr %9, align 8, !tbaa !5
  %217 = getelementptr inbounds %struct.perturbs, ptr %216, i32 0, i32 186
  %218 = getelementptr inbounds [2048 x i8], ptr %217, i64 0, i64 0
  %219 = getelementptr inbounds [2048 x i8], ptr %39, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %218, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1477, ptr noundef %219)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %39) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %38) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

220:                                              ; preds = %183
  %221 = load ptr, ptr %8, align 8, !tbaa !5
  %222 = getelementptr inbounds %struct.thermo, ptr %221, i32 0, i32 67
  %223 = load double, ptr %222, align 8, !tbaa !30
  store double %223, ptr %22, align 8, !tbaa !21
  %224 = load ptr, ptr %7, align 8, !tbaa !5
  %225 = load double, ptr %22, align 8, !tbaa !21
  %226 = load ptr, ptr %7, align 8, !tbaa !5
  %227 = getelementptr inbounds %struct.background, ptr %226, i32 0, i32 141
  %228 = load i16, ptr %227, align 8, !tbaa !22
  %229 = load ptr, ptr %7, align 8, !tbaa !5
  %230 = getelementptr inbounds %struct.background, ptr %229, i32 0, i32 144
  %231 = load i16, ptr %230, align 2, !tbaa !23
  %232 = load ptr, ptr %29, align 8, !tbaa !5
  %233 = call i32 @background_at_tau(ptr noundef %224, double noundef %225, i16 noundef signext %228, i16 noundef signext %231, ptr noundef %17, ptr noundef %232)
  %234 = icmp eq i32 %233, 1
  br i1 %234, label %235, label %244

235:                                              ; preds = %220
  call void @llvm.lifetime.start.p0(i64 2048, ptr %40) #6
  %236 = getelementptr inbounds [2048 x i8], ptr %40, i64 0, i64 0
  %237 = load ptr, ptr %7, align 8, !tbaa !5
  %238 = getelementptr inbounds %struct.background, ptr %237, i32 0, i32 149
  %239 = getelementptr inbounds [2048 x i8], ptr %238, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %236, ptr noundef @.str, ptr noundef @.str.134, ptr noundef %239)
  %240 = load ptr, ptr %9, align 8, !tbaa !5
  %241 = getelementptr inbounds %struct.perturbs, ptr %240, i32 0, i32 186
  %242 = getelementptr inbounds [2048 x i8], ptr %241, i64 0, i64 0
  %243 = getelementptr inbounds [2048 x i8], ptr %40, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %242, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1489, ptr noundef %243)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %40) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

244:                                              ; preds = %220
  %245 = load ptr, ptr %7, align 8, !tbaa !5
  %246 = load ptr, ptr %8, align 8, !tbaa !5
  %247 = load ptr, ptr %29, align 8, !tbaa !5
  %248 = load ptr, ptr %7, align 8, !tbaa !5
  %249 = getelementptr inbounds %struct.background, ptr %248, i32 0, i32 59
  %250 = load i32, ptr %249, align 8, !tbaa !24
  %251 = sext i32 %250 to i64
  %252 = getelementptr inbounds double, ptr %247, i64 %251
  %253 = load double, ptr %252, align 8, !tbaa !21
  %254 = fdiv double 1.000000e+00, %253
  %255 = fsub double %254, 1.000000e+00
  %256 = load ptr, ptr %8, align 8, !tbaa !5
  %257 = getelementptr inbounds %struct.thermo, ptr %256, i32 0, i32 90
  %258 = load i16, ptr %257, align 8, !tbaa !25
  %259 = load ptr, ptr %29, align 8, !tbaa !5
  %260 = load ptr, ptr %30, align 8, !tbaa !5
  %261 = call i32 @thermodynamics_at_z(ptr noundef %245, ptr noundef %246, double noundef %255, i16 noundef signext %258, ptr noundef %18, ptr noundef %259, ptr noundef %260)
  %262 = icmp eq i32 %261, 1
  br i1 %262, label %263, label %272

263:                                              ; preds = %244
  call void @llvm.lifetime.start.p0(i64 2048, ptr %41) #6
  %264 = getelementptr inbounds [2048 x i8], ptr %41, i64 0, i64 0
  %265 = load ptr, ptr %8, align 8, !tbaa !5
  %266 = getelementptr inbounds %struct.thermo, ptr %265, i32 0, i32 93
  %267 = getelementptr inbounds [2048 x i8], ptr %266, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %264, ptr noundef @.str, ptr noundef @.str.131, ptr noundef %267)
  %268 = load ptr, ptr %9, align 8, !tbaa !5
  %269 = getelementptr inbounds %struct.perturbs, ptr %268, i32 0, i32 186
  %270 = getelementptr inbounds [2048 x i8], ptr %269, i64 0, i64 0
  %271 = getelementptr inbounds [2048 x i8], ptr %41, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %270, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1499, ptr noundef %271)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %41) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

272:                                              ; preds = %244
  %273 = load ptr, ptr %29, align 8, !tbaa !5
  %274 = load ptr, ptr %7, align 8, !tbaa !5
  %275 = getelementptr inbounds %struct.background, ptr %274, i32 0, i32 59
  %276 = load i32, ptr %275, align 8, !tbaa !24
  %277 = sext i32 %276 to i64
  %278 = getelementptr inbounds double, ptr %273, i64 %277
  %279 = load double, ptr %278, align 8, !tbaa !21
  %280 = load ptr, ptr %29, align 8, !tbaa !5
  %281 = load ptr, ptr %7, align 8, !tbaa !5
  %282 = getelementptr inbounds %struct.background, ptr %281, i32 0, i32 60
  %283 = load i32, ptr %282, align 4, !tbaa !26
  %284 = sext i32 %283 to i64
  %285 = getelementptr inbounds double, ptr %280, i64 %284
  %286 = load double, ptr %285, align 8, !tbaa !21
  %287 = fmul double %279, %286
  %288 = load ptr, ptr %30, align 8, !tbaa !5
  %289 = load ptr, ptr %8, align 8, !tbaa !5
  %290 = getelementptr inbounds %struct.thermo, ptr %289, i32 0, i32 37
  %291 = load i32, ptr %290, align 4, !tbaa !27
  %292 = sext i32 %291 to i64
  %293 = getelementptr inbounds double, ptr %288, i64 %292
  %294 = load double, ptr %293, align 8, !tbaa !21
  %295 = fdiv double %287, %294
  %296 = load ptr, ptr %6, align 8, !tbaa !5
  %297 = getelementptr inbounds %struct.precision, ptr %296, i32 0, i32 67
  %298 = load double, ptr %297, align 8, !tbaa !28
  %299 = fcmp olt double %295, %298
  br i1 %299, label %300, label %308

300:                                              ; preds = %272
  call void @llvm.lifetime.start.p0(i64 2048, ptr %42) #6
  %301 = getelementptr inbounds [2048 x i8], ptr %42, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %301, ptr noundef @.str.135)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %43) #6
  %302 = getelementptr inbounds [2048 x i8], ptr %43, i64 0, i64 0
  %303 = getelementptr inbounds [2048 x i8], ptr %42, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %302, ptr noundef @.str.7, ptr noundef @.str.136, ptr noundef %303)
  %304 = load ptr, ptr %9, align 8, !tbaa !5
  %305 = getelementptr inbounds %struct.perturbs, ptr %304, i32 0, i32 186
  %306 = getelementptr inbounds [2048 x i8], ptr %305, i64 0, i64 0
  %307 = getelementptr inbounds [2048 x i8], ptr %43, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %306, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1506, ptr noundef %307)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %43) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %42) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

308:                                              ; preds = %272
  %309 = load double, ptr %21, align 8, !tbaa !21
  %310 = load double, ptr %22, align 8, !tbaa !21
  %311 = fadd double %309, %310
  %312 = fmul double 5.000000e-01, %311
  store double %312, ptr %23, align 8, !tbaa !21
  br label %313

313:                                              ; preds = %402, %308
  %314 = load double, ptr %22, align 8, !tbaa !21
  %315 = load double, ptr %21, align 8, !tbaa !21
  %316 = fsub double %314, %315
  %317 = load ptr, ptr %6, align 8, !tbaa !5
  %318 = getelementptr inbounds %struct.precision, ptr %317, i32 0, i32 86
  %319 = load double, ptr %318, align 8, !tbaa !31
  %320 = fcmp ogt double %316, %319
  br i1 %320, label %321, label %407

321:                                              ; preds = %313
  %322 = load ptr, ptr %7, align 8, !tbaa !5
  %323 = load double, ptr %23, align 8, !tbaa !21
  %324 = load ptr, ptr %7, align 8, !tbaa !5
  %325 = getelementptr inbounds %struct.background, ptr %324, i32 0, i32 141
  %326 = load i16, ptr %325, align 8, !tbaa !22
  %327 = load ptr, ptr %7, align 8, !tbaa !5
  %328 = getelementptr inbounds %struct.background, ptr %327, i32 0, i32 144
  %329 = load i16, ptr %328, align 2, !tbaa !23
  %330 = load ptr, ptr %29, align 8, !tbaa !5
  %331 = call i32 @background_at_tau(ptr noundef %322, double noundef %323, i16 noundef signext %326, i16 noundef signext %329, ptr noundef %17, ptr noundef %330)
  %332 = icmp eq i32 %331, 1
  br i1 %332, label %333, label %342

333:                                              ; preds = %321
  call void @llvm.lifetime.start.p0(i64 2048, ptr %44) #6
  %334 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  %335 = load ptr, ptr %7, align 8, !tbaa !5
  %336 = getelementptr inbounds %struct.background, ptr %335, i32 0, i32 149
  %337 = getelementptr inbounds [2048 x i8], ptr %336, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %334, ptr noundef @.str, ptr noundef @.str.137, ptr noundef %337)
  %338 = load ptr, ptr %9, align 8, !tbaa !5
  %339 = getelementptr inbounds %struct.perturbs, ptr %338, i32 0, i32 186
  %340 = getelementptr inbounds [2048 x i8], ptr %339, i64 0, i64 0
  %341 = getelementptr inbounds [2048 x i8], ptr %44, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %340, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1519, ptr noundef %341)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %44) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

342:                                              ; preds = %321
  %343 = load ptr, ptr %7, align 8, !tbaa !5
  %344 = load ptr, ptr %8, align 8, !tbaa !5
  %345 = load ptr, ptr %29, align 8, !tbaa !5
  %346 = load ptr, ptr %7, align 8, !tbaa !5
  %347 = getelementptr inbounds %struct.background, ptr %346, i32 0, i32 59
  %348 = load i32, ptr %347, align 8, !tbaa !24
  %349 = sext i32 %348 to i64
  %350 = getelementptr inbounds double, ptr %345, i64 %349
  %351 = load double, ptr %350, align 8, !tbaa !21
  %352 = fdiv double 1.000000e+00, %351
  %353 = fsub double %352, 1.000000e+00
  %354 = load ptr, ptr %8, align 8, !tbaa !5
  %355 = getelementptr inbounds %struct.thermo, ptr %354, i32 0, i32 90
  %356 = load i16, ptr %355, align 8, !tbaa !25
  %357 = load ptr, ptr %29, align 8, !tbaa !5
  %358 = load ptr, ptr %30, align 8, !tbaa !5
  %359 = call i32 @thermodynamics_at_z(ptr noundef %343, ptr noundef %344, double noundef %353, i16 noundef signext %356, ptr noundef %18, ptr noundef %357, ptr noundef %358)
  %360 = icmp eq i32 %359, 1
  br i1 %360, label %361, label %370

361:                                              ; preds = %342
  call void @llvm.lifetime.start.p0(i64 2048, ptr %45) #6
  %362 = getelementptr inbounds [2048 x i8], ptr %45, i64 0, i64 0
  %363 = load ptr, ptr %8, align 8, !tbaa !5
  %364 = getelementptr inbounds %struct.thermo, ptr %363, i32 0, i32 93
  %365 = getelementptr inbounds [2048 x i8], ptr %364, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %362, ptr noundef @.str, ptr noundef @.str.131, ptr noundef %365)
  %366 = load ptr, ptr %9, align 8, !tbaa !5
  %367 = getelementptr inbounds %struct.perturbs, ptr %366, i32 0, i32 186
  %368 = getelementptr inbounds [2048 x i8], ptr %367, i64 0, i64 0
  %369 = getelementptr inbounds [2048 x i8], ptr %45, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %368, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1529, ptr noundef %369)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %45) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

370:                                              ; preds = %342
  %371 = load ptr, ptr %29, align 8, !tbaa !5
  %372 = load ptr, ptr %7, align 8, !tbaa !5
  %373 = getelementptr inbounds %struct.background, ptr %372, i32 0, i32 59
  %374 = load i32, ptr %373, align 8, !tbaa !24
  %375 = sext i32 %374 to i64
  %376 = getelementptr inbounds double, ptr %371, i64 %375
  %377 = load double, ptr %376, align 8, !tbaa !21
  %378 = load ptr, ptr %29, align 8, !tbaa !5
  %379 = load ptr, ptr %7, align 8, !tbaa !5
  %380 = getelementptr inbounds %struct.background, ptr %379, i32 0, i32 60
  %381 = load i32, ptr %380, align 4, !tbaa !26
  %382 = sext i32 %381 to i64
  %383 = getelementptr inbounds double, ptr %378, i64 %382
  %384 = load double, ptr %383, align 8, !tbaa !21
  %385 = fmul double %377, %384
  %386 = load ptr, ptr %30, align 8, !tbaa !5
  %387 = load ptr, ptr %8, align 8, !tbaa !5
  %388 = getelementptr inbounds %struct.thermo, ptr %387, i32 0, i32 37
  %389 = load i32, ptr %388, align 4, !tbaa !27
  %390 = sext i32 %389 to i64
  %391 = getelementptr inbounds double, ptr %386, i64 %390
  %392 = load double, ptr %391, align 8, !tbaa !21
  %393 = fdiv double %385, %392
  %394 = load ptr, ptr %6, align 8, !tbaa !5
  %395 = getelementptr inbounds %struct.precision, ptr %394, i32 0, i32 67
  %396 = load double, ptr %395, align 8, !tbaa !28
  %397 = fcmp ogt double %393, %396
  br i1 %397, label %398, label %400

398:                                              ; preds = %370
  %399 = load double, ptr %23, align 8, !tbaa !21
  store double %399, ptr %22, align 8, !tbaa !21
  br label %402

400:                                              ; preds = %370
  %401 = load double, ptr %23, align 8, !tbaa !21
  store double %401, ptr %21, align 8, !tbaa !21
  br label %402

402:                                              ; preds = %400, %398
  %403 = load double, ptr %21, align 8, !tbaa !21
  %404 = load double, ptr %22, align 8, !tbaa !21
  %405 = fadd double %403, %404
  %406 = fmul double 5.000000e-01, %405
  store double %406, ptr %23, align 8, !tbaa !21
  br label %313, !llvm.loop !32

407:                                              ; preds = %313
  %408 = load double, ptr %23, align 8, !tbaa !21
  store double %408, ptr %20, align 8, !tbaa !21
  br label %476

409:                                              ; preds = %125
  %410 = load ptr, ptr %7, align 8, !tbaa !5
  %411 = load ptr, ptr %9, align 8, !tbaa !5
  %412 = getelementptr inbounds %struct.perturbs, ptr %411, i32 0, i32 49
  %413 = load double, ptr %412, align 8, !tbaa !34
  %414 = fadd double %413, 1.000000e+00
  %415 = call i32 @background_tau_of_z(ptr noundef %410, double noundef %414, ptr noundef %20)
  %416 = icmp eq i32 %415, 1
  br i1 %416, label %417, label %426

417:                                              ; preds = %409
  call void @llvm.lifetime.start.p0(i64 2048, ptr %46) #6
  %418 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  %419 = load ptr, ptr %7, align 8, !tbaa !5
  %420 = getelementptr inbounds %struct.background, ptr %419, i32 0, i32 149
  %421 = getelementptr inbounds [2048 x i8], ptr %420, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %418, ptr noundef @.str, ptr noundef @.str.138, ptr noundef %421)
  %422 = load ptr, ptr %9, align 8, !tbaa !5
  %423 = getelementptr inbounds %struct.perturbs, ptr %422, i32 0, i32 186
  %424 = getelementptr inbounds [2048 x i8], ptr %423, i64 0, i64 0
  %425 = getelementptr inbounds [2048 x i8], ptr %46, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %424, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1555, ptr noundef %425)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %46) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

426:                                              ; preds = %409
  %427 = load ptr, ptr %7, align 8, !tbaa !5
  %428 = load double, ptr %20, align 8, !tbaa !21
  %429 = load ptr, ptr %7, align 8, !tbaa !5
  %430 = getelementptr inbounds %struct.background, ptr %429, i32 0, i32 141
  %431 = load i16, ptr %430, align 8, !tbaa !22
  %432 = load ptr, ptr %7, align 8, !tbaa !5
  %433 = getelementptr inbounds %struct.background, ptr %432, i32 0, i32 144
  %434 = load i16, ptr %433, align 2, !tbaa !23
  %435 = load ptr, ptr %29, align 8, !tbaa !5
  %436 = call i32 @background_at_tau(ptr noundef %427, double noundef %428, i16 noundef signext %431, i16 noundef signext %434, ptr noundef %17, ptr noundef %435)
  %437 = icmp eq i32 %436, 1
  br i1 %437, label %438, label %447

438:                                              ; preds = %426
  call void @llvm.lifetime.start.p0(i64 2048, ptr %47) #6
  %439 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  %440 = load ptr, ptr %7, align 8, !tbaa !5
  %441 = getelementptr inbounds %struct.background, ptr %440, i32 0, i32 149
  %442 = getelementptr inbounds [2048 x i8], ptr %441, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %439, ptr noundef @.str, ptr noundef @.str.139, ptr noundef %442)
  %443 = load ptr, ptr %9, align 8, !tbaa !5
  %444 = getelementptr inbounds %struct.perturbs, ptr %443, i32 0, i32 186
  %445 = getelementptr inbounds [2048 x i8], ptr %444, i64 0, i64 0
  %446 = getelementptr inbounds [2048 x i8], ptr %47, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %445, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1568, ptr noundef %446)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %47) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

447:                                              ; preds = %426
  %448 = load ptr, ptr %7, align 8, !tbaa !5
  %449 = load ptr, ptr %8, align 8, !tbaa !5
  %450 = load ptr, ptr %29, align 8, !tbaa !5
  %451 = load ptr, ptr %7, align 8, !tbaa !5
  %452 = getelementptr inbounds %struct.background, ptr %451, i32 0, i32 59
  %453 = load i32, ptr %452, align 8, !tbaa !24
  %454 = sext i32 %453 to i64
  %455 = getelementptr inbounds double, ptr %450, i64 %454
  %456 = load double, ptr %455, align 8, !tbaa !21
  %457 = fdiv double 1.000000e+00, %456
  %458 = fsub double %457, 1.000000e+00
  %459 = load ptr, ptr %8, align 8, !tbaa !5
  %460 = getelementptr inbounds %struct.thermo, ptr %459, i32 0, i32 90
  %461 = load i16, ptr %460, align 8, !tbaa !25
  %462 = load ptr, ptr %29, align 8, !tbaa !5
  %463 = load ptr, ptr %30, align 8, !tbaa !5
  %464 = call i32 @thermodynamics_at_z(ptr noundef %448, ptr noundef %449, double noundef %458, i16 noundef signext %461, ptr noundef %18, ptr noundef %462, ptr noundef %463)
  %465 = icmp eq i32 %464, 1
  br i1 %465, label %466, label %475

466:                                              ; preds = %447
  call void @llvm.lifetime.start.p0(i64 2048, ptr %48) #6
  %467 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  %468 = load ptr, ptr %8, align 8, !tbaa !5
  %469 = getelementptr inbounds %struct.thermo, ptr %468, i32 0, i32 93
  %470 = getelementptr inbounds [2048 x i8], ptr %469, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %467, ptr noundef @.str, ptr noundef @.str.131, ptr noundef %470)
  %471 = load ptr, ptr %9, align 8, !tbaa !5
  %472 = getelementptr inbounds %struct.perturbs, ptr %471, i32 0, i32 186
  %473 = getelementptr inbounds [2048 x i8], ptr %472, i64 0, i64 0
  %474 = getelementptr inbounds [2048 x i8], ptr %48, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %473, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1578, ptr noundef %474)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %48) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

475:                                              ; preds = %447
  br label %476

476:                                              ; preds = %475, %407
  store i32 1, ptr %10, align 4, !tbaa !14
  %477 = load i32, ptr %17, align 4, !tbaa !14
  store i32 %477, ptr %15, align 4, !tbaa !14
  %478 = load i32, ptr %18, align 4, !tbaa !14
  store i32 %478, ptr %16, align 4, !tbaa !14
  %479 = load double, ptr %20, align 8, !tbaa !21
  store double %479, ptr %19, align 8, !tbaa !21
  br label %480

480:                                              ; preds = %651, %476
  %481 = load double, ptr %19, align 8, !tbaa !21
  %482 = load ptr, ptr %7, align 8, !tbaa !5
  %483 = getelementptr inbounds %struct.background, ptr %482, i32 0, i32 44
  %484 = load double, ptr %483, align 8, !tbaa !35
  %485 = fcmp olt double %481, %484
  br i1 %485, label %486, label %660

486:                                              ; preds = %480
  %487 = load ptr, ptr %7, align 8, !tbaa !5
  %488 = load double, ptr %19, align 8, !tbaa !21
  %489 = load ptr, ptr %7, align 8, !tbaa !5
  %490 = getelementptr inbounds %struct.background, ptr %489, i32 0, i32 141
  %491 = load i16, ptr %490, align 8, !tbaa !22
  %492 = load ptr, ptr %7, align 8, !tbaa !5
  %493 = getelementptr inbounds %struct.background, ptr %492, i32 0, i32 145
  %494 = load i16, ptr %493, align 8, !tbaa !36
  %495 = load ptr, ptr %29, align 8, !tbaa !5
  %496 = call i32 @background_at_tau(ptr noundef %487, double noundef %488, i16 noundef signext %491, i16 noundef signext %494, ptr noundef %15, ptr noundef %495)
  %497 = icmp eq i32 %496, 1
  br i1 %497, label %498, label %507

498:                                              ; preds = %486
  call void @llvm.lifetime.start.p0(i64 2048, ptr %49) #6
  %499 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  %500 = load ptr, ptr %7, align 8, !tbaa !5
  %501 = getelementptr inbounds %struct.background, ptr %500, i32 0, i32 149
  %502 = getelementptr inbounds [2048 x i8], ptr %501, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %499, ptr noundef @.str, ptr noundef @.str.140, ptr noundef %502)
  %503 = load ptr, ptr %9, align 8, !tbaa !5
  %504 = getelementptr inbounds %struct.perturbs, ptr %503, i32 0, i32 186
  %505 = getelementptr inbounds [2048 x i8], ptr %504, i64 0, i64 0
  %506 = getelementptr inbounds [2048 x i8], ptr %49, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %505, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1604, ptr noundef %506)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %49) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

507:                                              ; preds = %486
  %508 = load ptr, ptr %7, align 8, !tbaa !5
  %509 = load ptr, ptr %8, align 8, !tbaa !5
  %510 = load ptr, ptr %29, align 8, !tbaa !5
  %511 = load ptr, ptr %7, align 8, !tbaa !5
  %512 = getelementptr inbounds %struct.background, ptr %511, i32 0, i32 59
  %513 = load i32, ptr %512, align 8, !tbaa !24
  %514 = sext i32 %513 to i64
  %515 = getelementptr inbounds double, ptr %510, i64 %514
  %516 = load double, ptr %515, align 8, !tbaa !21
  %517 = fdiv double 1.000000e+00, %516
  %518 = fsub double %517, 1.000000e+00
  %519 = load ptr, ptr %8, align 8, !tbaa !5
  %520 = getelementptr inbounds %struct.thermo, ptr %519, i32 0, i32 91
  %521 = load i16, ptr %520, align 2, !tbaa !37
  %522 = load ptr, ptr %29, align 8, !tbaa !5
  %523 = load ptr, ptr %30, align 8, !tbaa !5
  %524 = call i32 @thermodynamics_at_z(ptr noundef %508, ptr noundef %509, double noundef %518, i16 noundef signext %521, ptr noundef %16, ptr noundef %522, ptr noundef %523)
  %525 = icmp eq i32 %524, 1
  br i1 %525, label %526, label %535

526:                                              ; preds = %507
  call void @llvm.lifetime.start.p0(i64 2048, ptr %50) #6
  %527 = getelementptr inbounds [2048 x i8], ptr %50, i64 0, i64 0
  %528 = load ptr, ptr %8, align 8, !tbaa !5
  %529 = getelementptr inbounds %struct.thermo, ptr %528, i32 0, i32 93
  %530 = getelementptr inbounds [2048 x i8], ptr %529, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %527, ptr noundef @.str, ptr noundef @.str.141, ptr noundef %530)
  %531 = load ptr, ptr %9, align 8, !tbaa !5
  %532 = getelementptr inbounds %struct.perturbs, ptr %531, i32 0, i32 186
  %533 = getelementptr inbounds [2048 x i8], ptr %532, i64 0, i64 0
  %534 = getelementptr inbounds [2048 x i8], ptr %50, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %533, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1614, ptr noundef %534)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %50) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

535:                                              ; preds = %507
  %536 = load ptr, ptr %9, align 8, !tbaa !5
  %537 = getelementptr inbounds %struct.perturbs, ptr %536, i32 0, i32 53
  %538 = load i16, ptr %537, align 4, !tbaa !17
  %539 = sext i16 %538 to i32
  %540 = icmp eq i32 %539, 1
  br i1 %540, label %541, label %595

541:                                              ; preds = %535
  %542 = load ptr, ptr %30, align 8, !tbaa !5
  %543 = load ptr, ptr %8, align 8, !tbaa !5
  %544 = getelementptr inbounds %struct.thermo, ptr %543, i32 0, i32 59
  %545 = load i32, ptr %544, align 4, !tbaa !38
  %546 = sext i32 %545 to i64
  %547 = getelementptr inbounds double, ptr %542, i64 %546
  %548 = load double, ptr %547, align 8, !tbaa !21
  store double %548, ptr %25, align 8, !tbaa !21
  %549 = load ptr, ptr %29, align 8, !tbaa !5
  %550 = load ptr, ptr %7, align 8, !tbaa !5
  %551 = getelementptr inbounds %struct.background, ptr %550, i32 0, i32 60
  %552 = load i32, ptr %551, align 4, !tbaa !26
  %553 = sext i32 %552 to i64
  %554 = getelementptr inbounds double, ptr %549, i64 %553
  %555 = load double, ptr %554, align 8, !tbaa !21
  %556 = load ptr, ptr %29, align 8, !tbaa !5
  %557 = load ptr, ptr %7, align 8, !tbaa !5
  %558 = getelementptr inbounds %struct.background, ptr %557, i32 0, i32 59
  %559 = load i32, ptr %558, align 8, !tbaa !24
  %560 = sext i32 %559 to i64
  %561 = getelementptr inbounds double, ptr %556, i64 %560
  %562 = load double, ptr %561, align 8, !tbaa !21
  %563 = fmul double %555, %562
  store double %563, ptr %27, align 8, !tbaa !21
  %564 = load ptr, ptr %29, align 8, !tbaa !5
  %565 = load ptr, ptr %7, align 8, !tbaa !5
  %566 = getelementptr inbounds %struct.background, ptr %565, i32 0, i32 61
  %567 = load i32, ptr %566, align 8, !tbaa !39
  %568 = sext i32 %567 to i64
  %569 = getelementptr inbounds double, ptr %564, i64 %568
  %570 = load double, ptr %569, align 8, !tbaa !21
  %571 = load ptr, ptr %29, align 8, !tbaa !5
  %572 = load ptr, ptr %7, align 8, !tbaa !5
  %573 = getelementptr inbounds %struct.background, ptr %572, i32 0, i32 59
  %574 = load i32, ptr %573, align 8, !tbaa !24
  %575 = sext i32 %574 to i64
  %576 = getelementptr inbounds double, ptr %571, i64 %575
  %577 = load double, ptr %576, align 8, !tbaa !21
  %578 = load double, ptr %27, align 8, !tbaa !21
  %579 = fmul double 2.000000e+00, %578
  %580 = load double, ptr %27, align 8, !tbaa !21
  %581 = fmul double %579, %580
  %582 = call double @llvm.fmuladd.f64(double %570, double %577, double %581)
  store double %582, ptr %28, align 8, !tbaa !21
  %583 = load double, ptr %28, align 8, !tbaa !21
  %584 = load double, ptr %27, align 8, !tbaa !21
  %585 = load double, ptr %27, align 8, !tbaa !21
  %586 = fmul double %584, %585
  %587 = fneg double %586
  %588 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %583, double %587)
  %589 = call double @llvm.fabs.f64(double %588)
  store double %589, ptr %26, align 8, !tbaa !21
  %590 = load double, ptr %25, align 8, !tbaa !21
  %591 = load double, ptr %25, align 8, !tbaa !21
  %592 = load double, ptr %26, align 8, !tbaa !21
  %593 = call double @llvm.fmuladd.f64(double %590, double %591, double %592)
  %594 = call double @sqrt(double noundef %593) #6
  store double %594, ptr %24, align 8, !tbaa !21
  br label %612

595:                                              ; preds = %535
  %596 = load ptr, ptr %29, align 8, !tbaa !5
  %597 = load ptr, ptr %7, align 8, !tbaa !5
  %598 = getelementptr inbounds %struct.background, ptr %597, i32 0, i32 60
  %599 = load i32, ptr %598, align 4, !tbaa !26
  %600 = sext i32 %599 to i64
  %601 = getelementptr inbounds double, ptr %596, i64 %600
  %602 = load double, ptr %601, align 8, !tbaa !21
  %603 = load ptr, ptr %29, align 8, !tbaa !5
  %604 = load ptr, ptr %7, align 8, !tbaa !5
  %605 = getelementptr inbounds %struct.background, ptr %604, i32 0, i32 59
  %606 = load i32, ptr %605, align 8, !tbaa !24
  %607 = sext i32 %606 to i64
  %608 = getelementptr inbounds double, ptr %603, i64 %607
  %609 = load double, ptr %608, align 8, !tbaa !21
  %610 = fmul double %602, %609
  store double %610, ptr %27, align 8, !tbaa !21
  %611 = load double, ptr %27, align 8, !tbaa !21
  store double %611, ptr %24, align 8, !tbaa !21
  br label %612

612:                                              ; preds = %595, %541
  %613 = load double, ptr %24, align 8, !tbaa !21
  %614 = fcmp oeq double %613, 0.000000e+00
  br i1 %614, label %615, label %623

615:                                              ; preds = %612
  call void @llvm.lifetime.start.p0(i64 2048, ptr %51) #6
  %616 = getelementptr inbounds [2048 x i8], ptr %51, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %616, ptr noundef @.str.142)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %52) #6
  %617 = getelementptr inbounds [2048 x i8], ptr %52, i64 0, i64 0
  %618 = getelementptr inbounds [2048 x i8], ptr %51, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %617, ptr noundef @.str.7, ptr noundef @.str.143, ptr noundef %618)
  %619 = load ptr, ptr %9, align 8, !tbaa !5
  %620 = getelementptr inbounds %struct.perturbs, ptr %619, i32 0, i32 186
  %621 = getelementptr inbounds [2048 x i8], ptr %620, i64 0, i64 0
  %622 = getelementptr inbounds [2048 x i8], ptr %52, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %621, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1640, ptr noundef %622)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %52) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %51) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

623:                                              ; preds = %612
  %624 = load double, ptr %24, align 8, !tbaa !21
  %625 = fdiv double 1.000000e+00, %624
  store double %625, ptr %24, align 8, !tbaa !21
  %626 = load ptr, ptr %6, align 8, !tbaa !5
  %627 = getelementptr inbounds %struct.precision, ptr %626, i32 0, i32 83
  %628 = load double, ptr %627, align 8, !tbaa !40
  %629 = load double, ptr %24, align 8, !tbaa !21
  %630 = fmul double %628, %629
  %631 = load double, ptr %19, align 8, !tbaa !21
  %632 = fdiv double %630, %631
  %633 = call double @llvm.fabs.f64(double %632)
  %634 = load ptr, ptr %6, align 8, !tbaa !5
  %635 = getelementptr inbounds %struct.precision, ptr %634, i32 0, i32 175
  %636 = load double, ptr %635, align 8, !tbaa !41
  %637 = fcmp olt double %633, %636
  br i1 %637, label %638, label %651

638:                                              ; preds = %623
  call void @llvm.lifetime.start.p0(i64 2048, ptr %53) #6
  %639 = getelementptr inbounds [2048 x i8], ptr %53, i64 0, i64 0
  %640 = load ptr, ptr %6, align 8, !tbaa !5
  %641 = getelementptr inbounds %struct.precision, ptr %640, i32 0, i32 83
  %642 = load double, ptr %641, align 8, !tbaa !40
  %643 = load double, ptr %24, align 8, !tbaa !21
  %644 = fmul double %642, %643
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %639, ptr noundef @.str.144, double noundef %644)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %54) #6
  %645 = getelementptr inbounds [2048 x i8], ptr %54, i64 0, i64 0
  %646 = getelementptr inbounds [2048 x i8], ptr %53, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %645, ptr noundef @.str.7, ptr noundef @.str.145, ptr noundef %646)
  %647 = load ptr, ptr %9, align 8, !tbaa !5
  %648 = getelementptr inbounds %struct.perturbs, ptr %647, i32 0, i32 186
  %649 = getelementptr inbounds [2048 x i8], ptr %648, i64 0, i64 0
  %650 = getelementptr inbounds [2048 x i8], ptr %54, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %649, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1647, ptr noundef %650)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %54) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %53) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

651:                                              ; preds = %623
  %652 = load double, ptr %19, align 8, !tbaa !21
  %653 = load ptr, ptr %6, align 8, !tbaa !5
  %654 = getelementptr inbounds %struct.precision, ptr %653, i32 0, i32 83
  %655 = load double, ptr %654, align 8, !tbaa !40
  %656 = load double, ptr %24, align 8, !tbaa !21
  %657 = call double @llvm.fmuladd.f64(double %655, double %656, double %652)
  store double %657, ptr %19, align 8, !tbaa !21
  %658 = load i32, ptr %10, align 4, !tbaa !14
  %659 = add nsw i32 %658, 1
  store i32 %659, ptr %10, align 4, !tbaa !14
  br label %480, !llvm.loop !42

660:                                              ; preds = %480
  %661 = load i32, ptr %10, align 4, !tbaa !14
  %662 = load ptr, ptr %9, align 8, !tbaa !5
  %663 = getelementptr inbounds %struct.perturbs, ptr %662, i32 0, i32 159
  store i32 %661, ptr %663, align 8, !tbaa !43
  %664 = load ptr, ptr %9, align 8, !tbaa !5
  %665 = getelementptr inbounds %struct.perturbs, ptr %664, i32 0, i32 159
  %666 = load i32, ptr %665, align 8, !tbaa !43
  %667 = sext i32 %666 to i64
  %668 = mul i64 %667, 8
  %669 = call noalias ptr @malloc(i64 noundef %668) #7
  %670 = load ptr, ptr %9, align 8, !tbaa !5
  %671 = getelementptr inbounds %struct.perturbs, ptr %670, i32 0, i32 158
  store ptr %669, ptr %671, align 8, !tbaa !44
  %672 = load ptr, ptr %9, align 8, !tbaa !5
  %673 = getelementptr inbounds %struct.perturbs, ptr %672, i32 0, i32 158
  %674 = load ptr, ptr %673, align 8, !tbaa !44
  %675 = icmp eq ptr %674, null
  br i1 %675, label %676, label %689

676:                                              ; preds = %660
  call void @llvm.lifetime.start.p0(i64 4, ptr %55) #6
  %677 = load ptr, ptr %9, align 8, !tbaa !5
  %678 = getelementptr inbounds %struct.perturbs, ptr %677, i32 0, i32 159
  %679 = load i32, ptr %678, align 8, !tbaa !43
  %680 = sext i32 %679 to i64
  %681 = mul i64 %680, 8
  %682 = trunc i64 %681 to i32
  store i32 %682, ptr %55, align 4, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 2048, ptr %56) #6
  %683 = getelementptr inbounds [2048 x i8], ptr %56, i64 0, i64 0
  %684 = load i32, ptr %55, align 4, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %683, ptr noundef @.str.4, ptr noundef @.str.146, i32 noundef %684)
  %685 = load ptr, ptr %9, align 8, !tbaa !5
  %686 = getelementptr inbounds %struct.perturbs, ptr %685, i32 0, i32 186
  %687 = getelementptr inbounds [2048 x i8], ptr %686, i64 0, i64 0
  %688 = getelementptr inbounds [2048 x i8], ptr %56, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %687, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1658, ptr noundef %688)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %56) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %55) #6
  br label %1345

689:                                              ; preds = %660
  store i32 0, ptr %10, align 4, !tbaa !14
  %690 = load double, ptr %20, align 8, !tbaa !21
  %691 = load ptr, ptr %9, align 8, !tbaa !5
  %692 = getelementptr inbounds %struct.perturbs, ptr %691, i32 0, i32 158
  %693 = load ptr, ptr %692, align 8, !tbaa !44
  %694 = load i32, ptr %10, align 4, !tbaa !14
  %695 = sext i32 %694 to i64
  %696 = getelementptr inbounds double, ptr %693, i64 %695
  store double %690, ptr %696, align 8, !tbaa !21
  %697 = load i32, ptr %17, align 4, !tbaa !14
  store i32 %697, ptr %15, align 4, !tbaa !14
  %698 = load i32, ptr %18, align 4, !tbaa !14
  store i32 %698, ptr %16, align 4, !tbaa !14
  %699 = load double, ptr %20, align 8, !tbaa !21
  store double %699, ptr %19, align 8, !tbaa !21
  br label %700

700:                                              ; preds = %871, %689
  %701 = load double, ptr %19, align 8, !tbaa !21
  %702 = load ptr, ptr %7, align 8, !tbaa !5
  %703 = getelementptr inbounds %struct.background, ptr %702, i32 0, i32 44
  %704 = load double, ptr %703, align 8, !tbaa !35
  %705 = fcmp olt double %701, %704
  br i1 %705, label %706, label %887

706:                                              ; preds = %700
  %707 = load ptr, ptr %7, align 8, !tbaa !5
  %708 = load double, ptr %19, align 8, !tbaa !21
  %709 = load ptr, ptr %7, align 8, !tbaa !5
  %710 = getelementptr inbounds %struct.background, ptr %709, i32 0, i32 141
  %711 = load i16, ptr %710, align 8, !tbaa !22
  %712 = load ptr, ptr %7, align 8, !tbaa !5
  %713 = getelementptr inbounds %struct.background, ptr %712, i32 0, i32 145
  %714 = load i16, ptr %713, align 8, !tbaa !36
  %715 = load ptr, ptr %29, align 8, !tbaa !5
  %716 = call i32 @background_at_tau(ptr noundef %707, double noundef %708, i16 noundef signext %711, i16 noundef signext %714, ptr noundef %15, ptr noundef %715)
  %717 = icmp eq i32 %716, 1
  br i1 %717, label %718, label %727

718:                                              ; preds = %706
  call void @llvm.lifetime.start.p0(i64 2048, ptr %57) #6
  %719 = getelementptr inbounds [2048 x i8], ptr %57, i64 0, i64 0
  %720 = load ptr, ptr %7, align 8, !tbaa !5
  %721 = getelementptr inbounds %struct.background, ptr %720, i32 0, i32 149
  %722 = getelementptr inbounds [2048 x i8], ptr %721, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %719, ptr noundef @.str, ptr noundef @.str.140, ptr noundef %722)
  %723 = load ptr, ptr %9, align 8, !tbaa !5
  %724 = getelementptr inbounds %struct.perturbs, ptr %723, i32 0, i32 186
  %725 = getelementptr inbounds [2048 x i8], ptr %724, i64 0, i64 0
  %726 = getelementptr inbounds [2048 x i8], ptr %57, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %725, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1687, ptr noundef %726)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %57) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

727:                                              ; preds = %706
  %728 = load ptr, ptr %7, align 8, !tbaa !5
  %729 = load ptr, ptr %8, align 8, !tbaa !5
  %730 = load ptr, ptr %29, align 8, !tbaa !5
  %731 = load ptr, ptr %7, align 8, !tbaa !5
  %732 = getelementptr inbounds %struct.background, ptr %731, i32 0, i32 59
  %733 = load i32, ptr %732, align 8, !tbaa !24
  %734 = sext i32 %733 to i64
  %735 = getelementptr inbounds double, ptr %730, i64 %734
  %736 = load double, ptr %735, align 8, !tbaa !21
  %737 = fdiv double 1.000000e+00, %736
  %738 = fsub double %737, 1.000000e+00
  %739 = load ptr, ptr %8, align 8, !tbaa !5
  %740 = getelementptr inbounds %struct.thermo, ptr %739, i32 0, i32 91
  %741 = load i16, ptr %740, align 2, !tbaa !37
  %742 = load ptr, ptr %29, align 8, !tbaa !5
  %743 = load ptr, ptr %30, align 8, !tbaa !5
  %744 = call i32 @thermodynamics_at_z(ptr noundef %728, ptr noundef %729, double noundef %738, i16 noundef signext %741, ptr noundef %16, ptr noundef %742, ptr noundef %743)
  %745 = icmp eq i32 %744, 1
  br i1 %745, label %746, label %755

746:                                              ; preds = %727
  call void @llvm.lifetime.start.p0(i64 2048, ptr %58) #6
  %747 = getelementptr inbounds [2048 x i8], ptr %58, i64 0, i64 0
  %748 = load ptr, ptr %8, align 8, !tbaa !5
  %749 = getelementptr inbounds %struct.thermo, ptr %748, i32 0, i32 93
  %750 = getelementptr inbounds [2048 x i8], ptr %749, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %747, ptr noundef @.str, ptr noundef @.str.141, ptr noundef %750)
  %751 = load ptr, ptr %9, align 8, !tbaa !5
  %752 = getelementptr inbounds %struct.perturbs, ptr %751, i32 0, i32 186
  %753 = getelementptr inbounds [2048 x i8], ptr %752, i64 0, i64 0
  %754 = getelementptr inbounds [2048 x i8], ptr %58, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %753, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1697, ptr noundef %754)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %58) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

755:                                              ; preds = %727
  %756 = load ptr, ptr %9, align 8, !tbaa !5
  %757 = getelementptr inbounds %struct.perturbs, ptr %756, i32 0, i32 53
  %758 = load i16, ptr %757, align 4, !tbaa !17
  %759 = sext i16 %758 to i32
  %760 = icmp eq i32 %759, 1
  br i1 %760, label %761, label %815

761:                                              ; preds = %755
  %762 = load ptr, ptr %30, align 8, !tbaa !5
  %763 = load ptr, ptr %8, align 8, !tbaa !5
  %764 = getelementptr inbounds %struct.thermo, ptr %763, i32 0, i32 59
  %765 = load i32, ptr %764, align 4, !tbaa !38
  %766 = sext i32 %765 to i64
  %767 = getelementptr inbounds double, ptr %762, i64 %766
  %768 = load double, ptr %767, align 8, !tbaa !21
  store double %768, ptr %25, align 8, !tbaa !21
  %769 = load ptr, ptr %29, align 8, !tbaa !5
  %770 = load ptr, ptr %7, align 8, !tbaa !5
  %771 = getelementptr inbounds %struct.background, ptr %770, i32 0, i32 60
  %772 = load i32, ptr %771, align 4, !tbaa !26
  %773 = sext i32 %772 to i64
  %774 = getelementptr inbounds double, ptr %769, i64 %773
  %775 = load double, ptr %774, align 8, !tbaa !21
  %776 = load ptr, ptr %29, align 8, !tbaa !5
  %777 = load ptr, ptr %7, align 8, !tbaa !5
  %778 = getelementptr inbounds %struct.background, ptr %777, i32 0, i32 59
  %779 = load i32, ptr %778, align 8, !tbaa !24
  %780 = sext i32 %779 to i64
  %781 = getelementptr inbounds double, ptr %776, i64 %780
  %782 = load double, ptr %781, align 8, !tbaa !21
  %783 = fmul double %775, %782
  store double %783, ptr %27, align 8, !tbaa !21
  %784 = load ptr, ptr %29, align 8, !tbaa !5
  %785 = load ptr, ptr %7, align 8, !tbaa !5
  %786 = getelementptr inbounds %struct.background, ptr %785, i32 0, i32 61
  %787 = load i32, ptr %786, align 8, !tbaa !39
  %788 = sext i32 %787 to i64
  %789 = getelementptr inbounds double, ptr %784, i64 %788
  %790 = load double, ptr %789, align 8, !tbaa !21
  %791 = load ptr, ptr %29, align 8, !tbaa !5
  %792 = load ptr, ptr %7, align 8, !tbaa !5
  %793 = getelementptr inbounds %struct.background, ptr %792, i32 0, i32 59
  %794 = load i32, ptr %793, align 8, !tbaa !24
  %795 = sext i32 %794 to i64
  %796 = getelementptr inbounds double, ptr %791, i64 %795
  %797 = load double, ptr %796, align 8, !tbaa !21
  %798 = load double, ptr %27, align 8, !tbaa !21
  %799 = fmul double 2.000000e+00, %798
  %800 = load double, ptr %27, align 8, !tbaa !21
  %801 = fmul double %799, %800
  %802 = call double @llvm.fmuladd.f64(double %790, double %797, double %801)
  store double %802, ptr %28, align 8, !tbaa !21
  %803 = load double, ptr %28, align 8, !tbaa !21
  %804 = load double, ptr %27, align 8, !tbaa !21
  %805 = load double, ptr %27, align 8, !tbaa !21
  %806 = fmul double %804, %805
  %807 = fneg double %806
  %808 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %803, double %807)
  %809 = call double @llvm.fabs.f64(double %808)
  store double %809, ptr %26, align 8, !tbaa !21
  %810 = load double, ptr %25, align 8, !tbaa !21
  %811 = load double, ptr %25, align 8, !tbaa !21
  %812 = load double, ptr %26, align 8, !tbaa !21
  %813 = call double @llvm.fmuladd.f64(double %810, double %811, double %812)
  %814 = call double @sqrt(double noundef %813) #6
  store double %814, ptr %24, align 8, !tbaa !21
  br label %832

815:                                              ; preds = %755
  %816 = load ptr, ptr %29, align 8, !tbaa !5
  %817 = load ptr, ptr %7, align 8, !tbaa !5
  %818 = getelementptr inbounds %struct.background, ptr %817, i32 0, i32 60
  %819 = load i32, ptr %818, align 4, !tbaa !26
  %820 = sext i32 %819 to i64
  %821 = getelementptr inbounds double, ptr %816, i64 %820
  %822 = load double, ptr %821, align 8, !tbaa !21
  %823 = load ptr, ptr %29, align 8, !tbaa !5
  %824 = load ptr, ptr %7, align 8, !tbaa !5
  %825 = getelementptr inbounds %struct.background, ptr %824, i32 0, i32 59
  %826 = load i32, ptr %825, align 8, !tbaa !24
  %827 = sext i32 %826 to i64
  %828 = getelementptr inbounds double, ptr %823, i64 %827
  %829 = load double, ptr %828, align 8, !tbaa !21
  %830 = fmul double %822, %829
  store double %830, ptr %27, align 8, !tbaa !21
  %831 = load double, ptr %27, align 8, !tbaa !21
  store double %831, ptr %24, align 8, !tbaa !21
  br label %832

832:                                              ; preds = %815, %761
  %833 = load double, ptr %24, align 8, !tbaa !21
  %834 = fcmp oeq double %833, 0.000000e+00
  br i1 %834, label %835, label %843

835:                                              ; preds = %832
  call void @llvm.lifetime.start.p0(i64 2048, ptr %59) #6
  %836 = getelementptr inbounds [2048 x i8], ptr %59, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %836, ptr noundef @.str.142)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %60) #6
  %837 = getelementptr inbounds [2048 x i8], ptr %60, i64 0, i64 0
  %838 = getelementptr inbounds [2048 x i8], ptr %59, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %837, ptr noundef @.str.7, ptr noundef @.str.143, ptr noundef %838)
  %839 = load ptr, ptr %9, align 8, !tbaa !5
  %840 = getelementptr inbounds %struct.perturbs, ptr %839, i32 0, i32 186
  %841 = getelementptr inbounds [2048 x i8], ptr %840, i64 0, i64 0
  %842 = getelementptr inbounds [2048 x i8], ptr %60, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %841, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1721, ptr noundef %842)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %60) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %59) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

843:                                              ; preds = %832
  %844 = load double, ptr %24, align 8, !tbaa !21
  %845 = fdiv double 1.000000e+00, %844
  store double %845, ptr %24, align 8, !tbaa !21
  %846 = load ptr, ptr %6, align 8, !tbaa !5
  %847 = getelementptr inbounds %struct.precision, ptr %846, i32 0, i32 83
  %848 = load double, ptr %847, align 8, !tbaa !40
  %849 = load double, ptr %24, align 8, !tbaa !21
  %850 = fmul double %848, %849
  %851 = load double, ptr %19, align 8, !tbaa !21
  %852 = fdiv double %850, %851
  %853 = call double @llvm.fabs.f64(double %852)
  %854 = load ptr, ptr %6, align 8, !tbaa !5
  %855 = getelementptr inbounds %struct.precision, ptr %854, i32 0, i32 175
  %856 = load double, ptr %855, align 8, !tbaa !41
  %857 = fcmp olt double %853, %856
  br i1 %857, label %858, label %871

858:                                              ; preds = %843
  call void @llvm.lifetime.start.p0(i64 2048, ptr %61) #6
  %859 = getelementptr inbounds [2048 x i8], ptr %61, i64 0, i64 0
  %860 = load ptr, ptr %6, align 8, !tbaa !5
  %861 = getelementptr inbounds %struct.precision, ptr %860, i32 0, i32 83
  %862 = load double, ptr %861, align 8, !tbaa !40
  %863 = load double, ptr %24, align 8, !tbaa !21
  %864 = fmul double %862, %863
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %859, ptr noundef @.str.144, double noundef %864)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %62) #6
  %865 = getelementptr inbounds [2048 x i8], ptr %62, i64 0, i64 0
  %866 = getelementptr inbounds [2048 x i8], ptr %61, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %865, ptr noundef @.str.7, ptr noundef @.str.145, ptr noundef %866)
  %867 = load ptr, ptr %9, align 8, !tbaa !5
  %868 = getelementptr inbounds %struct.perturbs, ptr %867, i32 0, i32 186
  %869 = getelementptr inbounds [2048 x i8], ptr %868, i64 0, i64 0
  %870 = getelementptr inbounds [2048 x i8], ptr %62, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %869, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1728, ptr noundef %870)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %62) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %61) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

871:                                              ; preds = %843
  %872 = load double, ptr %19, align 8, !tbaa !21
  %873 = load ptr, ptr %6, align 8, !tbaa !5
  %874 = getelementptr inbounds %struct.precision, ptr %873, i32 0, i32 83
  %875 = load double, ptr %874, align 8, !tbaa !40
  %876 = load double, ptr %24, align 8, !tbaa !21
  %877 = call double @llvm.fmuladd.f64(double %875, double %876, double %872)
  store double %877, ptr %19, align 8, !tbaa !21
  %878 = load i32, ptr %10, align 4, !tbaa !14
  %879 = add nsw i32 %878, 1
  store i32 %879, ptr %10, align 4, !tbaa !14
  %880 = load double, ptr %19, align 8, !tbaa !21
  %881 = load ptr, ptr %9, align 8, !tbaa !5
  %882 = getelementptr inbounds %struct.perturbs, ptr %881, i32 0, i32 158
  %883 = load ptr, ptr %882, align 8, !tbaa !44
  %884 = load i32, ptr %10, align 4, !tbaa !14
  %885 = sext i32 %884 to i64
  %886 = getelementptr inbounds double, ptr %883, i64 %885
  store double %880, ptr %886, align 8, !tbaa !21
  br label %700, !llvm.loop !45

887:                                              ; preds = %700
  %888 = load ptr, ptr %7, align 8, !tbaa !5
  %889 = getelementptr inbounds %struct.background, ptr %888, i32 0, i32 44
  %890 = load double, ptr %889, align 8, !tbaa !35
  %891 = load ptr, ptr %9, align 8, !tbaa !5
  %892 = getelementptr inbounds %struct.perturbs, ptr %891, i32 0, i32 158
  %893 = load ptr, ptr %892, align 8, !tbaa !44
  %894 = load i32, ptr %10, align 4, !tbaa !14
  %895 = sext i32 %894 to i64
  %896 = getelementptr inbounds double, ptr %893, i64 %895
  store double %890, ptr %896, align 8, !tbaa !21
  %897 = load ptr, ptr %29, align 8, !tbaa !5
  call void @free(ptr noundef %897) #6
  %898 = load ptr, ptr %30, align 8, !tbaa !5
  call void @free(ptr noundef %898) #6
  %899 = load ptr, ptr %9, align 8, !tbaa !5
  %900 = getelementptr inbounds %struct.perturbs, ptr %899, i32 0, i32 49
  %901 = load double, ptr %900, align 8, !tbaa !34
  %902 = fcmp olt double %901, 0.000000e+00
  br i1 %902, label %903, label %914

903:                                              ; preds = %887
  call void @llvm.lifetime.start.p0(i64 2048, ptr %63) #6
  %904 = getelementptr inbounds [2048 x i8], ptr %63, i64 0, i64 0
  %905 = load ptr, ptr %9, align 8, !tbaa !5
  %906 = getelementptr inbounds %struct.perturbs, ptr %905, i32 0, i32 49
  %907 = load double, ptr %906, align 8, !tbaa !34
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %904, ptr noundef @.str.147, double noundef %907)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %64) #6
  %908 = getelementptr inbounds [2048 x i8], ptr %64, i64 0, i64 0
  %909 = getelementptr inbounds [2048 x i8], ptr %63, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %908, ptr noundef @.str.7, ptr noundef @.str.148, ptr noundef %909)
  %910 = load ptr, ptr %9, align 8, !tbaa !5
  %911 = getelementptr inbounds %struct.perturbs, ptr %910, i32 0, i32 186
  %912 = getelementptr inbounds [2048 x i8], ptr %911, i64 0, i64 0
  %913 = getelementptr inbounds [2048 x i8], ptr %64, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %912, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1752, ptr noundef %913)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %64) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %63) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

914:                                              ; preds = %887
  %915 = load ptr, ptr %9, align 8, !tbaa !5
  %916 = getelementptr inbounds %struct.perturbs, ptr %915, i32 0, i32 49
  %917 = load double, ptr %916, align 8, !tbaa !34
  %918 = fcmp oeq double %917, 0.000000e+00
  br i1 %918, label %919, label %922

919:                                              ; preds = %914
  %920 = load ptr, ptr %9, align 8, !tbaa !5
  %921 = getelementptr inbounds %struct.perturbs, ptr %920, i32 0, i32 169
  store i32 1, ptr %921, align 8, !tbaa !46
  br label %1083

922:                                              ; preds = %914
  %923 = load ptr, ptr %7, align 8, !tbaa !5
  %924 = load ptr, ptr %9, align 8, !tbaa !5
  %925 = getelementptr inbounds %struct.perturbs, ptr %924, i32 0, i32 49
  %926 = load double, ptr %925, align 8, !tbaa !34
  %927 = call i32 @background_tau_of_z(ptr noundef %923, double noundef %926, ptr noundef %21)
  %928 = icmp eq i32 %927, 1
  br i1 %928, label %929, label %938

929:                                              ; preds = %922
  call void @llvm.lifetime.start.p0(i64 2048, ptr %65) #6
  %930 = getelementptr inbounds [2048 x i8], ptr %65, i64 0, i64 0
  %931 = load ptr, ptr %7, align 8, !tbaa !5
  %932 = getelementptr inbounds %struct.background, ptr %931, i32 0, i32 149
  %933 = getelementptr inbounds [2048 x i8], ptr %932, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %930, ptr noundef @.str, ptr noundef @.str.149, ptr noundef %933)
  %934 = load ptr, ptr %9, align 8, !tbaa !5
  %935 = getelementptr inbounds %struct.perturbs, ptr %934, i32 0, i32 186
  %936 = getelementptr inbounds [2048 x i8], ptr %935, i64 0, i64 0
  %937 = getelementptr inbounds [2048 x i8], ptr %65, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %936, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1765, ptr noundef %937)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %65) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

938:                                              ; preds = %922
  store i32 0, ptr %14, align 4, !tbaa !14
  %939 = load double, ptr %21, align 8, !tbaa !21
  %940 = load ptr, ptr %9, align 8, !tbaa !5
  %941 = getelementptr inbounds %struct.perturbs, ptr %940, i32 0, i32 158
  %942 = load ptr, ptr %941, align 8, !tbaa !44
  %943 = load i32, ptr %14, align 4, !tbaa !14
  %944 = sext i32 %943 to i64
  %945 = getelementptr inbounds double, ptr %942, i64 %944
  %946 = load double, ptr %945, align 8, !tbaa !21
  %947 = fcmp ole double %939, %946
  br i1 %947, label %948, label %965

948:                                              ; preds = %938
  call void @llvm.lifetime.start.p0(i64 2048, ptr %66) #6
  %949 = getelementptr inbounds [2048 x i8], ptr %66, i64 0, i64 0
  %950 = load ptr, ptr %9, align 8, !tbaa !5
  %951 = getelementptr inbounds %struct.perturbs, ptr %950, i32 0, i32 49
  %952 = load double, ptr %951, align 8, !tbaa !34
  %953 = load double, ptr %21, align 8, !tbaa !21
  %954 = load ptr, ptr %9, align 8, !tbaa !5
  %955 = getelementptr inbounds %struct.perturbs, ptr %954, i32 0, i32 158
  %956 = load ptr, ptr %955, align 8, !tbaa !44
  %957 = getelementptr inbounds double, ptr %956, i64 0
  %958 = load double, ptr %957, align 8, !tbaa !21
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %949, ptr noundef @.str.150, double noundef %952, double noundef %953, double noundef %958)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %67) #6
  %959 = getelementptr inbounds [2048 x i8], ptr %67, i64 0, i64 0
  %960 = getelementptr inbounds [2048 x i8], ptr %66, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %959, ptr noundef @.str.7, ptr noundef @.str.151, ptr noundef %960)
  %961 = load ptr, ptr %9, align 8, !tbaa !5
  %962 = getelementptr inbounds %struct.perturbs, ptr %961, i32 0, i32 186
  %963 = getelementptr inbounds [2048 x i8], ptr %962, i64 0, i64 0
  %964 = getelementptr inbounds [2048 x i8], ptr %67, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %963, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1770, ptr noundef %964)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %67) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %66) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

965:                                              ; preds = %938
  br label %966

966:                                              ; preds = %976, %965
  %967 = load ptr, ptr %9, align 8, !tbaa !5
  %968 = getelementptr inbounds %struct.perturbs, ptr %967, i32 0, i32 158
  %969 = load ptr, ptr %968, align 8, !tbaa !44
  %970 = load i32, ptr %14, align 4, !tbaa !14
  %971 = sext i32 %970 to i64
  %972 = getelementptr inbounds double, ptr %969, i64 %971
  %973 = load double, ptr %972, align 8, !tbaa !21
  %974 = load double, ptr %21, align 8, !tbaa !21
  %975 = fcmp olt double %973, %974
  br i1 %975, label %976, label %979

976:                                              ; preds = %966
  %977 = load i32, ptr %14, align 4, !tbaa !14
  %978 = add nsw i32 %977, 1
  store i32 %978, ptr %14, align 4, !tbaa !14
  br label %966, !llvm.loop !47

979:                                              ; preds = %966
  %980 = load i32, ptr %14, align 4, !tbaa !14
  %981 = add nsw i32 %980, -1
  store i32 %981, ptr %14, align 4, !tbaa !14
  %982 = load i32, ptr %14, align 4, !tbaa !14
  %983 = icmp slt i32 %982, 0
  br i1 %983, label %984, label %992

984:                                              ; preds = %979
  call void @llvm.lifetime.start.p0(i64 2048, ptr %68) #6
  %985 = getelementptr inbounds [2048 x i8], ptr %68, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %985, ptr noundef @.str.152)
  call void @llvm.lifetime.start.p0(i64 2048, ptr %69) #6
  %986 = getelementptr inbounds [2048 x i8], ptr %69, i64 0, i64 0
  %987 = getelementptr inbounds [2048 x i8], ptr %68, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %986, ptr noundef @.str.7, ptr noundef @.str.153, ptr noundef %987)
  %988 = load ptr, ptr %9, align 8, !tbaa !5
  %989 = getelementptr inbounds %struct.perturbs, ptr %988, i32 0, i32 186
  %990 = getelementptr inbounds [2048 x i8], ptr %989, i64 0, i64 0
  %991 = getelementptr inbounds [2048 x i8], ptr %69, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %990, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1778, ptr noundef %991)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %69) #6
  call void @llvm.lifetime.end.p0(i64 2048, ptr %68) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

992:                                              ; preds = %979
  %993 = load i32, ptr %14, align 4, !tbaa !14
  %994 = icmp sgt i32 %993, 0
  br i1 %994, label %995, label %998

995:                                              ; preds = %992
  %996 = load i32, ptr %14, align 4, !tbaa !14
  %997 = add nsw i32 %996, -1
  store i32 %997, ptr %14, align 4, !tbaa !14
  br label %998

998:                                              ; preds = %995, %992
  %999 = load i32, ptr %14, align 4, !tbaa !14
  %1000 = icmp sgt i32 %999, 0
  br i1 %1000, label %1001, label %1004

1001:                                             ; preds = %998
  %1002 = load i32, ptr %14, align 4, !tbaa !14
  %1003 = add nsw i32 %1002, -1
  store i32 %1003, ptr %14, align 4, !tbaa !14
  br label %1004

1004:                                             ; preds = %1001, %998
  %1005 = load i32, ptr %14, align 4, !tbaa !14
  %1006 = icmp sgt i32 %1005, 0
  br i1 %1006, label %1007, label %1010

1007:                                             ; preds = %1004
  %1008 = load i32, ptr %14, align 4, !tbaa !14
  %1009 = add nsw i32 %1008, -1
  store i32 %1009, ptr %14, align 4, !tbaa !14
  br label %1010

1010:                                             ; preds = %1007, %1004
  %1011 = load i32, ptr %14, align 4, !tbaa !14
  %1012 = icmp sgt i32 %1011, 0
  br i1 %1012, label %1013, label %1016

1013:                                             ; preds = %1010
  %1014 = load i32, ptr %14, align 4, !tbaa !14
  %1015 = add nsw i32 %1014, -1
  store i32 %1015, ptr %14, align 4, !tbaa !14
  br label %1016

1016:                                             ; preds = %1013, %1010
  %1017 = load ptr, ptr %9, align 8, !tbaa !5
  %1018 = getelementptr inbounds %struct.perturbs, ptr %1017, i32 0, i32 159
  %1019 = load i32, ptr %1018, align 8, !tbaa !43
  %1020 = load i32, ptr %14, align 4, !tbaa !14
  %1021 = sub nsw i32 %1019, %1020
  %1022 = load ptr, ptr %9, align 8, !tbaa !5
  %1023 = getelementptr inbounds %struct.perturbs, ptr %1022, i32 0, i32 169
  store i32 %1021, ptr %1023, align 8, !tbaa !46
  %1024 = load ptr, ptr %9, align 8, !tbaa !5
  %1025 = getelementptr inbounds %struct.perturbs, ptr %1024, i32 0, i32 169
  %1026 = load i32, ptr %1025, align 8, !tbaa !46
  %1027 = sext i32 %1026 to i64
  %1028 = mul i64 %1027, 8
  %1029 = call noalias ptr @malloc(i64 noundef %1028) #7
  %1030 = load ptr, ptr %9, align 8, !tbaa !5
  %1031 = getelementptr inbounds %struct.perturbs, ptr %1030, i32 0, i32 168
  store ptr %1029, ptr %1031, align 8, !tbaa !48
  %1032 = load ptr, ptr %9, align 8, !tbaa !5
  %1033 = getelementptr inbounds %struct.perturbs, ptr %1032, i32 0, i32 168
  %1034 = load ptr, ptr %1033, align 8, !tbaa !48
  %1035 = icmp eq ptr %1034, null
  br i1 %1035, label %1036, label %1049

1036:                                             ; preds = %1016
  call void @llvm.lifetime.start.p0(i64 4, ptr %70) #6
  %1037 = load ptr, ptr %9, align 8, !tbaa !5
  %1038 = getelementptr inbounds %struct.perturbs, ptr %1037, i32 0, i32 169
  %1039 = load i32, ptr %1038, align 8, !tbaa !46
  %1040 = sext i32 %1039 to i64
  %1041 = mul i64 %1040, 8
  %1042 = trunc i64 %1041 to i32
  store i32 %1042, ptr %70, align 4, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 2048, ptr %71) #6
  %1043 = getelementptr inbounds [2048 x i8], ptr %71, i64 0, i64 0
  %1044 = load i32, ptr %70, align 4, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1043, ptr noundef @.str.4, ptr noundef @.str.154, i32 noundef %1044)
  %1045 = load ptr, ptr %9, align 8, !tbaa !5
  %1046 = getelementptr inbounds %struct.perturbs, ptr %1045, i32 0, i32 186
  %1047 = getelementptr inbounds [2048 x i8], ptr %1046, i64 0, i64 0
  %1048 = getelementptr inbounds [2048 x i8], ptr %71, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1047, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1788, ptr noundef %1048)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %71) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %70) #6
  br label %1345

1049:                                             ; preds = %1016
  store i32 0, ptr %14, align 4, !tbaa !14
  br label %1050

1050:                                             ; preds = %1079, %1049
  %1051 = load i32, ptr %14, align 4, !tbaa !14
  %1052 = load ptr, ptr %9, align 8, !tbaa !5
  %1053 = getelementptr inbounds %struct.perturbs, ptr %1052, i32 0, i32 169
  %1054 = load i32, ptr %1053, align 8, !tbaa !46
  %1055 = icmp slt i32 %1051, %1054
  br i1 %1055, label %1056, label %1082

1056:                                             ; preds = %1050
  %1057 = load ptr, ptr %9, align 8, !tbaa !5
  %1058 = getelementptr inbounds %struct.perturbs, ptr %1057, i32 0, i32 158
  %1059 = load ptr, ptr %1058, align 8, !tbaa !44
  %1060 = load i32, ptr %14, align 4, !tbaa !14
  %1061 = load ptr, ptr %9, align 8, !tbaa !5
  %1062 = getelementptr inbounds %struct.perturbs, ptr %1061, i32 0, i32 169
  %1063 = load i32, ptr %1062, align 8, !tbaa !46
  %1064 = sub nsw i32 %1060, %1063
  %1065 = load ptr, ptr %9, align 8, !tbaa !5
  %1066 = getelementptr inbounds %struct.perturbs, ptr %1065, i32 0, i32 159
  %1067 = load i32, ptr %1066, align 8, !tbaa !43
  %1068 = add nsw i32 %1064, %1067
  %1069 = sext i32 %1068 to i64
  %1070 = getelementptr inbounds double, ptr %1059, i64 %1069
  %1071 = load double, ptr %1070, align 8, !tbaa !21
  %1072 = call double @log(double noundef %1071) #6
  %1073 = load ptr, ptr %9, align 8, !tbaa !5
  %1074 = getelementptr inbounds %struct.perturbs, ptr %1073, i32 0, i32 168
  %1075 = load ptr, ptr %1074, align 8, !tbaa !48
  %1076 = load i32, ptr %14, align 4, !tbaa !14
  %1077 = sext i32 %1076 to i64
  %1078 = getelementptr inbounds double, ptr %1075, i64 %1077
  store double %1072, ptr %1078, align 8, !tbaa !21
  br label %1079

1079:                                             ; preds = %1056
  %1080 = load i32, ptr %14, align 4, !tbaa !14
  %1081 = add nsw i32 %1080, 1
  store i32 %1081, ptr %14, align 4, !tbaa !14
  br label %1050, !llvm.loop !49

1082:                                             ; preds = %1050
  br label %1083

1083:                                             ; preds = %1082, %919
  store i32 0, ptr %11, align 4, !tbaa !14
  br label %1084

1084:                                             ; preds = %1341, %1083
  %1085 = load i32, ptr %11, align 4, !tbaa !14
  %1086 = load ptr, ptr %9, align 8, !tbaa !5
  %1087 = getelementptr inbounds %struct.perturbs, ptr %1086, i32 0, i32 59
  %1088 = load i32, ptr %1087, align 8, !tbaa !50
  %1089 = icmp slt i32 %1085, %1088
  br i1 %1089, label %1090, label %1344

1090:                                             ; preds = %1084
  store i32 0, ptr %13, align 4, !tbaa !14
  br label %1091

1091:                                             ; preds = %1337, %1090
  %1092 = load i32, ptr %13, align 4, !tbaa !14
  %1093 = load ptr, ptr %9, align 8, !tbaa !5
  %1094 = getelementptr inbounds %struct.perturbs, ptr %1093, i32 0, i32 66
  %1095 = load ptr, ptr %1094, align 8, !tbaa !51
  %1096 = load i32, ptr %11, align 4, !tbaa !14
  %1097 = sext i32 %1096 to i64
  %1098 = getelementptr inbounds i32, ptr %1095, i64 %1097
  %1099 = load i32, ptr %1098, align 4, !tbaa !14
  %1100 = icmp slt i32 %1092, %1099
  br i1 %1100, label %1101, label %1340

1101:                                             ; preds = %1091
  store i32 0, ptr %12, align 4, !tbaa !14
  br label %1102

1102:                                             ; preds = %1333, %1101
  %1103 = load i32, ptr %12, align 4, !tbaa !14
  %1104 = load ptr, ptr %9, align 8, !tbaa !5
  %1105 = getelementptr inbounds %struct.perturbs, ptr %1104, i32 0, i32 151
  %1106 = load ptr, ptr %1105, align 8, !tbaa !52
  %1107 = load i32, ptr %11, align 4, !tbaa !14
  %1108 = sext i32 %1107 to i64
  %1109 = getelementptr inbounds i32, ptr %1106, i64 %1108
  %1110 = load i32, ptr %1109, align 4, !tbaa !14
  %1111 = icmp slt i32 %1103, %1110
  br i1 %1111, label %1112, label %1336

1112:                                             ; preds = %1102
  %1113 = load ptr, ptr %9, align 8, !tbaa !5
  %1114 = getelementptr inbounds %struct.perturbs, ptr %1113, i32 0, i32 154
  %1115 = load ptr, ptr %1114, align 8, !tbaa !53
  %1116 = load i32, ptr %11, align 4, !tbaa !14
  %1117 = sext i32 %1116 to i64
  %1118 = getelementptr inbounds i32, ptr %1115, i64 %1117
  %1119 = load i32, ptr %1118, align 4, !tbaa !14
  %1120 = load ptr, ptr %9, align 8, !tbaa !5
  %1121 = getelementptr inbounds %struct.perturbs, ptr %1120, i32 0, i32 159
  %1122 = load i32, ptr %1121, align 8, !tbaa !43
  %1123 = mul nsw i32 %1119, %1122
  %1124 = sext i32 %1123 to i64
  %1125 = mul i64 %1124, 8
  %1126 = call noalias ptr @malloc(i64 noundef %1125) #7
  %1127 = load ptr, ptr %9, align 8, !tbaa !5
  %1128 = getelementptr inbounds %struct.perturbs, ptr %1127, i32 0, i32 167
  %1129 = load ptr, ptr %1128, align 8, !tbaa !54
  %1130 = load i32, ptr %11, align 4, !tbaa !14
  %1131 = sext i32 %1130 to i64
  %1132 = getelementptr inbounds ptr, ptr %1129, i64 %1131
  %1133 = load ptr, ptr %1132, align 8, !tbaa !5
  %1134 = load i32, ptr %13, align 4, !tbaa !14
  %1135 = load ptr, ptr %9, align 8, !tbaa !5
  %1136 = getelementptr inbounds %struct.perturbs, ptr %1135, i32 0, i32 151
  %1137 = load ptr, ptr %1136, align 8, !tbaa !52
  %1138 = load i32, ptr %11, align 4, !tbaa !14
  %1139 = sext i32 %1138 to i64
  %1140 = getelementptr inbounds i32, ptr %1137, i64 %1139
  %1141 = load i32, ptr %1140, align 4, !tbaa !14
  %1142 = mul nsw i32 %1134, %1141
  %1143 = load i32, ptr %12, align 4, !tbaa !14
  %1144 = add nsw i32 %1142, %1143
  %1145 = sext i32 %1144 to i64
  %1146 = getelementptr inbounds ptr, ptr %1133, i64 %1145
  store ptr %1126, ptr %1146, align 8, !tbaa !5
  %1147 = load ptr, ptr %9, align 8, !tbaa !5
  %1148 = getelementptr inbounds %struct.perturbs, ptr %1147, i32 0, i32 167
  %1149 = load ptr, ptr %1148, align 8, !tbaa !54
  %1150 = load i32, ptr %11, align 4, !tbaa !14
  %1151 = sext i32 %1150 to i64
  %1152 = getelementptr inbounds ptr, ptr %1149, i64 %1151
  %1153 = load ptr, ptr %1152, align 8, !tbaa !5
  %1154 = load i32, ptr %13, align 4, !tbaa !14
  %1155 = load ptr, ptr %9, align 8, !tbaa !5
  %1156 = getelementptr inbounds %struct.perturbs, ptr %1155, i32 0, i32 151
  %1157 = load ptr, ptr %1156, align 8, !tbaa !52
  %1158 = load i32, ptr %11, align 4, !tbaa !14
  %1159 = sext i32 %1158 to i64
  %1160 = getelementptr inbounds i32, ptr %1157, i64 %1159
  %1161 = load i32, ptr %1160, align 4, !tbaa !14
  %1162 = mul nsw i32 %1154, %1161
  %1163 = load i32, ptr %12, align 4, !tbaa !14
  %1164 = add nsw i32 %1162, %1163
  %1165 = sext i32 %1164 to i64
  %1166 = getelementptr inbounds ptr, ptr %1153, i64 %1165
  %1167 = load ptr, ptr %1166, align 8, !tbaa !5
  %1168 = icmp eq ptr %1167, null
  br i1 %1168, label %1169, label %1190

1169:                                             ; preds = %1112
  call void @llvm.lifetime.start.p0(i64 4, ptr %72) #6
  %1170 = load ptr, ptr %9, align 8, !tbaa !5
  %1171 = getelementptr inbounds %struct.perturbs, ptr %1170, i32 0, i32 154
  %1172 = load ptr, ptr %1171, align 8, !tbaa !53
  %1173 = load i32, ptr %11, align 4, !tbaa !14
  %1174 = sext i32 %1173 to i64
  %1175 = getelementptr inbounds i32, ptr %1172, i64 %1174
  %1176 = load i32, ptr %1175, align 4, !tbaa !14
  %1177 = load ptr, ptr %9, align 8, !tbaa !5
  %1178 = getelementptr inbounds %struct.perturbs, ptr %1177, i32 0, i32 159
  %1179 = load i32, ptr %1178, align 8, !tbaa !43
  %1180 = mul nsw i32 %1176, %1179
  %1181 = sext i32 %1180 to i64
  %1182 = mul i64 %1181, 8
  %1183 = trunc i64 %1182 to i32
  store i32 %1183, ptr %72, align 4, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 2048, ptr %73) #6
  %1184 = getelementptr inbounds [2048 x i8], ptr %73, i64 0, i64 0
  %1185 = load i32, ptr %72, align 4, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1184, ptr noundef @.str.4, ptr noundef @.str.155, i32 noundef %1185)
  %1186 = load ptr, ptr %9, align 8, !tbaa !5
  %1187 = getelementptr inbounds %struct.perturbs, ptr %1186, i32 0, i32 186
  %1188 = getelementptr inbounds [2048 x i8], ptr %1187, i64 0, i64 0
  %1189 = getelementptr inbounds [2048 x i8], ptr %73, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1188, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1804, ptr noundef %1189)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %73) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %72) #6
  br label %1345

1190:                                             ; preds = %1112
  %1191 = load ptr, ptr %9, align 8, !tbaa !5
  %1192 = getelementptr inbounds %struct.perturbs, ptr %1191, i32 0, i32 169
  %1193 = load i32, ptr %1192, align 8, !tbaa !46
  %1194 = icmp sgt i32 %1193, 1
  br i1 %1194, label %1195, label %1332

1195:                                             ; preds = %1190
  %1196 = load ptr, ptr %9, align 8, !tbaa !5
  %1197 = getelementptr inbounds %struct.perturbs, ptr %1196, i32 0, i32 167
  %1198 = load ptr, ptr %1197, align 8, !tbaa !54
  %1199 = load i32, ptr %11, align 4, !tbaa !14
  %1200 = sext i32 %1199 to i64
  %1201 = getelementptr inbounds ptr, ptr %1198, i64 %1200
  %1202 = load ptr, ptr %1201, align 8, !tbaa !5
  %1203 = load i32, ptr %13, align 4, !tbaa !14
  %1204 = load ptr, ptr %9, align 8, !tbaa !5
  %1205 = getelementptr inbounds %struct.perturbs, ptr %1204, i32 0, i32 151
  %1206 = load ptr, ptr %1205, align 8, !tbaa !52
  %1207 = load i32, ptr %11, align 4, !tbaa !14
  %1208 = sext i32 %1207 to i64
  %1209 = getelementptr inbounds i32, ptr %1206, i64 %1208
  %1210 = load i32, ptr %1209, align 4, !tbaa !14
  %1211 = mul nsw i32 %1203, %1210
  %1212 = load i32, ptr %12, align 4, !tbaa !14
  %1213 = add nsw i32 %1211, %1212
  %1214 = sext i32 %1213 to i64
  %1215 = getelementptr inbounds ptr, ptr %1202, i64 %1214
  %1216 = load ptr, ptr %1215, align 8, !tbaa !5
  %1217 = load ptr, ptr %9, align 8, !tbaa !5
  %1218 = getelementptr inbounds %struct.perturbs, ptr %1217, i32 0, i32 159
  %1219 = load i32, ptr %1218, align 8, !tbaa !43
  %1220 = load ptr, ptr %9, align 8, !tbaa !5
  %1221 = getelementptr inbounds %struct.perturbs, ptr %1220, i32 0, i32 169
  %1222 = load i32, ptr %1221, align 8, !tbaa !46
  %1223 = sub nsw i32 %1219, %1222
  %1224 = load ptr, ptr %9, align 8, !tbaa !5
  %1225 = getelementptr inbounds %struct.perturbs, ptr %1224, i32 0, i32 154
  %1226 = load ptr, ptr %1225, align 8, !tbaa !53
  %1227 = load i32, ptr %11, align 4, !tbaa !14
  %1228 = sext i32 %1227 to i64
  %1229 = getelementptr inbounds i32, ptr %1226, i64 %1228
  %1230 = load i32, ptr %1229, align 4, !tbaa !14
  %1231 = mul nsw i32 %1223, %1230
  %1232 = sext i32 %1231 to i64
  %1233 = getelementptr inbounds double, ptr %1216, i64 %1232
  %1234 = load ptr, ptr %9, align 8, !tbaa !5
  %1235 = getelementptr inbounds %struct.perturbs, ptr %1234, i32 0, i32 170
  %1236 = load ptr, ptr %1235, align 8, !tbaa !55
  %1237 = load i32, ptr %11, align 4, !tbaa !14
  %1238 = sext i32 %1237 to i64
  %1239 = getelementptr inbounds ptr, ptr %1236, i64 %1238
  %1240 = load ptr, ptr %1239, align 8, !tbaa !5
  %1241 = load i32, ptr %13, align 4, !tbaa !14
  %1242 = load ptr, ptr %9, align 8, !tbaa !5
  %1243 = getelementptr inbounds %struct.perturbs, ptr %1242, i32 0, i32 151
  %1244 = load ptr, ptr %1243, align 8, !tbaa !52
  %1245 = load i32, ptr %11, align 4, !tbaa !14
  %1246 = sext i32 %1245 to i64
  %1247 = getelementptr inbounds i32, ptr %1244, i64 %1246
  %1248 = load i32, ptr %1247, align 4, !tbaa !14
  %1249 = mul nsw i32 %1241, %1248
  %1250 = load i32, ptr %12, align 4, !tbaa !14
  %1251 = add nsw i32 %1249, %1250
  %1252 = sext i32 %1251 to i64
  %1253 = getelementptr inbounds ptr, ptr %1240, i64 %1252
  store ptr %1233, ptr %1253, align 8, !tbaa !5
  %1254 = load ptr, ptr %9, align 8, !tbaa !5
  %1255 = getelementptr inbounds %struct.perturbs, ptr %1254, i32 0, i32 154
  %1256 = load ptr, ptr %1255, align 8, !tbaa !53
  %1257 = load i32, ptr %11, align 4, !tbaa !14
  %1258 = sext i32 %1257 to i64
  %1259 = getelementptr inbounds i32, ptr %1256, i64 %1258
  %1260 = load i32, ptr %1259, align 4, !tbaa !14
  %1261 = load ptr, ptr %9, align 8, !tbaa !5
  %1262 = getelementptr inbounds %struct.perturbs, ptr %1261, i32 0, i32 169
  %1263 = load i32, ptr %1262, align 8, !tbaa !46
  %1264 = mul nsw i32 %1260, %1263
  %1265 = sext i32 %1264 to i64
  %1266 = mul i64 %1265, 8
  %1267 = call noalias ptr @malloc(i64 noundef %1266) #7
  %1268 = load ptr, ptr %9, align 8, !tbaa !5
  %1269 = getelementptr inbounds %struct.perturbs, ptr %1268, i32 0, i32 171
  %1270 = load ptr, ptr %1269, align 8, !tbaa !56
  %1271 = load i32, ptr %11, align 4, !tbaa !14
  %1272 = sext i32 %1271 to i64
  %1273 = getelementptr inbounds ptr, ptr %1270, i64 %1272
  %1274 = load ptr, ptr %1273, align 8, !tbaa !5
  %1275 = load i32, ptr %13, align 4, !tbaa !14
  %1276 = load ptr, ptr %9, align 8, !tbaa !5
  %1277 = getelementptr inbounds %struct.perturbs, ptr %1276, i32 0, i32 151
  %1278 = load ptr, ptr %1277, align 8, !tbaa !52
  %1279 = load i32, ptr %11, align 4, !tbaa !14
  %1280 = sext i32 %1279 to i64
  %1281 = getelementptr inbounds i32, ptr %1278, i64 %1280
  %1282 = load i32, ptr %1281, align 4, !tbaa !14
  %1283 = mul nsw i32 %1275, %1282
  %1284 = load i32, ptr %12, align 4, !tbaa !14
  %1285 = add nsw i32 %1283, %1284
  %1286 = sext i32 %1285 to i64
  %1287 = getelementptr inbounds ptr, ptr %1274, i64 %1286
  store ptr %1267, ptr %1287, align 8, !tbaa !5
  %1288 = load ptr, ptr %9, align 8, !tbaa !5
  %1289 = getelementptr inbounds %struct.perturbs, ptr %1288, i32 0, i32 171
  %1290 = load ptr, ptr %1289, align 8, !tbaa !56
  %1291 = load i32, ptr %11, align 4, !tbaa !14
  %1292 = sext i32 %1291 to i64
  %1293 = getelementptr inbounds ptr, ptr %1290, i64 %1292
  %1294 = load ptr, ptr %1293, align 8, !tbaa !5
  %1295 = load i32, ptr %13, align 4, !tbaa !14
  %1296 = load ptr, ptr %9, align 8, !tbaa !5
  %1297 = getelementptr inbounds %struct.perturbs, ptr %1296, i32 0, i32 151
  %1298 = load ptr, ptr %1297, align 8, !tbaa !52
  %1299 = load i32, ptr %11, align 4, !tbaa !14
  %1300 = sext i32 %1299 to i64
  %1301 = getelementptr inbounds i32, ptr %1298, i64 %1300
  %1302 = load i32, ptr %1301, align 4, !tbaa !14
  %1303 = mul nsw i32 %1295, %1302
  %1304 = load i32, ptr %12, align 4, !tbaa !14
  %1305 = add nsw i32 %1303, %1304
  %1306 = sext i32 %1305 to i64
  %1307 = getelementptr inbounds ptr, ptr %1294, i64 %1306
  %1308 = load ptr, ptr %1307, align 8, !tbaa !5
  %1309 = icmp eq ptr %1308, null
  br i1 %1309, label %1310, label %1331

1310:                                             ; preds = %1195
  call void @llvm.lifetime.start.p0(i64 4, ptr %74) #6
  %1311 = load ptr, ptr %9, align 8, !tbaa !5
  %1312 = getelementptr inbounds %struct.perturbs, ptr %1311, i32 0, i32 154
  %1313 = load ptr, ptr %1312, align 8, !tbaa !53
  %1314 = load i32, ptr %11, align 4, !tbaa !14
  %1315 = sext i32 %1314 to i64
  %1316 = getelementptr inbounds i32, ptr %1313, i64 %1315
  %1317 = load i32, ptr %1316, align 4, !tbaa !14
  %1318 = load ptr, ptr %9, align 8, !tbaa !5
  %1319 = getelementptr inbounds %struct.perturbs, ptr %1318, i32 0, i32 169
  %1320 = load i32, ptr %1319, align 8, !tbaa !46
  %1321 = mul nsw i32 %1317, %1320
  %1322 = sext i32 %1321 to i64
  %1323 = mul i64 %1322, 8
  %1324 = trunc i64 %1323 to i32
  store i32 %1324, ptr %74, align 4, !tbaa !14
  call void @llvm.lifetime.start.p0(i64 2048, ptr %75) #6
  %1325 = getelementptr inbounds [2048 x i8], ptr %75, i64 0, i64 0
  %1326 = load i32, ptr %74, align 4, !tbaa !14
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1325, ptr noundef @.str.4, ptr noundef @.str.156, i32 noundef %1326)
  %1327 = load ptr, ptr %9, align 8, !tbaa !5
  %1328 = getelementptr inbounds %struct.perturbs, ptr %1327, i32 0, i32 186
  %1329 = getelementptr inbounds [2048 x i8], ptr %1328, i64 0, i64 0
  %1330 = getelementptr inbounds [2048 x i8], ptr %75, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %1329, ptr noundef @.str.2, ptr noundef @__func__.perturb_timesampling_for_sources, i32 noundef 1814, ptr noundef %1330)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %75) #6
  store i32 1, ptr %5, align 4
  store i32 1, ptr %33, align 4
  call void @llvm.lifetime.end.p0(i64 4, ptr %74) #6
  br label %1345

1331:                                             ; preds = %1195
  br label %1332

1332:                                             ; preds = %1331, %1190
  br label %1333

1333:                                             ; preds = %1332
  %1334 = load i32, ptr %12, align 4, !tbaa !14
  %1335 = add nsw i32 %1334, 1
  store i32 %1335, ptr %12, align 4, !tbaa !14
  br label %1102, !llvm.loop !57

1336:                                             ; preds = %1102
  br label %1337

1337:                                             ; preds = %1336
  %1338 = load i32, ptr %13, align 4, !tbaa !14
  %1339 = add nsw i32 %1338, 1
  store i32 %1339, ptr %13, align 4, !tbaa !14
  br label %1091, !llvm.loop !58

1340:                                             ; preds = %1091
  br label %1341

1341:                                             ; preds = %1340
  %1342 = load i32, ptr %11, align 4, !tbaa !14
  %1343 = add nsw i32 %1342, 1
  store i32 %1343, ptr %11, align 4, !tbaa !14
  br label %1084, !llvm.loop !59

1344:                                             ; preds = %1084
  store i32 0, ptr %5, align 4
  store i32 1, ptr %33, align 4
  br label %1345

1345:                                             ; preds = %1344, %1310, %1169, %1036, %984, %948, %929, %903, %858, %835, %746, %718, %676, %638, %615, %526, %498, %466, %438, %417, %361, %333, %300, %263, %235, %211, %174, %146, %106, %84
  call void @llvm.lifetime.end.p0(i64 8, ptr %30) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %29) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %26) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %22) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #6
  call void @llvm.lifetime.end.p0(i64 8, ptr %19) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %18) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %17) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %16) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %15) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %14) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %13) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %12) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %11) #6
  call void @llvm.lifetime.end.p0(i64 4, ptr %10) #6
  %1346 = load i32, ptr %5, align 4
  ret i32 %1346
}

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

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

declare i32 @background_at_tau(ptr noundef, double noundef, i16 noundef signext, i16 noundef signext, ptr noundef, ptr noundef) #2

declare i32 @thermodynamics_at_z(ptr noundef, ptr noundef, double noundef, i16 noundef signext, ptr noundef, ptr noundef, ptr noundef) #2

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

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

!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, !12, i64 616}
!10 = !{!"background", !11, i64 0, !11, i64 8, !11, i64 16, !11, i64 24, !11, i64 32, !11, i64 40, !11, i64 48, !12, i64 56, !11, i64 64, !11, i64 72, !11, i64 80, !11, i64 88, !13, i64 96, !11, i64 104, !11, i64 112, !11, i64 120, !11, i64 128, !11, i64 136, !11, i64 144, !11, i64 152, !11, i64 160, !11, i64 168, !13, i64 176, !11, i64 184, !11, i64 192, !6, i64 200, !12, i64 208, !12, i64 212, !11, i64 216, !12, i64 224, !6, i64 232, !6, i64 240, !11, i64 248, !6, i64 256, !11, i64 264, !6, i64 272, !11, i64 280, !6, i64 288, !11, i64 296, !6, i64 304, !6, i64 312, !6, i64 320, !11, i64 328, !11, i64 336, !11, i64 344, !11, i64 352, !12, i64 360, !6, i64 368, !11, i64 376, !11, i64 384, !11, i64 392, !11, i64 400, !11, i64 408, !11, i64 416, !11, i64 424, !11, i64 432, !11, i64 440, !11, i64 448, !11, i64 456, !12, i64 464, !12, i64 468, !12, i64 472, !12, i64 476, !12, i64 480, !12, i64 484, !12, i64 488, !12, i64 492, !12, i64 496, !12, i64 500, !12, i64 504, !12, i64 508, !12, i64 512, !12, i64 516, !12, i64 520, !12, i64 524, !12, i64 528, !12, i64 532, !12, i64 536, !12, i64 540, !12, i64 544, !12, i64 548, !12, i64 552, !12, i64 556, !12, i64 560, !12, i64 564, !12, i64 568, !12, i64 572, !12, i64 576, !12, i64 580, !12, i64 584, !12, i64 588, !12, i64 592, !12, i64 596, !12, i64 600, !12, i64 604, !12, i64 608, !12, i64 612, !12, i64 616, !12, i64 620, !12, i64 624, !12, i64 628, !6, i64 632, !6, i64 640, !6, i64 648, !6, i64 656, !6, i64 664, !12, i64 672, !12, i64 676, !12, i64 680, !12, i64 684, !12, i64 688, !12, i64 692, !12, i64 696, !12, i64 700, !12, i64 704, !12, i64 708, !12, i64 712, !12, i64 716, !12, i64 720, !13, i64 724, !13, i64 726, !13, i64 728, !13, i64 730, !13, i64 732, !13, i64 734, !13, i64 736, !13, i64 738, !13, i64 740, !13, i64 742, !13, i64 744, !6, i64 752, !6, i64 760, !6, i64 768, !6, i64 776, !6, i64 784, !6, i64 792, !6, i64 800, !6, i64 808, !6, i64 816, !6, i64 824, !6, i64 832, !13, i64 840, !13, i64 842, !13, i64 844, !13, i64 846, !13, i64 848, !13, i64 850, !7, i64 852, !13, i64 2900, !7, i64 2902}
!11 = !{!"double", !7, i64 0}
!12 = !{!"int", !7, i64 0}
!13 = !{!"short", !7, i64 0}
!14 = !{!12, !12, i64 0}
!15 = !{!16, !12, i64 372}
!16 = !{!"thermo", !11, i64 0, !12, i64 8, !12, i64 12, !12, i64 16, !11, i64 24, !11, i64 32, !13, i64 40, !13, i64 42, !11, i64 48, !11, i64 56, !11, i64 64, !11, i64 72, !12, i64 80, !6, i64 88, !6, i64 96, !11, i64 104, !12, i64 112, !6, i64 120, !6, i64 128, !11, i64 136, !12, i64 144, !6, i64 152, !6, i64 160, !11, i64 168, !13, i64 176, !11, i64 184, !11, i64 192, !11, i64 200, !11, i64 208, !11, i64 216, !11, i64 224, !11, i64 232, !11, i64 240, !11, i64 248, !11, i64 256, !11, i64 264, !12, i64 272, !12, i64 276, !12, i64 280, !12, i64 284, !12, i64 288, !12, i64 292, !12, i64 296, !12, i64 300, !12, i64 304, !12, i64 308, !12, i64 312, !12, i64 316, !12, i64 320, !12, i64 324, !12, i64 328, !12, i64 332, !12, i64 336, !12, i64 340, !12, i64 344, !12, i64 348, !12, i64 352, !12, i64 356, !12, i64 360, !12, i64 364, !12, i64 368, !12, i64 372, !12, i64 376, !6, i64 384, !6, i64 392, !6, i64 400, !11, i64 408, !11, i64 416, !11, i64 424, !11, i64 432, !11, i64 440, !11, i64 448, !11, i64 456, !11, i64 464, !11, i64 472, !11, i64 480, !11, i64 488, !11, i64 496, !11, i64 504, !11, i64 512, !11, i64 520, !11, i64 528, !11, i64 536, !11, i64 544, !11, i64 552, !11, i64 560, !11, i64 568, !11, i64 576, !11, i64 584, !11, i64 592, !13, i64 600, !13, i64 602, !13, i64 604, !7, i64 606}
!17 = !{!18, !13, i64 2020}
!18 = !{!"perturbs", !13, i64 0, !13, i64 2, !13, i64 4, !13, i64 6, !13, i64 8, !13, i64 10, !13, i64 12, !13, i64 14, !13, i64 16, !13, i64 18, !13, i64 20, !12, i64 24, !13, i64 28, !13, i64 30, !13, i64 32, !13, i64 34, !13, i64 36, !13, i64 38, !13, i64 40, !13, i64 42, !13, i64 44, !13, i64 46, !13, i64 48, !13, i64 50, !13, i64 52, !13, i64 54, !13, i64 56, !13, i64 58, !13, i64 60, !12, i64 64, !12, i64 68, !12, i64 72, !12, i64 76, !11, i64 80, !12, i64 88, !12, i64 92, !7, i64 96, !7, i64 896, !12, i64 1696, !12, i64 1700, !12, i64 1704, !12, i64 1708, !12, i64 1712, !11, i64 1720, !12, i64 1728, !12, i64 1732, !7, i64 1736, !11, i64 1976, !11, i64 1984, !11, i64 1992, !6, i64 2000, !6, i64 2008, !12, i64 2016, !13, i64 2020, !13, i64 2022, !12, i64 2024, !12, i64 2028, !12, i64 2032, !12, i64 2036, !12, i64 2040, !12, i64 2044, !12, i64 2048, !12, i64 2052, !12, i64 2056, !12, i64 2060, !12, i64 2064, !6, i64 2072, !13, i64 2080, !13, i64 2082, !13, i64 2084, !13, i64 2086, !13, i64 2088, !13, i64 2090, !13, i64 2092, !13, i64 2094, !13, i64 2096, !13, i64 2098, !13, i64 2100, !13, i64 2102, !13, i64 2104, !13, i64 2106, !13, i64 2108, !13, i64 2110, !13, i64 2112, !13, i64 2114, !13, i64 2116, !13, i64 2118, !13, i64 2120, !13, i64 2122, !13, i64 2124, !13, i64 2126, !13, i64 2128, !13, i64 2130, !13, i64 2132, !13, i64 2134, !13, i64 2136, !13, i64 2138, !13, i64 2140, !13, i64 2142, !13, i64 2144, !13, i64 2146, !13, i64 2148, !13, i64 2150, !13, i64 2152, !13, i64 2154, !13, i64 2156, !13, i64 2158, !12, i64 2160, !12, i64 2164, !12, i64 2168, !12, i64 2172, !12, i64 2176, !12, i64 2180, !12, i64 2184, !12, i64 2188, !12, i64 2192, !12, i64 2196, !12, i64 2200, !12, i64 2204, !12, i64 2208, !12, i64 2212, !12, i64 2216, !12, i64 2220, !12, i64 2224, !12, i64 2228, !12, i64 2232, !12, i64 2236, !12, i64 2240, !12, i64 2244, !12, i64 2248, !12, i64 2252, !12, i64 2256, !12, i64 2260, !12, i64 2264, !12, i64 2268, !12, i64 2272, !12, i64 2276, !12, i64 2280, !12, i64 2284, !12, i64 2288, !12, i64 2292, !12, i64 2296, !12, i64 2300, !12, i64 2304, !12, i64 2308, !12, i64 2312, !12, i64 2316, !12, i64 2320, !12, i64 2324, !12, i64 2328, !12, i64 2332, !6, i64 2336, !6, i64 2344, !6, i64 2352, !6, i64 2360, !6, i64 2368, !11, i64 2376, !11, i64 2384, !6, i64 2392, !12, i64 2400, !11, i64 2408, !11, i64 2416, !11, i64 2424, !6, i64 2432, !6, i64 2440, !6, i64 2448, !6, i64 2456, !6, i64 2464, !6, i64 2472, !12, i64 2480, !6, i64 2488, !6, i64 2496, !6, i64 2504, !7, i64 2512, !7, i64 10512, !7, i64 18512, !12, i64 26512, !12, i64 26516, !12, i64 26520, !7, i64 26528, !7, i64 26768, !7, i64 27008, !7, i64 27248, !7, i64 27368, !7, i64 27488, !13, i64 27608, !7, i64 27610}
!19 = !{!18, !13, i64 20}
!20 = !{!16, !11, i64 584}
!21 = !{!11, !11, i64 0}
!22 = !{!10, !13, i64 840}
!23 = !{!10, !13, i64 846}
!24 = !{!10, !12, i64 464}
!25 = !{!16, !13, i64 600}
!26 = !{!10, !12, i64 468}
!27 = !{!16, !12, i64 276}
!28 = !{!29, !11, i64 1520}
!29 = !{!"precision", !11, i64 0, !11, i64 8, !11, i64 16, !11, i64 24, !11, i64 32, !11, i64 40, !11, i64 48, !11, i64 56, !11, i64 64, !11, i64 72, !11, i64 80, !11, i64 88, !11, i64 96, !7, i64 104, !11, i64 360, !12, i64 368, !11, i64 376, !12, i64 384, !12, i64 388, !11, i64 392, !12, i64 400, !11, i64 408, !12, i64 416, !11, i64 424, !11, i64 432, !11, i64 440, !11, i64 448, !11, i64 456, !11, i64 464, !11, i64 472, !11, i64 480, !11, i64 488, !11, i64 496, !11, i64 504, !11, i64 512, !11, i64 520, !11, i64 528, !11, i64 536, !11, i64 544, !11, i64 552, !11, i64 560, !11, i64 568, !11, i64 576, !11, i64 584, !11, i64 592, !11, i64 600, !11, i64 608, !11, i64 616, !12, i64 624, !7, i64 628, !7, i64 884, !7, i64 1140, !11, i64 1400, !11, i64 1408, !11, i64 1416, !11, i64 1424, !11, i64 1432, !11, i64 1440, !11, i64 1448, !11, i64 1456, !11, i64 1464, !11, i64 1472, !11, i64 1480, !11, i64 1488, !11, i64 1496, !11, i64 1504, !11, i64 1512, !11, i64 1520, !12, i64 1528, !11, i64 1536, !11, i64 1544, !12, i64 1552, !12, i64 1556, !12, i64 1560, !12, i64 1564, !12, i64 1568, !12, i64 1572, !12, i64 1576, !12, i64 1580, !11, i64 1584, !11, i64 1592, !11, i64 1600, !11, i64 1608, !11, i64 1616, !11, i64 1624, !11, i64 1632, !11, i64 1640, !12, i64 1648, !11, i64 1656, !11, i64 1664, !12, i64 1672, !11, i64 1680, !11, i64 1688, !12, i64 1696, !11, i64 1704, !12, i64 1712, !11, i64 1720, !11, i64 1728, !12, i64 1736, !11, i64 1744, !11, i64 1752, !11, i64 1760, !12, i64 1768, !11, i64 1776, !11, i64 1784, !11, i64 1792, !11, i64 1800, !11, i64 1808, !12, i64 1816, !11, i64 1824, !11, i64 1832, !11, i64 1840, !11, i64 1848, !11, i64 1856, !11, i64 1864, !11, i64 1872, !12, i64 1880, !11, i64 1888, !11, i64 1896, !11, i64 1904, !11, i64 1912, !11, i64 1920, !11, i64 1928, !11, i64 1936, !11, i64 1944, !11, i64 1952, !11, i64 1960, !11, i64 1968, !11, i64 1976, !11, i64 1984, !11, i64 1992, !11, i64 2000, !11, i64 2008, !11, i64 2016, !11, i64 2024, !11, i64 2032, !11, i64 2040, !11, i64 2048, !11, i64 2056, !11, i64 2064, !11, i64 2072, !11, i64 2080, !11, i64 2088, !11, i64 2096, !11, i64 2104, !11, i64 2112, !11, i64 2120, !11, i64 2128, !11, i64 2136, !11, i64 2144, !11, i64 2152, !11, i64 2160, !11, i64 2168, !11, i64 2176, !11, i64 2184, !11, i64 2192, !11, i64 2200, !11, i64 2208, !11, i64 2216, !11, i64 2224, !11, i64 2232, !11, i64 2240, !11, i64 2248, !12, i64 2256, !11, i64 2264, !11, i64 2272, !11, i64 2280, !11, i64 2288, !12, i64 2296, !11, i64 2304, !11, i64 2312, !12, i64 2320, !12, i64 2324, !12, i64 2328, !11, i64 2336, !11, i64 2344, !7, i64 2352}
!30 = !{!16, !11, i64 416}
!31 = !{!29, !11, i64 1640}
!32 = distinct !{!32, !33}
!33 = !{!"llvm.loop.mustprogress"}
!34 = !{!18, !11, i64 1992}
!35 = !{!10, !11, i64 344}
!36 = !{!10, !13, i64 848}
!37 = !{!16, !13, i64 602}
!38 = !{!16, !12, i64 364}
!39 = !{!10, !12, i64 472}
!40 = !{!29, !11, i64 1616}
!41 = !{!29, !11, i64 2344}
!42 = distinct !{!42, !33}
!43 = !{!18, !12, i64 2400}
!44 = !{!18, !6, i64 2392}
!45 = distinct !{!45, !33}
!46 = !{!18, !12, i64 2480}
!47 = distinct !{!47, !33}
!48 = !{!18, !6, i64 2472}
!49 = distinct !{!49, !33}
!50 = !{!18, !12, i64 2040}
!51 = !{!18, !6, i64 2072}
!52 = !{!18, !6, i64 2336}
!53 = !{!18, !6, i64 2360}
!54 = !{!18, !6, i64 2464}
!55 = !{!18, !6, i64 2488}
!56 = !{!18, !6, i64 2496}
!57 = distinct !{!57, !33}
!58 = distinct !{!58, !33}
!59 = distinct !{!59, !33}


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

@.str = external hidden unnamed_addr constant [18 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@.str.4 = external hidden unnamed_addr constant [35 x i8], align 1
@.str.7 = external hidden unnamed_addr constant [27 x i8], align 1
@.str.128 = external hidden unnamed_addr constant [9 x i8], align 1
@__func__.perturb_timesampling_for_sources = external hidden unnamed_addr constant [33 x i8], align 1
@.str.129 = external hidden unnamed_addr constant [11 x i8], align 1
@.str.130 = external hidden unnamed_addr constant [99 x i8], align 1
@.str.131 = external hidden unnamed_addr constant [125 x i8], align 1
@.str.132 = external hidden unnamed_addr constant [274 x i8], align 1
@.str.133 = external hidden unnamed_addr constant [128 x i8], align 1
@.str.134 = external hidden unnamed_addr constant [99 x i8], align 1
@.str.135 = external hidden unnamed_addr constant [171 x i8], align 1
@.str.136 = external hidden unnamed_addr constant [128 x i8], align 1
@.str.137 = external hidden unnamed_addr constant [97 x i8], align 1
@.str.138 = external hidden unnamed_addr constant [52 x i8], align 1
@.str.139 = external hidden unnamed_addr constant [97 x i8], align 1
@.str.140 = external hidden unnamed_addr constant [93 x i8], align 1
@.str.141 = external hidden unnamed_addr constant [125 x i8], align 1
@.str.142 = external hidden unnamed_addr constant [46 x i8], align 1
@.str.143 = external hidden unnamed_addr constant [23 x i8], align 1
@.str.144 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.145 = external hidden unnamed_addr constant [92 x i8], align 1
@.str.146 = external hidden unnamed_addr constant [18 x i8], align 1
@.str.147 = external hidden unnamed_addr constant [33 x i8], align 1
@.str.148 = external hidden unnamed_addr constant [18 x i8], align 1
@.str.149 = external hidden unnamed_addr constant [50 x i8], align 1
@.str.150 = external hidden unnamed_addr constant [155 x i8], align 1
@.str.151 = external hidden unnamed_addr constant [44 x i8], align 1
@.str.152 = external hidden unnamed_addr constant [85 x i8], align 1
@.str.153 = external hidden unnamed_addr constant [12 x i8], align 1
@.str.154 = external hidden unnamed_addr constant [12 x i8], align 1
@.str.155 = external hidden unnamed_addr constant [65 x i8], align 1
@.str.156 = external hidden unnamed_addr constant [72 x i8], align 1

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

declare void @class_protect_sprintf(ptr noundef, ptr noundef, ...) local_unnamed_addr #2

; 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 allockind("alloc,uninitialized") allocsize(0) memory(inaccessiblemem: readwrite)
declare noalias noundef ptr @malloc(i64 noundef) local_unnamed_addr #3

declare i32 @background_tau_of_z(ptr noundef, double noundef, ptr noundef) local_unnamed_addr #2

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

; Function Attrs: nounwind uwtable
define dso_local range(i32 0, 2) i32 @perturb_timesampling_for_sources(ptr nocapture noundef readonly %0, ptr noundef %1, ptr noundef %2, ptr noundef %3) local_unnamed_addr #5 {
  %5 = alloca i32, align 4
  %6 = alloca i32, align 4
  %7 = alloca i32, align 4
  %8 = alloca i32, align 4
  %9 = alloca double, align 8
  %10 = alloca double, align 8
  %11 = alloca [2048 x i8], align 16
  %12 = alloca [2048 x i8], align 16
  %13 = alloca [2048 x i8], align 16
  %14 = alloca [2048 x i8], align 16
  %15 = alloca [2048 x i8], align 16
  %16 = alloca [2048 x i8], align 16
  %17 = alloca [2048 x i8], align 16
  %18 = alloca [2048 x i8], align 16
  %19 = alloca [2048 x i8], align 16
  %20 = alloca [2048 x i8], align 16
  %21 = alloca [2048 x i8], align 16
  %22 = alloca [2048 x i8], align 16
  %23 = alloca [2048 x i8], align 16
  %24 = alloca [2048 x i8], align 16
  %25 = alloca [2048 x i8], align 16
  %26 = alloca [2048 x i8], align 16
  %27 = alloca [2048 x i8], align 16
  %28 = alloca [2048 x i8], align 16
  %29 = alloca [2048 x i8], align 16
  %30 = alloca [2048 x i8], align 16
  %31 = alloca [2048 x i8], align 16
  %32 = alloca [2048 x i8], align 16
  %33 = alloca [2048 x i8], align 16
  %34 = alloca [2048 x i8], align 16
  %35 = alloca [2048 x i8], align 16
  %36 = alloca [2048 x i8], align 16
  %37 = alloca [2048 x i8], align 16
  %38 = alloca [2048 x i8], align 16
  %39 = alloca [2048 x i8], align 16
  %40 = alloca [2048 x i8], align 16
  %41 = alloca [2048 x i8], align 16
  %42 = alloca [2048 x i8], align 16
  %43 = alloca [2048 x i8], align 16
  %44 = alloca [2048 x i8], align 16
  %45 = alloca [2048 x i8], align 16
  %46 = alloca [2048 x i8], align 16
  %47 = alloca [2048 x i8], align 16
  %48 = alloca [2048 x i8], align 16
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %5) #8
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %6) #8
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %7) #8
  call void @llvm.lifetime.start.p0(i64 4, ptr nonnull %8) #8
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %9) #8
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %10) #8
  %49 = getelementptr inbounds i8, ptr %1, i64 616
  %50 = load i32, ptr %49, align 8, !tbaa !5
  %51 = sext i32 %50 to i64
  %52 = shl nsw i64 %51, 3
  %53 = tail call noalias ptr @malloc(i64 noundef %52) #9
  %54 = icmp eq ptr %53, null
  br i1 %54, label %55, label %58

55:                                               ; preds = %4
  %56 = shl i32 %50, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %11) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %11, ptr noundef nonnull @.str.4, ptr noundef nonnull @.str.128, i32 noundef %56) #8
  %57 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %57, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1431, ptr noundef nonnull %11) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %11) #8
  br label %.loopexit

58:                                               ; preds = %4
  %59 = getelementptr inbounds i8, ptr %2, i64 372
  %60 = load i32, ptr %59, align 4, !tbaa !13
  %61 = sext i32 %60 to i64
  %62 = shl nsw i64 %61, 3
  %63 = tail call noalias ptr @malloc(i64 noundef %62) #9
  %64 = icmp eq ptr %63, null
  br i1 %64, label %65, label %68

65:                                               ; preds = %58
  %66 = shl i32 %60, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %12) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %12, ptr noundef nonnull @.str.4, ptr noundef nonnull @.str.129, i32 noundef %66) #8
  %67 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %67, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1432, ptr noundef nonnull %12) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %12) #8
  br label %.loopexit

68:                                               ; preds = %58
  %69 = getelementptr inbounds i8, ptr %3, i64 2020
  %70 = load i16, ptr %69, align 4, !tbaa !15
  %71 = icmp eq i16 %70, 1
  br i1 %71, label %76, label %72

72:                                               ; preds = %68
  %73 = getelementptr inbounds i8, ptr %3, i64 20
  %74 = load i16, ptr %73, align 4, !tbaa !17
  %75 = icmp eq i16 %74, 1
  br i1 %75, label %76, label %219

76:                                               ; preds = %72, %68
  %77 = getelementptr inbounds i8, ptr %2, i64 584
  %78 = load double, ptr %77, align 8, !tbaa !18
  store double %78, ptr %10, align 8, !tbaa !19
  %79 = getelementptr inbounds i8, ptr %1, i64 840
  %80 = load i16, ptr %79, align 8, !tbaa !20
  %81 = getelementptr inbounds i8, ptr %1, i64 846
  %82 = load i16, ptr %81, align 2, !tbaa !21
  %83 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %78, i16 noundef signext %80, i16 noundef signext %82, ptr noundef nonnull %7, ptr noundef nonnull %53) #8
  %84 = icmp eq i32 %83, 1
  br i1 %84, label %85, label %88

85:                                               ; preds = %76
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %13) #8
  %86 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %13, ptr noundef nonnull @.str, ptr noundef nonnull @.str.130, ptr noundef nonnull %86) #8
  %87 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %87, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1459, ptr noundef nonnull %13) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %13) #8
  br label %.loopexit

88:                                               ; preds = %76
  %89 = getelementptr inbounds i8, ptr %1, i64 464
  %90 = load i32, ptr %89, align 8, !tbaa !22
  %91 = sext i32 %90 to i64
  %92 = getelementptr inbounds double, ptr %53, i64 %91
  %93 = load double, ptr %92, align 8, !tbaa !19
  %94 = fdiv double 1.000000e+00, %93
  %95 = fadd double %94, -1.000000e+00
  %96 = getelementptr inbounds i8, ptr %2, i64 600
  %97 = load i16, ptr %96, align 8, !tbaa !23
  %98 = call i32 @thermodynamics_at_z(ptr noundef nonnull %1, ptr noundef nonnull %2, double noundef %95, i16 noundef signext %97, ptr noundef nonnull %8, ptr noundef nonnull %53, ptr noundef nonnull %63) #8
  %99 = icmp eq i32 %98, 1
  br i1 %99, label %100, label %103

100:                                              ; preds = %88
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %14) #8
  %101 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %14, ptr noundef nonnull @.str, ptr noundef nonnull @.str.131, ptr noundef nonnull %101) #8
  %102 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %102, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1469, ptr noundef nonnull %14) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %14) #8
  br label %.loopexit

103:                                              ; preds = %88
  %104 = load i32, ptr %89, align 8, !tbaa !22
  %105 = sext i32 %104 to i64
  %106 = getelementptr inbounds double, ptr %53, i64 %105
  %107 = load double, ptr %106, align 8, !tbaa !19
  %108 = getelementptr inbounds i8, ptr %1, i64 468
  %109 = load i32, ptr %108, align 4, !tbaa !24
  %110 = sext i32 %109 to i64
  %111 = getelementptr inbounds double, ptr %53, i64 %110
  %112 = load double, ptr %111, align 8, !tbaa !19
  %113 = fmul double %107, %112
  %114 = getelementptr inbounds i8, ptr %2, i64 276
  %115 = load i32, ptr %114, align 4, !tbaa !25
  %116 = sext i32 %115 to i64
  %117 = getelementptr inbounds double, ptr %63, i64 %116
  %118 = load double, ptr %117, align 8, !tbaa !19
  %119 = fdiv double %113, %118
  %120 = getelementptr inbounds i8, ptr %0, i64 1520
  %121 = load double, ptr %120, align 8, !tbaa !26
  %122 = fcmp ogt double %119, %121
  br i1 %122, label %123, label %126

123:                                              ; preds = %103
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %15) #8
  %124 = load double, ptr %10, align 8, !tbaa !19
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %15, ptr noundef nonnull @.str.132, double noundef %124) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %16) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %16, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.133, ptr noundef nonnull %15) #8
  %125 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %125, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1477, ptr noundef nonnull %16) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %16) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %15) #8
  br label %.loopexit

126:                                              ; preds = %103
  %127 = getelementptr inbounds i8, ptr %2, i64 416
  %128 = load double, ptr %127, align 8, !tbaa !28
  %129 = load i16, ptr %79, align 8, !tbaa !20
  %130 = load i16, ptr %81, align 2, !tbaa !21
  %131 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %128, i16 noundef signext %129, i16 noundef signext %130, ptr noundef nonnull %7, ptr noundef nonnull %53) #8
  %132 = icmp eq i32 %131, 1
  br i1 %132, label %133, label %136

133:                                              ; preds = %126
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %17) #8
  %134 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %17, ptr noundef nonnull @.str, ptr noundef nonnull @.str.134, ptr noundef nonnull %134) #8
  %135 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %135, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1489, ptr noundef nonnull %17) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %17) #8
  br label %.loopexit

136:                                              ; preds = %126
  %137 = load i32, ptr %89, align 8, !tbaa !22
  %138 = sext i32 %137 to i64
  %139 = getelementptr inbounds double, ptr %53, i64 %138
  %140 = load double, ptr %139, align 8, !tbaa !19
  %141 = fdiv double 1.000000e+00, %140
  %142 = fadd double %141, -1.000000e+00
  %143 = load i16, ptr %96, align 8, !tbaa !23
  %144 = call i32 @thermodynamics_at_z(ptr noundef nonnull %1, ptr noundef nonnull %2, double noundef %142, i16 noundef signext %143, ptr noundef nonnull %8, ptr noundef nonnull %53, ptr noundef nonnull %63) #8
  %145 = icmp eq i32 %144, 1
  br i1 %145, label %146, label %149

146:                                              ; preds = %136
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %18) #8
  %147 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %18, ptr noundef nonnull @.str, ptr noundef nonnull @.str.131, ptr noundef nonnull %147) #8
  %148 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %148, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1499, ptr noundef nonnull %18) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %18) #8
  br label %.loopexit

149:                                              ; preds = %136
  %150 = load i32, ptr %89, align 8, !tbaa !22
  %151 = sext i32 %150 to i64
  %152 = getelementptr inbounds double, ptr %53, i64 %151
  %153 = load double, ptr %152, align 8, !tbaa !19
  %154 = load i32, ptr %108, align 4, !tbaa !24
  %155 = sext i32 %154 to i64
  %156 = getelementptr inbounds double, ptr %53, i64 %155
  %157 = load double, ptr %156, align 8, !tbaa !19
  %158 = fmul double %153, %157
  %159 = load i32, ptr %114, align 4, !tbaa !25
  %160 = sext i32 %159 to i64
  %161 = getelementptr inbounds double, ptr %63, i64 %160
  %162 = load double, ptr %161, align 8, !tbaa !19
  %163 = fdiv double %158, %162
  %164 = load double, ptr %120, align 8, !tbaa !26
  %165 = fcmp olt double %163, %164
  br i1 %165, label %166, label %168

166:                                              ; preds = %149
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %19) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %19, ptr noundef nonnull @.str.135) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %20) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %20, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.136, ptr noundef nonnull %19) #8
  %167 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %167, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1506, ptr noundef nonnull %20) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %20) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %19) #8
  br label %.loopexit

168:                                              ; preds = %149
  %169 = load double, ptr %10, align 8, !tbaa !19
  %170 = fadd double %128, %169
  %171 = getelementptr inbounds i8, ptr %0, i64 1640
  %.0366419 = fmul double %170, 5.000000e-01
  %172 = fsub double %128, %169
  %173 = load double, ptr %171, align 8, !tbaa !29
  %174 = fcmp ogt double %172, %173
  br i1 %174, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %168, %213
  %.0366421 = phi double [ %.0366, %213 ], [ %.0366419, %168 ]
  %.0362420 = phi double [ %.1363, %213 ], [ %128, %168 ]
  %175 = load i16, ptr %79, align 8, !tbaa !20
  %176 = load i16, ptr %81, align 2, !tbaa !21
  %177 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %.0366421, i16 noundef signext %175, i16 noundef signext %176, ptr noundef nonnull %7, ptr noundef nonnull %53) #8
  %178 = icmp eq i32 %177, 1
  br i1 %178, label %179, label %182

179:                                              ; preds = %.lr.ph
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %21) #8
  %180 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %21, ptr noundef nonnull @.str, ptr noundef nonnull @.str.137, ptr noundef nonnull %180) #8
  %181 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %181, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1519, ptr noundef nonnull %21) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %21) #8
  br label %.loopexit

182:                                              ; preds = %.lr.ph
  %183 = load i32, ptr %89, align 8, !tbaa !22
  %184 = sext i32 %183 to i64
  %185 = getelementptr inbounds double, ptr %53, i64 %184
  %186 = load double, ptr %185, align 8, !tbaa !19
  %187 = fdiv double 1.000000e+00, %186
  %188 = fadd double %187, -1.000000e+00
  %189 = load i16, ptr %96, align 8, !tbaa !23
  %190 = call i32 @thermodynamics_at_z(ptr noundef nonnull %1, ptr noundef nonnull %2, double noundef %188, i16 noundef signext %189, ptr noundef nonnull %8, ptr noundef nonnull %53, ptr noundef nonnull %63) #8
  %191 = icmp eq i32 %190, 1
  br i1 %191, label %192, label %195

192:                                              ; preds = %182
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %22) #8
  %193 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %22, ptr noundef nonnull @.str, ptr noundef nonnull @.str.131, ptr noundef nonnull %193) #8
  %194 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %194, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1529, ptr noundef nonnull %22) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %22) #8
  br label %.loopexit

195:                                              ; preds = %182
  %196 = load i32, ptr %89, align 8, !tbaa !22
  %197 = sext i32 %196 to i64
  %198 = getelementptr inbounds double, ptr %53, i64 %197
  %199 = load double, ptr %198, align 8, !tbaa !19
  %200 = load i32, ptr %108, align 4, !tbaa !24
  %201 = sext i32 %200 to i64
  %202 = getelementptr inbounds double, ptr %53, i64 %201
  %203 = load double, ptr %202, align 8, !tbaa !19
  %204 = fmul double %199, %203
  %205 = load i32, ptr %114, align 4, !tbaa !25
  %206 = sext i32 %205 to i64
  %207 = getelementptr inbounds double, ptr %63, i64 %206
  %208 = load double, ptr %207, align 8, !tbaa !19
  %209 = fdiv double %204, %208
  %210 = load double, ptr %120, align 8, !tbaa !26
  %211 = fcmp ogt double %209, %210
  br i1 %211, label %._crit_edge468, label %212

._crit_edge468:                                   ; preds = %195
  %.pre = load double, ptr %10, align 8, !tbaa !19
  br label %213

212:                                              ; preds = %195
  store double %.0366421, ptr %10, align 8, !tbaa !19
  br label %213

213:                                              ; preds = %._crit_edge468, %212
  %214 = phi double [ %.0366421, %212 ], [ %.pre, %._crit_edge468 ]
  %.1363 = phi double [ %.0362420, %212 ], [ %.0366421, %._crit_edge468 ]
  %215 = fadd double %.1363, %214
  %.0366 = fmul double %215, 5.000000e-01
  %216 = fsub double %.1363, %214
  %217 = load double, ptr %171, align 8, !tbaa !29
  %218 = fcmp ogt double %216, %217
  br i1 %218, label %.lr.ph, label %._crit_edge, !llvm.loop !30

._crit_edge:                                      ; preds = %213, %168
  %.0366.lcssa = phi double [ %.0366419, %168 ], [ %.0366, %213 ]
  store double %.0366.lcssa, ptr %9, align 8, !tbaa !19
  br label %254

219:                                              ; preds = %72
  %220 = getelementptr inbounds i8, ptr %3, i64 1992
  %221 = load double, ptr %220, align 8, !tbaa !32
  %222 = fadd double %221, 1.000000e+00
  %223 = call i32 @background_tau_of_z(ptr noundef nonnull %1, double noundef %222, ptr noundef nonnull %9) #8
  %224 = icmp eq i32 %223, 1
  br i1 %224, label %225, label %228

225:                                              ; preds = %219
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %23) #8
  %226 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %23, ptr noundef nonnull @.str, ptr noundef nonnull @.str.138, ptr noundef nonnull %226) #8
  %227 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %227, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1555, ptr noundef nonnull %23) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %23) #8
  br label %.loopexit

228:                                              ; preds = %219
  %229 = load double, ptr %9, align 8, !tbaa !19
  %230 = getelementptr inbounds i8, ptr %1, i64 840
  %231 = load i16, ptr %230, align 8, !tbaa !20
  %232 = getelementptr inbounds i8, ptr %1, i64 846
  %233 = load i16, ptr %232, align 2, !tbaa !21
  %234 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %229, i16 noundef signext %231, i16 noundef signext %233, ptr noundef nonnull %7, ptr noundef nonnull %53) #8
  %235 = icmp eq i32 %234, 1
  br i1 %235, label %236, label %239

236:                                              ; preds = %228
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %24) #8
  %237 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %24, ptr noundef nonnull @.str, ptr noundef nonnull @.str.139, ptr noundef nonnull %237) #8
  %238 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %238, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1568, ptr noundef nonnull %24) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %24) #8
  br label %.loopexit

239:                                              ; preds = %228
  %240 = getelementptr inbounds i8, ptr %1, i64 464
  %241 = load i32, ptr %240, align 8, !tbaa !22
  %242 = sext i32 %241 to i64
  %243 = getelementptr inbounds double, ptr %53, i64 %242
  %244 = load double, ptr %243, align 8, !tbaa !19
  %245 = fdiv double 1.000000e+00, %244
  %246 = fadd double %245, -1.000000e+00
  %247 = getelementptr inbounds i8, ptr %2, i64 600
  %248 = load i16, ptr %247, align 8, !tbaa !23
  %249 = call i32 @thermodynamics_at_z(ptr noundef nonnull %1, ptr noundef nonnull %2, double noundef %246, i16 noundef signext %248, ptr noundef nonnull %8, ptr noundef nonnull %53, ptr noundef nonnull %63) #8
  %250 = icmp eq i32 %249, 1
  br i1 %250, label %251, label %._crit_edge469

._crit_edge469:                                   ; preds = %239
  %.pre470 = load double, ptr %9, align 8, !tbaa !19
  br label %254

251:                                              ; preds = %239
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %25) #8
  %252 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %25, ptr noundef nonnull @.str, ptr noundef nonnull @.str.131, ptr noundef nonnull %252) #8
  %253 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %253, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1578, ptr noundef nonnull %25) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %25) #8
  br label %.loopexit

254:                                              ; preds = %._crit_edge469, %._crit_edge
  %255 = phi double [ %.pre470, %._crit_edge469 ], [ %.0366.lcssa, %._crit_edge ]
  %256 = load i32, ptr %7, align 4, !tbaa !33
  store i32 %256, ptr %5, align 4, !tbaa !33
  %257 = load i32, ptr %8, align 4, !tbaa !33
  store i32 %257, ptr %6, align 4, !tbaa !33
  %258 = getelementptr inbounds i8, ptr %1, i64 344
  %259 = load double, ptr %258, align 8, !tbaa !34
  %260 = fcmp olt double %255, %259
  br i1 %260, label %.lr.ph425, label %._crit_edge426

.lr.ph425:                                        ; preds = %254
  %261 = getelementptr inbounds i8, ptr %1, i64 840
  %262 = getelementptr inbounds i8, ptr %1, i64 848
  %263 = getelementptr inbounds i8, ptr %1, i64 464
  %264 = getelementptr inbounds i8, ptr %2, i64 602
  %265 = getelementptr inbounds i8, ptr %1, i64 468
  %266 = getelementptr inbounds i8, ptr %2, i64 364
  %267 = getelementptr inbounds i8, ptr %1, i64 472
  %268 = getelementptr inbounds i8, ptr %0, i64 1616
  %269 = getelementptr inbounds i8, ptr %0, i64 2344
  br label %270

270:                                              ; preds = %.lr.ph425, %344
  %.0354423 = phi i32 [ 1, %.lr.ph425 ], [ %346, %344 ]
  %.0360422 = phi double [ %255, %.lr.ph425 ], [ %345, %344 ]
  %271 = load i16, ptr %261, align 8, !tbaa !20
  %272 = load i16, ptr %262, align 8, !tbaa !35
  %273 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %.0360422, i16 noundef signext %271, i16 noundef signext %272, ptr noundef nonnull %5, ptr noundef nonnull %53) #8
  %274 = icmp eq i32 %273, 1
  br i1 %274, label %275, label %278

275:                                              ; preds = %270
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %26) #8
  %276 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %26, ptr noundef nonnull @.str, ptr noundef nonnull @.str.140, ptr noundef nonnull %276) #8
  %277 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %277, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1604, ptr noundef nonnull %26) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %26) #8
  br label %.loopexit

278:                                              ; preds = %270
  %279 = load i32, ptr %263, align 8, !tbaa !22
  %280 = sext i32 %279 to i64
  %281 = getelementptr inbounds double, ptr %53, i64 %280
  %282 = load double, ptr %281, align 8, !tbaa !19
  %283 = fdiv double 1.000000e+00, %282
  %284 = fadd double %283, -1.000000e+00
  %285 = load i16, ptr %264, align 2, !tbaa !36
  %286 = call i32 @thermodynamics_at_z(ptr noundef nonnull %1, ptr noundef nonnull %2, double noundef %284, i16 noundef signext %285, ptr noundef nonnull %6, ptr noundef nonnull %53, ptr noundef nonnull %63) #8
  %287 = icmp eq i32 %286, 1
  br i1 %287, label %288, label %291

288:                                              ; preds = %278
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %27) #8
  %289 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %27, ptr noundef nonnull @.str, ptr noundef nonnull @.str.141, ptr noundef nonnull %289) #8
  %290 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %290, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1614, ptr noundef nonnull %27) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %27) #8
  br label %.loopexit

291:                                              ; preds = %278
  %292 = load i16, ptr %69, align 4, !tbaa !15
  %293 = icmp eq i16 %292, 1
  br i1 %293, label %294, label %320

294:                                              ; preds = %291
  %295 = load i32, ptr %266, align 4, !tbaa !37
  %296 = sext i32 %295 to i64
  %297 = getelementptr inbounds double, ptr %63, i64 %296
  %298 = load double, ptr %297, align 8, !tbaa !19
  %299 = load i32, ptr %265, align 4, !tbaa !24
  %300 = sext i32 %299 to i64
  %301 = getelementptr inbounds double, ptr %53, i64 %300
  %302 = load double, ptr %301, align 8, !tbaa !19
  %303 = load i32, ptr %263, align 8, !tbaa !22
  %304 = sext i32 %303 to i64
  %305 = getelementptr inbounds double, ptr %53, i64 %304
  %306 = load double, ptr %305, align 8, !tbaa !19
  %307 = fmul double %302, %306
  %308 = load i32, ptr %267, align 8, !tbaa !38
  %309 = sext i32 %308 to i64
  %310 = getelementptr inbounds double, ptr %53, i64 %309
  %311 = load double, ptr %310, align 8, !tbaa !19
  %312 = fmul double %307, 2.000000e+00
  %313 = fmul double %307, %312
  %314 = call double @llvm.fmuladd.f64(double %311, double %306, double %313)
  %315 = fneg double %307
  %316 = fmul double %307, %315
  %317 = call double @llvm.fmuladd.f64(double %314, double 2.000000e+00, double %316)
  %318 = call double @llvm.fabs.f64(double %317)
  %319 = call double @llvm.fmuladd.f64(double %298, double %298, double %318)
  %sqrt = call double @llvm.sqrt.f64(double %319)
  br label %330

320:                                              ; preds = %291
  %321 = load i32, ptr %265, align 4, !tbaa !24
  %322 = sext i32 %321 to i64
  %323 = getelementptr inbounds double, ptr %53, i64 %322
  %324 = load double, ptr %323, align 8, !tbaa !19
  %325 = load i32, ptr %263, align 8, !tbaa !22
  %326 = sext i32 %325 to i64
  %327 = getelementptr inbounds double, ptr %53, i64 %326
  %328 = load double, ptr %327, align 8, !tbaa !19
  %329 = fmul double %324, %328
  br label %330

330:                                              ; preds = %320, %294
  %.0364 = phi double [ %sqrt, %294 ], [ %329, %320 ]
  %331 = fcmp oeq double %.0364, 0.000000e+00
  br i1 %331, label %332, label %334

332:                                              ; preds = %330
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %28) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %28, ptr noundef nonnull @.str.142) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %29) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %29, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.143, ptr noundef nonnull %28) #8
  %333 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %333, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1640, ptr noundef nonnull %29) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %29) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %28) #8
  br label %.loopexit

334:                                              ; preds = %330
  %335 = fdiv double 1.000000e+00, %.0364
  %336 = load double, ptr %268, align 8, !tbaa !39
  %337 = fmul double %335, %336
  %338 = fdiv double %337, %.0360422
  %339 = call double @llvm.fabs.f64(double %338)
  %340 = load double, ptr %269, align 8, !tbaa !40
  %341 = fcmp olt double %339, %340
  br i1 %341, label %342, label %344

342:                                              ; preds = %334
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %30) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %30, ptr noundef nonnull @.str.144, double noundef %337) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %31) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %31, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.145, ptr noundef nonnull %30) #8
  %343 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %343, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1647, ptr noundef nonnull %31) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %31) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %30) #8
  br label %.loopexit

344:                                              ; preds = %334
  %345 = call double @llvm.fmuladd.f64(double %336, double %335, double %.0360422)
  %346 = add nuw nsw i32 %.0354423, 1
  %347 = load double, ptr %258, align 8, !tbaa !34
  %348 = fcmp olt double %345, %347
  br i1 %348, label %270, label %._crit_edge426, !llvm.loop !41

._crit_edge426:                                   ; preds = %344, %254
  %349 = phi double [ %259, %254 ], [ %347, %344 ]
  %.0354.lcssa = phi i32 [ 1, %254 ], [ %346, %344 ]
  %350 = getelementptr inbounds i8, ptr %3, i64 2400
  store i32 %.0354.lcssa, ptr %350, align 8, !tbaa !42
  %351 = zext nneg i32 %.0354.lcssa to i64
  %352 = shl nuw nsw i64 %351, 3
  %353 = call noalias ptr @malloc(i64 noundef %352) #9
  %354 = getelementptr inbounds i8, ptr %3, i64 2392
  store ptr %353, ptr %354, align 8, !tbaa !43
  %355 = icmp eq ptr %353, null
  br i1 %355, label %356, label %359

356:                                              ; preds = %._crit_edge426
  %357 = shl i32 %.0354.lcssa, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %32) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %32, ptr noundef nonnull @.str.4, ptr noundef nonnull @.str.146, i32 noundef %357) #8
  %358 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %358, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1658, ptr noundef nonnull %32) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %32) #8
  br label %.loopexit

359:                                              ; preds = %._crit_edge426
  %360 = load double, ptr %9, align 8, !tbaa !19
  store double %360, ptr %353, align 8, !tbaa !19
  %361 = load i32, ptr %7, align 4, !tbaa !33
  store i32 %361, ptr %5, align 4, !tbaa !33
  %362 = load i32, ptr %8, align 4, !tbaa !33
  store i32 %362, ptr %6, align 4, !tbaa !33
  %363 = fcmp olt double %360, %349
  br i1 %363, label %.lr.ph431, label %._crit_edge432

.lr.ph431:                                        ; preds = %359
  %364 = getelementptr inbounds i8, ptr %1, i64 840
  %365 = getelementptr inbounds i8, ptr %1, i64 848
  %366 = getelementptr inbounds i8, ptr %1, i64 464
  %367 = getelementptr inbounds i8, ptr %2, i64 602
  %368 = getelementptr inbounds i8, ptr %1, i64 468
  %369 = getelementptr inbounds i8, ptr %2, i64 364
  %370 = getelementptr inbounds i8, ptr %1, i64 472
  %371 = getelementptr inbounds i8, ptr %0, i64 1616
  %372 = getelementptr inbounds i8, ptr %0, i64 2344
  br label %373

373:                                              ; preds = %.lr.ph431, %447
  %indvars.iv = phi i64 [ 0, %.lr.ph431 ], [ %indvars.iv.next, %447 ]
  %.1361428 = phi double [ %360, %.lr.ph431 ], [ %448, %447 ]
  %374 = load i16, ptr %364, align 8, !tbaa !20
  %375 = load i16, ptr %365, align 8, !tbaa !35
  %376 = call i32 @background_at_tau(ptr noundef nonnull %1, double noundef %.1361428, i16 noundef signext %374, i16 noundef signext %375, ptr noundef nonnull %5, ptr noundef nonnull %53) #8
  %377 = icmp eq i32 %376, 1
  br i1 %377, label %378, label %381

378:                                              ; preds = %373
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %33) #8
  %379 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %33, ptr noundef nonnull @.str, ptr noundef nonnull @.str.140, ptr noundef nonnull %379) #8
  %380 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %380, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1687, ptr noundef nonnull %33) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %33) #8
  br label %.loopexit

381:                                              ; preds = %373
  %382 = load i32, ptr %366, align 8, !tbaa !22
  %383 = sext i32 %382 to i64
  %384 = getelementptr inbounds double, ptr %53, i64 %383
  %385 = load double, ptr %384, align 8, !tbaa !19
  %386 = fdiv double 1.000000e+00, %385
  %387 = fadd double %386, -1.000000e+00
  %388 = load i16, ptr %367, align 2, !tbaa !36
  %389 = call i32 @thermodynamics_at_z(ptr noundef nonnull %1, ptr noundef nonnull %2, double noundef %387, i16 noundef signext %388, ptr noundef nonnull %6, ptr noundef nonnull %53, ptr noundef nonnull %63) #8
  %390 = icmp eq i32 %389, 1
  br i1 %390, label %391, label %394

391:                                              ; preds = %381
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %34) #8
  %392 = getelementptr inbounds i8, ptr %2, i64 606
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %34, ptr noundef nonnull @.str, ptr noundef nonnull @.str.141, ptr noundef nonnull %392) #8
  %393 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %393, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1697, ptr noundef nonnull %34) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %34) #8
  br label %.loopexit

394:                                              ; preds = %381
  %395 = load i16, ptr %69, align 4, !tbaa !15
  %396 = icmp eq i16 %395, 1
  br i1 %396, label %397, label %423

397:                                              ; preds = %394
  %398 = load i32, ptr %369, align 4, !tbaa !37
  %399 = sext i32 %398 to i64
  %400 = getelementptr inbounds double, ptr %63, i64 %399
  %401 = load double, ptr %400, align 8, !tbaa !19
  %402 = load i32, ptr %368, align 4, !tbaa !24
  %403 = sext i32 %402 to i64
  %404 = getelementptr inbounds double, ptr %53, i64 %403
  %405 = load double, ptr %404, align 8, !tbaa !19
  %406 = load i32, ptr %366, align 8, !tbaa !22
  %407 = sext i32 %406 to i64
  %408 = getelementptr inbounds double, ptr %53, i64 %407
  %409 = load double, ptr %408, align 8, !tbaa !19
  %410 = fmul double %405, %409
  %411 = load i32, ptr %370, align 8, !tbaa !38
  %412 = sext i32 %411 to i64
  %413 = getelementptr inbounds double, ptr %53, i64 %412
  %414 = load double, ptr %413, align 8, !tbaa !19
  %415 = fmul double %410, 2.000000e+00
  %416 = fmul double %410, %415
  %417 = call double @llvm.fmuladd.f64(double %414, double %409, double %416)
  %418 = fneg double %410
  %419 = fmul double %410, %418
  %420 = call double @llvm.fmuladd.f64(double %417, double 2.000000e+00, double %419)
  %421 = call double @llvm.fabs.f64(double %420)
  %422 = call double @llvm.fmuladd.f64(double %401, double %401, double %421)
  %sqrt385 = call double @llvm.sqrt.f64(double %422)
  br label %433

423:                                              ; preds = %394
  %424 = load i32, ptr %368, align 4, !tbaa !24
  %425 = sext i32 %424 to i64
  %426 = getelementptr inbounds double, ptr %53, i64 %425
  %427 = load double, ptr %426, align 8, !tbaa !19
  %428 = load i32, ptr %366, align 8, !tbaa !22
  %429 = sext i32 %428 to i64
  %430 = getelementptr inbounds double, ptr %53, i64 %429
  %431 = load double, ptr %430, align 8, !tbaa !19
  %432 = fmul double %427, %431
  br label %433

433:                                              ; preds = %423, %397
  %.1365 = phi double [ %sqrt385, %397 ], [ %432, %423 ]
  %434 = fcmp oeq double %.1365, 0.000000e+00
  br i1 %434, label %435, label %437

435:                                              ; preds = %433
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %35) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %35, ptr noundef nonnull @.str.142) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %36) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %36, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.143, ptr noundef nonnull %35) #8
  %436 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %436, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1721, ptr noundef nonnull %36) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %36) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %35) #8
  br label %.loopexit

437:                                              ; preds = %433
  %438 = fdiv double 1.000000e+00, %.1365
  %439 = load double, ptr %371, align 8, !tbaa !39
  %440 = fmul double %438, %439
  %441 = fdiv double %440, %.1361428
  %442 = call double @llvm.fabs.f64(double %441)
  %443 = load double, ptr %372, align 8, !tbaa !40
  %444 = fcmp olt double %442, %443
  br i1 %444, label %445, label %447

445:                                              ; preds = %437
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %37) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %37, ptr noundef nonnull @.str.144, double noundef %440) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %38) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %38, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.145, ptr noundef nonnull %37) #8
  %446 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %446, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1728, ptr noundef nonnull %38) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %38) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %37) #8
  br label %.loopexit

447:                                              ; preds = %437
  %448 = call double @llvm.fmuladd.f64(double %439, double %438, double %.1361428)
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %449 = load ptr, ptr %354, align 8, !tbaa !43
  %450 = getelementptr inbounds double, ptr %449, i64 %indvars.iv.next
  store double %448, ptr %450, align 8, !tbaa !19
  %451 = load double, ptr %258, align 8, !tbaa !34
  %452 = fcmp olt double %448, %451
  br i1 %452, label %373, label %._crit_edge432.loopexit, !llvm.loop !44

._crit_edge432.loopexit:                          ; preds = %447
  %453 = and i64 %indvars.iv.next, 4294967295
  br label %._crit_edge432

._crit_edge432:                                   ; preds = %._crit_edge432.loopexit, %359
  %454 = phi ptr [ %353, %359 ], [ %449, %._crit_edge432.loopexit ]
  %.1.lcssa = phi i64 [ 0, %359 ], [ %453, %._crit_edge432.loopexit ]
  %.lcssa397 = phi double [ %349, %359 ], [ %451, %._crit_edge432.loopexit ]
  %455 = getelementptr inbounds double, ptr %454, i64 %.1.lcssa
  store double %.lcssa397, ptr %455, align 8, !tbaa !19
  call void @free(ptr noundef nonnull %53) #8
  call void @free(ptr noundef %63) #8
  %456 = getelementptr inbounds i8, ptr %3, i64 1992
  %457 = load double, ptr %456, align 8, !tbaa !32
  %458 = fcmp olt double %457, 0.000000e+00
  br i1 %458, label %459, label %461

459:                                              ; preds = %._crit_edge432
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %39) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %39, ptr noundef nonnull @.str.147, double noundef %457) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %40) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %40, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.148, ptr noundef nonnull %39) #8
  %460 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %460, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1752, ptr noundef nonnull %40) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %40) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %39) #8
  br label %.loopexit

461:                                              ; preds = %._crit_edge432
  %462 = fcmp oeq double %457, 0.000000e+00
  br i1 %462, label %463, label %465

463:                                              ; preds = %461
  %464 = getelementptr inbounds i8, ptr %3, i64 2480
  store i32 1, ptr %464, align 8, !tbaa !45
  br label %.loopexit388

465:                                              ; preds = %461
  %466 = call i32 @background_tau_of_z(ptr noundef nonnull %1, double noundef %457, ptr noundef nonnull %10) #8
  %467 = icmp eq i32 %466, 1
  br i1 %467, label %468, label %471

468:                                              ; preds = %465
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %41) #8
  %469 = getelementptr inbounds i8, ptr %1, i64 2902
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %41, ptr noundef nonnull @.str, ptr noundef nonnull @.str.149, ptr noundef nonnull %469) #8
  %470 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %470, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1765, ptr noundef nonnull %41) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %41) #8
  br label %.loopexit

471:                                              ; preds = %465
  %472 = load double, ptr %10, align 8, !tbaa !19
  %473 = load ptr, ptr %354, align 8, !tbaa !43
  %474 = load double, ptr %473, align 8, !tbaa !19
  %475 = fcmp ugt double %472, %474
  br i1 %475, label %.preheader389, label %476

476:                                              ; preds = %471
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %42) #8
  %477 = load double, ptr %456, align 8, !tbaa !32
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %42, ptr noundef nonnull @.str.150, double noundef %477, double noundef %472, double noundef %474) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %43) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %43, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.151, ptr noundef nonnull %42) #8
  %478 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %478, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1770, ptr noundef nonnull %43) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %43) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %42) #8
  br label %.loopexit

.preheader389:                                    ; preds = %471, %.preheader389
  %indvars.iv459 = phi i64 [ %indvars.iv.next460, %.preheader389 ], [ 0, %471 ]
  %479 = getelementptr inbounds double, ptr %473, i64 %indvars.iv459
  %480 = load double, ptr %479, align 8, !tbaa !19
  %481 = fcmp olt double %480, %472
  %indvars.iv.next460 = add nuw nsw i64 %indvars.iv459, 1
  br i1 %481, label %.preheader389, label %482, !llvm.loop !46

482:                                              ; preds = %.preheader389
  %483 = trunc nuw nsw i64 %indvars.iv459 to i32
  switch i32 %483, label %486 [
    i32 0, label %484
    i32 1, label %.thread383
    i32 3, label %.thread383
    i32 2, label %.thread383
  ]

484:                                              ; preds = %482
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %44) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %44, ptr noundef nonnull @.str.152) #8
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %45) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %45, ptr noundef nonnull @.str.7, ptr noundef nonnull @.str.153, ptr noundef nonnull %44) #8
  %485 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %485, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1778, ptr noundef nonnull %45) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %45) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %44) #8
  br label %.loopexit

486:                                              ; preds = %482
  %487 = icmp ugt i64 %indvars.iv459, 4
  %488 = add nsw i32 %483, -5
  %spec.select = select i1 %487, i32 %488, i32 0
  br label %.thread383

.thread383:                                       ; preds = %482, %482, %482, %486
  %489 = phi i32 [ 0, %482 ], [ %spec.select, %486 ], [ 0, %482 ], [ 0, %482 ]
  %490 = load i32, ptr %350, align 8, !tbaa !42
  %491 = sub nsw i32 %490, %489
  %492 = getelementptr inbounds i8, ptr %3, i64 2480
  store i32 %491, ptr %492, align 8, !tbaa !45
  %493 = sext i32 %491 to i64
  %494 = shl nsw i64 %493, 3
  %495 = call noalias ptr @malloc(i64 noundef %494) #9
  %496 = getelementptr inbounds i8, ptr %3, i64 2472
  store ptr %495, ptr %496, align 8, !tbaa !47
  %497 = icmp eq ptr %495, null
  br i1 %497, label %499, label %.preheader387

.preheader387:                                    ; preds = %.thread383
  %498 = icmp sgt i32 %491, 0
  br i1 %498, label %.lr.ph436, label %.loopexit388

499:                                              ; preds = %.thread383
  %500 = shl i32 %491, 3
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %46) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %46, ptr noundef nonnull @.str.4, ptr noundef nonnull @.str.154, i32 noundef %500) #8
  %501 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %501, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1788, ptr noundef nonnull %46) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %46) #8
  br label %.loopexit

.lr.ph436:                                        ; preds = %.preheader387, %.lr.ph436
  %indvars.iv462 = phi i64 [ %indvars.iv.next463, %.lr.ph436 ], [ 0, %.preheader387 ]
  %502 = phi i32 [ %514, %.lr.ph436 ], [ %491, %.preheader387 ]
  %503 = load ptr, ptr %354, align 8, !tbaa !43
  %504 = trunc nuw nsw i64 %indvars.iv462 to i32
  %505 = sub i32 %504, %502
  %506 = load i32, ptr %350, align 8, !tbaa !42
  %507 = add nsw i32 %505, %506
  %508 = sext i32 %507 to i64
  %509 = getelementptr inbounds double, ptr %503, i64 %508
  %510 = load double, ptr %509, align 8, !tbaa !19
  %511 = call double @log(double noundef %510) #8
  %512 = load ptr, ptr %496, align 8, !tbaa !47
  %513 = getelementptr inbounds double, ptr %512, i64 %indvars.iv462
  store double %511, ptr %513, align 8, !tbaa !19
  %indvars.iv.next463 = add nuw nsw i64 %indvars.iv462, 1
  %514 = load i32, ptr %492, align 8, !tbaa !45
  %515 = sext i32 %514 to i64
  %516 = icmp slt i64 %indvars.iv.next463, %515
  br i1 %516, label %.lr.ph436, label %.loopexit388, !llvm.loop !48

.loopexit388:                                     ; preds = %.lr.ph436, %.preheader387, %463
  %517 = phi i32 [ %491, %.preheader387 ], [ 1, %463 ], [ %514, %.lr.ph436 ]
  %518 = getelementptr inbounds i8, ptr %3, i64 2040
  %519 = load i32, ptr %518, align 8, !tbaa !49
  %520 = icmp sgt i32 %519, 0
  br i1 %520, label %.preheader386.lr.ph, label %.loopexit

.preheader386.lr.ph:                              ; preds = %.loopexit388
  %521 = getelementptr inbounds i8, ptr %3, i64 2072
  %522 = getelementptr inbounds i8, ptr %3, i64 2336
  %523 = getelementptr inbounds i8, ptr %3, i64 2360
  %524 = getelementptr inbounds i8, ptr %3, i64 2464
  %525 = getelementptr inbounds i8, ptr %3, i64 2488
  %526 = getelementptr inbounds i8, ptr %3, i64 2496
  %wide.trip.count = zext nneg i32 %519 to i64
  %.pre471 = load ptr, ptr %521, align 8, !tbaa !50
  %527 = icmp sgt i32 %517, 1
  br label %.preheader386

.preheader386:                                    ; preds = %.preheader386.lr.ph, %._crit_edge441
  %528 = phi ptr [ %.pre471, %.preheader386.lr.ph ], [ %638, %._crit_edge441 ]
  %529 = phi ptr [ %.pre471, %.preheader386.lr.ph ], [ %639, %._crit_edge441 ]
  %indvars.iv465 = phi i64 [ 0, %.preheader386.lr.ph ], [ %indvars.iv.next466, %._crit_edge441 ]
  %530 = getelementptr inbounds i32, ptr %529, i64 %indvars.iv465
  %531 = load i32, ptr %530, align 4, !tbaa !33
  %532 = icmp sgt i32 %531, 0
  br i1 %532, label %.preheader.lr.ph, label %._crit_edge441

.preheader.lr.ph:                                 ; preds = %.preheader386
  %533 = load ptr, ptr %522, align 8, !tbaa !51
  %534 = getelementptr inbounds i32, ptr %533, i64 %indvars.iv465
  %535 = load i32, ptr %534, align 4, !tbaa !33
  %536 = icmp sgt i32 %535, 0
  br i1 %536, label %.preheader, label %._crit_edge441

.preheader:                                       ; preds = %.preheader.lr.ph, %._crit_edge439
  %537 = phi ptr [ %632, %._crit_edge439 ], [ %528, %.preheader.lr.ph ]
  %538 = phi ptr [ %633, %._crit_edge439 ], [ %533, %.preheader.lr.ph ]
  %.0357440 = phi i32 [ %634, %._crit_edge439 ], [ 0, %.preheader.lr.ph ]
  %539 = getelementptr inbounds i32, ptr %538, i64 %indvars.iv465
  %540 = load i32, ptr %539, align 4, !tbaa !33
  %541 = icmp sgt i32 %540, 0
  br i1 %541, label %.lr.ph438, label %._crit_edge439

.lr.ph438:                                        ; preds = %.preheader
  %542 = load i32, ptr %350, align 8, !tbaa !42
  %543 = sub nsw i32 %542, %517
  br label %544

544:                                              ; preds = %.lr.ph438, %627
  %545 = phi i32 [ %540, %.lr.ph438 ], [ %628, %627 ]
  %.0356437 = phi i32 [ 0, %.lr.ph438 ], [ %630, %627 ]
  %546 = load ptr, ptr %523, align 8, !tbaa !52
  %547 = getelementptr inbounds i32, ptr %546, i64 %indvars.iv465
  %548 = load i32, ptr %547, align 4, !tbaa !33
  %549 = mul nsw i32 %542, %548
  %550 = sext i32 %549 to i64
  %551 = shl nsw i64 %550, 3
  %552 = call noalias ptr @malloc(i64 noundef %551) #9
  %553 = load ptr, ptr %524, align 8, !tbaa !53
  %554 = getelementptr inbounds ptr, ptr %553, i64 %indvars.iv465
  %555 = load ptr, ptr %554, align 8, !tbaa !54
  %556 = mul nsw i32 %545, %.0357440
  %557 = add nsw i32 %556, %.0356437
  %558 = sext i32 %557 to i64
  %559 = getelementptr inbounds ptr, ptr %555, i64 %558
  store ptr %552, ptr %559, align 8, !tbaa !54
  %560 = load ptr, ptr %524, align 8, !tbaa !53
  %561 = getelementptr inbounds ptr, ptr %560, i64 %indvars.iv465
  %562 = load ptr, ptr %561, align 8, !tbaa !54
  %563 = load ptr, ptr %522, align 8, !tbaa !51
  %564 = getelementptr inbounds i32, ptr %563, i64 %indvars.iv465
  %565 = load i32, ptr %564, align 4, !tbaa !33
  %566 = mul nsw i32 %565, %.0357440
  %567 = add nsw i32 %566, %.0356437
  %568 = sext i32 %567 to i64
  %569 = getelementptr inbounds ptr, ptr %562, i64 %568
  %570 = load ptr, ptr %569, align 8, !tbaa !54
  %571 = icmp eq ptr %570, null
  br i1 %571, label %572, label %579

572:                                              ; preds = %544
  %573 = load ptr, ptr %523, align 8, !tbaa !52
  %574 = getelementptr inbounds i32, ptr %573, i64 %indvars.iv465
  %575 = load i32, ptr %574, align 4, !tbaa !33
  %576 = shl i32 %542, 3
  %577 = mul i32 %576, %575
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %47) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %47, ptr noundef nonnull @.str.4, ptr noundef nonnull @.str.155, i32 noundef %577) #8
  %578 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %578, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1804, ptr noundef nonnull %47) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %47) #8
  br label %.loopexit

579:                                              ; preds = %544
  br i1 %527, label %580, label %627

580:                                              ; preds = %579
  %581 = load ptr, ptr %523, align 8, !tbaa !52
  %582 = getelementptr inbounds i32, ptr %581, i64 %indvars.iv465
  %583 = load i32, ptr %582, align 4, !tbaa !33
  %584 = mul nsw i32 %583, %543
  %585 = sext i32 %584 to i64
  %586 = getelementptr inbounds double, ptr %570, i64 %585
  %587 = load ptr, ptr %525, align 8, !tbaa !55
  %588 = getelementptr inbounds ptr, ptr %587, i64 %indvars.iv465
  %589 = load ptr, ptr %588, align 8, !tbaa !54
  %590 = getelementptr inbounds ptr, ptr %589, i64 %568
  store ptr %586, ptr %590, align 8, !tbaa !54
  %591 = load ptr, ptr %523, align 8, !tbaa !52
  %592 = getelementptr inbounds i32, ptr %591, i64 %indvars.iv465
  %593 = load i32, ptr %592, align 4, !tbaa !33
  %594 = mul nsw i32 %593, %517
  %595 = sext i32 %594 to i64
  %596 = shl nsw i64 %595, 3
  %597 = call noalias ptr @malloc(i64 noundef %596) #9
  %598 = load ptr, ptr %526, align 8, !tbaa !56
  %599 = getelementptr inbounds ptr, ptr %598, i64 %indvars.iv465
  %600 = load ptr, ptr %599, align 8, !tbaa !54
  %601 = load ptr, ptr %522, align 8, !tbaa !51
  %602 = getelementptr inbounds i32, ptr %601, i64 %indvars.iv465
  %603 = load i32, ptr %602, align 4, !tbaa !33
  %604 = mul nsw i32 %603, %.0357440
  %605 = add nsw i32 %604, %.0356437
  %606 = sext i32 %605 to i64
  %607 = getelementptr inbounds ptr, ptr %600, i64 %606
  store ptr %597, ptr %607, align 8, !tbaa !54
  %608 = load ptr, ptr %526, align 8, !tbaa !56
  %609 = getelementptr inbounds ptr, ptr %608, i64 %indvars.iv465
  %610 = load ptr, ptr %609, align 8, !tbaa !54
  %611 = load ptr, ptr %522, align 8, !tbaa !51
  %612 = getelementptr inbounds i32, ptr %611, i64 %indvars.iv465
  %613 = load i32, ptr %612, align 4, !tbaa !33
  %614 = mul nsw i32 %613, %.0357440
  %615 = add nsw i32 %614, %.0356437
  %616 = sext i32 %615 to i64
  %617 = getelementptr inbounds ptr, ptr %610, i64 %616
  %618 = load ptr, ptr %617, align 8, !tbaa !54
  %619 = icmp eq ptr %618, null
  br i1 %619, label %620, label %627

620:                                              ; preds = %580
  %621 = load ptr, ptr %523, align 8, !tbaa !52
  %622 = getelementptr inbounds i32, ptr %621, i64 %indvars.iv465
  %623 = load i32, ptr %622, align 4, !tbaa !33
  %624 = shl i32 %517, 3
  %625 = mul i32 %624, %623
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %48) #8
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %48, ptr noundef nonnull @.str.4, ptr noundef nonnull @.str.156, i32 noundef %625) #8
  %626 = getelementptr inbounds i8, ptr %3, i64 27610
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %626, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.perturb_timesampling_for_sources, i32 noundef 1814, ptr noundef nonnull %48) #8
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %48) #8
  br label %.loopexit

627:                                              ; preds = %579, %580
  %628 = phi i32 [ %565, %579 ], [ %613, %580 ]
  %629 = phi ptr [ %563, %579 ], [ %611, %580 ]
  %630 = add nuw nsw i32 %.0356437, 1
  %631 = icmp slt i32 %630, %628
  br i1 %631, label %544, label %._crit_edge439.loopexit, !llvm.loop !57

._crit_edge439.loopexit:                          ; preds = %627
  %.pre472 = load ptr, ptr %521, align 8, !tbaa !50
  br label %._crit_edge439

._crit_edge439:                                   ; preds = %._crit_edge439.loopexit, %.preheader
  %632 = phi ptr [ %.pre472, %._crit_edge439.loopexit ], [ %537, %.preheader ]
  %633 = phi ptr [ %629, %._crit_edge439.loopexit ], [ %538, %.preheader ]
  %634 = add nuw nsw i32 %.0357440, 1
  %635 = getelementptr inbounds i32, ptr %632, i64 %indvars.iv465
  %636 = load i32, ptr %635, align 4, !tbaa !33
  %637 = icmp slt i32 %634, %636
  br i1 %637, label %.preheader, label %._crit_edge441, !llvm.loop !58

._crit_edge441:                                   ; preds = %._crit_edge439, %.preheader.lr.ph, %.preheader386
  %638 = phi ptr [ %528, %.preheader386 ], [ %528, %.preheader.lr.ph ], [ %632, %._crit_edge439 ]
  %639 = phi ptr [ %529, %.preheader386 ], [ %529, %.preheader.lr.ph ], [ %632, %._crit_edge439 ]
  %indvars.iv.next466 = add nuw nsw i64 %indvars.iv465, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next466, %wide.trip.count
  br i1 %exitcond.not, label %.loopexit, label %.preheader386, !llvm.loop !60

.loopexit:                                        ; preds = %._crit_edge441, %.loopexit388, %620, %572, %499, %484, %476, %468, %459, %445, %435, %391, %378, %356, %342, %332, %288, %275, %251, %236, %225, %192, %179, %166, %146, %133, %123, %100, %85, %65, %55
  %.0 = phi i32 [ 1, %55 ], [ 1, %65 ], [ 1, %85 ], [ 1, %100 ], [ 1, %123 ], [ 1, %133 ], [ 1, %146 ], [ 1, %166 ], [ 1, %179 ], [ 1, %192 ], [ 1, %275 ], [ 1, %288 ], [ 1, %332 ], [ 1, %342 ], [ 1, %356 ], [ 1, %378 ], [ 1, %391 ], [ 1, %435 ], [ 1, %445 ], [ 1, %459 ], [ 1, %572 ], [ 1, %620 ], [ 1, %468 ], [ 1, %476 ], [ 1, %484 ], [ 1, %499 ], [ 1, %225 ], [ 1, %236 ], [ 1, %251 ], [ 0, %.loopexit388 ], [ 0, %._crit_edge441 ]
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %10) #8
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %9) #8
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %8) #8
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %7) #8
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %6) #8
  call void @llvm.lifetime.end.p0(i64 4, ptr nonnull %5) #8
  ret i32 %.0
}

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

declare i32 @background_at_tau(ptr noundef, double noundef, i16 noundef signext, i16 noundef signext, ptr noundef, ptr noundef) local_unnamed_addr #2

declare i32 @thermodynamics_at_z(ptr noundef, ptr noundef, double noundef, i16 noundef signext, ptr noundef, ptr noundef, ptr noundef) local_unnamed_addr #2

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

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

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

!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, !10, i64 616}
!6 = !{!"background", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !7, i64 40, !7, i64 48, !10, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !11, i64 96, !7, i64 104, !7, i64 112, !7, i64 120, !7, i64 128, !7, i64 136, !7, i64 144, !7, i64 152, !7, i64 160, !7, i64 168, !11, i64 176, !7, i64 184, !7, i64 192, !12, i64 200, !10, i64 208, !10, i64 212, !7, i64 216, !10, i64 224, !12, i64 232, !12, i64 240, !7, i64 248, !12, i64 256, !7, i64 264, !12, i64 272, !7, i64 280, !12, i64 288, !7, i64 296, !12, i64 304, !12, i64 312, !12, i64 320, !7, i64 328, !7, i64 336, !7, i64 344, !7, i64 352, !10, i64 360, !12, i64 368, !7, i64 376, !7, i64 384, !7, i64 392, !7, i64 400, !7, i64 408, !7, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !10, i64 464, !10, i64 468, !10, i64 472, !10, i64 476, !10, i64 480, !10, i64 484, !10, i64 488, !10, i64 492, !10, i64 496, !10, i64 500, !10, i64 504, !10, i64 508, !10, i64 512, !10, i64 516, !10, i64 520, !10, i64 524, !10, i64 528, !10, i64 532, !10, i64 536, !10, i64 540, !10, i64 544, !10, i64 548, !10, i64 552, !10, i64 556, !10, i64 560, !10, i64 564, !10, i64 568, !10, i64 572, !10, i64 576, !10, i64 580, !10, i64 584, !10, i64 588, !10, i64 592, !10, i64 596, !10, i64 600, !10, i64 604, !10, i64 608, !10, i64 612, !10, i64 616, !10, i64 620, !10, i64 624, !10, i64 628, !12, i64 632, !12, i64 640, !12, i64 648, !12, i64 656, !12, i64 664, !10, i64 672, !10, i64 676, !10, i64 680, !10, i64 684, !10, i64 688, !10, i64 692, !10, i64 696, !10, i64 700, !10, i64 704, !10, i64 708, !10, i64 712, !10, i64 716, !10, i64 720, !11, i64 724, !11, i64 726, !11, i64 728, !11, i64 730, !11, i64 732, !11, i64 734, !11, i64 736, !11, i64 738, !11, i64 740, !11, i64 742, !11, i64 744, !12, i64 752, !12, i64 760, !12, i64 768, !12, i64 776, !12, i64 784, !12, i64 792, !12, i64 800, !12, i64 808, !12, i64 816, !12, i64 824, !12, i64 832, !11, i64 840, !11, i64 842, !11, i64 844, !11, i64 846, !11, i64 848, !11, i64 850, !8, i64 852, !11, i64 2900, !8, i64 2902}
!7 = !{!"double", !8, i64 0}
!8 = !{!"omnipotent char", !9, i64 0}
!9 = !{!"Simple C/C++ TBAA"}
!10 = !{!"int", !8, i64 0}
!11 = !{!"short", !8, i64 0}
!12 = !{!"any pointer", !8, i64 0}
!13 = !{!14, !10, i64 372}
!14 = !{!"thermo", !7, i64 0, !10, i64 8, !10, i64 12, !10, i64 16, !7, i64 24, !7, i64 32, !11, i64 40, !11, i64 42, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !10, i64 80, !12, i64 88, !12, i64 96, !7, i64 104, !10, i64 112, !12, i64 120, !12, i64 128, !7, i64 136, !10, i64 144, !12, i64 152, !12, i64 160, !7, i64 168, !11, i64 176, !7, i64 184, !7, i64 192, !7, i64 200, !7, i64 208, !7, i64 216, !7, i64 224, !7, i64 232, !7, i64 240, !7, i64 248, !7, i64 256, !7, i64 264, !10, i64 272, !10, i64 276, !10, i64 280, !10, i64 284, !10, i64 288, !10, i64 292, !10, i64 296, !10, i64 300, !10, i64 304, !10, i64 308, !10, i64 312, !10, i64 316, !10, i64 320, !10, i64 324, !10, i64 328, !10, i64 332, !10, i64 336, !10, i64 340, !10, i64 344, !10, i64 348, !10, i64 352, !10, i64 356, !10, i64 360, !10, i64 364, !10, i64 368, !10, i64 372, !10, i64 376, !12, i64 384, !12, i64 392, !12, i64 400, !7, i64 408, !7, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !7, i64 464, !7, i64 472, !7, i64 480, !7, i64 488, !7, i64 496, !7, i64 504, !7, i64 512, !7, i64 520, !7, i64 528, !7, i64 536, !7, i64 544, !7, i64 552, !7, i64 560, !7, i64 568, !7, i64 576, !7, i64 584, !7, i64 592, !11, i64 600, !11, i64 602, !11, i64 604, !8, i64 606}
!15 = !{!16, !11, i64 2020}
!16 = !{!"perturbs", !11, i64 0, !11, i64 2, !11, i64 4, !11, i64 6, !11, i64 8, !11, i64 10, !11, i64 12, !11, i64 14, !11, i64 16, !11, i64 18, !11, i64 20, !10, i64 24, !11, i64 28, !11, i64 30, !11, i64 32, !11, i64 34, !11, i64 36, !11, i64 38, !11, i64 40, !11, i64 42, !11, i64 44, !11, i64 46, !11, i64 48, !11, i64 50, !11, i64 52, !11, i64 54, !11, i64 56, !11, i64 58, !11, i64 60, !10, i64 64, !10, i64 68, !10, i64 72, !10, i64 76, !7, i64 80, !10, i64 88, !10, i64 92, !8, i64 96, !8, i64 896, !10, i64 1696, !10, i64 1700, !10, i64 1704, !10, i64 1708, !10, i64 1712, !7, i64 1720, !10, i64 1728, !10, i64 1732, !8, i64 1736, !7, i64 1976, !7, i64 1984, !7, i64 1992, !12, i64 2000, !12, i64 2008, !10, i64 2016, !11, i64 2020, !11, i64 2022, !10, i64 2024, !10, i64 2028, !10, i64 2032, !10, i64 2036, !10, i64 2040, !10, i64 2044, !10, i64 2048, !10, i64 2052, !10, i64 2056, !10, i64 2060, !10, i64 2064, !12, i64 2072, !11, i64 2080, !11, i64 2082, !11, i64 2084, !11, i64 2086, !11, i64 2088, !11, i64 2090, !11, i64 2092, !11, i64 2094, !11, i64 2096, !11, i64 2098, !11, i64 2100, !11, i64 2102, !11, i64 2104, !11, i64 2106, !11, i64 2108, !11, i64 2110, !11, i64 2112, !11, i64 2114, !11, i64 2116, !11, i64 2118, !11, i64 2120, !11, i64 2122, !11, i64 2124, !11, i64 2126, !11, i64 2128, !11, i64 2130, !11, i64 2132, !11, i64 2134, !11, i64 2136, !11, i64 2138, !11, i64 2140, !11, i64 2142, !11, i64 2144, !11, i64 2146, !11, i64 2148, !11, i64 2150, !11, i64 2152, !11, i64 2154, !11, i64 2156, !11, i64 2158, !10, i64 2160, !10, i64 2164, !10, i64 2168, !10, i64 2172, !10, i64 2176, !10, i64 2180, !10, i64 2184, !10, i64 2188, !10, i64 2192, !10, i64 2196, !10, i64 2200, !10, i64 2204, !10, i64 2208, !10, i64 2212, !10, i64 2216, !10, i64 2220, !10, i64 2224, !10, i64 2228, !10, i64 2232, !10, i64 2236, !10, i64 2240, !10, i64 2244, !10, i64 2248, !10, i64 2252, !10, i64 2256, !10, i64 2260, !10, i64 2264, !10, i64 2268, !10, i64 2272, !10, i64 2276, !10, i64 2280, !10, i64 2284, !10, i64 2288, !10, i64 2292, !10, i64 2296, !10, i64 2300, !10, i64 2304, !10, i64 2308, !10, i64 2312, !10, i64 2316, !10, i64 2320, !10, i64 2324, !10, i64 2328, !10, i64 2332, !12, i64 2336, !12, i64 2344, !12, i64 2352, !12, i64 2360, !12, i64 2368, !7, i64 2376, !7, i64 2384, !12, i64 2392, !10, i64 2400, !7, i64 2408, !7, i64 2416, !7, i64 2424, !12, i64 2432, !12, i64 2440, !12, i64 2448, !12, i64 2456, !12, i64 2464, !12, i64 2472, !10, i64 2480, !12, i64 2488, !12, i64 2496, !12, i64 2504, !8, i64 2512, !8, i64 10512, !8, i64 18512, !10, i64 26512, !10, i64 26516, !10, i64 26520, !8, i64 26528, !8, i64 26768, !8, i64 27008, !8, i64 27248, !8, i64 27368, !8, i64 27488, !11, i64 27608, !8, i64 27610}
!17 = !{!16, !11, i64 20}
!18 = !{!14, !7, i64 584}
!19 = !{!7, !7, i64 0}
!20 = !{!6, !11, i64 840}
!21 = !{!6, !11, i64 846}
!22 = !{!6, !10, i64 464}
!23 = !{!14, !11, i64 600}
!24 = !{!6, !10, i64 468}
!25 = !{!14, !10, i64 276}
!26 = !{!27, !7, i64 1520}
!27 = !{!"precision", !7, i64 0, !7, i64 8, !7, i64 16, !7, i64 24, !7, i64 32, !7, i64 40, !7, i64 48, !7, i64 56, !7, i64 64, !7, i64 72, !7, i64 80, !7, i64 88, !7, i64 96, !8, i64 104, !7, i64 360, !10, i64 368, !7, i64 376, !10, i64 384, !10, i64 388, !7, i64 392, !10, i64 400, !7, i64 408, !10, i64 416, !7, i64 424, !7, i64 432, !7, i64 440, !7, i64 448, !7, i64 456, !7, i64 464, !7, i64 472, !7, i64 480, !7, i64 488, !7, i64 496, !7, i64 504, !7, i64 512, !7, i64 520, !7, i64 528, !7, i64 536, !7, i64 544, !7, i64 552, !7, i64 560, !7, i64 568, !7, i64 576, !7, i64 584, !7, i64 592, !7, i64 600, !7, i64 608, !7, i64 616, !10, i64 624, !8, i64 628, !8, i64 884, !8, i64 1140, !7, i64 1400, !7, i64 1408, !7, i64 1416, !7, i64 1424, !7, i64 1432, !7, i64 1440, !7, i64 1448, !7, i64 1456, !7, i64 1464, !7, i64 1472, !7, i64 1480, !7, i64 1488, !7, i64 1496, !7, i64 1504, !7, i64 1512, !7, i64 1520, !10, i64 1528, !7, i64 1536, !7, i64 1544, !10, i64 1552, !10, i64 1556, !10, i64 1560, !10, i64 1564, !10, i64 1568, !10, i64 1572, !10, i64 1576, !10, i64 1580, !7, i64 1584, !7, i64 1592, !7, i64 1600, !7, i64 1608, !7, i64 1616, !7, i64 1624, !7, i64 1632, !7, i64 1640, !10, i64 1648, !7, i64 1656, !7, i64 1664, !10, i64 1672, !7, i64 1680, !7, i64 1688, !10, i64 1696, !7, i64 1704, !10, i64 1712, !7, i64 1720, !7, i64 1728, !10, i64 1736, !7, i64 1744, !7, i64 1752, !7, i64 1760, !10, i64 1768, !7, i64 1776, !7, i64 1784, !7, i64 1792, !7, i64 1800, !7, i64 1808, !10, i64 1816, !7, i64 1824, !7, i64 1832, !7, i64 1840, !7, i64 1848, !7, i64 1856, !7, i64 1864, !7, i64 1872, !10, i64 1880, !7, i64 1888, !7, i64 1896, !7, i64 1904, !7, i64 1912, !7, i64 1920, !7, i64 1928, !7, i64 1936, !7, i64 1944, !7, i64 1952, !7, i64 1960, !7, i64 1968, !7, i64 1976, !7, i64 1984, !7, i64 1992, !7, i64 2000, !7, i64 2008, !7, i64 2016, !7, i64 2024, !7, i64 2032, !7, i64 2040, !7, i64 2048, !7, i64 2056, !7, i64 2064, !7, i64 2072, !7, i64 2080, !7, i64 2088, !7, i64 2096, !7, i64 2104, !7, i64 2112, !7, i64 2120, !7, i64 2128, !7, i64 2136, !7, i64 2144, !7, i64 2152, !7, i64 2160, !7, i64 2168, !7, i64 2176, !7, i64 2184, !7, i64 2192, !7, i64 2200, !7, i64 2208, !7, i64 2216, !7, i64 2224, !7, i64 2232, !7, i64 2240, !7, i64 2248, !10, i64 2256, !7, i64 2264, !7, i64 2272, !7, i64 2280, !7, i64 2288, !10, i64 2296, !7, i64 2304, !7, i64 2312, !10, i64 2320, !10, i64 2324, !10, i64 2328, !7, i64 2336, !7, i64 2344, !8, i64 2352}
!28 = !{!14, !7, i64 416}
!29 = !{!27, !7, i64 1640}
!30 = distinct !{!30, !31}
!31 = !{!"llvm.loop.mustprogress"}
!32 = !{!16, !7, i64 1992}
!33 = !{!10, !10, i64 0}
!34 = !{!6, !7, i64 344}
!35 = !{!6, !11, i64 848}
!36 = !{!14, !11, i64 602}
!37 = !{!14, !10, i64 364}
!38 = !{!6, !10, i64 472}
!39 = !{!27, !7, i64 1616}
!40 = !{!27, !7, i64 2344}
!41 = distinct !{!41, !31}
!42 = !{!16, !10, i64 2400}
!43 = !{!16, !12, i64 2392}
!44 = distinct !{!44, !31}
!45 = !{!16, !10, i64 2480}
!46 = distinct !{!46, !31}
!47 = !{!16, !12, i64 2472}
!48 = distinct !{!48, !31}
!49 = !{!16, !10, i64 2040}
!50 = !{!16, !12, i64 2072}
!51 = !{!16, !12, i64 2336}
!52 = !{!16, !12, i64 2360}
!53 = !{!16, !12, i64 2464}
!54 = !{!12, !12, i64 0}
!55 = !{!16, !12, i64 2488}
!56 = !{!16, !12, i64 2496}
!57 = distinct !{!57, !31}
!58 = distinct !{!58, !31, !59}
!59 = !{!"llvm.loop.unswitch.partial.disable"}
!60 = distinct !{!60, !31}


; Effective O3 passes: AggressiveInstCombinePass,BDCEPass,CorrelatedValuePropagationPass,EarlyCSEPass,GVNPass,GlobalOptPass,IPSCCPPass,IndVarSimplifyPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,JumpThreadingPass,LCSSAPass,LICMPass,LoopRotatePass,LoopSimplifyCFGPass,LoopSimplifyPass,LoopUnrollPass,PostOrderFunctionAttrsPass,ReassociatePass,SROAPass,SimpleLoopUnswitchPass,SimplifyCFGPass,TailCallElimPass