; LLVM IR from zhenrong-wang/co2 :: co2_property_calc_origeqns (phird)

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

@coeffr1 = external dso_local constant [7 x [3 x double]], align 16
@coeffr2 = external dso_local constant [27 x [4 x double]], align 16
@coeffr3 = external dso_local constant [5 x [7 x double]], align 16
@coeffr4 = external dso_local constant [3 x [8 x double]], align 16

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

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

; Function Attrs: nounwind
declare double @exp(double 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 uwtable
declare dso_local double @deltd(double noundef, double noundef, double noundef, double noundef, double noundef, double noundef, double noundef) #3

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

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

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

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

; Function Attrs: nounwind uwtable
define dso_local double @phird(double noundef %0, double noundef %1) #3 {
  %3 = alloca double, align 8
  %4 = alloca double, align 8
  %5 = alloca i32, align 4
  %6 = alloca double, align 8
  %7 = alloca double, align 8
  %8 = alloca double, align 8
  %9 = alloca double, align 8
  %10 = alloca double, align 8
  %11 = alloca double, align 8
  %12 = alloca double, align 8
  %13 = alloca double, align 8
  %14 = alloca double, align 8
  %15 = alloca double, align 8
  %16 = alloca double, align 8
  %17 = alloca double, align 8
  %18 = alloca double, align 8
  %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
  store double %0, ptr %3, align 8, !tbaa !5
  store double %1, ptr %4, align 8, !tbaa !5
  call void @llvm.lifetime.start.p0(i64 4, ptr %5) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %6) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %7) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %8) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %9) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %10) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %11) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %12) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %13) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %14) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %15) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %16) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %17) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %18) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %19) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %20) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %21) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %22) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %23) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %24) #4
  call void @llvm.lifetime.start.p0(i64 8, ptr %25) #4
  store double 0.000000e+00, ptr %25, align 8, !tbaa !5
  store i32 0, ptr %5, align 4, !tbaa !9
  br label %26

26:                                               ; preds = %58, %2
  %27 = load i32, ptr %5, align 4, !tbaa !9
  %28 = icmp slt i32 %27, 7
  br i1 %28, label %29, label %61

29:                                               ; preds = %26
  %30 = load i32, ptr %5, align 4, !tbaa !9
  %31 = sext i32 %30 to i64
  %32 = getelementptr inbounds [7 x [3 x double]], ptr @coeffr1, i64 0, i64 %31
  %33 = getelementptr inbounds [3 x double], ptr %32, i64 0, i64 0
  %34 = load double, ptr %33, align 8, !tbaa !5
  store double %34, ptr %6, align 8, !tbaa !5
  %35 = load i32, ptr %5, align 4, !tbaa !9
  %36 = sext i32 %35 to i64
  %37 = getelementptr inbounds [7 x [3 x double]], ptr @coeffr1, i64 0, i64 %36
  %38 = getelementptr inbounds [3 x double], ptr %37, i64 0, i64 1
  %39 = load double, ptr %38, align 8, !tbaa !5
  store double %39, ptr %7, align 8, !tbaa !5
  %40 = load i32, ptr %5, align 4, !tbaa !9
  %41 = sext i32 %40 to i64
  %42 = getelementptr inbounds [7 x [3 x double]], ptr @coeffr1, i64 0, i64 %41
  %43 = getelementptr inbounds [3 x double], ptr %42, i64 0, i64 2
  %44 = load double, ptr %43, align 8, !tbaa !5
  store double %44, ptr %8, align 8, !tbaa !5
  %45 = load double, ptr %6, align 8, !tbaa !5
  %46 = load double, ptr %7, align 8, !tbaa !5
  %47 = fmul double %45, %46
  %48 = load double, ptr %3, align 8, !tbaa !5
  %49 = load double, ptr %7, align 8, !tbaa !5
  %50 = fsub double %49, 1.000000e+00
  %51 = call double @pow(double noundef %48, double noundef %50) #4
  %52 = fmul double %47, %51
  %53 = load double, ptr %4, align 8, !tbaa !5
  %54 = load double, ptr %8, align 8, !tbaa !5
  %55 = call double @pow(double noundef %53, double noundef %54) #4
  %56 = load double, ptr %25, align 8, !tbaa !5
  %57 = call double @llvm.fmuladd.f64(double %52, double %55, double %56)
  store double %57, ptr %25, align 8, !tbaa !5
  br label %58

58:                                               ; preds = %29
  %59 = load i32, ptr %5, align 4, !tbaa !9
  %60 = add nsw i32 %59, 1
  store i32 %60, ptr %5, align 4, !tbaa !9
  br label %26, !llvm.loop !11

61:                                               ; preds = %26
  store i32 0, ptr %5, align 4, !tbaa !9
  br label %62

62:                                               ; preds = %111, %61
  %63 = load i32, ptr %5, align 4, !tbaa !9
  %64 = icmp slt i32 %63, 27
  br i1 %64, label %65, label %114

65:                                               ; preds = %62
  %66 = load i32, ptr %5, align 4, !tbaa !9
  %67 = sext i32 %66 to i64
  %68 = getelementptr inbounds [27 x [4 x double]], ptr @coeffr2, i64 0, i64 %67
  %69 = getelementptr inbounds [4 x double], ptr %68, i64 0, i64 0
  %70 = load double, ptr %69, align 16, !tbaa !5
  store double %70, ptr %6, align 8, !tbaa !5
  %71 = load i32, ptr %5, align 4, !tbaa !9
  %72 = sext i32 %71 to i64
  %73 = getelementptr inbounds [27 x [4 x double]], ptr @coeffr2, i64 0, i64 %72
  %74 = getelementptr inbounds [4 x double], ptr %73, i64 0, i64 1
  %75 = load double, ptr %74, align 8, !tbaa !5
  store double %75, ptr %7, align 8, !tbaa !5
  %76 = load i32, ptr %5, align 4, !tbaa !9
  %77 = sext i32 %76 to i64
  %78 = getelementptr inbounds [27 x [4 x double]], ptr @coeffr2, i64 0, i64 %77
  %79 = getelementptr inbounds [4 x double], ptr %78, i64 0, i64 2
  %80 = load double, ptr %79, align 16, !tbaa !5
  store double %80, ptr %8, align 8, !tbaa !5
  %81 = load i32, ptr %5, align 4, !tbaa !9
  %82 = sext i32 %81 to i64
  %83 = getelementptr inbounds [27 x [4 x double]], ptr @coeffr2, i64 0, i64 %82
  %84 = getelementptr inbounds [4 x double], ptr %83, i64 0, i64 3
  %85 = load double, ptr %84, align 8, !tbaa !5
  store double %85, ptr %9, align 8, !tbaa !5
  %86 = load double, ptr %6, align 8, !tbaa !5
  %87 = load double, ptr %4, align 8, !tbaa !5
  %88 = load double, ptr %8, align 8, !tbaa !5
  %89 = call double @pow(double noundef %87, double noundef %88) #4
  %90 = fmul double %86, %89
  %91 = load double, ptr %3, align 8, !tbaa !5
  %92 = load double, ptr %9, align 8, !tbaa !5
  %93 = call double @pow(double noundef %91, double noundef %92) #4
  %94 = fneg double %93
  %95 = call double @exp(double noundef %94) #4
  %96 = fmul double %90, %95
  %97 = load double, ptr %3, align 8, !tbaa !5
  %98 = load double, ptr %7, align 8, !tbaa !5
  %99 = fsub double %98, 1.000000e+00
  %100 = call double @pow(double noundef %97, double noundef %99) #4
  %101 = fmul double %96, %100
  %102 = load double, ptr %7, align 8, !tbaa !5
  %103 = load double, ptr %9, align 8, !tbaa !5
  %104 = load double, ptr %3, align 8, !tbaa !5
  %105 = load double, ptr %9, align 8, !tbaa !5
  %106 = call double @pow(double noundef %104, double noundef %105) #4
  %107 = fneg double %103
  %108 = call double @llvm.fmuladd.f64(double %107, double %106, double %102)
  %109 = load double, ptr %25, align 8, !tbaa !5
  %110 = call double @llvm.fmuladd.f64(double %101, double %108, double %109)
  store double %110, ptr %25, align 8, !tbaa !5
  br label %111

111:                                              ; preds = %65
  %112 = load i32, ptr %5, align 4, !tbaa !9
  %113 = add nsw i32 %112, 1
  store i32 %113, ptr %5, align 4, !tbaa !9
  br label %62, !llvm.loop !13

114:                                              ; preds = %62
  store i32 0, ptr %5, align 4, !tbaa !9
  br label %115

115:                                              ; preds = %197, %114
  %116 = load i32, ptr %5, align 4, !tbaa !9
  %117 = icmp slt i32 %116, 5
  br i1 %117, label %118, label %200

118:                                              ; preds = %115
  %119 = load i32, ptr %5, align 4, !tbaa !9
  %120 = sext i32 %119 to i64
  %121 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %120
  %122 = getelementptr inbounds [7 x double], ptr %121, i64 0, i64 0
  %123 = load double, ptr %122, align 8, !tbaa !5
  store double %123, ptr %6, align 8, !tbaa !5
  %124 = load i32, ptr %5, align 4, !tbaa !9
  %125 = sext i32 %124 to i64
  %126 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %125
  %127 = getelementptr inbounds [7 x double], ptr %126, i64 0, i64 1
  %128 = load double, ptr %127, align 8, !tbaa !5
  store double %128, ptr %7, align 8, !tbaa !5
  %129 = load i32, ptr %5, align 4, !tbaa !9
  %130 = sext i32 %129 to i64
  %131 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %130
  %132 = getelementptr inbounds [7 x double], ptr %131, i64 0, i64 2
  %133 = load double, ptr %132, align 8, !tbaa !5
  store double %133, ptr %8, align 8, !tbaa !5
  %134 = load i32, ptr %5, align 4, !tbaa !9
  %135 = sext i32 %134 to i64
  %136 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %135
  %137 = getelementptr inbounds [7 x double], ptr %136, i64 0, i64 3
  %138 = load double, ptr %137, align 8, !tbaa !5
  store double %138, ptr %10, align 8, !tbaa !5
  %139 = load i32, ptr %5, align 4, !tbaa !9
  %140 = sext i32 %139 to i64
  %141 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %140
  %142 = getelementptr inbounds [7 x double], ptr %141, i64 0, i64 4
  %143 = load double, ptr %142, align 8, !tbaa !5
  store double %143, ptr %11, align 8, !tbaa !5
  %144 = load i32, ptr %5, align 4, !tbaa !9
  %145 = sext i32 %144 to i64
  %146 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %145
  %147 = getelementptr inbounds [7 x double], ptr %146, i64 0, i64 5
  %148 = load double, ptr %147, align 8, !tbaa !5
  store double %148, ptr %12, align 8, !tbaa !5
  %149 = load i32, ptr %5, align 4, !tbaa !9
  %150 = sext i32 %149 to i64
  %151 = getelementptr inbounds [5 x [7 x double]], ptr @coeffr3, i64 0, i64 %150
  %152 = getelementptr inbounds [7 x double], ptr %151, i64 0, i64 6
  %153 = load double, ptr %152, align 8, !tbaa !5
  store double %153, ptr %13, align 8, !tbaa !5
  %154 = load double, ptr %6, align 8, !tbaa !5
  %155 = load double, ptr %3, align 8, !tbaa !5
  %156 = load double, ptr %7, align 8, !tbaa !5
  %157 = call double @pow(double noundef %155, double noundef %156) #4
  %158 = fmul double %154, %157
  %159 = load double, ptr %4, align 8, !tbaa !5
  %160 = load double, ptr %8, align 8, !tbaa !5
  %161 = call double @pow(double noundef %159, double noundef %160) #4
  %162 = fmul double %158, %161
  %163 = load double, ptr %10, align 8, !tbaa !5
  %164 = fneg double %163
  %165 = load double, ptr %3, align 8, !tbaa !5
  %166 = load double, ptr %13, align 8, !tbaa !5
  %167 = fsub double %165, %166
  %168 = fmul double %164, %167
  %169 = load double, ptr %3, align 8, !tbaa !5
  %170 = load double, ptr %13, align 8, !tbaa !5
  %171 = fsub double %169, %170
  %172 = load double, ptr %11, align 8, !tbaa !5
  %173 = load double, ptr %4, align 8, !tbaa !5
  %174 = load double, ptr %12, align 8, !tbaa !5
  %175 = fsub double %173, %174
  %176 = fmul double %172, %175
  %177 = load double, ptr %4, align 8, !tbaa !5
  %178 = load double, ptr %12, align 8, !tbaa !5
  %179 = fsub double %177, %178
  %180 = fmul double %176, %179
  %181 = fneg double %180
  %182 = call double @llvm.fmuladd.f64(double %168, double %171, double %181)
  %183 = call double @exp(double noundef %182) #4
  %184 = fmul double %162, %183
  %185 = load double, ptr %7, align 8, !tbaa !5
  %186 = load double, ptr %3, align 8, !tbaa !5
  %187 = fdiv double %185, %186
  %188 = load double, ptr %10, align 8, !tbaa !5
  %189 = fmul double 2.000000e+00, %188
  %190 = load double, ptr %3, align 8, !tbaa !5
  %191 = load double, ptr %13, align 8, !tbaa !5
  %192 = fsub double %190, %191
  %193 = fneg double %189
  %194 = call double @llvm.fmuladd.f64(double %193, double %192, double %187)
  %195 = load double, ptr %25, align 8, !tbaa !5
  %196 = call double @llvm.fmuladd.f64(double %184, double %194, double %195)
  store double %196, ptr %25, align 8, !tbaa !5
  br label %197

197:                                              ; preds = %118
  %198 = load i32, ptr %5, align 4, !tbaa !9
  %199 = add nsw i32 %198, 1
  store i32 %199, ptr %5, align 4, !tbaa !9
  br label %115, !llvm.loop !14

200:                                              ; preds = %115
  store i32 0, ptr %5, align 4, !tbaa !9
  br label %201

201:                                              ; preds = %325, %200
  %202 = load i32, ptr %5, align 4, !tbaa !9
  %203 = icmp slt i32 %202, 3
  br i1 %203, label %204, label %328

204:                                              ; preds = %201
  %205 = load i32, ptr %5, align 4, !tbaa !9
  %206 = sext i32 %205 to i64
  %207 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %206
  %208 = getelementptr inbounds [8 x double], ptr %207, i64 0, i64 0
  %209 = load double, ptr %208, align 16, !tbaa !5
  store double %209, ptr %6, align 8, !tbaa !5
  %210 = load i32, ptr %5, align 4, !tbaa !9
  %211 = sext i32 %210 to i64
  %212 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %211
  %213 = getelementptr inbounds [8 x double], ptr %212, i64 0, i64 1
  %214 = load double, ptr %213, align 8, !tbaa !5
  store double %214, ptr %14, align 8, !tbaa !5
  %215 = load i32, ptr %5, align 4, !tbaa !9
  %216 = sext i32 %215 to i64
  %217 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %216
  %218 = getelementptr inbounds [8 x double], ptr %217, i64 0, i64 2
  %219 = load double, ptr %218, align 16, !tbaa !5
  store double %219, ptr %15, align 8, !tbaa !5
  %220 = load i32, ptr %5, align 4, !tbaa !9
  %221 = sext i32 %220 to i64
  %222 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %221
  %223 = getelementptr inbounds [8 x double], ptr %222, i64 0, i64 3
  %224 = load double, ptr %223, align 8, !tbaa !5
  store double %224, ptr %11, align 8, !tbaa !5
  %225 = load i32, ptr %5, align 4, !tbaa !9
  %226 = sext i32 %225 to i64
  %227 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %226
  %228 = getelementptr inbounds [8 x double], ptr %227, i64 0, i64 4
  %229 = load double, ptr %228, align 16, !tbaa !5
  store double %229, ptr %16, align 8, !tbaa !5
  %230 = load i32, ptr %5, align 4, !tbaa !9
  %231 = sext i32 %230 to i64
  %232 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %231
  %233 = getelementptr inbounds [8 x double], ptr %232, i64 0, i64 5
  %234 = load double, ptr %233, align 8, !tbaa !5
  store double %234, ptr %17, align 8, !tbaa !5
  %235 = load i32, ptr %5, align 4, !tbaa !9
  %236 = sext i32 %235 to i64
  %237 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %236
  %238 = getelementptr inbounds [8 x double], ptr %237, i64 0, i64 6
  %239 = load double, ptr %238, align 16, !tbaa !5
  store double %239, ptr %18, align 8, !tbaa !5
  %240 = load i32, ptr %5, align 4, !tbaa !9
  %241 = sext i32 %240 to i64
  %242 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %241
  %243 = getelementptr inbounds [8 x double], ptr %242, i64 0, i64 7
  %244 = load double, ptr %243, align 8, !tbaa !5
  store double %244, ptr %19, align 8, !tbaa !5
  %245 = load double, ptr %4, align 8, !tbaa !5
  %246 = fsub double 1.000000e+00, %245
  %247 = load double, ptr %16, align 8, !tbaa !5
  %248 = load double, ptr %3, align 8, !tbaa !5
  %249 = fsub double %248, 1.000000e+00
  %250 = load double, ptr %3, align 8, !tbaa !5
  %251 = fsub double %250, 1.000000e+00
  %252 = fmul double %249, %251
  %253 = load double, ptr %11, align 8, !tbaa !5
  %254 = fdiv double 5.000000e-01, %253
  %255 = call double @pow(double noundef %252, double noundef %254) #4
  %256 = call double @llvm.fmuladd.f64(double %247, double %255, double %246)
  store double %256, ptr %22, align 8, !tbaa !5
  %257 = load double, ptr %22, align 8, !tbaa !5
  %258 = load double, ptr %22, align 8, !tbaa !5
  %259 = load double, ptr %17, align 8, !tbaa !5
  %260 = load double, ptr %3, align 8, !tbaa !5
  %261 = fsub double %260, 1.000000e+00
  %262 = load double, ptr %3, align 8, !tbaa !5
  %263 = fsub double %262, 1.000000e+00
  %264 = fmul double %261, %263
  %265 = load double, ptr %14, align 8, !tbaa !5
  %266 = call double @pow(double noundef %264, double noundef %265) #4
  %267 = fmul double %259, %266
  %268 = call double @llvm.fmuladd.f64(double %257, double %258, double %267)
  store double %268, ptr %20, align 8, !tbaa !5
  %269 = load double, ptr %18, align 8, !tbaa !5
  %270 = fneg double %269
  %271 = load double, ptr %3, align 8, !tbaa !5
  %272 = fsub double %271, 1.000000e+00
  %273 = fmul double %270, %272
  %274 = load double, ptr %3, align 8, !tbaa !5
  %275 = fsub double %274, 1.000000e+00
  %276 = load double, ptr %19, align 8, !tbaa !5
  %277 = load double, ptr %4, align 8, !tbaa !5
  %278 = fsub double %277, 1.000000e+00
  %279 = fmul double %276, %278
  %280 = load double, ptr %4, align 8, !tbaa !5
  %281 = fsub double %280, 1.000000e+00
  %282 = fmul double %279, %281
  %283 = fneg double %282
  %284 = call double @llvm.fmuladd.f64(double %273, double %275, double %283)
  %285 = call double @exp(double noundef %284) #4
  store double %285, ptr %21, align 8, !tbaa !5
  %286 = load double, ptr %3, align 8, !tbaa !5
  %287 = load double, ptr %20, align 8, !tbaa !5
  %288 = load double, ptr %16, align 8, !tbaa !5
  %289 = load double, ptr %22, align 8, !tbaa !5
  %290 = load double, ptr %11, align 8, !tbaa !5
  %291 = load double, ptr %17, align 8, !tbaa !5
  %292 = load double, ptr %14, align 8, !tbaa !5
  %293 = call double @deltd(double noundef %286, double noundef %287, double noundef %288, double noundef %289, double noundef %290, double noundef %291, double noundef %292)
  store double %293, ptr %23, align 8, !tbaa !5
  %294 = load double, ptr %23, align 8, !tbaa !5
  %295 = load double, ptr %3, align 8, !tbaa !5
  %296 = load double, ptr %20, align 8, !tbaa !5
  %297 = load double, ptr %16, align 8, !tbaa !5
  %298 = load double, ptr %22, align 8, !tbaa !5
  %299 = load double, ptr %11, align 8, !tbaa !5
  %300 = load double, ptr %17, align 8, !tbaa !5
  %301 = load double, ptr %14, align 8, !tbaa !5
  %302 = call double @deltdd(double noundef %294, double noundef %295, double noundef %296, double noundef %297, double noundef %298, double noundef %299, double noundef %300, double noundef %301)
  store double %302, ptr %24, align 8, !tbaa !5
  %303 = load double, ptr %6, align 8, !tbaa !5
  %304 = load double, ptr %20, align 8, !tbaa !5
  %305 = load double, ptr %15, align 8, !tbaa !5
  %306 = call double @pow(double noundef %304, double noundef %305) #4
  %307 = load double, ptr %21, align 8, !tbaa !5
  %308 = load double, ptr %3, align 8, !tbaa !5
  %309 = load double, ptr %3, align 8, !tbaa !5
  %310 = load double, ptr %21, align 8, !tbaa !5
  %311 = load double, ptr %18, align 8, !tbaa !5
  %312 = call double @psid(double noundef %309, double noundef %310, double noundef %311)
  %313 = call double @llvm.fmuladd.f64(double %308, double %312, double %307)
  %314 = load double, ptr %15, align 8, !tbaa !5
  %315 = load double, ptr %20, align 8, !tbaa !5
  %316 = load double, ptr %23, align 8, !tbaa !5
  %317 = call double @deltabid(double noundef %314, double noundef %315, double noundef %316)
  %318 = load double, ptr %3, align 8, !tbaa !5
  %319 = fmul double %317, %318
  %320 = load double, ptr %21, align 8, !tbaa !5
  %321 = fmul double %319, %320
  %322 = call double @llvm.fmuladd.f64(double %306, double %313, double %321)
  %323 = load double, ptr %25, align 8, !tbaa !5
  %324 = call double @llvm.fmuladd.f64(double %303, double %322, double %323)
  store double %324, ptr %25, align 8, !tbaa !5
  br label %325

325:                                              ; preds = %204
  %326 = load i32, ptr %5, align 4, !tbaa !9
  %327 = add nsw i32 %326, 1
  store i32 %327, ptr %5, align 4, !tbaa !9
  br label %201, !llvm.loop !15

328:                                              ; preds = %201
  %329 = load double, ptr %25, align 8, !tbaa !5
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %22) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %21) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %20) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %19) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %18) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %17) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %16) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %15) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %14) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %13) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %12) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %11) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %10) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %9) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %8) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %7) #4
  call void @llvm.lifetime.end.p0(i64 8, ptr %6) #4
  call void @llvm.lifetime.end.p0(i64 4, ptr %5) #4
  ret double %329
}

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

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

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


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

@coeffr1 = external dso_local local_unnamed_addr constant [7 x [3 x double]], align 16
@coeffr2 = external dso_local local_unnamed_addr constant [27 x [4 x double]], align 16
@coeffr3 = external dso_local local_unnamed_addr constant [5 x [7 x double]], align 16
@coeffr4 = external dso_local local_unnamed_addr constant [3 x [8 x double]], align 16

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

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

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

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

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

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

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

; Function Attrs: nounwind uwtable
define dso_local double @phird(double noundef %0, double noundef %1) local_unnamed_addr #2 {
.preheader121.preheader:
  %2 = load double, ptr @coeffr1, align 16, !tbaa !5
  %3 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 8), align 8, !tbaa !5
  %4 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 16), align 16, !tbaa !5
  %5 = fmul double %2, %3
  %6 = fadd double %3, -1.000000e+00
  %7 = tail call double @pow(double noundef %0, double noundef %6) #3
  %8 = fmul double %5, %7
  %9 = tail call double @pow(double noundef %1, double noundef %4) #3
  %10 = tail call double @llvm.fmuladd.f64(double %8, double %9, double 0.000000e+00)
  %11 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 24), align 8, !tbaa !5
  %12 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 32), align 16, !tbaa !5
  %13 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 40), align 8, !tbaa !5
  %14 = fmul double %11, %12
  %15 = fadd double %12, -1.000000e+00
  %16 = tail call double @pow(double noundef %0, double noundef %15) #3
  %17 = fmul double %14, %16
  %18 = tail call double @pow(double noundef %1, double noundef %13) #3
  %19 = tail call double @llvm.fmuladd.f64(double %17, double %18, double %10)
  %20 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 48), align 16, !tbaa !5
  %21 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 56), align 8, !tbaa !5
  %22 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 64), align 16, !tbaa !5
  %23 = fmul double %20, %21
  %24 = fadd double %21, -1.000000e+00
  %25 = tail call double @pow(double noundef %0, double noundef %24) #3
  %26 = fmul double %23, %25
  %27 = tail call double @pow(double noundef %1, double noundef %22) #3
  %28 = tail call double @llvm.fmuladd.f64(double %26, double %27, double %19)
  %29 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 72), align 8, !tbaa !5
  %30 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 80), align 16, !tbaa !5
  %31 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 88), align 8, !tbaa !5
  %32 = fmul double %29, %30
  %33 = fadd double %30, -1.000000e+00
  %34 = tail call double @pow(double noundef %0, double noundef %33) #3
  %35 = fmul double %32, %34
  %36 = tail call double @pow(double noundef %1, double noundef %31) #3
  %37 = tail call double @llvm.fmuladd.f64(double %35, double %36, double %28)
  %38 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 96), align 16, !tbaa !5
  %39 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 104), align 8, !tbaa !5
  %40 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 112), align 16, !tbaa !5
  %41 = fmul double %38, %39
  %42 = fadd double %39, -1.000000e+00
  %43 = tail call double @pow(double noundef %0, double noundef %42) #3
  %44 = fmul double %41, %43
  %45 = tail call double @pow(double noundef %1, double noundef %40) #3
  %46 = tail call double @llvm.fmuladd.f64(double %44, double %45, double %37)
  %47 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 120), align 8, !tbaa !5
  %48 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 128), align 16, !tbaa !5
  %49 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 136), align 8, !tbaa !5
  %50 = fmul double %47, %48
  %51 = fadd double %48, -1.000000e+00
  %52 = tail call double @pow(double noundef %0, double noundef %51) #3
  %53 = fmul double %50, %52
  %54 = tail call double @pow(double noundef %1, double noundef %49) #3
  %55 = tail call double @llvm.fmuladd.f64(double %53, double %54, double %46)
  %56 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 144), align 16, !tbaa !5
  %57 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 152), align 8, !tbaa !5
  %58 = load double, ptr getelementptr inbounds (i8, ptr @coeffr1, i64 160), align 16, !tbaa !5
  %59 = fmul double %56, %57
  %60 = fadd double %57, -1.000000e+00
  %61 = tail call double @pow(double noundef %0, double noundef %60) #3
  %62 = fmul double %59, %61
  %63 = tail call double @pow(double noundef %1, double noundef %58) #3
  %64 = tail call double @llvm.fmuladd.f64(double %62, double %63, double %55)
  br label %.preheader121

.preheader121:                                    ; preds = %.preheader121.preheader, %.preheader121
  %indvars.iv = phi i64 [ 0, %.preheader121.preheader ], [ %indvars.iv.next, %.preheader121 ]
  %.1125 = phi double [ %64, %.preheader121.preheader ], [ %85, %.preheader121 ]
  %65 = getelementptr inbounds [27 x [4 x double]], ptr @coeffr2, i64 0, i64 %indvars.iv
  %66 = load double, ptr %65, align 16, !tbaa !5
  %67 = getelementptr inbounds i8, ptr %65, i64 8
  %68 = load double, ptr %67, align 8, !tbaa !5
  %69 = getelementptr inbounds i8, ptr %65, i64 16
  %70 = load double, ptr %69, align 16, !tbaa !5
  %71 = getelementptr inbounds i8, ptr %65, i64 24
  %72 = load double, ptr %71, align 8, !tbaa !5
  %73 = tail call double @pow(double noundef %1, double noundef %70) #3
  %74 = fmul double %66, %73
  %75 = tail call double @pow(double noundef %0, double noundef %72) #3
  %76 = fneg double %75
  %77 = tail call double @exp(double noundef %76) #3
  %78 = fmul double %74, %77
  %79 = fadd double %68, -1.000000e+00
  %80 = tail call double @pow(double noundef %0, double noundef %79) #3
  %81 = fmul double %78, %80
  %82 = tail call double @pow(double noundef %0, double noundef %72) #3
  %83 = fneg double %72
  %84 = tail call double @llvm.fmuladd.f64(double %83, double %82, double %68)
  %85 = tail call double @llvm.fmuladd.f64(double %81, double %84, double %.1125)
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, 27
  br i1 %exitcond.not, label %.preheader120.preheader, label %.preheader121, !llvm.loop !9

.preheader120.preheader:                          ; preds = %.preheader121
  %86 = load double, ptr @coeffr3, align 16, !tbaa !5
  %87 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 8), align 8, !tbaa !5
  %88 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 16), align 16, !tbaa !5
  %89 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 24), align 8, !tbaa !5
  %90 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 32), align 16, !tbaa !5
  %91 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 40), align 8, !tbaa !5
  %92 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 48), align 16, !tbaa !5
  %93 = tail call double @pow(double noundef %0, double noundef %87) #3
  %94 = fmul double %86, %93
  %95 = tail call double @pow(double noundef %1, double noundef %88) #3
  %96 = fmul double %94, %95
  %97 = fneg double %89
  %98 = fsub double %0, %92
  %99 = fmul double %98, %97
  %100 = fsub double %1, %91
  %101 = fneg double %90
  %102 = fmul double %100, %101
  %103 = fmul double %100, %102
  %104 = tail call double @llvm.fmuladd.f64(double %99, double %98, double %103)
  %105 = tail call double @exp(double noundef %104) #3
  %106 = fmul double %96, %105
  %107 = fdiv double %87, %0
  %108 = fmul double %89, -2.000000e+00
  %109 = tail call double @llvm.fmuladd.f64(double %108, double %98, double %107)
  %110 = tail call double @llvm.fmuladd.f64(double %106, double %109, double %85)
  %111 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 56), align 8, !tbaa !5
  %112 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 64), align 16, !tbaa !5
  %113 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 72), align 8, !tbaa !5
  %114 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 80), align 16, !tbaa !5
  %115 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 88), align 8, !tbaa !5
  %116 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 96), align 16, !tbaa !5
  %117 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 104), align 8, !tbaa !5
  %118 = tail call double @pow(double noundef %0, double noundef %112) #3
  %119 = fmul double %111, %118
  %120 = tail call double @pow(double noundef %1, double noundef %113) #3
  %121 = fmul double %119, %120
  %122 = fneg double %114
  %123 = fsub double %0, %117
  %124 = fmul double %123, %122
  %125 = fsub double %1, %116
  %126 = fneg double %115
  %127 = fmul double %125, %126
  %128 = fmul double %125, %127
  %129 = tail call double @llvm.fmuladd.f64(double %124, double %123, double %128)
  %130 = tail call double @exp(double noundef %129) #3
  %131 = fmul double %121, %130
  %132 = fdiv double %112, %0
  %133 = fmul double %114, -2.000000e+00
  %134 = tail call double @llvm.fmuladd.f64(double %133, double %123, double %132)
  %135 = tail call double @llvm.fmuladd.f64(double %131, double %134, double %110)
  %136 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 112), align 16, !tbaa !5
  %137 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 120), align 8, !tbaa !5
  %138 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 128), align 16, !tbaa !5
  %139 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 136), align 8, !tbaa !5
  %140 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 144), align 16, !tbaa !5
  %141 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 152), align 8, !tbaa !5
  %142 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 160), align 16, !tbaa !5
  %143 = tail call double @pow(double noundef %0, double noundef %137) #3
  %144 = fmul double %136, %143
  %145 = tail call double @pow(double noundef %1, double noundef %138) #3
  %146 = fmul double %144, %145
  %147 = fneg double %139
  %148 = fsub double %0, %142
  %149 = fmul double %148, %147
  %150 = fsub double %1, %141
  %151 = fneg double %140
  %152 = fmul double %150, %151
  %153 = fmul double %150, %152
  %154 = tail call double @llvm.fmuladd.f64(double %149, double %148, double %153)
  %155 = tail call double @exp(double noundef %154) #3
  %156 = fmul double %146, %155
  %157 = fdiv double %137, %0
  %158 = fmul double %139, -2.000000e+00
  %159 = tail call double @llvm.fmuladd.f64(double %158, double %148, double %157)
  %160 = tail call double @llvm.fmuladd.f64(double %156, double %159, double %135)
  %161 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 168), align 8, !tbaa !5
  %162 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 176), align 16, !tbaa !5
  %163 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 184), align 8, !tbaa !5
  %164 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 192), align 16, !tbaa !5
  %165 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 200), align 8, !tbaa !5
  %166 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 208), align 16, !tbaa !5
  %167 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 216), align 8, !tbaa !5
  %168 = tail call double @pow(double noundef %0, double noundef %162) #3
  %169 = fmul double %161, %168
  %170 = tail call double @pow(double noundef %1, double noundef %163) #3
  %171 = fmul double %169, %170
  %172 = fneg double %164
  %173 = fsub double %0, %167
  %174 = fmul double %173, %172
  %175 = fsub double %1, %166
  %176 = fneg double %165
  %177 = fmul double %175, %176
  %178 = fmul double %175, %177
  %179 = tail call double @llvm.fmuladd.f64(double %174, double %173, double %178)
  %180 = tail call double @exp(double noundef %179) #3
  %181 = fmul double %171, %180
  %182 = fdiv double %162, %0
  %183 = fmul double %164, -2.000000e+00
  %184 = tail call double @llvm.fmuladd.f64(double %183, double %173, double %182)
  %185 = tail call double @llvm.fmuladd.f64(double %181, double %184, double %160)
  %186 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 224), align 16, !tbaa !5
  %187 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 232), align 8, !tbaa !5
  %188 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 240), align 16, !tbaa !5
  %189 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 248), align 8, !tbaa !5
  %190 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 256), align 16, !tbaa !5
  %191 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 264), align 8, !tbaa !5
  %192 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 272), align 16, !tbaa !5
  %193 = tail call double @pow(double noundef %0, double noundef %187) #3
  %194 = fmul double %186, %193
  %195 = tail call double @pow(double noundef %1, double noundef %188) #3
  %196 = fmul double %194, %195
  %197 = fneg double %189
  %198 = fsub double %0, %192
  %199 = fmul double %198, %197
  %200 = fsub double %1, %191
  %201 = fneg double %190
  %202 = fmul double %200, %201
  %203 = fmul double %200, %202
  %204 = tail call double @llvm.fmuladd.f64(double %199, double %198, double %203)
  %205 = tail call double @exp(double noundef %204) #3
  %206 = fmul double %196, %205
  %207 = fdiv double %187, %0
  %208 = fmul double %189, -2.000000e+00
  %209 = tail call double @llvm.fmuladd.f64(double %208, double %198, double %207)
  %210 = tail call double @llvm.fmuladd.f64(double %206, double %209, double %185)
  %211 = fsub double 1.000000e+00, %1
  %212 = fadd double %0, -1.000000e+00
  %213 = fmul double %212, %212
  %214 = fadd double %1, -1.000000e+00
  %215 = load double, ptr @coeffr4, align 16, !tbaa !5
  %216 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 8), align 8, !tbaa !5
  %217 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 16), align 16, !tbaa !5
  %218 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 24), align 8, !tbaa !5
  %219 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 32), align 16, !tbaa !5
  %220 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 40), align 8, !tbaa !5
  %221 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 48), align 16, !tbaa !5
  %222 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 56), align 8, !tbaa !5
  %223 = fdiv double 5.000000e-01, %218
  %224 = tail call double @pow(double noundef %213, double noundef %223) #3
  %225 = tail call double @llvm.fmuladd.f64(double %219, double %224, double %211)
  %226 = tail call double @pow(double noundef %213, double noundef %216) #3
  %227 = fmul double %220, %226
  %228 = tail call double @llvm.fmuladd.f64(double %225, double %225, double %227)
  %229 = fneg double %221
  %230 = fmul double %212, %229
  %231 = fneg double %222
  %232 = fmul double %214, %231
  %233 = fmul double %214, %232
  %234 = tail call double @llvm.fmuladd.f64(double %230, double %212, double %233)
  %235 = tail call double @exp(double noundef %234) #3
  %236 = tail call double @deltd(double noundef %0, double noundef %228, double noundef %219, double noundef %225, double noundef %218, double noundef %220, double noundef %216)
  %237 = tail call double @deltdd(double noundef %236, double noundef %0, double noundef %228, double noundef %219, double noundef %225, double noundef %218, double noundef %220, double noundef %216)
  %238 = tail call double @pow(double noundef %228, double noundef %217) #3
  %239 = tail call double @psid(double noundef %0, double noundef %235, double noundef %221)
  %240 = tail call double @llvm.fmuladd.f64(double %0, double %239, double %235)
  %241 = tail call double @deltabid(double noundef %217, double noundef %228, double noundef %236)
  %242 = fmul double %241, %0
  %243 = fmul double %235, %242
  %244 = tail call double @llvm.fmuladd.f64(double %238, double %240, double %243)
  %245 = tail call double @llvm.fmuladd.f64(double %215, double %244, double %210)
  %246 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 64), align 16, !tbaa !5
  %247 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 72), align 8, !tbaa !5
  %248 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 80), align 16, !tbaa !5
  %249 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 88), align 8, !tbaa !5
  %250 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 96), align 16, !tbaa !5
  %251 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 104), align 8, !tbaa !5
  %252 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 112), align 16, !tbaa !5
  %253 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 120), align 8, !tbaa !5
  %254 = fdiv double 5.000000e-01, %249
  %255 = tail call double @pow(double noundef %213, double noundef %254) #3
  %256 = tail call double @llvm.fmuladd.f64(double %250, double %255, double %211)
  %257 = tail call double @pow(double noundef %213, double noundef %247) #3
  %258 = fmul double %251, %257
  %259 = tail call double @llvm.fmuladd.f64(double %256, double %256, double %258)
  %260 = fneg double %252
  %261 = fmul double %212, %260
  %262 = fneg double %253
  %263 = fmul double %214, %262
  %264 = fmul double %214, %263
  %265 = tail call double @llvm.fmuladd.f64(double %261, double %212, double %264)
  %266 = tail call double @exp(double noundef %265) #3
  %267 = tail call double @deltd(double noundef %0, double noundef %259, double noundef %250, double noundef %256, double noundef %249, double noundef %251, double noundef %247)
  %268 = tail call double @deltdd(double noundef %267, double noundef %0, double noundef %259, double noundef %250, double noundef %256, double noundef %249, double noundef %251, double noundef %247)
  %269 = tail call double @pow(double noundef %259, double noundef %248) #3
  %270 = tail call double @psid(double noundef %0, double noundef %266, double noundef %252)
  %271 = tail call double @llvm.fmuladd.f64(double %0, double %270, double %266)
  %272 = tail call double @deltabid(double noundef %248, double noundef %259, double noundef %267)
  %273 = fmul double %272, %0
  %274 = fmul double %266, %273
  %275 = tail call double @llvm.fmuladd.f64(double %269, double %271, double %274)
  %276 = tail call double @llvm.fmuladd.f64(double %246, double %275, double %245)
  %277 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 128), align 16, !tbaa !5
  %278 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 136), align 8, !tbaa !5
  %279 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 144), align 16, !tbaa !5
  %280 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 152), align 8, !tbaa !5
  %281 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 160), align 16, !tbaa !5
  %282 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 168), align 8, !tbaa !5
  %283 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 176), align 16, !tbaa !5
  %284 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 184), align 8, !tbaa !5
  %285 = fdiv double 5.000000e-01, %280
  %286 = tail call double @pow(double noundef %213, double noundef %285) #3
  %287 = tail call double @llvm.fmuladd.f64(double %281, double %286, double %211)
  %288 = tail call double @pow(double noundef %213, double noundef %278) #3
  %289 = fmul double %282, %288
  %290 = tail call double @llvm.fmuladd.f64(double %287, double %287, double %289)
  %291 = fneg double %283
  %292 = fmul double %212, %291
  %293 = fneg double %284
  %294 = fmul double %214, %293
  %295 = fmul double %214, %294
  %296 = tail call double @llvm.fmuladd.f64(double %292, double %212, double %295)
  %297 = tail call double @exp(double noundef %296) #3
  %298 = tail call double @deltd(double noundef %0, double noundef %290, double noundef %281, double noundef %287, double noundef %280, double noundef %282, double noundef %278)
  %299 = tail call double @deltdd(double noundef %298, double noundef %0, double noundef %290, double noundef %281, double noundef %287, double noundef %280, double noundef %282, double noundef %278)
  %300 = tail call double @pow(double noundef %290, double noundef %279) #3
  %301 = tail call double @psid(double noundef %0, double noundef %297, double noundef %283)
  %302 = tail call double @llvm.fmuladd.f64(double %0, double %301, double %297)
  %303 = tail call double @deltabid(double noundef %279, double noundef %290, double noundef %298)
  %304 = fmul double %303, %0
  %305 = fmul double %297, %304
  %306 = tail call double @llvm.fmuladd.f64(double %300, double %302, double %305)
  %307 = tail call double @llvm.fmuladd.f64(double %277, double %306, double %276)
  ret double %307
}

attributes #0 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #1 = { mustprogress nofree nounwind willreturn memory(write) "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #2 = { nounwind uwtable "min-legal-vector-width"="0" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #3 = { nounwind }

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

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


; Effective O3 passes: EarlyCSEPass,GVNPass,GlobalOptPass,IndVarSimplifyPass,InferAlignmentPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,LCSSAPass,LICMPass,LoopRotatePass,LoopSimplifyPass,LoopUnrollPass,ReassociatePass,SROAPass,SimplifyCFGPass,TailCallElimPass