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

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

@.str = external hidden unnamed_addr constant [18 x i8], align 1
@.str.1 = external hidden unnamed_addr constant [156 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@__func__.thermodynamics_at_z = external hidden unnamed_addr constant [20 x i8], align 1
@.str.3 = external hidden unnamed_addr constant [189 x i8], align 1
@.str.4 = external hidden unnamed_addr constant [205 x i8], align 1

; Function Attrs: nounwind uwtable
define dso_local i32 @thermodynamics_at_z(ptr noundef %0, ptr noundef %1, double noundef %2, i16 noundef signext %3, ptr noundef %4, ptr noundef %5, ptr noundef %6) #0 {
  %8 = alloca i32, align 4
  %9 = alloca ptr, align 8
  %10 = alloca ptr, align 8
  %11 = alloca double, align 8
  %12 = alloca i16, align 2
  %13 = alloca ptr, align 8
  %14 = alloca ptr, align 8
  %15 = alloca ptr, align 8
  %16 = alloca double, align 8
  %17 = alloca [2048 x i8], align 16
  %18 = alloca i32, align 4
  %19 = alloca [2048 x i8], align 16
  %20 = alloca [2048 x i8], align 16
  store ptr %0, ptr %9, align 8, !tbaa !5
  store ptr %1, ptr %10, align 8, !tbaa !5
  store double %2, ptr %11, align 8, !tbaa !9
  store i16 %3, ptr %12, align 2, !tbaa !11
  store ptr %4, ptr %13, align 8, !tbaa !5
  store ptr %5, ptr %14, align 8, !tbaa !5
  store ptr %6, ptr %15, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #5
  %21 = load double, ptr %11, align 8, !tbaa !9
  %22 = load ptr, ptr %10, align 8, !tbaa !5
  %23 = getelementptr inbounds %struct.thermo, ptr %22, i32 0, i32 63
  %24 = load ptr, ptr %23, align 8, !tbaa !13
  %25 = load ptr, ptr %10, align 8, !tbaa !5
  %26 = getelementptr inbounds %struct.thermo, ptr %25, i32 0, i32 62
  %27 = load i32, ptr %26, align 8, !tbaa !16
  %28 = sub nsw i32 %27, 1
  %29 = sext i32 %28 to i64
  %30 = getelementptr inbounds double, ptr %24, i64 %29
  %31 = load double, ptr %30, align 8, !tbaa !9
  %32 = fcmp oge double %21, %31
  br i1 %32, label %33, label %774

33:                                               ; preds = %7
  %34 = load ptr, ptr %10, align 8, !tbaa !5
  %35 = getelementptr inbounds %struct.thermo, ptr %34, i32 0, i32 64
  %36 = load ptr, ptr %35, align 8, !tbaa !17
  %37 = load ptr, ptr %10, align 8, !tbaa !5
  %38 = getelementptr inbounds %struct.thermo, ptr %37, i32 0, i32 62
  %39 = load i32, ptr %38, align 8, !tbaa !16
  %40 = sub nsw i32 %39, 1
  %41 = load ptr, ptr %10, align 8, !tbaa !5
  %42 = getelementptr inbounds %struct.thermo, ptr %41, i32 0, i32 61
  %43 = load i32, ptr %42, align 4, !tbaa !18
  %44 = mul nsw i32 %40, %43
  %45 = load ptr, ptr %10, align 8, !tbaa !5
  %46 = getelementptr inbounds %struct.thermo, ptr %45, i32 0, i32 36
  %47 = load i32, ptr %46, align 8, !tbaa !19
  %48 = add nsw i32 %44, %47
  %49 = sext i32 %48 to i64
  %50 = getelementptr inbounds double, ptr %36, i64 %49
  %51 = load double, ptr %50, align 8, !tbaa !9
  store double %51, ptr %16, align 8, !tbaa !9
  %52 = load double, ptr %16, align 8, !tbaa !9
  %53 = load ptr, ptr %15, align 8, !tbaa !5
  %54 = load ptr, ptr %10, align 8, !tbaa !5
  %55 = getelementptr inbounds %struct.thermo, ptr %54, i32 0, i32 36
  %56 = load i32, ptr %55, align 8, !tbaa !19
  %57 = sext i32 %56 to i64
  %58 = getelementptr inbounds double, ptr %53, i64 %57
  store double %52, ptr %58, align 8, !tbaa !9
  %59 = load double, ptr %11, align 8, !tbaa !9
  %60 = fadd double 1.000000e+00, %59
  %61 = load double, ptr %11, align 8, !tbaa !9
  %62 = fadd double 1.000000e+00, %61
  %63 = fmul double %60, %62
  %64 = load ptr, ptr %10, align 8, !tbaa !5
  %65 = getelementptr inbounds %struct.thermo, ptr %64, i32 0, i32 89
  %66 = load double, ptr %65, align 8, !tbaa !20
  %67 = fmul double %63, %66
  %68 = load double, ptr %16, align 8, !tbaa !9
  %69 = fmul double %67, %68
  %70 = fmul double %69, 0x3A15151E37FBF236
  %71 = fmul double %70, 0x449A22FF260E3826
  %72 = load ptr, ptr %15, align 8, !tbaa !5
  %73 = load ptr, ptr %10, align 8, !tbaa !5
  %74 = getelementptr inbounds %struct.thermo, ptr %73, i32 0, i32 37
  %75 = load i32, ptr %74, align 4, !tbaa !21
  %76 = sext i32 %75 to i64
  %77 = getelementptr inbounds double, ptr %72, i64 %76
  store double %71, ptr %77, align 8, !tbaa !9
  %78 = load ptr, ptr %10, align 8, !tbaa !5
  %79 = getelementptr inbounds %struct.thermo, ptr %78, i32 0, i32 64
  %80 = load ptr, ptr %79, align 8, !tbaa !17
  %81 = load ptr, ptr %10, align 8, !tbaa !5
  %82 = getelementptr inbounds %struct.thermo, ptr %81, i32 0, i32 62
  %83 = load i32, ptr %82, align 8, !tbaa !16
  %84 = sub nsw i32 %83, 1
  %85 = load ptr, ptr %10, align 8, !tbaa !5
  %86 = getelementptr inbounds %struct.thermo, ptr %85, i32 0, i32 61
  %87 = load i32, ptr %86, align 4, !tbaa !18
  %88 = mul nsw i32 %84, %87
  %89 = load ptr, ptr %10, align 8, !tbaa !5
  %90 = getelementptr inbounds %struct.thermo, ptr %89, i32 0, i32 38
  %91 = load i32, ptr %90, align 8, !tbaa !22
  %92 = add nsw i32 %88, %91
  %93 = sext i32 %92 to i64
  %94 = getelementptr inbounds double, ptr %80, i64 %93
  %95 = load double, ptr %94, align 8, !tbaa !9
  %96 = load double, ptr %11, align 8, !tbaa !9
  %97 = fadd double 1.000000e+00, %96
  %98 = load ptr, ptr %10, align 8, !tbaa !5
  %99 = getelementptr inbounds %struct.thermo, ptr %98, i32 0, i32 63
  %100 = load ptr, ptr %99, align 8, !tbaa !13
  %101 = load ptr, ptr %10, align 8, !tbaa !5
  %102 = getelementptr inbounds %struct.thermo, ptr %101, i32 0, i32 62
  %103 = load i32, ptr %102, align 8, !tbaa !16
  %104 = sub nsw i32 %103, 1
  %105 = sext i32 %104 to i64
  %106 = getelementptr inbounds double, ptr %100, i64 %105
  %107 = load double, ptr %106, align 8, !tbaa !9
  %108 = fadd double 1.000000e+00, %107
  %109 = fdiv double %97, %108
  %110 = call double @pow(double noundef %109, double noundef 2.000000e+00) #5
  %111 = fmul double %95, %110
  %112 = load ptr, ptr %15, align 8, !tbaa !5
  %113 = load ptr, ptr %10, align 8, !tbaa !5
  %114 = getelementptr inbounds %struct.thermo, ptr %113, i32 0, i32 38
  %115 = load i32, ptr %114, align 8, !tbaa !22
  %116 = sext i32 %115 to i64
  %117 = getelementptr inbounds double, ptr %112, i64 %116
  store double %111, ptr %117, align 8, !tbaa !9
  %118 = load ptr, ptr %10, align 8, !tbaa !5
  %119 = getelementptr inbounds %struct.thermo, ptr %118, i32 0, i32 7
  %120 = load i16, ptr %119, align 2, !tbaa !23
  %121 = sext i16 %120 to i32
  %122 = icmp eq i32 %121, 1
  br i1 %122, label %123, label %164

123:                                              ; preds = %33
  %124 = load ptr, ptr %10, align 8, !tbaa !5
  %125 = getelementptr inbounds %struct.thermo, ptr %124, i32 0, i32 64
  %126 = load ptr, ptr %125, align 8, !tbaa !17
  %127 = load ptr, ptr %10, align 8, !tbaa !5
  %128 = getelementptr inbounds %struct.thermo, ptr %127, i32 0, i32 62
  %129 = load i32, ptr %128, align 8, !tbaa !16
  %130 = sub nsw i32 %129, 1
  %131 = load ptr, ptr %10, align 8, !tbaa !5
  %132 = getelementptr inbounds %struct.thermo, ptr %131, i32 0, i32 61
  %133 = load i32, ptr %132, align 4, !tbaa !18
  %134 = mul nsw i32 %130, %133
  %135 = load ptr, ptr %10, align 8, !tbaa !5
  %136 = getelementptr inbounds %struct.thermo, ptr %135, i32 0, i32 60
  %137 = load i32, ptr %136, align 8, !tbaa !24
  %138 = add nsw i32 %134, %137
  %139 = sext i32 %138 to i64
  %140 = getelementptr inbounds double, ptr %126, i64 %139
  %141 = load double, ptr %140, align 8, !tbaa !9
  %142 = load double, ptr %11, align 8, !tbaa !9
  %143 = fadd double 1.000000e+00, %142
  %144 = load ptr, ptr %10, align 8, !tbaa !5
  %145 = getelementptr inbounds %struct.thermo, ptr %144, i32 0, i32 63
  %146 = load ptr, ptr %145, align 8, !tbaa !13
  %147 = load ptr, ptr %10, align 8, !tbaa !5
  %148 = getelementptr inbounds %struct.thermo, ptr %147, i32 0, i32 62
  %149 = load i32, ptr %148, align 8, !tbaa !16
  %150 = sub nsw i32 %149, 1
  %151 = sext i32 %150 to i64
  %152 = getelementptr inbounds double, ptr %146, i64 %151
  %153 = load double, ptr %152, align 8, !tbaa !9
  %154 = fadd double 1.000000e+00, %153
  %155 = fdiv double %143, %154
  %156 = call double @pow(double noundef %155, double noundef -1.500000e+00) #5
  %157 = fmul double %141, %156
  %158 = load ptr, ptr %15, align 8, !tbaa !5
  %159 = load ptr, ptr %10, align 8, !tbaa !5
  %160 = getelementptr inbounds %struct.thermo, ptr %159, i32 0, i32 60
  %161 = load i32, ptr %160, align 8, !tbaa !24
  %162 = sext i32 %161 to i64
  %163 = getelementptr inbounds double, ptr %158, i64 %162
  store double %157, ptr %163, align 8, !tbaa !9
  br label %164

164:                                              ; preds = %123, %33
  %165 = load ptr, ptr %14, align 8, !tbaa !5
  %166 = load ptr, ptr %9, align 8, !tbaa !5
  %167 = getelementptr inbounds %struct.background, ptr %166, i32 0, i32 60
  %168 = load i32, ptr %167, align 4, !tbaa !25
  %169 = sext i32 %168 to i64
  %170 = getelementptr inbounds double, ptr %165, i64 %169
  %171 = load double, ptr %170, align 8, !tbaa !9
  %172 = fneg double %171
  %173 = fmul double %172, 2.000000e+00
  %174 = load double, ptr %11, align 8, !tbaa !9
  %175 = fadd double 1.000000e+00, %174
  %176 = fdiv double %173, %175
  %177 = load ptr, ptr %15, align 8, !tbaa !5
  %178 = load ptr, ptr %10, align 8, !tbaa !5
  %179 = getelementptr inbounds %struct.thermo, ptr %178, i32 0, i32 37
  %180 = load i32, ptr %179, align 4, !tbaa !21
  %181 = sext i32 %180 to i64
  %182 = getelementptr inbounds double, ptr %177, i64 %181
  %183 = load double, ptr %182, align 8, !tbaa !9
  %184 = fmul double %176, %183
  %185 = load ptr, ptr %15, align 8, !tbaa !5
  %186 = load ptr, ptr %10, align 8, !tbaa !5
  %187 = getelementptr inbounds %struct.thermo, ptr %186, i32 0, i32 39
  %188 = load i32, ptr %187, align 4, !tbaa !27
  %189 = sext i32 %188 to i64
  %190 = getelementptr inbounds double, ptr %185, i64 %189
  store double %184, ptr %190, align 8, !tbaa !9
  %191 = load ptr, ptr %14, align 8, !tbaa !5
  %192 = load ptr, ptr %9, align 8, !tbaa !5
  %193 = getelementptr inbounds %struct.background, ptr %192, i32 0, i32 60
  %194 = load i32, ptr %193, align 4, !tbaa !25
  %195 = sext i32 %194 to i64
  %196 = getelementptr inbounds double, ptr %191, i64 %195
  %197 = load double, ptr %196, align 8, !tbaa !9
  %198 = load ptr, ptr %14, align 8, !tbaa !5
  %199 = load ptr, ptr %9, align 8, !tbaa !5
  %200 = getelementptr inbounds %struct.background, ptr %199, i32 0, i32 60
  %201 = load i32, ptr %200, align 4, !tbaa !25
  %202 = sext i32 %201 to i64
  %203 = getelementptr inbounds double, ptr %198, i64 %202
  %204 = load double, ptr %203, align 8, !tbaa !9
  %205 = fmul double %197, %204
  %206 = load double, ptr %11, align 8, !tbaa !9
  %207 = fadd double 1.000000e+00, %206
  %208 = fdiv double %205, %207
  %209 = load ptr, ptr %14, align 8, !tbaa !5
  %210 = load ptr, ptr %9, align 8, !tbaa !5
  %211 = getelementptr inbounds %struct.background, ptr %210, i32 0, i32 61
  %212 = load i32, ptr %211, align 8, !tbaa !28
  %213 = sext i32 %212 to i64
  %214 = getelementptr inbounds double, ptr %209, i64 %213
  %215 = load double, ptr %214, align 8, !tbaa !9
  %216 = fsub double %208, %215
  %217 = fmul double %216, 2.000000e+00
  %218 = load double, ptr %11, align 8, !tbaa !9
  %219 = fadd double 1.000000e+00, %218
  %220 = fdiv double %217, %219
  %221 = load ptr, ptr %15, align 8, !tbaa !5
  %222 = load ptr, ptr %10, align 8, !tbaa !5
  %223 = getelementptr inbounds %struct.thermo, ptr %222, i32 0, i32 37
  %224 = load i32, ptr %223, align 4, !tbaa !21
  %225 = sext i32 %224 to i64
  %226 = getelementptr inbounds double, ptr %221, i64 %225
  %227 = load double, ptr %226, align 8, !tbaa !9
  %228 = fmul double %220, %227
  %229 = load ptr, ptr %15, align 8, !tbaa !5
  %230 = load ptr, ptr %10, align 8, !tbaa !5
  %231 = getelementptr inbounds %struct.thermo, ptr %230, i32 0, i32 40
  %232 = load i32, ptr %231, align 8, !tbaa !29
  %233 = sext i32 %232 to i64
  %234 = getelementptr inbounds double, ptr %229, i64 %233
  store double %228, ptr %234, align 8, !tbaa !9
  %235 = load ptr, ptr %15, align 8, !tbaa !5
  %236 = load ptr, ptr %10, align 8, !tbaa !5
  %237 = getelementptr inbounds %struct.thermo, ptr %236, i32 0, i32 41
  %238 = load i32, ptr %237, align 4, !tbaa !30
  %239 = sext i32 %238 to i64
  %240 = getelementptr inbounds double, ptr %235, i64 %239
  store double 0.000000e+00, ptr %240, align 8, !tbaa !9
  %241 = load ptr, ptr %15, align 8, !tbaa !5
  %242 = load ptr, ptr %10, align 8, !tbaa !5
  %243 = getelementptr inbounds %struct.thermo, ptr %242, i32 0, i32 42
  %244 = load i32, ptr %243, align 8, !tbaa !31
  %245 = sext i32 %244 to i64
  %246 = getelementptr inbounds double, ptr %241, i64 %245
  store double 0.000000e+00, ptr %246, align 8, !tbaa !9
  %247 = load ptr, ptr %15, align 8, !tbaa !5
  %248 = load ptr, ptr %10, align 8, !tbaa !5
  %249 = getelementptr inbounds %struct.thermo, ptr %248, i32 0, i32 43
  %250 = load i32, ptr %249, align 4, !tbaa !32
  %251 = sext i32 %250 to i64
  %252 = getelementptr inbounds double, ptr %247, i64 %251
  store double 0.000000e+00, ptr %252, align 8, !tbaa !9
  %253 = load ptr, ptr %15, align 8, !tbaa !5
  %254 = load ptr, ptr %10, align 8, !tbaa !5
  %255 = getelementptr inbounds %struct.thermo, ptr %254, i32 0, i32 44
  %256 = load i32, ptr %255, align 8, !tbaa !33
  %257 = sext i32 %256 to i64
  %258 = getelementptr inbounds double, ptr %253, i64 %257
  store double 0.000000e+00, ptr %258, align 8, !tbaa !9
  %259 = load ptr, ptr %9, align 8, !tbaa !5
  %260 = getelementptr inbounds %struct.background, ptr %259, i32 0, i32 2
  %261 = load double, ptr %260, align 8, !tbaa !34
  %262 = load double, ptr %11, align 8, !tbaa !9
  %263 = fadd double 1.000000e+00, %262
  %264 = fmul double %261, %263
  %265 = load ptr, ptr %15, align 8, !tbaa !5
  %266 = load ptr, ptr %10, align 8, !tbaa !5
  %267 = getelementptr inbounds %struct.thermo, ptr %266, i32 0, i32 54
  %268 = load i32, ptr %267, align 8, !tbaa !35
  %269 = sext i32 %268 to i64
  %270 = getelementptr inbounds double, ptr %265, i64 %269
  store double %264, ptr %270, align 8, !tbaa !9
  %271 = load ptr, ptr %10, align 8, !tbaa !5
  %272 = getelementptr inbounds %struct.thermo, ptr %271, i32 0, i32 0
  %273 = load double, ptr %272, align 8, !tbaa !36
  %274 = call double @llvm.fmuladd.f64(double 0xBFE7F14DA42541E8, double %273, double 1.000000e+00)
  %275 = load double, ptr %16, align 8, !tbaa !9
  %276 = load ptr, ptr %10, align 8, !tbaa !5
  %277 = getelementptr inbounds %struct.thermo, ptr %276, i32 0, i32 0
  %278 = load double, ptr %277, align 8, !tbaa !36
  %279 = fsub double 1.000000e+00, %278
  %280 = call double @llvm.fmuladd.f64(double %275, double %279, double %274)
  %281 = fmul double 0x3D39D63187BB40C4, %280
  %282 = load ptr, ptr %9, align 8, !tbaa !5
  %283 = getelementptr inbounds %struct.background, ptr %282, i32 0, i32 2
  %284 = load double, ptr %283, align 8, !tbaa !34
  %285 = fmul double %281, %284
  %286 = load double, ptr %11, align 8, !tbaa !9
  %287 = fadd double 1.000000e+00, %286
  %288 = fmul double %285, %287
  %289 = load ptr, ptr %15, align 8, !tbaa !5
  %290 = load ptr, ptr %10, align 8, !tbaa !5
  %291 = getelementptr inbounds %struct.thermo, ptr %290, i32 0, i32 55
  %292 = load i32, ptr %291, align 4, !tbaa !37
  %293 = sext i32 %292 to i64
  %294 = getelementptr inbounds double, ptr %289, i64 %293
  store double %288, ptr %294, align 8, !tbaa !9
  %295 = load ptr, ptr %15, align 8, !tbaa !5
  %296 = load ptr, ptr %10, align 8, !tbaa !5
  %297 = getelementptr inbounds %struct.thermo, ptr %296, i32 0, i32 55
  %298 = load i32, ptr %297, align 4, !tbaa !37
  %299 = sext i32 %298 to i64
  %300 = getelementptr inbounds double, ptr %295, i64 %299
  %301 = load double, ptr %300, align 8, !tbaa !9
  %302 = fmul double %301, 4.000000e+00
  %303 = fdiv double %302, 3.000000e+00
  %304 = load ptr, ptr %15, align 8, !tbaa !5
  %305 = load ptr, ptr %10, align 8, !tbaa !5
  %306 = getelementptr inbounds %struct.thermo, ptr %305, i32 0, i32 56
  %307 = load i32, ptr %306, align 8, !tbaa !38
  %308 = sext i32 %307 to i64
  %309 = getelementptr inbounds double, ptr %304, i64 %308
  store double %303, ptr %309, align 8, !tbaa !9
  %310 = load ptr, ptr %10, align 8, !tbaa !5
  %311 = getelementptr inbounds %struct.thermo, ptr %310, i32 0, i32 6
  %312 = load i16, ptr %311, align 8, !tbaa !39
  %313 = sext i16 %312 to i32
  %314 = icmp eq i32 %313, 1
  br i1 %314, label %315, label %376

315:                                              ; preds = %164
  %316 = load ptr, ptr %14, align 8, !tbaa !5
  %317 = load ptr, ptr %9, align 8, !tbaa !5
  %318 = getelementptr inbounds %struct.background, ptr %317, i32 0, i32 60
  %319 = load i32, ptr %318, align 4, !tbaa !25
  %320 = sext i32 %319 to i64
  %321 = getelementptr inbounds double, ptr %316, i64 %320
  %322 = load double, ptr %321, align 8, !tbaa !9
  %323 = fneg double %322
  %324 = load ptr, ptr %14, align 8, !tbaa !5
  %325 = load ptr, ptr %9, align 8, !tbaa !5
  %326 = getelementptr inbounds %struct.background, ptr %325, i32 0, i32 59
  %327 = load i32, ptr %326, align 8, !tbaa !40
  %328 = sext i32 %327 to i64
  %329 = getelementptr inbounds double, ptr %324, i64 %328
  %330 = load double, ptr %329, align 8, !tbaa !9
  %331 = fmul double %323, %330
  %332 = load ptr, ptr %15, align 8, !tbaa !5
  %333 = load ptr, ptr %10, align 8, !tbaa !5
  %334 = getelementptr inbounds %struct.thermo, ptr %333, i32 0, i32 56
  %335 = load i32, ptr %334, align 8, !tbaa !38
  %336 = sext i32 %335 to i64
  %337 = getelementptr inbounds double, ptr %332, i64 %336
  %338 = load double, ptr %337, align 8, !tbaa !9
  %339 = fmul double %331, %338
  %340 = load ptr, ptr %15, align 8, !tbaa !5
  %341 = load ptr, ptr %10, align 8, !tbaa !5
  %342 = getelementptr inbounds %struct.thermo, ptr %341, i32 0, i32 57
  %343 = load i32, ptr %342, align 4, !tbaa !41
  %344 = sext i32 %343 to i64
  %345 = getelementptr inbounds double, ptr %340, i64 %344
  store double %339, ptr %345, align 8, !tbaa !9
  %346 = load ptr, ptr %14, align 8, !tbaa !5
  %347 = load ptr, ptr %9, align 8, !tbaa !5
  %348 = getelementptr inbounds %struct.background, ptr %347, i32 0, i32 61
  %349 = load i32, ptr %348, align 8, !tbaa !28
  %350 = sext i32 %349 to i64
  %351 = getelementptr inbounds double, ptr %346, i64 %350
  %352 = load double, ptr %351, align 8, !tbaa !9
  %353 = fneg double %352
  %354 = load ptr, ptr %14, align 8, !tbaa !5
  %355 = load ptr, ptr %9, align 8, !tbaa !5
  %356 = getelementptr inbounds %struct.background, ptr %355, i32 0, i32 59
  %357 = load i32, ptr %356, align 8, !tbaa !40
  %358 = sext i32 %357 to i64
  %359 = getelementptr inbounds double, ptr %354, i64 %358
  %360 = load double, ptr %359, align 8, !tbaa !9
  %361 = fmul double %353, %360
  %362 = load ptr, ptr %15, align 8, !tbaa !5
  %363 = load ptr, ptr %10, align 8, !tbaa !5
  %364 = getelementptr inbounds %struct.thermo, ptr %363, i32 0, i32 56
  %365 = load i32, ptr %364, align 8, !tbaa !38
  %366 = sext i32 %365 to i64
  %367 = getelementptr inbounds double, ptr %362, i64 %366
  %368 = load double, ptr %367, align 8, !tbaa !9
  %369 = fmul double %361, %368
  %370 = load ptr, ptr %15, align 8, !tbaa !5
  %371 = load ptr, ptr %10, align 8, !tbaa !5
  %372 = getelementptr inbounds %struct.thermo, ptr %371, i32 0, i32 58
  %373 = load i32, ptr %372, align 8, !tbaa !42
  %374 = sext i32 %373 to i64
  %375 = getelementptr inbounds double, ptr %370, i64 %374
  store double %369, ptr %375, align 8, !tbaa !9
  br label %376

376:                                              ; preds = %315, %164
  %377 = load ptr, ptr %15, align 8, !tbaa !5
  %378 = load ptr, ptr %10, align 8, !tbaa !5
  %379 = getelementptr inbounds %struct.thermo, ptr %378, i32 0, i32 37
  %380 = load i32, ptr %379, align 4, !tbaa !21
  %381 = sext i32 %380 to i64
  %382 = getelementptr inbounds double, ptr %377, i64 %381
  %383 = load double, ptr %382, align 8, !tbaa !9
  %384 = load ptr, ptr %15, align 8, !tbaa !5
  %385 = load ptr, ptr %10, align 8, !tbaa !5
  %386 = getelementptr inbounds %struct.thermo, ptr %385, i32 0, i32 59
  %387 = load i32, ptr %386, align 4, !tbaa !43
  %388 = sext i32 %387 to i64
  %389 = getelementptr inbounds double, ptr %384, i64 %388
  store double %383, ptr %389, align 8, !tbaa !9
  %390 = load ptr, ptr %9, align 8, !tbaa !5
  %391 = getelementptr inbounds %struct.background, ptr %390, i32 0, i32 128
  %392 = load i16, ptr %391, align 2, !tbaa !44
  %393 = sext i16 %392 to i32
  %394 = icmp eq i32 %393, 1
  br i1 %394, label %395, label %773

395:                                              ; preds = %376
  %396 = load ptr, ptr %10, align 8, !tbaa !5
  %397 = getelementptr inbounds %struct.thermo, ptr %396, i32 0, i32 32
  %398 = load double, ptr %397, align 8, !tbaa !45
  %399 = load double, ptr %11, align 8, !tbaa !9
  %400 = fadd double 1.000000e+00, %399
  %401 = fdiv double %400, 1.000000e+07
  %402 = load ptr, ptr %10, align 8, !tbaa !5
  %403 = getelementptr inbounds %struct.thermo, ptr %402, i32 0, i32 34
  %404 = load double, ptr %403, align 8, !tbaa !46
  %405 = call double @pow(double noundef %401, double noundef %404) #5
  %406 = fmul double %398, %405
  %407 = load ptr, ptr %9, align 8, !tbaa !5
  %408 = getelementptr inbounds %struct.background, ptr %407, i32 0, i32 17
  %409 = load double, ptr %408, align 8, !tbaa !47
  %410 = fmul double %406, %409
  %411 = load ptr, ptr %9, align 8, !tbaa !5
  %412 = getelementptr inbounds %struct.background, ptr %411, i32 0, i32 42
  %413 = load double, ptr %412, align 8, !tbaa !48
  %414 = call double @pow(double noundef %413, double noundef 2.000000e+00) #5
  %415 = fmul double %410, %414
  %416 = load ptr, ptr %15, align 8, !tbaa !5
  %417 = load ptr, ptr %10, align 8, !tbaa !5
  %418 = getelementptr inbounds %struct.thermo, ptr %417, i32 0, i32 45
  %419 = load i32, ptr %418, align 4, !tbaa !49
  %420 = sext i32 %419 to i64
  %421 = getelementptr inbounds double, ptr %416, i64 %420
  store double %415, ptr %421, align 8, !tbaa !9
  %422 = load ptr, ptr %14, align 8, !tbaa !5
  %423 = load ptr, ptr %9, align 8, !tbaa !5
  %424 = getelementptr inbounds %struct.background, ptr %423, i32 0, i32 60
  %425 = load i32, ptr %424, align 4, !tbaa !25
  %426 = sext i32 %425 to i64
  %427 = getelementptr inbounds double, ptr %422, i64 %426
  %428 = load double, ptr %427, align 8, !tbaa !9
  %429 = fneg double %428
  %430 = load ptr, ptr %10, align 8, !tbaa !5
  %431 = getelementptr inbounds %struct.thermo, ptr %430, i32 0, i32 34
  %432 = load double, ptr %431, align 8, !tbaa !46
  %433 = fmul double %429, %432
  %434 = load double, ptr %11, align 8, !tbaa !9
  %435 = fadd double 1.000000e+00, %434
  %436 = fdiv double %433, %435
  %437 = load ptr, ptr %15, align 8, !tbaa !5
  %438 = load ptr, ptr %10, align 8, !tbaa !5
  %439 = getelementptr inbounds %struct.thermo, ptr %438, i32 0, i32 45
  %440 = load i32, ptr %439, align 4, !tbaa !49
  %441 = sext i32 %440 to i64
  %442 = getelementptr inbounds double, ptr %437, i64 %441
  %443 = load double, ptr %442, align 8, !tbaa !9
  %444 = fmul double %436, %443
  %445 = load ptr, ptr %15, align 8, !tbaa !5
  %446 = load ptr, ptr %10, align 8, !tbaa !5
  %447 = getelementptr inbounds %struct.thermo, ptr %446, i32 0, i32 46
  %448 = load i32, ptr %447, align 8, !tbaa !50
  %449 = sext i32 %448 to i64
  %450 = getelementptr inbounds double, ptr %445, i64 %449
  store double %444, ptr %450, align 8, !tbaa !9
  %451 = load ptr, ptr %14, align 8, !tbaa !5
  %452 = load ptr, ptr %9, align 8, !tbaa !5
  %453 = getelementptr inbounds %struct.background, ptr %452, i32 0, i32 60
  %454 = load i32, ptr %453, align 4, !tbaa !25
  %455 = sext i32 %454 to i64
  %456 = getelementptr inbounds double, ptr %451, i64 %455
  %457 = load double, ptr %456, align 8, !tbaa !9
  %458 = load ptr, ptr %14, align 8, !tbaa !5
  %459 = load ptr, ptr %9, align 8, !tbaa !5
  %460 = getelementptr inbounds %struct.background, ptr %459, i32 0, i32 60
  %461 = load i32, ptr %460, align 4, !tbaa !25
  %462 = sext i32 %461 to i64
  %463 = getelementptr inbounds double, ptr %458, i64 %462
  %464 = load double, ptr %463, align 8, !tbaa !9
  %465 = fmul double %457, %464
  %466 = load double, ptr %11, align 8, !tbaa !9
  %467 = fadd double 1.000000e+00, %466
  %468 = fdiv double %465, %467
  %469 = load ptr, ptr %14, align 8, !tbaa !5
  %470 = load ptr, ptr %9, align 8, !tbaa !5
  %471 = getelementptr inbounds %struct.background, ptr %470, i32 0, i32 61
  %472 = load i32, ptr %471, align 8, !tbaa !28
  %473 = sext i32 %472 to i64
  %474 = getelementptr inbounds double, ptr %469, i64 %473
  %475 = load double, ptr %474, align 8, !tbaa !9
  %476 = fsub double %468, %475
  %477 = load ptr, ptr %10, align 8, !tbaa !5
  %478 = getelementptr inbounds %struct.thermo, ptr %477, i32 0, i32 34
  %479 = load double, ptr %478, align 8, !tbaa !46
  %480 = fmul double %476, %479
  %481 = load double, ptr %11, align 8, !tbaa !9
  %482 = fadd double 1.000000e+00, %481
  %483 = fdiv double %480, %482
  %484 = load ptr, ptr %15, align 8, !tbaa !5
  %485 = load ptr, ptr %10, align 8, !tbaa !5
  %486 = getelementptr inbounds %struct.thermo, ptr %485, i32 0, i32 45
  %487 = load i32, ptr %486, align 4, !tbaa !49
  %488 = sext i32 %487 to i64
  %489 = getelementptr inbounds double, ptr %484, i64 %488
  %490 = load double, ptr %489, align 8, !tbaa !9
  %491 = fmul double %483, %490
  %492 = load ptr, ptr %15, align 8, !tbaa !5
  %493 = load ptr, ptr %10, align 8, !tbaa !5
  %494 = getelementptr inbounds %struct.thermo, ptr %493, i32 0, i32 47
  %495 = load i32, ptr %494, align 4, !tbaa !51
  %496 = sext i32 %495 to i64
  %497 = getelementptr inbounds double, ptr %492, i64 %496
  store double %491, ptr %497, align 8, !tbaa !9
  %498 = load ptr, ptr %10, align 8, !tbaa !5
  %499 = getelementptr inbounds %struct.thermo, ptr %498, i32 0, i32 33
  %500 = load double, ptr %499, align 8, !tbaa !52
  %501 = load double, ptr %11, align 8, !tbaa !9
  %502 = fadd double 1.000000e+00, %501
  %503 = fdiv double %502, 1.000000e+07
  %504 = load ptr, ptr %10, align 8, !tbaa !5
  %505 = getelementptr inbounds %struct.thermo, ptr %504, i32 0, i32 34
  %506 = load double, ptr %505, align 8, !tbaa !46
  %507 = call double @pow(double noundef %503, double noundef %506) #5
  %508 = fmul double %500, %507
  %509 = load ptr, ptr %9, align 8, !tbaa !5
  %510 = getelementptr inbounds %struct.background, ptr %509, i32 0, i32 15
  %511 = load double, ptr %510, align 8, !tbaa !53
  %512 = fmul double %508, %511
  %513 = load ptr, ptr %9, align 8, !tbaa !5
  %514 = getelementptr inbounds %struct.background, ptr %513, i32 0, i32 42
  %515 = load double, ptr %514, align 8, !tbaa !48
  %516 = call double @pow(double noundef %515, double noundef 2.000000e+00) #5
  %517 = fmul double %512, %516
  %518 = load ptr, ptr %15, align 8, !tbaa !5
  %519 = load ptr, ptr %10, align 8, !tbaa !5
  %520 = getelementptr inbounds %struct.thermo, ptr %519, i32 0, i32 48
  %521 = load i32, ptr %520, align 8, !tbaa !54
  %522 = sext i32 %521 to i64
  %523 = getelementptr inbounds double, ptr %518, i64 %522
  store double %517, ptr %523, align 8, !tbaa !9
  %524 = load ptr, ptr %10, align 8, !tbaa !5
  %525 = getelementptr inbounds %struct.thermo, ptr %524, i32 0, i32 64
  %526 = load ptr, ptr %525, align 8, !tbaa !17
  %527 = load ptr, ptr %10, align 8, !tbaa !5
  %528 = getelementptr inbounds %struct.thermo, ptr %527, i32 0, i32 62
  %529 = load i32, ptr %528, align 8, !tbaa !16
  %530 = sub nsw i32 %529, 1
  %531 = load ptr, ptr %10, align 8, !tbaa !5
  %532 = getelementptr inbounds %struct.thermo, ptr %531, i32 0, i32 61
  %533 = load i32, ptr %532, align 4, !tbaa !18
  %534 = mul nsw i32 %530, %533
  %535 = load ptr, ptr %10, align 8, !tbaa !5
  %536 = getelementptr inbounds %struct.thermo, ptr %535, i32 0, i32 49
  %537 = load i32, ptr %536, align 4, !tbaa !55
  %538 = add nsw i32 %534, %537
  %539 = sext i32 %538 to i64
  %540 = getelementptr inbounds double, ptr %526, i64 %539
  %541 = load double, ptr %540, align 8, !tbaa !9
  %542 = load ptr, ptr %10, align 8, !tbaa !5
  %543 = getelementptr inbounds %struct.thermo, ptr %542, i32 0, i32 64
  %544 = load ptr, ptr %543, align 8, !tbaa !17
  %545 = load ptr, ptr %10, align 8, !tbaa !5
  %546 = getelementptr inbounds %struct.thermo, ptr %545, i32 0, i32 62
  %547 = load i32, ptr %546, align 8, !tbaa !16
  %548 = sub nsw i32 %547, 1
  %549 = load ptr, ptr %10, align 8, !tbaa !5
  %550 = getelementptr inbounds %struct.thermo, ptr %549, i32 0, i32 61
  %551 = load i32, ptr %550, align 4, !tbaa !18
  %552 = mul nsw i32 %548, %551
  %553 = load ptr, ptr %10, align 8, !tbaa !5
  %554 = getelementptr inbounds %struct.thermo, ptr %553, i32 0, i32 49
  %555 = load i32, ptr %554, align 4, !tbaa !55
  %556 = add nsw i32 %552, %555
  %557 = sext i32 %556 to i64
  %558 = getelementptr inbounds double, ptr %544, i64 %557
  %559 = load double, ptr %558, align 8, !tbaa !9
  %560 = load ptr, ptr %10, align 8, !tbaa !5
  %561 = getelementptr inbounds %struct.thermo, ptr %560, i32 0, i32 64
  %562 = load ptr, ptr %561, align 8, !tbaa !17
  %563 = load ptr, ptr %10, align 8, !tbaa !5
  %564 = getelementptr inbounds %struct.thermo, ptr %563, i32 0, i32 62
  %565 = load i32, ptr %564, align 8, !tbaa !16
  %566 = sub nsw i32 %565, 2
  %567 = load ptr, ptr %10, align 8, !tbaa !5
  %568 = getelementptr inbounds %struct.thermo, ptr %567, i32 0, i32 61
  %569 = load i32, ptr %568, align 4, !tbaa !18
  %570 = mul nsw i32 %566, %569
  %571 = load ptr, ptr %10, align 8, !tbaa !5
  %572 = getelementptr inbounds %struct.thermo, ptr %571, i32 0, i32 49
  %573 = load i32, ptr %572, align 4, !tbaa !55
  %574 = add nsw i32 %570, %573
  %575 = sext i32 %574 to i64
  %576 = getelementptr inbounds double, ptr %562, i64 %575
  %577 = load double, ptr %576, align 8, !tbaa !9
  %578 = fsub double %559, %577
  %579 = load double, ptr %11, align 8, !tbaa !9
  %580 = load ptr, ptr %10, align 8, !tbaa !5
  %581 = getelementptr inbounds %struct.thermo, ptr %580, i32 0, i32 63
  %582 = load ptr, ptr %581, align 8, !tbaa !13
  %583 = load ptr, ptr %10, align 8, !tbaa !5
  %584 = getelementptr inbounds %struct.thermo, ptr %583, i32 0, i32 62
  %585 = load i32, ptr %584, align 8, !tbaa !16
  %586 = sub nsw i32 %585, 1
  %587 = sext i32 %586 to i64
  %588 = getelementptr inbounds double, ptr %582, i64 %587
  %589 = load double, ptr %588, align 8, !tbaa !9
  %590 = fsub double %579, %589
  %591 = fmul double %578, %590
  %592 = load ptr, ptr %10, align 8, !tbaa !5
  %593 = getelementptr inbounds %struct.thermo, ptr %592, i32 0, i32 63
  %594 = load ptr, ptr %593, align 8, !tbaa !13
  %595 = load ptr, ptr %10, align 8, !tbaa !5
  %596 = getelementptr inbounds %struct.thermo, ptr %595, i32 0, i32 62
  %597 = load i32, ptr %596, align 8, !tbaa !16
  %598 = sub nsw i32 %597, 1
  %599 = sext i32 %598 to i64
  %600 = getelementptr inbounds double, ptr %594, i64 %599
  %601 = load double, ptr %600, align 8, !tbaa !9
  %602 = load ptr, ptr %10, align 8, !tbaa !5
  %603 = getelementptr inbounds %struct.thermo, ptr %602, i32 0, i32 63
  %604 = load ptr, ptr %603, align 8, !tbaa !13
  %605 = load ptr, ptr %10, align 8, !tbaa !5
  %606 = getelementptr inbounds %struct.thermo, ptr %605, i32 0, i32 62
  %607 = load i32, ptr %606, align 8, !tbaa !16
  %608 = sub nsw i32 %607, 2
  %609 = sext i32 %608 to i64
  %610 = getelementptr inbounds double, ptr %604, i64 %609
  %611 = load double, ptr %610, align 8, !tbaa !9
  %612 = fsub double %601, %611
  %613 = fdiv double %591, %612
  %614 = fadd double %541, %613
  %615 = load ptr, ptr %15, align 8, !tbaa !5
  %616 = load ptr, ptr %10, align 8, !tbaa !5
  %617 = getelementptr inbounds %struct.thermo, ptr %616, i32 0, i32 49
  %618 = load i32, ptr %617, align 4, !tbaa !55
  %619 = sext i32 %618 to i64
  %620 = getelementptr inbounds double, ptr %615, i64 %619
  store double %614, ptr %620, align 8, !tbaa !9
  %621 = load ptr, ptr %10, align 8, !tbaa !5
  %622 = getelementptr inbounds %struct.thermo, ptr %621, i32 0, i32 64
  %623 = load ptr, ptr %622, align 8, !tbaa !17
  %624 = load ptr, ptr %10, align 8, !tbaa !5
  %625 = getelementptr inbounds %struct.thermo, ptr %624, i32 0, i32 62
  %626 = load i32, ptr %625, align 8, !tbaa !16
  %627 = sub nsw i32 %626, 1
  %628 = load ptr, ptr %10, align 8, !tbaa !5
  %629 = getelementptr inbounds %struct.thermo, ptr %628, i32 0, i32 61
  %630 = load i32, ptr %629, align 4, !tbaa !18
  %631 = mul nsw i32 %627, %630
  %632 = load ptr, ptr %10, align 8, !tbaa !5
  %633 = getelementptr inbounds %struct.thermo, ptr %632, i32 0, i32 50
  %634 = load i32, ptr %633, align 8, !tbaa !56
  %635 = add nsw i32 %631, %634
  %636 = sext i32 %635 to i64
  %637 = getelementptr inbounds double, ptr %623, i64 %636
  %638 = load double, ptr %637, align 8, !tbaa !9
  %639 = load ptr, ptr %10, align 8, !tbaa !5
  %640 = getelementptr inbounds %struct.thermo, ptr %639, i32 0, i32 64
  %641 = load ptr, ptr %640, align 8, !tbaa !17
  %642 = load ptr, ptr %10, align 8, !tbaa !5
  %643 = getelementptr inbounds %struct.thermo, ptr %642, i32 0, i32 62
  %644 = load i32, ptr %643, align 8, !tbaa !16
  %645 = sub nsw i32 %644, 1
  %646 = load ptr, ptr %10, align 8, !tbaa !5
  %647 = getelementptr inbounds %struct.thermo, ptr %646, i32 0, i32 61
  %648 = load i32, ptr %647, align 4, !tbaa !18
  %649 = mul nsw i32 %645, %648
  %650 = load ptr, ptr %10, align 8, !tbaa !5
  %651 = getelementptr inbounds %struct.thermo, ptr %650, i32 0, i32 50
  %652 = load i32, ptr %651, align 8, !tbaa !56
  %653 = add nsw i32 %649, %652
  %654 = sext i32 %653 to i64
  %655 = getelementptr inbounds double, ptr %641, i64 %654
  %656 = load double, ptr %655, align 8, !tbaa !9
  %657 = load ptr, ptr %10, align 8, !tbaa !5
  %658 = getelementptr inbounds %struct.thermo, ptr %657, i32 0, i32 64
  %659 = load ptr, ptr %658, align 8, !tbaa !17
  %660 = load ptr, ptr %10, align 8, !tbaa !5
  %661 = getelementptr inbounds %struct.thermo, ptr %660, i32 0, i32 62
  %662 = load i32, ptr %661, align 8, !tbaa !16
  %663 = sub nsw i32 %662, 2
  %664 = load ptr, ptr %10, align 8, !tbaa !5
  %665 = getelementptr inbounds %struct.thermo, ptr %664, i32 0, i32 61
  %666 = load i32, ptr %665, align 4, !tbaa !18
  %667 = mul nsw i32 %663, %666
  %668 = load ptr, ptr %10, align 8, !tbaa !5
  %669 = getelementptr inbounds %struct.thermo, ptr %668, i32 0, i32 50
  %670 = load i32, ptr %669, align 8, !tbaa !56
  %671 = add nsw i32 %667, %670
  %672 = sext i32 %671 to i64
  %673 = getelementptr inbounds double, ptr %659, i64 %672
  %674 = load double, ptr %673, align 8, !tbaa !9
  %675 = fsub double %656, %674
  %676 = load double, ptr %11, align 8, !tbaa !9
  %677 = load ptr, ptr %10, align 8, !tbaa !5
  %678 = getelementptr inbounds %struct.thermo, ptr %677, i32 0, i32 63
  %679 = load ptr, ptr %678, align 8, !tbaa !13
  %680 = load ptr, ptr %10, align 8, !tbaa !5
  %681 = getelementptr inbounds %struct.thermo, ptr %680, i32 0, i32 62
  %682 = load i32, ptr %681, align 8, !tbaa !16
  %683 = sub nsw i32 %682, 1
  %684 = sext i32 %683 to i64
  %685 = getelementptr inbounds double, ptr %679, i64 %684
  %686 = load double, ptr %685, align 8, !tbaa !9
  %687 = fsub double %676, %686
  %688 = fmul double %675, %687
  %689 = load ptr, ptr %10, align 8, !tbaa !5
  %690 = getelementptr inbounds %struct.thermo, ptr %689, i32 0, i32 63
  %691 = load ptr, ptr %690, align 8, !tbaa !13
  %692 = load ptr, ptr %10, align 8, !tbaa !5
  %693 = getelementptr inbounds %struct.thermo, ptr %692, i32 0, i32 62
  %694 = load i32, ptr %693, align 8, !tbaa !16
  %695 = sub nsw i32 %694, 1
  %696 = sext i32 %695 to i64
  %697 = getelementptr inbounds double, ptr %691, i64 %696
  %698 = load double, ptr %697, align 8, !tbaa !9
  %699 = load ptr, ptr %10, align 8, !tbaa !5
  %700 = getelementptr inbounds %struct.thermo, ptr %699, i32 0, i32 63
  %701 = load ptr, ptr %700, align 8, !tbaa !13
  %702 = load ptr, ptr %10, align 8, !tbaa !5
  %703 = getelementptr inbounds %struct.thermo, ptr %702, i32 0, i32 62
  %704 = load i32, ptr %703, align 8, !tbaa !16
  %705 = sub nsw i32 %704, 2
  %706 = sext i32 %705 to i64
  %707 = getelementptr inbounds double, ptr %701, i64 %706
  %708 = load double, ptr %707, align 8, !tbaa !9
  %709 = fsub double %698, %708
  %710 = fdiv double %688, %709
  %711 = fadd double %638, %710
  %712 = load ptr, ptr %15, align 8, !tbaa !5
  %713 = load ptr, ptr %10, align 8, !tbaa !5
  %714 = getelementptr inbounds %struct.thermo, ptr %713, i32 0, i32 50
  %715 = load i32, ptr %714, align 8, !tbaa !56
  %716 = sext i32 %715 to i64
  %717 = getelementptr inbounds double, ptr %712, i64 %716
  store double %711, ptr %717, align 8, !tbaa !9
  %718 = load ptr, ptr %10, align 8, !tbaa !5
  %719 = getelementptr inbounds %struct.thermo, ptr %718, i32 0, i32 64
  %720 = load ptr, ptr %719, align 8, !tbaa !17
  %721 = load ptr, ptr %10, align 8, !tbaa !5
  %722 = getelementptr inbounds %struct.thermo, ptr %721, i32 0, i32 62
  %723 = load i32, ptr %722, align 8, !tbaa !16
  %724 = sub nsw i32 %723, 1
  %725 = load ptr, ptr %10, align 8, !tbaa !5
  %726 = getelementptr inbounds %struct.thermo, ptr %725, i32 0, i32 61
  %727 = load i32, ptr %726, align 4, !tbaa !18
  %728 = mul nsw i32 %724, %727
  %729 = load ptr, ptr %10, align 8, !tbaa !5
  %730 = getelementptr inbounds %struct.thermo, ptr %729, i32 0, i32 51
  %731 = load i32, ptr %730, align 4, !tbaa !57
  %732 = add nsw i32 %728, %731
  %733 = sext i32 %732 to i64
  %734 = getelementptr inbounds double, ptr %720, i64 %733
  %735 = load double, ptr %734, align 8, !tbaa !9
  %736 = load ptr, ptr %15, align 8, !tbaa !5
  %737 = load ptr, ptr %10, align 8, !tbaa !5
  %738 = getelementptr inbounds %struct.thermo, ptr %737, i32 0, i32 51
  %739 = load i32, ptr %738, align 4, !tbaa !57
  %740 = sext i32 %739 to i64
  %741 = getelementptr inbounds double, ptr %736, i64 %740
  store double %735, ptr %741, align 8, !tbaa !9
  %742 = load ptr, ptr %9, align 8, !tbaa !5
  %743 = getelementptr inbounds %struct.background, ptr %742, i32 0, i32 16
  %744 = load double, ptr %743, align 8, !tbaa !58
  %745 = fmul double 0x3B50B0E7F1528C6D, %744
  %746 = load double, ptr %11, align 8, !tbaa !9
  %747 = fadd double 1.000000e+00, %746
  %748 = fmul double %745, %747
  %749 = fdiv double %748, 0x3C07A4DA04A6DAF5
  %750 = fdiv double %749, 3.000000e+00
  %751 = load ptr, ptr %10, align 8, !tbaa !5
  %752 = getelementptr inbounds %struct.thermo, ptr %751, i32 0, i32 35
  %753 = load double, ptr %752, align 8, !tbaa !59
  %754 = fdiv double %750, %753
  %755 = load ptr, ptr %15, align 8, !tbaa !5
  %756 = load ptr, ptr %10, align 8, !tbaa !5
  %757 = getelementptr inbounds %struct.thermo, ptr %756, i32 0, i32 52
  %758 = load i32, ptr %757, align 8, !tbaa !60
  %759 = sext i32 %758 to i64
  %760 = getelementptr inbounds double, ptr %755, i64 %759
  store double %754, ptr %760, align 8, !tbaa !9
  %761 = load ptr, ptr %9, align 8, !tbaa !5
  %762 = getelementptr inbounds %struct.background, ptr %761, i32 0, i32 16
  %763 = load double, ptr %762, align 8, !tbaa !58
  %764 = load double, ptr %11, align 8, !tbaa !9
  %765 = fadd double 1.000000e+00, %764
  %766 = fmul double %763, %765
  %767 = load ptr, ptr %15, align 8, !tbaa !5
  %768 = load ptr, ptr %10, align 8, !tbaa !5
  %769 = getelementptr inbounds %struct.thermo, ptr %768, i32 0, i32 53
  %770 = load i32, ptr %769, align 4, !tbaa !61
  %771 = sext i32 %770 to i64
  %772 = getelementptr inbounds double, ptr %767, i64 %771
  store double %766, ptr %772, align 8, !tbaa !9
  br label %773

773:                                              ; preds = %395, %376
  br label %920

774:                                              ; preds = %7
  %775 = load ptr, ptr %10, align 8, !tbaa !5
  %776 = getelementptr inbounds %struct.thermo, ptr %775, i32 0, i32 2
  %777 = load i32, ptr %776, align 4, !tbaa !62
  %778 = icmp eq i32 %777, 3
  br i1 %778, label %779, label %786

779:                                              ; preds = %774
  %780 = load double, ptr %11, align 8, !tbaa !9
  %781 = load ptr, ptr %10, align 8, !tbaa !5
  %782 = getelementptr inbounds %struct.thermo, ptr %781, i32 0, i32 5
  %783 = load double, ptr %782, align 8, !tbaa !63
  %784 = fmul double 2.000000e+00, %783
  %785 = fcmp olt double %780, %784
  br i1 %785, label %794, label %786

786:                                              ; preds = %779, %774
  %787 = load ptr, ptr %10, align 8, !tbaa !5
  %788 = getelementptr inbounds %struct.thermo, ptr %787, i32 0, i32 2
  %789 = load i32, ptr %788, align 4, !tbaa !62
  %790 = icmp eq i32 %789, 5
  br i1 %790, label %791, label %828

791:                                              ; preds = %786
  %792 = load double, ptr %11, align 8, !tbaa !9
  %793 = fcmp olt double %792, 5.000000e+01
  br i1 %793, label %794, label %828

794:                                              ; preds = %791, %779
  %795 = load ptr, ptr %10, align 8, !tbaa !5
  %796 = getelementptr inbounds %struct.thermo, ptr %795, i32 0, i32 63
  %797 = load ptr, ptr %796, align 8, !tbaa !13
  %798 = load ptr, ptr %10, align 8, !tbaa !5
  %799 = getelementptr inbounds %struct.thermo, ptr %798, i32 0, i32 62
  %800 = load i32, ptr %799, align 8, !tbaa !16
  %801 = load ptr, ptr %10, align 8, !tbaa !5
  %802 = getelementptr inbounds %struct.thermo, ptr %801, i32 0, i32 64
  %803 = load ptr, ptr %802, align 8, !tbaa !17
  %804 = load ptr, ptr %10, align 8, !tbaa !5
  %805 = getelementptr inbounds %struct.thermo, ptr %804, i32 0, i32 61
  %806 = load i32, ptr %805, align 4, !tbaa !18
  %807 = load double, ptr %11, align 8, !tbaa !9
  %808 = load ptr, ptr %13, align 8, !tbaa !5
  %809 = load ptr, ptr %15, align 8, !tbaa !5
  %810 = load ptr, ptr %10, align 8, !tbaa !5
  %811 = getelementptr inbounds %struct.thermo, ptr %810, i32 0, i32 61
  %812 = load i32, ptr %811, align 4, !tbaa !18
  %813 = load ptr, ptr %10, align 8, !tbaa !5
  %814 = getelementptr inbounds %struct.thermo, ptr %813, i32 0, i32 93
  %815 = getelementptr inbounds [2048 x i8], ptr %814, i64 0, i64 0
  %816 = call i32 @array_interpolate_linear(ptr noundef %797, i32 noundef %800, ptr noundef %803, i32 noundef %806, double noundef %807, ptr noundef %808, ptr noundef %809, i32 noundef %812, ptr noundef %815)
  %817 = icmp eq i32 %816, 1
  br i1 %817, label %818, label %827

818:                                              ; preds = %794
  call void @llvm.lifetime.start.p0(i64 2048, ptr %17) #5
  %819 = getelementptr inbounds [2048 x i8], ptr %17, i64 0, i64 0
  %820 = load ptr, ptr %10, align 8, !tbaa !5
  %821 = getelementptr inbounds %struct.thermo, ptr %820, i32 0, i32 93
  %822 = getelementptr inbounds [2048 x i8], ptr %821, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %819, ptr noundef @.str, ptr noundef @.str.1, ptr noundef %822)
  %823 = load ptr, ptr %10, align 8, !tbaa !5
  %824 = getelementptr inbounds %struct.thermo, ptr %823, i32 0, i32 93
  %825 = getelementptr inbounds [2048 x i8], ptr %824, i64 0, i64 0
  %826 = getelementptr inbounds [2048 x i8], ptr %17, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %825, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_at_z, i32 noundef 233, ptr noundef %826)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %17) #5
  store i32 1, ptr %8, align 4
  store i32 1, ptr %18, align 4
  br label %921

827:                                              ; preds = %794
  br label %919

828:                                              ; preds = %791, %786
  %829 = load i16, ptr %12, align 2, !tbaa !11
  %830 = sext i16 %829 to i32
  %831 = load ptr, ptr %10, align 8, !tbaa !5
  %832 = getelementptr inbounds %struct.thermo, ptr %831, i32 0, i32 90
  %833 = load i16, ptr %832, align 8, !tbaa !64
  %834 = sext i16 %833 to i32
  %835 = icmp eq i32 %830, %834
  br i1 %835, label %836, label %873

836:                                              ; preds = %828
  %837 = load ptr, ptr %10, align 8, !tbaa !5
  %838 = getelementptr inbounds %struct.thermo, ptr %837, i32 0, i32 63
  %839 = load ptr, ptr %838, align 8, !tbaa !13
  %840 = load ptr, ptr %10, align 8, !tbaa !5
  %841 = getelementptr inbounds %struct.thermo, ptr %840, i32 0, i32 62
  %842 = load i32, ptr %841, align 8, !tbaa !16
  %843 = load ptr, ptr %10, align 8, !tbaa !5
  %844 = getelementptr inbounds %struct.thermo, ptr %843, i32 0, i32 64
  %845 = load ptr, ptr %844, align 8, !tbaa !17
  %846 = load ptr, ptr %10, align 8, !tbaa !5
  %847 = getelementptr inbounds %struct.thermo, ptr %846, i32 0, i32 65
  %848 = load ptr, ptr %847, align 8, !tbaa !65
  %849 = load ptr, ptr %10, align 8, !tbaa !5
  %850 = getelementptr inbounds %struct.thermo, ptr %849, i32 0, i32 61
  %851 = load i32, ptr %850, align 4, !tbaa !18
  %852 = load double, ptr %11, align 8, !tbaa !9
  %853 = load ptr, ptr %13, align 8, !tbaa !5
  %854 = load ptr, ptr %15, align 8, !tbaa !5
  %855 = load ptr, ptr %10, align 8, !tbaa !5
  %856 = getelementptr inbounds %struct.thermo, ptr %855, i32 0, i32 61
  %857 = load i32, ptr %856, align 4, !tbaa !18
  %858 = load ptr, ptr %10, align 8, !tbaa !5
  %859 = getelementptr inbounds %struct.thermo, ptr %858, i32 0, i32 93
  %860 = getelementptr inbounds [2048 x i8], ptr %859, i64 0, i64 0
  %861 = call i32 @array_interpolate_spline(ptr noundef %839, i32 noundef %842, ptr noundef %845, ptr noundef %848, i32 noundef %851, double noundef %852, ptr noundef %853, ptr noundef %854, i32 noundef %857, ptr noundef %860)
  %862 = icmp eq i32 %861, 1
  br i1 %862, label %863, label %872

863:                                              ; preds = %836
  call void @llvm.lifetime.start.p0(i64 2048, ptr %19) #5
  %864 = getelementptr inbounds [2048 x i8], ptr %19, i64 0, i64 0
  %865 = load ptr, ptr %10, align 8, !tbaa !5
  %866 = getelementptr inbounds %struct.thermo, ptr %865, i32 0, i32 93
  %867 = getelementptr inbounds [2048 x i8], ptr %866, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %864, ptr noundef @.str, ptr noundef @.str.3, ptr noundef %867)
  %868 = load ptr, ptr %10, align 8, !tbaa !5
  %869 = getelementptr inbounds %struct.thermo, ptr %868, i32 0, i32 93
  %870 = getelementptr inbounds [2048 x i8], ptr %869, i64 0, i64 0
  %871 = getelementptr inbounds [2048 x i8], ptr %19, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %870, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_at_z, i32 noundef 253, ptr noundef %871)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %19) #5
  store i32 1, ptr %8, align 4
  store i32 1, ptr %18, align 4
  br label %921

872:                                              ; preds = %836
  br label %873

873:                                              ; preds = %872, %828
  %874 = load i16, ptr %12, align 2, !tbaa !11
  %875 = sext i16 %874 to i32
  %876 = load ptr, ptr %10, align 8, !tbaa !5
  %877 = getelementptr inbounds %struct.thermo, ptr %876, i32 0, i32 91
  %878 = load i16, ptr %877, align 2, !tbaa !66
  %879 = sext i16 %878 to i32
  %880 = icmp eq i32 %875, %879
  br i1 %880, label %881, label %918

881:                                              ; preds = %873
  %882 = load ptr, ptr %10, align 8, !tbaa !5
  %883 = getelementptr inbounds %struct.thermo, ptr %882, i32 0, i32 63
  %884 = load ptr, ptr %883, align 8, !tbaa !13
  %885 = load ptr, ptr %10, align 8, !tbaa !5
  %886 = getelementptr inbounds %struct.thermo, ptr %885, i32 0, i32 62
  %887 = load i32, ptr %886, align 8, !tbaa !16
  %888 = load ptr, ptr %10, align 8, !tbaa !5
  %889 = getelementptr inbounds %struct.thermo, ptr %888, i32 0, i32 64
  %890 = load ptr, ptr %889, align 8, !tbaa !17
  %891 = load ptr, ptr %10, align 8, !tbaa !5
  %892 = getelementptr inbounds %struct.thermo, ptr %891, i32 0, i32 65
  %893 = load ptr, ptr %892, align 8, !tbaa !65
  %894 = load ptr, ptr %10, align 8, !tbaa !5
  %895 = getelementptr inbounds %struct.thermo, ptr %894, i32 0, i32 61
  %896 = load i32, ptr %895, align 4, !tbaa !18
  %897 = load double, ptr %11, align 8, !tbaa !9
  %898 = load ptr, ptr %13, align 8, !tbaa !5
  %899 = load ptr, ptr %15, align 8, !tbaa !5
  %900 = load ptr, ptr %10, align 8, !tbaa !5
  %901 = getelementptr inbounds %struct.thermo, ptr %900, i32 0, i32 61
  %902 = load i32, ptr %901, align 4, !tbaa !18
  %903 = load ptr, ptr %10, align 8, !tbaa !5
  %904 = getelementptr inbounds %struct.thermo, ptr %903, i32 0, i32 93
  %905 = getelementptr inbounds [2048 x i8], ptr %904, i64 0, i64 0
  %906 = call i32 @array_interpolate_spline_growing_closeby(ptr noundef %884, i32 noundef %887, ptr noundef %890, ptr noundef %893, i32 noundef %896, double noundef %897, ptr noundef %898, ptr noundef %899, i32 noundef %902, ptr noundef %905)
  %907 = icmp eq i32 %906, 1
  br i1 %907, label %908, label %917

908:                                              ; preds = %881
  call void @llvm.lifetime.start.p0(i64 2048, ptr %20) #5
  %909 = getelementptr inbounds [2048 x i8], ptr %20, i64 0, i64 0
  %910 = load ptr, ptr %10, align 8, !tbaa !5
  %911 = getelementptr inbounds %struct.thermo, ptr %910, i32 0, i32 93
  %912 = getelementptr inbounds [2048 x i8], ptr %911, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %909, ptr noundef @.str, ptr noundef @.str.4, ptr noundef %912)
  %913 = load ptr, ptr %10, align 8, !tbaa !5
  %914 = getelementptr inbounds %struct.thermo, ptr %913, i32 0, i32 93
  %915 = getelementptr inbounds [2048 x i8], ptr %914, i64 0, i64 0
  %916 = getelementptr inbounds [2048 x i8], ptr %20, i64 0, i64 0
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef %915, ptr noundef @.str.2, ptr noundef @__func__.thermodynamics_at_z, i32 noundef 270, ptr noundef %916)
  call void @llvm.lifetime.end.p0(i64 2048, ptr %20) #5
  store i32 1, ptr %8, align 4
  store i32 1, ptr %18, align 4
  br label %921

917:                                              ; preds = %881
  br label %918

918:                                              ; preds = %917, %873
  br label %919

919:                                              ; preds = %918, %827
  br label %920

920:                                              ; preds = %919, %773
  store i32 0, ptr %8, align 4
  store i32 1, ptr %18, align 4
  br label %921

921:                                              ; preds = %920, %908, %863, %818
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #5
  %922 = load i32, ptr %8, align 4
  ret i32 %922
}

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

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

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

declare i32 @array_interpolate_linear(ptr noundef, i32 noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) #4

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

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

declare i32 @array_interpolate_spline(ptr noundef, i32 noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) #4

declare i32 @array_interpolate_spline_growing_closeby(ptr noundef, i32 noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) #4

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

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !6, i64 0}
!6 = !{!"any pointer", !7, i64 0}
!7 = !{!"omnipotent char", !8, i64 0}
!8 = !{!"Simple C/C++ TBAA"}
!9 = !{!10, !10, i64 0}
!10 = !{!"double", !7, i64 0}
!11 = !{!12, !12, i64 0}
!12 = !{!"short", !7, i64 0}
!13 = !{!14, !6, i64 384}
!14 = !{!"thermo", !10, i64 0, !15, i64 8, !15, i64 12, !15, i64 16, !10, i64 24, !10, i64 32, !12, i64 40, !12, i64 42, !10, i64 48, !10, i64 56, !10, i64 64, !10, i64 72, !15, i64 80, !6, i64 88, !6, i64 96, !10, i64 104, !15, i64 112, !6, i64 120, !6, i64 128, !10, i64 136, !15, i64 144, !6, i64 152, !6, i64 160, !10, i64 168, !12, i64 176, !10, i64 184, !10, i64 192, !10, i64 200, !10, i64 208, !10, i64 216, !10, i64 224, !10, i64 232, !10, i64 240, !10, i64 248, !10, i64 256, !10, i64 264, !15, i64 272, !15, i64 276, !15, i64 280, !15, i64 284, !15, i64 288, !15, i64 292, !15, i64 296, !15, i64 300, !15, i64 304, !15, i64 308, !15, i64 312, !15, i64 316, !15, i64 320, !15, i64 324, !15, i64 328, !15, i64 332, !15, i64 336, !15, i64 340, !15, i64 344, !15, i64 348, !15, i64 352, !15, i64 356, !15, i64 360, !15, i64 364, !15, i64 368, !15, i64 372, !15, i64 376, !6, i64 384, !6, i64 392, !6, i64 400, !10, i64 408, !10, i64 416, !10, i64 424, !10, i64 432, !10, i64 440, !10, i64 448, !10, i64 456, !10, i64 464, !10, i64 472, !10, i64 480, !10, i64 488, !10, i64 496, !10, i64 504, !10, i64 512, !10, i64 520, !10, i64 528, !10, i64 536, !10, i64 544, !10, i64 552, !10, i64 560, !10, i64 568, !10, i64 576, !10, i64 584, !10, i64 592, !12, i64 600, !12, i64 602, !12, i64 604, !7, i64 606}
!15 = !{!"int", !7, i64 0}
!16 = !{!14, !15, i64 376}
!17 = !{!14, !6, i64 392}
!18 = !{!14, !15, i64 372}
!19 = !{!14, !15, i64 272}
!20 = !{!14, !10, i64 592}
!21 = !{!14, !15, i64 276}
!22 = !{!14, !15, i64 280}
!23 = !{!14, !12, i64 42}
!24 = !{!14, !15, i64 368}
!25 = !{!26, !15, i64 468}
!26 = !{!"background", !10, i64 0, !10, i64 8, !10, i64 16, !10, i64 24, !10, i64 32, !10, i64 40, !10, i64 48, !15, i64 56, !10, i64 64, !10, i64 72, !10, i64 80, !10, i64 88, !12, i64 96, !10, i64 104, !10, i64 112, !10, i64 120, !10, i64 128, !10, i64 136, !10, i64 144, !10, i64 152, !10, i64 160, !10, i64 168, !12, i64 176, !10, i64 184, !10, i64 192, !6, i64 200, !15, i64 208, !15, i64 212, !10, i64 216, !15, i64 224, !6, i64 232, !6, i64 240, !10, i64 248, !6, i64 256, !10, i64 264, !6, i64 272, !10, i64 280, !6, i64 288, !10, i64 296, !6, i64 304, !6, i64 312, !6, i64 320, !10, i64 328, !10, i64 336, !10, i64 344, !10, i64 352, !15, i64 360, !6, i64 368, !10, i64 376, !10, i64 384, !10, i64 392, !10, i64 400, !10, i64 408, !10, i64 416, !10, i64 424, !10, i64 432, !10, i64 440, !10, i64 448, !10, i64 456, !15, i64 464, !15, i64 468, !15, i64 472, !15, i64 476, !15, i64 480, !15, i64 484, !15, i64 488, !15, i64 492, !15, i64 496, !15, i64 500, !15, i64 504, !15, i64 508, !15, i64 512, !15, i64 516, !15, i64 520, !15, i64 524, !15, i64 528, !15, i64 532, !15, i64 536, !15, i64 540, !15, i64 544, !15, i64 548, !15, i64 552, !15, i64 556, !15, i64 560, !15, i64 564, !15, i64 568, !15, i64 572, !15, i64 576, !15, i64 580, !15, i64 584, !15, i64 588, !15, i64 592, !15, i64 596, !15, i64 600, !15, i64 604, !15, i64 608, !15, i64 612, !15, i64 616, !15, i64 620, !15, i64 624, !15, i64 628, !6, i64 632, !6, i64 640, !6, i64 648, !6, i64 656, !6, i64 664, !15, i64 672, !15, i64 676, !15, i64 680, !15, i64 684, !15, i64 688, !15, i64 692, !15, i64 696, !15, i64 700, !15, i64 704, !15, i64 708, !15, i64 712, !15, i64 716, !15, i64 720, !12, i64 724, !12, i64 726, !12, i64 728, !12, i64 730, !12, i64 732, !12, i64 734, !12, i64 736, !12, i64 738, !12, i64 740, !12, i64 742, !12, 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, !12, i64 840, !12, i64 842, !12, i64 844, !12, i64 846, !12, i64 848, !12, i64 850, !7, i64 852, !12, i64 2900, !7, i64 2902}
!27 = !{!14, !15, i64 284}
!28 = !{!26, !15, i64 472}
!29 = !{!14, !15, i64 288}
!30 = !{!14, !15, i64 292}
!31 = !{!14, !15, i64 296}
!32 = !{!14, !15, i64 300}
!33 = !{!14, !15, i64 304}
!34 = !{!26, !10, i64 16}
!35 = !{!14, !15, i64 344}
!36 = !{!14, !10, i64 0}
!37 = !{!14, !15, i64 348}
!38 = !{!14, !15, i64 352}
!39 = !{!14, !12, i64 40}
!40 = !{!26, !15, i64 464}
!41 = !{!14, !15, i64 356}
!42 = !{!14, !15, i64 360}
!43 = !{!14, !15, i64 364}
!44 = !{!26, !12, i64 742}
!45 = !{!14, !10, i64 240}
!46 = !{!14, !10, i64 256}
!47 = !{!26, !10, i64 136}
!48 = !{!26, !10, i64 328}
!49 = !{!14, !15, i64 308}
!50 = !{!14, !15, i64 312}
!51 = !{!14, !15, i64 316}
!52 = !{!14, !10, i64 248}
!53 = !{!26, !10, i64 120}
!54 = !{!14, !15, i64 320}
!55 = !{!14, !15, i64 324}
!56 = !{!14, !15, i64 328}
!57 = !{!14, !15, i64 332}
!58 = !{!26, !10, i64 128}
!59 = !{!14, !10, i64 264}
!60 = !{!14, !15, i64 336}
!61 = !{!14, !15, i64 340}
!62 = !{!14, !15, i64 12}
!63 = !{!14, !10, i64 32}
!64 = !{!14, !12, i64 600}
!65 = !{!14, !6, i64 400}
!66 = !{!14, !12, i64 602}


; 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.1 = external hidden unnamed_addr constant [156 x i8], align 1
@.str.2 = external hidden unnamed_addr constant [13 x i8], align 1
@__func__.thermodynamics_at_z = external hidden unnamed_addr constant [20 x i8], align 1
@.str.3 = external hidden unnamed_addr constant [189 x i8], align 1
@.str.4 = external hidden unnamed_addr constant [205 x i8], align 1

; Function Attrs: nounwind uwtable
define dso_local range(i32 0, 2) i32 @thermodynamics_at_z(ptr nocapture noundef readonly %0, ptr noundef %1, double noundef %2, i16 noundef signext %3, ptr noundef %4, ptr nocapture noundef readonly %5, ptr noundef %6) local_unnamed_addr #0 {
  %8 = alloca [2048 x i8], align 16
  %9 = alloca [2048 x i8], align 16
  %10 = alloca [2048 x i8], align 16
  %11 = getelementptr inbounds i8, ptr %1, i64 384
  %12 = load ptr, ptr %11, align 8, !tbaa !5
  %13 = getelementptr inbounds i8, ptr %1, i64 376
  %14 = load i32, ptr %13, align 8, !tbaa !13
  %15 = add nsw i32 %14, -1
  %16 = sext i32 %15 to i64
  %17 = getelementptr inbounds double, ptr %12, i64 %16
  %18 = load double, ptr %17, align 8, !tbaa !14
  %19 = fcmp ugt double %18, %2
  br i1 %19, label %326, label %20

20:                                               ; preds = %7
  %21 = getelementptr inbounds i8, ptr %1, i64 392
  %22 = load ptr, ptr %21, align 8, !tbaa !15
  %23 = getelementptr inbounds i8, ptr %1, i64 372
  %24 = load i32, ptr %23, align 4, !tbaa !16
  %25 = mul nsw i32 %24, %15
  %26 = getelementptr inbounds i8, ptr %1, i64 272
  %27 = load i32, ptr %26, align 8, !tbaa !17
  %28 = add nsw i32 %25, %27
  %29 = sext i32 %28 to i64
  %30 = getelementptr inbounds double, ptr %22, i64 %29
  %31 = load double, ptr %30, align 8, !tbaa !14
  %32 = sext i32 %27 to i64
  %33 = getelementptr inbounds double, ptr %6, i64 %32
  store double %31, ptr %33, align 8, !tbaa !14
  %34 = fadd double %2, 1.000000e+00
  %35 = fmul double %34, %34
  %36 = getelementptr inbounds i8, ptr %1, i64 592
  %37 = load double, ptr %36, align 8, !tbaa !18
  %38 = fmul double %35, %37
  %39 = fmul double %31, %38
  %40 = fmul double %39, 0x3A15151E37FBF236
  %41 = fmul double %40, 0x449A22FF260E3826
  %42 = getelementptr inbounds i8, ptr %1, i64 276
  %43 = load i32, ptr %42, align 4, !tbaa !19
  %44 = sext i32 %43 to i64
  %45 = getelementptr inbounds double, ptr %6, i64 %44
  store double %41, ptr %45, align 8, !tbaa !14
  %46 = getelementptr inbounds i8, ptr %1, i64 280
  %47 = load i32, ptr %46, align 8, !tbaa !20
  %48 = add nsw i32 %47, %25
  %49 = sext i32 %48 to i64
  %50 = getelementptr inbounds double, ptr %22, i64 %49
  %51 = load double, ptr %50, align 8, !tbaa !14
  %52 = load double, ptr %17, align 8, !tbaa !14
  %53 = fadd double %52, 1.000000e+00
  %54 = fdiv double %34, %53
  %square = fmul double %54, %54
  %55 = fmul double %51, %square
  %56 = sext i32 %47 to i64
  %57 = getelementptr inbounds double, ptr %6, i64 %56
  store double %55, ptr %57, align 8, !tbaa !14
  %58 = getelementptr inbounds i8, ptr %1, i64 42
  %59 = load i16, ptr %58, align 2, !tbaa !21
  %60 = icmp eq i16 %59, 1
  br i1 %60, label %61, label %76

61:                                               ; preds = %20
  %62 = getelementptr inbounds i8, ptr %1, i64 368
  %63 = load i32, ptr %62, align 8, !tbaa !22
  %64 = add nsw i32 %63, %25
  %65 = sext i32 %64 to i64
  %66 = getelementptr inbounds double, ptr %22, i64 %65
  %67 = load double, ptr %66, align 8, !tbaa !14
  %68 = load double, ptr %17, align 8, !tbaa !14
  %69 = fadd double %68, 1.000000e+00
  %70 = fdiv double %34, %69
  %71 = tail call double @pow(double noundef %70, double noundef -1.500000e+00) #5
  %72 = fmul double %67, %71
  %73 = load i32, ptr %62, align 8, !tbaa !22
  %74 = sext i32 %73 to i64
  %75 = getelementptr inbounds double, ptr %6, i64 %74
  store double %72, ptr %75, align 8, !tbaa !14
  %.pre = load i32, ptr %42, align 4, !tbaa !19
  %.pre240 = sext i32 %.pre to i64
  br label %76

76:                                               ; preds = %61, %20
  %.pre-phi = phi i64 [ %.pre240, %61 ], [ %44, %20 ]
  %77 = getelementptr inbounds i8, ptr %0, i64 468
  %78 = load i32, ptr %77, align 4, !tbaa !23
  %79 = sext i32 %78 to i64
  %80 = getelementptr inbounds double, ptr %5, i64 %79
  %81 = load double, ptr %80, align 8, !tbaa !14
  %82 = fmul double %81, -2.000000e+00
  %83 = fdiv double %82, %34
  %84 = getelementptr inbounds double, ptr %6, i64 %.pre-phi
  %85 = load double, ptr %84, align 8, !tbaa !14
  %86 = fmul double %83, %85
  %87 = getelementptr inbounds i8, ptr %1, i64 284
  %88 = load i32, ptr %87, align 4, !tbaa !25
  %89 = sext i32 %88 to i64
  %90 = getelementptr inbounds double, ptr %6, i64 %89
  store double %86, ptr %90, align 8, !tbaa !14
  %91 = load double, ptr %80, align 8, !tbaa !14
  %92 = fmul double %91, %91
  %93 = fdiv double %92, %34
  %94 = getelementptr inbounds i8, ptr %0, i64 472
  %95 = load i32, ptr %94, align 8, !tbaa !26
  %96 = sext i32 %95 to i64
  %97 = getelementptr inbounds double, ptr %5, i64 %96
  %98 = load double, ptr %97, align 8, !tbaa !14
  %99 = fsub double %93, %98
  %100 = fmul double %99, 2.000000e+00
  %101 = fdiv double %100, %34
  %102 = load double, ptr %84, align 8, !tbaa !14
  %103 = fmul double %102, %101
  %104 = getelementptr inbounds i8, ptr %1, i64 288
  %105 = load i32, ptr %104, align 8, !tbaa !27
  %106 = sext i32 %105 to i64
  %107 = getelementptr inbounds double, ptr %6, i64 %106
  store double %103, ptr %107, align 8, !tbaa !14
  %108 = getelementptr inbounds i8, ptr %1, i64 292
  %109 = load i32, ptr %108, align 4, !tbaa !28
  %110 = sext i32 %109 to i64
  %111 = getelementptr inbounds double, ptr %6, i64 %110
  store double 0.000000e+00, ptr %111, align 8, !tbaa !14
  %112 = getelementptr inbounds i8, ptr %1, i64 296
  %113 = load i32, ptr %112, align 8, !tbaa !29
  %114 = sext i32 %113 to i64
  %115 = getelementptr inbounds double, ptr %6, i64 %114
  store double 0.000000e+00, ptr %115, align 8, !tbaa !14
  %116 = getelementptr inbounds i8, ptr %1, i64 300
  %117 = load i32, ptr %116, align 4, !tbaa !30
  %118 = sext i32 %117 to i64
  %119 = getelementptr inbounds double, ptr %6, i64 %118
  store double 0.000000e+00, ptr %119, align 8, !tbaa !14
  %120 = getelementptr inbounds i8, ptr %1, i64 304
  %121 = load i32, ptr %120, align 8, !tbaa !31
  %122 = sext i32 %121 to i64
  %123 = getelementptr inbounds double, ptr %6, i64 %122
  store double 0.000000e+00, ptr %123, align 8, !tbaa !14
  %124 = getelementptr inbounds i8, ptr %0, i64 16
  %125 = load double, ptr %124, align 8, !tbaa !32
  %126 = fmul double %34, %125
  %127 = getelementptr inbounds i8, ptr %1, i64 344
  %128 = load i32, ptr %127, align 8, !tbaa !33
  %129 = sext i32 %128 to i64
  %130 = getelementptr inbounds double, ptr %6, i64 %129
  store double %126, ptr %130, align 8, !tbaa !14
  %131 = load double, ptr %1, align 8, !tbaa !34
  %132 = tail call double @llvm.fmuladd.f64(double %131, double 0xBFE7F14DA42541E8, double 1.000000e+00)
  %133 = fsub double 1.000000e+00, %131
  %134 = tail call double @llvm.fmuladd.f64(double %31, double %133, double %132)
  %135 = fmul double %134, 0x3D39D63187BB40C4
  %136 = load double, ptr %124, align 8, !tbaa !32
  %137 = fmul double %136, %135
  %138 = fmul double %34, %137
  %139 = getelementptr inbounds i8, ptr %1, i64 348
  %140 = load i32, ptr %139, align 4, !tbaa !35
  %141 = sext i32 %140 to i64
  %142 = getelementptr inbounds double, ptr %6, i64 %141
  store double %138, ptr %142, align 8, !tbaa !14
  %143 = fmul double %138, 4.000000e+00
  %144 = fdiv double %143, 3.000000e+00
  %145 = getelementptr inbounds i8, ptr %1, i64 352
  %146 = load i32, ptr %145, align 8, !tbaa !36
  %147 = sext i32 %146 to i64
  %148 = getelementptr inbounds double, ptr %6, i64 %147
  store double %144, ptr %148, align 8, !tbaa !14
  %149 = getelementptr inbounds i8, ptr %1, i64 40
  %150 = load i16, ptr %149, align 8, !tbaa !37
  %151 = icmp eq i16 %150, 1
  br i1 %151, label %152, label %176

152:                                              ; preds = %76
  %153 = load double, ptr %80, align 8, !tbaa !14
  %154 = fneg double %153
  %155 = getelementptr inbounds i8, ptr %0, i64 464
  %156 = load i32, ptr %155, align 8, !tbaa !38
  %157 = sext i32 %156 to i64
  %158 = getelementptr inbounds double, ptr %5, i64 %157
  %159 = load double, ptr %158, align 8, !tbaa !14
  %160 = fmul double %159, %154
  %161 = fmul double %144, %160
  %162 = getelementptr inbounds i8, ptr %1, i64 356
  %163 = load i32, ptr %162, align 4, !tbaa !39
  %164 = sext i32 %163 to i64
  %165 = getelementptr inbounds double, ptr %6, i64 %164
  store double %161, ptr %165, align 8, !tbaa !14
  %166 = load double, ptr %97, align 8, !tbaa !14
  %167 = fneg double %166
  %168 = load double, ptr %158, align 8, !tbaa !14
  %169 = fmul double %168, %167
  %170 = load double, ptr %148, align 8, !tbaa !14
  %171 = fmul double %169, %170
  %172 = getelementptr inbounds i8, ptr %1, i64 360
  %173 = load i32, ptr %172, align 8, !tbaa !40
  %174 = sext i32 %173 to i64
  %175 = getelementptr inbounds double, ptr %6, i64 %174
  store double %171, ptr %175, align 8, !tbaa !14
  br label %176

176:                                              ; preds = %152, %76
  %177 = load double, ptr %84, align 8, !tbaa !14
  %178 = getelementptr inbounds i8, ptr %1, i64 364
  %179 = load i32, ptr %178, align 4, !tbaa !41
  %180 = sext i32 %179 to i64
  %181 = getelementptr inbounds double, ptr %6, i64 %180
  store double %177, ptr %181, align 8, !tbaa !14
  %182 = getelementptr inbounds i8, ptr %0, i64 742
  %183 = load i16, ptr %182, align 2, !tbaa !42
  %184 = icmp eq i16 %183, 1
  br i1 %184, label %185, label %378

185:                                              ; preds = %176
  %186 = getelementptr inbounds i8, ptr %1, i64 240
  %187 = load double, ptr %186, align 8, !tbaa !43
  %188 = fdiv double %34, 1.000000e+07
  %189 = getelementptr inbounds i8, ptr %1, i64 256
  %190 = load double, ptr %189, align 8, !tbaa !44
  %191 = tail call double @pow(double noundef %188, double noundef %190) #5
  %192 = fmul double %187, %191
  %193 = getelementptr inbounds i8, ptr %0, i64 136
  %194 = load double, ptr %193, align 8, !tbaa !45
  %195 = fmul double %192, %194
  %196 = getelementptr inbounds i8, ptr %0, i64 328
  %197 = load double, ptr %196, align 8, !tbaa !46
  %square237 = fmul double %197, %197
  %198 = fmul double %195, %square237
  %199 = getelementptr inbounds i8, ptr %1, i64 308
  %200 = load i32, ptr %199, align 4, !tbaa !47
  %201 = sext i32 %200 to i64
  %202 = getelementptr inbounds double, ptr %6, i64 %201
  store double %198, ptr %202, align 8, !tbaa !14
  %203 = load i32, ptr %77, align 4, !tbaa !23
  %204 = sext i32 %203 to i64
  %205 = getelementptr inbounds double, ptr %5, i64 %204
  %206 = load double, ptr %205, align 8, !tbaa !14
  %207 = fneg double %206
  %208 = load double, ptr %189, align 8, !tbaa !44
  %209 = fmul double %208, %207
  %210 = fdiv double %209, %34
  %211 = fmul double %198, %210
  %212 = getelementptr inbounds i8, ptr %1, i64 312
  %213 = load i32, ptr %212, align 8, !tbaa !48
  %214 = sext i32 %213 to i64
  %215 = getelementptr inbounds double, ptr %6, i64 %214
  store double %211, ptr %215, align 8, !tbaa !14
  %216 = load double, ptr %205, align 8, !tbaa !14
  %217 = fmul double %216, %216
  %218 = fdiv double %217, %34
  %219 = load i32, ptr %94, align 8, !tbaa !26
  %220 = sext i32 %219 to i64
  %221 = getelementptr inbounds double, ptr %5, i64 %220
  %222 = load double, ptr %221, align 8, !tbaa !14
  %223 = fsub double %218, %222
  %224 = load double, ptr %189, align 8, !tbaa !44
  %225 = fmul double %223, %224
  %226 = fdiv double %225, %34
  %227 = load double, ptr %202, align 8, !tbaa !14
  %228 = fmul double %227, %226
  %229 = getelementptr inbounds i8, ptr %1, i64 316
  %230 = load i32, ptr %229, align 4, !tbaa !49
  %231 = sext i32 %230 to i64
  %232 = getelementptr inbounds double, ptr %6, i64 %231
  store double %228, ptr %232, align 8, !tbaa !14
  %233 = getelementptr inbounds i8, ptr %1, i64 248
  %234 = load double, ptr %233, align 8, !tbaa !50
  %235 = load double, ptr %189, align 8, !tbaa !44
  %236 = tail call double @pow(double noundef %188, double noundef %235) #5
  %237 = fmul double %234, %236
  %238 = getelementptr inbounds i8, ptr %0, i64 120
  %239 = load double, ptr %238, align 8, !tbaa !51
  %240 = fmul double %237, %239
  %241 = load double, ptr %196, align 8, !tbaa !46
  %square238 = fmul double %241, %241
  %242 = fmul double %240, %square238
  %243 = getelementptr inbounds i8, ptr %1, i64 320
  %244 = load i32, ptr %243, align 8, !tbaa !52
  %245 = sext i32 %244 to i64
  %246 = getelementptr inbounds double, ptr %6, i64 %245
  store double %242, ptr %246, align 8, !tbaa !14
  %247 = load ptr, ptr %21, align 8, !tbaa !15
  %248 = load i32, ptr %13, align 8, !tbaa !13
  %249 = add nsw i32 %248, -1
  %250 = load i32, ptr %23, align 4, !tbaa !16
  %251 = mul nsw i32 %249, %250
  %252 = getelementptr inbounds i8, ptr %1, i64 324
  %253 = load i32, ptr %252, align 4, !tbaa !53
  %254 = add nsw i32 %251, %253
  %255 = sext i32 %254 to i64
  %256 = getelementptr inbounds double, ptr %247, i64 %255
  %257 = load double, ptr %256, align 8, !tbaa !14
  %258 = add nsw i32 %248, -2
  %259 = mul nsw i32 %258, %250
  %260 = add nsw i32 %259, %253
  %261 = sext i32 %260 to i64
  %262 = getelementptr inbounds double, ptr %247, i64 %261
  %263 = load double, ptr %262, align 8, !tbaa !14
  %264 = fsub double %257, %263
  %265 = load ptr, ptr %11, align 8, !tbaa !5
  %266 = sext i32 %249 to i64
  %267 = getelementptr inbounds double, ptr %265, i64 %266
  %268 = load double, ptr %267, align 8, !tbaa !14
  %269 = fsub double %2, %268
  %270 = fmul double %264, %269
  %271 = sext i32 %258 to i64
  %272 = getelementptr inbounds double, ptr %265, i64 %271
  %273 = load double, ptr %272, align 8, !tbaa !14
  %274 = fsub double %268, %273
  %275 = fdiv double %270, %274
  %276 = fadd double %257, %275
  %277 = sext i32 %253 to i64
  %278 = getelementptr inbounds double, ptr %6, i64 %277
  store double %276, ptr %278, align 8, !tbaa !14
  %279 = getelementptr inbounds i8, ptr %1, i64 328
  %280 = load i32, ptr %279, align 8, !tbaa !54
  %281 = add nsw i32 %280, %251
  %282 = sext i32 %281 to i64
  %283 = getelementptr inbounds double, ptr %247, i64 %282
  %284 = load double, ptr %283, align 8, !tbaa !14
  %285 = add nsw i32 %280, %259
  %286 = sext i32 %285 to i64
  %287 = getelementptr inbounds double, ptr %247, i64 %286
  %288 = load double, ptr %287, align 8, !tbaa !14
  %289 = fsub double %284, %288
  %290 = load double, ptr %267, align 8, !tbaa !14
  %291 = fsub double %2, %290
  %292 = fmul double %289, %291
  %293 = load double, ptr %272, align 8, !tbaa !14
  %294 = fsub double %290, %293
  %295 = fdiv double %292, %294
  %296 = fadd double %284, %295
  %297 = sext i32 %280 to i64
  %298 = getelementptr inbounds double, ptr %6, i64 %297
  store double %296, ptr %298, align 8, !tbaa !14
  %299 = getelementptr inbounds i8, ptr %1, i64 332
  %300 = load i32, ptr %299, align 4, !tbaa !55
  %301 = add nsw i32 %300, %251
  %302 = sext i32 %301 to i64
  %303 = getelementptr inbounds double, ptr %247, i64 %302
  %304 = load double, ptr %303, align 8, !tbaa !14
  %305 = sext i32 %300 to i64
  %306 = getelementptr inbounds double, ptr %6, i64 %305
  store double %304, ptr %306, align 8, !tbaa !14
  %307 = getelementptr inbounds i8, ptr %0, i64 128
  %308 = load double, ptr %307, align 8, !tbaa !56
  %309 = fmul double %308, 0x3B50B0E7F1528C6D
  %310 = fmul double %34, %309
  %311 = fdiv double %310, 0x3C07A4DA04A6DAF5
  %312 = fdiv double %311, 3.000000e+00
  %313 = getelementptr inbounds i8, ptr %1, i64 264
  %314 = load double, ptr %313, align 8, !tbaa !57
  %315 = fdiv double %312, %314
  %316 = getelementptr inbounds i8, ptr %1, i64 336
  %317 = load i32, ptr %316, align 8, !tbaa !58
  %318 = sext i32 %317 to i64
  %319 = getelementptr inbounds double, ptr %6, i64 %318
  store double %315, ptr %319, align 8, !tbaa !14
  %320 = load double, ptr %307, align 8, !tbaa !56
  %321 = fmul double %34, %320
  %322 = getelementptr inbounds i8, ptr %1, i64 340
  %323 = load i32, ptr %322, align 4, !tbaa !59
  %324 = sext i32 %323 to i64
  %325 = getelementptr inbounds double, ptr %6, i64 %324
  store double %321, ptr %325, align 8, !tbaa !14
  br label %378

326:                                              ; preds = %7
  %327 = getelementptr inbounds i8, ptr %1, i64 12
  %328 = load i32, ptr %327, align 4, !tbaa !60
  %329 = icmp eq i32 %328, 3
  br i1 %329, label %330, label %335

330:                                              ; preds = %326
  %331 = getelementptr inbounds i8, ptr %1, i64 32
  %332 = load double, ptr %331, align 8, !tbaa !61
  %333 = fmul double %332, 2.000000e+00
  %334 = fcmp ogt double %333, %2
  br i1 %334, label %338, label %.thread

335:                                              ; preds = %326
  %336 = icmp eq i32 %328, 5
  %337 = fcmp olt double %2, 5.000000e+01
  %or.cond = and i1 %337, %336
  br i1 %or.cond, label %338, label %.thread

338:                                              ; preds = %335, %330
  %339 = getelementptr inbounds i8, ptr %1, i64 392
  %340 = load ptr, ptr %339, align 8, !tbaa !15
  %341 = getelementptr inbounds i8, ptr %1, i64 372
  %342 = load i32, ptr %341, align 4, !tbaa !16
  %343 = getelementptr inbounds i8, ptr %1, i64 606
  %344 = tail call i32 @array_interpolate_linear(ptr noundef nonnull %12, i32 noundef %14, ptr noundef %340, i32 noundef %342, double noundef %2, ptr noundef %4, ptr noundef %6, i32 noundef %342, ptr noundef nonnull %343) #5
  %345 = icmp eq i32 %344, 1
  br i1 %345, label %346, label %378

346:                                              ; preds = %338
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %8) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %8, ptr noundef nonnull @.str, ptr noundef nonnull @.str.1, ptr noundef nonnull %343) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %343, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_at_z, i32 noundef 233, ptr noundef nonnull %8) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %8) #5
  br label %378

.thread:                                          ; preds = %330, %335
  %347 = getelementptr inbounds i8, ptr %1, i64 600
  %348 = load i16, ptr %347, align 8, !tbaa !62
  %349 = icmp eq i16 %348, %3
  br i1 %349, label %350, label %361

350:                                              ; preds = %.thread
  %351 = getelementptr inbounds i8, ptr %1, i64 392
  %352 = load ptr, ptr %351, align 8, !tbaa !15
  %353 = getelementptr inbounds i8, ptr %1, i64 400
  %354 = load ptr, ptr %353, align 8, !tbaa !63
  %355 = getelementptr inbounds i8, ptr %1, i64 372
  %356 = load i32, ptr %355, align 4, !tbaa !16
  %357 = getelementptr inbounds i8, ptr %1, i64 606
  %358 = tail call i32 @array_interpolate_spline(ptr noundef nonnull %12, i32 noundef %14, ptr noundef %352, ptr noundef %354, i32 noundef %356, double noundef %2, ptr noundef %4, ptr noundef %6, i32 noundef %356, ptr noundef nonnull %357) #5
  %359 = icmp eq i32 %358, 1
  br i1 %359, label %360, label %361

360:                                              ; preds = %350
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %9) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %9, ptr noundef nonnull @.str, ptr noundef nonnull @.str.3, ptr noundef nonnull %357) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %357, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_at_z, i32 noundef 253, ptr noundef nonnull %9) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %9) #5
  br label %378

361:                                              ; preds = %350, %.thread
  %362 = getelementptr inbounds i8, ptr %1, i64 602
  %363 = load i16, ptr %362, align 2, !tbaa !64
  %364 = icmp eq i16 %363, %3
  br i1 %364, label %365, label %378

365:                                              ; preds = %361
  %366 = load ptr, ptr %11, align 8, !tbaa !5
  %367 = load i32, ptr %13, align 8, !tbaa !13
  %368 = getelementptr inbounds i8, ptr %1, i64 392
  %369 = load ptr, ptr %368, align 8, !tbaa !15
  %370 = getelementptr inbounds i8, ptr %1, i64 400
  %371 = load ptr, ptr %370, align 8, !tbaa !63
  %372 = getelementptr inbounds i8, ptr %1, i64 372
  %373 = load i32, ptr %372, align 4, !tbaa !16
  %374 = getelementptr inbounds i8, ptr %1, i64 606
  %375 = tail call i32 @array_interpolate_spline_growing_closeby(ptr noundef %366, i32 noundef %367, ptr noundef %369, ptr noundef %371, i32 noundef %373, double noundef %2, ptr noundef %4, ptr noundef %6, i32 noundef %373, ptr noundef nonnull %374) #5
  %376 = icmp eq i32 %375, 1
  br i1 %376, label %377, label %378

377:                                              ; preds = %365
  call void @llvm.lifetime.start.p0(i64 2048, ptr nonnull %10) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %10, ptr noundef nonnull @.str, ptr noundef nonnull @.str.4, ptr noundef nonnull %374) #5
  call void (ptr, ptr, ...) @class_protect_sprintf(ptr noundef nonnull %374, ptr noundef nonnull @.str.2, ptr noundef nonnull @__func__.thermodynamics_at_z, i32 noundef 270, ptr noundef nonnull %10) #5
  call void @llvm.lifetime.end.p0(i64 2048, ptr nonnull %10) #5
  br label %378

378:                                              ; preds = %185, %176, %361, %365, %338, %377, %360, %346
  %.0 = phi i32 [ 1, %346 ], [ 1, %360 ], [ 1, %377 ], [ 0, %338 ], [ 0, %365 ], [ 0, %361 ], [ 0, %176 ], [ 0, %185 ]
  ret i32 %.0
}

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

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

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

declare i32 @array_interpolate_linear(ptr noundef, i32 noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) local_unnamed_addr #4

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

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

declare i32 @array_interpolate_spline(ptr noundef, i32 noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) local_unnamed_addr #4

declare i32 @array_interpolate_spline_growing_closeby(ptr noundef, i32 noundef, ptr noundef, ptr noundef, i32 noundef, double noundef, ptr noundef, ptr noundef, i32 noundef, ptr noundef) local_unnamed_addr #4

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

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

!0 = !{i32 1, !"wchar_size", i32 4}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{i32 7, !"uwtable", i32 2}
!4 = !{!"clang version 19.1.0 (git@github.com:llvm/llvm-project.git a4bf6cd7cfb1a1421ba92bca9d017b49936c55e4)"}
!5 = !{!6, !12, i64 384}
!6 = !{!"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}
!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 = !{!6, !10, i64 376}
!14 = !{!7, !7, i64 0}
!15 = !{!6, !12, i64 392}
!16 = !{!6, !10, i64 372}
!17 = !{!6, !10, i64 272}
!18 = !{!6, !7, i64 592}
!19 = !{!6, !10, i64 276}
!20 = !{!6, !10, i64 280}
!21 = !{!6, !11, i64 42}
!22 = !{!6, !10, i64 368}
!23 = !{!24, !10, i64 468}
!24 = !{!"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}
!25 = !{!6, !10, i64 284}
!26 = !{!24, !10, i64 472}
!27 = !{!6, !10, i64 288}
!28 = !{!6, !10, i64 292}
!29 = !{!6, !10, i64 296}
!30 = !{!6, !10, i64 300}
!31 = !{!6, !10, i64 304}
!32 = !{!24, !7, i64 16}
!33 = !{!6, !10, i64 344}
!34 = !{!6, !7, i64 0}
!35 = !{!6, !10, i64 348}
!36 = !{!6, !10, i64 352}
!37 = !{!6, !11, i64 40}
!38 = !{!24, !10, i64 464}
!39 = !{!6, !10, i64 356}
!40 = !{!6, !10, i64 360}
!41 = !{!6, !10, i64 364}
!42 = !{!24, !11, i64 742}
!43 = !{!6, !7, i64 240}
!44 = !{!6, !7, i64 256}
!45 = !{!24, !7, i64 136}
!46 = !{!24, !7, i64 328}
!47 = !{!6, !10, i64 308}
!48 = !{!6, !10, i64 312}
!49 = !{!6, !10, i64 316}
!50 = !{!6, !7, i64 248}
!51 = !{!24, !7, i64 120}
!52 = !{!6, !10, i64 320}
!53 = !{!6, !10, i64 324}
!54 = !{!6, !10, i64 328}
!55 = !{!6, !10, i64 332}
!56 = !{!24, !7, i64 128}
!57 = !{!6, !7, i64 264}
!58 = !{!6, !10, i64 336}
!59 = !{!6, !10, i64 340}
!60 = !{!6, !10, i64 12}
!61 = !{!6, !7, i64 32}
!62 = !{!6, !11, i64 600}
!63 = !{!6, !12, i64 400}
!64 = !{!6, !11, i64 602}


; Effective O3 passes: BDCEPass,CorrelatedValuePropagationPass,EarlyCSEPass,GVNPass,GlobalOptPass,IPSCCPPass,InferFunctionAttrsPass,InstCombinePass,JumpThreadingPass,PostOrderFunctionAttrsPass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass