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

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

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

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

; Function Attrs: nounwind uwtable
define dso_local double @phirdt(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 = %61, %2
  %27 = load i32, ptr %5, align 4, !tbaa !9
  %28 = icmp slt i32 %27, 7
  br i1 %28, label %29, label %64

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 %8, align 8, !tbaa !5
  %47 = fmul double %45, %46
  %48 = load double, ptr %7, align 8, !tbaa !5
  %49 = fmul double %47, %48
  %50 = load double, ptr %4, align 8, !tbaa !5
  %51 = load double, ptr %8, align 8, !tbaa !5
  %52 = fsub double %51, 1.000000e+00
  %53 = call double @pow(double noundef %50, double noundef %52) #4
  %54 = fmul double %49, %53
  %55 = load double, ptr %3, align 8, !tbaa !5
  %56 = load double, ptr %7, align 8, !tbaa !5
  %57 = fsub double %56, 1.000000e+00
  %58 = call double @pow(double noundef %55, double noundef %57) #4
  %59 = load double, ptr %25, align 8, !tbaa !5
  %60 = call double @llvm.fmuladd.f64(double %54, double %58, double %59)
  store double %60, ptr %25, align 8, !tbaa !5
  br label %61

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

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

65:                                               ; preds = %117, %64
  %66 = load i32, ptr %5, align 4, !tbaa !9
  %67 = icmp slt i32 %66, 27
  br i1 %67, label %68, label %120

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

117:                                              ; preds = %68
  %118 = load i32, ptr %5, align 4, !tbaa !9
  %119 = add nsw i32 %118, 1
  store i32 %119, ptr %5, align 4, !tbaa !9
  br label %65, !llvm.loop !13

120:                                              ; preds = %65
  store i32 0, ptr %5, align 4, !tbaa !9
  br label %121

121:                                              ; preds = %214, %120
  %122 = load i32, ptr %5, align 4, !tbaa !9
  %123 = icmp slt i32 %122, 5
  br i1 %123, label %124, label %217

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

214:                                              ; preds = %124
  %215 = load i32, ptr %5, align 4, !tbaa !9
  %216 = add nsw i32 %215, 1
  store i32 %216, ptr %5, align 4, !tbaa !9
  br label %121, !llvm.loop !14

217:                                              ; preds = %121
  store i32 0, ptr %5, align 4, !tbaa !9
  br label %218

218:                                              ; preds = %374, %217
  %219 = load i32, ptr %5, align 4, !tbaa !9
  %220 = icmp slt i32 %219, 3
  br i1 %220, label %221, label %377

221:                                              ; preds = %218
  %222 = load i32, ptr %5, align 4, !tbaa !9
  %223 = sext i32 %222 to i64
  %224 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %223
  %225 = getelementptr inbounds [8 x double], ptr %224, i64 0, i64 0
  %226 = load double, ptr %225, align 16, !tbaa !5
  store double %226, ptr %6, align 8, !tbaa !5
  %227 = load i32, ptr %5, align 4, !tbaa !9
  %228 = sext i32 %227 to i64
  %229 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %228
  %230 = getelementptr inbounds [8 x double], ptr %229, i64 0, i64 1
  %231 = load double, ptr %230, align 8, !tbaa !5
  store double %231, ptr %14, align 8, !tbaa !5
  %232 = load i32, ptr %5, align 4, !tbaa !9
  %233 = sext i32 %232 to i64
  %234 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %233
  %235 = getelementptr inbounds [8 x double], ptr %234, i64 0, i64 2
  %236 = load double, ptr %235, align 16, !tbaa !5
  store double %236, ptr %15, align 8, !tbaa !5
  %237 = load i32, ptr %5, align 4, !tbaa !9
  %238 = sext i32 %237 to i64
  %239 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %238
  %240 = getelementptr inbounds [8 x double], ptr %239, i64 0, i64 3
  %241 = load double, ptr %240, align 8, !tbaa !5
  store double %241, ptr %11, align 8, !tbaa !5
  %242 = load i32, ptr %5, align 4, !tbaa !9
  %243 = sext i32 %242 to i64
  %244 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %243
  %245 = getelementptr inbounds [8 x double], ptr %244, i64 0, i64 4
  %246 = load double, ptr %245, align 16, !tbaa !5
  store double %246, ptr %16, align 8, !tbaa !5
  %247 = load i32, ptr %5, align 4, !tbaa !9
  %248 = sext i32 %247 to i64
  %249 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %248
  %250 = getelementptr inbounds [8 x double], ptr %249, i64 0, i64 5
  %251 = load double, ptr %250, align 8, !tbaa !5
  store double %251, ptr %17, align 8, !tbaa !5
  %252 = load i32, ptr %5, align 4, !tbaa !9
  %253 = sext i32 %252 to i64
  %254 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %253
  %255 = getelementptr inbounds [8 x double], ptr %254, i64 0, i64 6
  %256 = load double, ptr %255, align 16, !tbaa !5
  store double %256, ptr %18, align 8, !tbaa !5
  %257 = load i32, ptr %5, align 4, !tbaa !9
  %258 = sext i32 %257 to i64
  %259 = getelementptr inbounds [3 x [8 x double]], ptr @coeffr4, i64 0, i64 %258
  %260 = getelementptr inbounds [8 x double], ptr %259, i64 0, i64 7
  %261 = load double, ptr %260, align 8, !tbaa !5
  store double %261, ptr %19, align 8, !tbaa !5
  %262 = load double, ptr %4, align 8, !tbaa !5
  %263 = fsub double 1.000000e+00, %262
  %264 = load double, ptr %16, align 8, !tbaa !5
  %265 = load double, ptr %3, align 8, !tbaa !5
  %266 = fsub double %265, 1.000000e+00
  %267 = load double, ptr %3, align 8, !tbaa !5
  %268 = fsub double %267, 1.000000e+00
  %269 = fmul double %266, %268
  %270 = load double, ptr %11, align 8, !tbaa !5
  %271 = fdiv double 5.000000e-01, %270
  %272 = call double @pow(double noundef %269, double noundef %271) #4
  %273 = call double @llvm.fmuladd.f64(double %264, double %272, double %263)
  store double %273, ptr %22, align 8, !tbaa !5
  %274 = load double, ptr %22, align 8, !tbaa !5
  %275 = load double, ptr %22, align 8, !tbaa !5
  %276 = load double, ptr %17, align 8, !tbaa !5
  %277 = load double, ptr %3, align 8, !tbaa !5
  %278 = fsub double %277, 1.000000e+00
  %279 = load double, ptr %3, align 8, !tbaa !5
  %280 = fsub double %279, 1.000000e+00
  %281 = fmul double %278, %280
  %282 = load double, ptr %14, align 8, !tbaa !5
  %283 = call double @pow(double noundef %281, double noundef %282) #4
  %284 = fmul double %276, %283
  %285 = call double @llvm.fmuladd.f64(double %274, double %275, double %284)
  store double %285, ptr %20, align 8, !tbaa !5
  %286 = load double, ptr %18, align 8, !tbaa !5
  %287 = fneg double %286
  %288 = load double, ptr %3, align 8, !tbaa !5
  %289 = fsub double %288, 1.000000e+00
  %290 = fmul double %287, %289
  %291 = load double, ptr %3, align 8, !tbaa !5
  %292 = fsub double %291, 1.000000e+00
  %293 = load double, ptr %19, align 8, !tbaa !5
  %294 = load double, ptr %4, align 8, !tbaa !5
  %295 = fsub double %294, 1.000000e+00
  %296 = fmul double %293, %295
  %297 = load double, ptr %4, align 8, !tbaa !5
  %298 = fsub double %297, 1.000000e+00
  %299 = fmul double %296, %298
  %300 = fneg double %299
  %301 = call double @llvm.fmuladd.f64(double %290, double %292, double %300)
  %302 = call double @exp(double noundef %301) #4
  store double %302, ptr %21, align 8, !tbaa !5
  %303 = load double, ptr %3, align 8, !tbaa !5
  %304 = load double, ptr %20, align 8, !tbaa !5
  %305 = load double, ptr %16, align 8, !tbaa !5
  %306 = load double, ptr %22, align 8, !tbaa !5
  %307 = load double, ptr %11, align 8, !tbaa !5
  %308 = load double, ptr %17, align 8, !tbaa !5
  %309 = load double, ptr %14, align 8, !tbaa !5
  %310 = call double @deltd(double noundef %303, double noundef %304, double noundef %305, double noundef %306, double noundef %307, double noundef %308, double noundef %309)
  store double %310, ptr %23, align 8, !tbaa !5
  %311 = load double, ptr %23, align 8, !tbaa !5
  %312 = load double, ptr %3, align 8, !tbaa !5
  %313 = load double, ptr %20, align 8, !tbaa !5
  %314 = load double, ptr %16, align 8, !tbaa !5
  %315 = load double, ptr %22, align 8, !tbaa !5
  %316 = load double, ptr %11, align 8, !tbaa !5
  %317 = load double, ptr %17, align 8, !tbaa !5
  %318 = load double, ptr %14, align 8, !tbaa !5
  %319 = call double @deltdd(double noundef %311, double noundef %312, double noundef %313, double noundef %314, double noundef %315, double noundef %316, double noundef %317, double noundef %318)
  store double %319, ptr %24, align 8, !tbaa !5
  %320 = load double, ptr %6, align 8, !tbaa !5
  %321 = load double, ptr %20, align 8, !tbaa !5
  %322 = load double, ptr %15, align 8, !tbaa !5
  %323 = call double @pow(double noundef %321, double noundef %322) #4
  %324 = load double, ptr %4, align 8, !tbaa !5
  %325 = load double, ptr %21, align 8, !tbaa !5
  %326 = load double, ptr %19, align 8, !tbaa !5
  %327 = call double @psit(double noundef %324, double noundef %325, double noundef %326)
  %328 = load double, ptr %3, align 8, !tbaa !5
  %329 = load double, ptr %4, align 8, !tbaa !5
  %330 = load double, ptr %3, align 8, !tbaa !5
  %331 = load double, ptr %21, align 8, !tbaa !5
  %332 = load double, ptr %18, align 8, !tbaa !5
  %333 = load double, ptr %19, align 8, !tbaa !5
  %334 = call double @psidt(double noundef %329, double noundef %330, double noundef %331, double noundef %332, double noundef %333)
  %335 = call double @llvm.fmuladd.f64(double %328, double %334, double %327)
  %336 = load double, ptr %3, align 8, !tbaa !5
  %337 = load double, ptr %15, align 8, !tbaa !5
  %338 = load double, ptr %20, align 8, !tbaa !5
  %339 = load double, ptr %23, align 8, !tbaa !5
  %340 = call double @deltabid(double noundef %337, double noundef %338, double noundef %339)
  %341 = fmul double %336, %340
  %342 = load double, ptr %4, align 8, !tbaa !5
  %343 = load double, ptr %21, align 8, !tbaa !5
  %344 = load double, ptr %19, align 8, !tbaa !5
  %345 = call double @psit(double noundef %342, double noundef %343, double noundef %344)
  %346 = fmul double %341, %345
  %347 = call double @llvm.fmuladd.f64(double %323, double %335, double %346)
  %348 = load double, ptr %15, align 8, !tbaa !5
  %349 = load double, ptr %22, align 8, !tbaa !5
  %350 = load double, ptr %20, align 8, !tbaa !5
  %351 = call double @deltabit(double noundef %348, double noundef %349, double noundef %350)
  %352 = load double, ptr %21, align 8, !tbaa !5
  %353 = load double, ptr %3, align 8, !tbaa !5
  %354 = load double, ptr %3, align 8, !tbaa !5
  %355 = load double, ptr %21, align 8, !tbaa !5
  %356 = load double, ptr %18, align 8, !tbaa !5
  %357 = call double @psid(double noundef %354, double noundef %355, double noundef %356)
  %358 = call double @llvm.fmuladd.f64(double %353, double %357, double %352)
  %359 = call double @llvm.fmuladd.f64(double %351, double %358, double %347)
  %360 = load double, ptr %16, align 8, !tbaa !5
  %361 = load double, ptr %15, align 8, !tbaa !5
  %362 = load double, ptr %11, align 8, !tbaa !5
  %363 = load double, ptr %20, align 8, !tbaa !5
  %364 = load double, ptr %3, align 8, !tbaa !5
  %365 = load double, ptr %22, align 8, !tbaa !5
  %366 = load double, ptr %23, align 8, !tbaa !5
  %367 = call double @deltabidt(double noundef %360, double noundef %361, double noundef %362, double noundef %363, double noundef %364, double noundef %365, double noundef %366)
  %368 = load double, ptr %3, align 8, !tbaa !5
  %369 = fmul double %367, %368
  %370 = load double, ptr %21, align 8, !tbaa !5
  %371 = call double @llvm.fmuladd.f64(double %369, double %370, double %359)
  %372 = load double, ptr %25, align 8, !tbaa !5
  %373 = call double @llvm.fmuladd.f64(double %320, double %371, double %372)
  store double %373, ptr %25, align 8, !tbaa !5
  br label %374

374:                                              ; preds = %221
  %375 = load i32, ptr %5, align 4, !tbaa !9
  %376 = add nsw i32 %375, 1
  store i32 %376, ptr %5, align 4, !tbaa !9
  br label %218, !llvm.loop !15

377:                                              ; preds = %218
  %378 = 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 %378
}

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 @deltabit(double noundef, double noundef, double noundef) local_unnamed_addr #2

; Function Attrs: nounwind uwtable
declare dso_local double @deltabidt(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 @psid(double noundef, double noundef, double noundef) local_unnamed_addr #2

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

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

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

.preheader151:                                    ; preds = %.preheader151.preheader, %.preheader151
  %indvars.iv = phi i64 [ 0, %.preheader151.preheader ], [ %indvars.iv.next, %.preheader151 ]
  %.1155 = phi double [ %78, %.preheader151.preheader ], [ %101, %.preheader151 ]
  %79 = getelementptr inbounds [27 x [4 x double]], ptr @coeffr2, i64 0, i64 %indvars.iv
  %80 = load double, ptr %79, align 16, !tbaa !5
  %81 = getelementptr inbounds i8, ptr %79, i64 8
  %82 = load double, ptr %81, align 8, !tbaa !5
  %83 = getelementptr inbounds i8, ptr %79, i64 16
  %84 = load double, ptr %83, align 16, !tbaa !5
  %85 = getelementptr inbounds i8, ptr %79, i64 24
  %86 = load double, ptr %85, align 8, !tbaa !5
  %87 = tail call double @pow(double noundef %0, double noundef %86) #3
  %88 = fneg double %87
  %89 = tail call double @exp(double noundef %88) #3
  %90 = fmul double %80, %89
  %91 = fadd double %82, -1.000000e+00
  %92 = tail call double @pow(double noundef %0, double noundef %91) #3
  %93 = fmul double %90, %92
  %94 = fmul double %84, %93
  %95 = fadd double %84, -1.000000e+00
  %96 = tail call double @pow(double noundef %1, double noundef %95) #3
  %97 = fmul double %96, %94
  %98 = tail call double @pow(double noundef %0, double noundef %86) #3
  %99 = fneg double %86
  %100 = tail call double @llvm.fmuladd.f64(double %99, double %98, double %82)
  %101 = tail call double @llvm.fmuladd.f64(double %97, double %100, double %.1155)
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, 27
  br i1 %exitcond.not, label %.preheader150.preheader, label %.preheader151, !llvm.loop !9

.preheader150.preheader:                          ; preds = %.preheader151
  %102 = load double, ptr @coeffr3, align 16, !tbaa !5
  %103 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 8), align 8, !tbaa !5
  %104 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 16), align 16, !tbaa !5
  %105 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 24), align 8, !tbaa !5
  %106 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 32), align 16, !tbaa !5
  %107 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 40), align 8, !tbaa !5
  %108 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 48), align 16, !tbaa !5
  %109 = tail call double @pow(double noundef %0, double noundef %103) #3
  %110 = fmul double %102, %109
  %111 = tail call double @pow(double noundef %1, double noundef %104) #3
  %112 = fmul double %110, %111
  %113 = fneg double %105
  %114 = fsub double %0, %108
  %115 = fmul double %114, %113
  %116 = fsub double %1, %107
  %117 = fneg double %106
  %118 = fmul double %116, %117
  %119 = fmul double %116, %118
  %120 = tail call double @llvm.fmuladd.f64(double %115, double %114, double %119)
  %121 = tail call double @exp(double noundef %120) #3
  %122 = fmul double %112, %121
  %123 = fdiv double %103, %0
  %124 = fmul double %105, -2.000000e+00
  %125 = tail call double @llvm.fmuladd.f64(double %124, double %114, double %123)
  %126 = fmul double %125, %122
  %127 = fdiv double %104, %1
  %128 = fmul double %106, -2.000000e+00
  %129 = tail call double @llvm.fmuladd.f64(double %128, double %116, double %127)
  %130 = tail call double @llvm.fmuladd.f64(double %126, double %129, double %101)
  %131 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 56), align 8, !tbaa !5
  %132 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 64), align 16, !tbaa !5
  %133 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 72), align 8, !tbaa !5
  %134 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 80), align 16, !tbaa !5
  %135 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 88), align 8, !tbaa !5
  %136 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 96), align 16, !tbaa !5
  %137 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 104), align 8, !tbaa !5
  %138 = tail call double @pow(double noundef %0, double noundef %132) #3
  %139 = fmul double %131, %138
  %140 = tail call double @pow(double noundef %1, double noundef %133) #3
  %141 = fmul double %139, %140
  %142 = fneg double %134
  %143 = fsub double %0, %137
  %144 = fmul double %143, %142
  %145 = fsub double %1, %136
  %146 = fneg double %135
  %147 = fmul double %145, %146
  %148 = fmul double %145, %147
  %149 = tail call double @llvm.fmuladd.f64(double %144, double %143, double %148)
  %150 = tail call double @exp(double noundef %149) #3
  %151 = fmul double %141, %150
  %152 = fdiv double %132, %0
  %153 = fmul double %134, -2.000000e+00
  %154 = tail call double @llvm.fmuladd.f64(double %153, double %143, double %152)
  %155 = fmul double %154, %151
  %156 = fdiv double %133, %1
  %157 = fmul double %135, -2.000000e+00
  %158 = tail call double @llvm.fmuladd.f64(double %157, double %145, double %156)
  %159 = tail call double @llvm.fmuladd.f64(double %155, double %158, double %130)
  %160 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 112), align 16, !tbaa !5
  %161 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 120), align 8, !tbaa !5
  %162 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 128), align 16, !tbaa !5
  %163 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 136), align 8, !tbaa !5
  %164 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 144), align 16, !tbaa !5
  %165 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 152), align 8, !tbaa !5
  %166 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 160), align 16, !tbaa !5
  %167 = tail call double @pow(double noundef %0, double noundef %161) #3
  %168 = fmul double %160, %167
  %169 = tail call double @pow(double noundef %1, double noundef %162) #3
  %170 = fmul double %168, %169
  %171 = fneg double %163
  %172 = fsub double %0, %166
  %173 = fmul double %172, %171
  %174 = fsub double %1, %165
  %175 = fneg double %164
  %176 = fmul double %174, %175
  %177 = fmul double %174, %176
  %178 = tail call double @llvm.fmuladd.f64(double %173, double %172, double %177)
  %179 = tail call double @exp(double noundef %178) #3
  %180 = fmul double %170, %179
  %181 = fdiv double %161, %0
  %182 = fmul double %163, -2.000000e+00
  %183 = tail call double @llvm.fmuladd.f64(double %182, double %172, double %181)
  %184 = fmul double %183, %180
  %185 = fdiv double %162, %1
  %186 = fmul double %164, -2.000000e+00
  %187 = tail call double @llvm.fmuladd.f64(double %186, double %174, double %185)
  %188 = tail call double @llvm.fmuladd.f64(double %184, double %187, double %159)
  %189 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 168), align 8, !tbaa !5
  %190 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 176), align 16, !tbaa !5
  %191 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 184), align 8, !tbaa !5
  %192 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 192), align 16, !tbaa !5
  %193 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 200), align 8, !tbaa !5
  %194 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 208), align 16, !tbaa !5
  %195 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 216), align 8, !tbaa !5
  %196 = tail call double @pow(double noundef %0, double noundef %190) #3
  %197 = fmul double %189, %196
  %198 = tail call double @pow(double noundef %1, double noundef %191) #3
  %199 = fmul double %197, %198
  %200 = fneg double %192
  %201 = fsub double %0, %195
  %202 = fmul double %201, %200
  %203 = fsub double %1, %194
  %204 = fneg double %193
  %205 = fmul double %203, %204
  %206 = fmul double %203, %205
  %207 = tail call double @llvm.fmuladd.f64(double %202, double %201, double %206)
  %208 = tail call double @exp(double noundef %207) #3
  %209 = fmul double %199, %208
  %210 = fdiv double %190, %0
  %211 = fmul double %192, -2.000000e+00
  %212 = tail call double @llvm.fmuladd.f64(double %211, double %201, double %210)
  %213 = fmul double %212, %209
  %214 = fdiv double %191, %1
  %215 = fmul double %193, -2.000000e+00
  %216 = tail call double @llvm.fmuladd.f64(double %215, double %203, double %214)
  %217 = tail call double @llvm.fmuladd.f64(double %213, double %216, double %188)
  %218 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 224), align 16, !tbaa !5
  %219 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 232), align 8, !tbaa !5
  %220 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 240), align 16, !tbaa !5
  %221 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 248), align 8, !tbaa !5
  %222 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 256), align 16, !tbaa !5
  %223 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 264), align 8, !tbaa !5
  %224 = load double, ptr getelementptr inbounds (i8, ptr @coeffr3, i64 272), align 16, !tbaa !5
  %225 = tail call double @pow(double noundef %0, double noundef %219) #3
  %226 = fmul double %218, %225
  %227 = tail call double @pow(double noundef %1, double noundef %220) #3
  %228 = fmul double %226, %227
  %229 = fneg double %221
  %230 = fsub double %0, %224
  %231 = fmul double %230, %229
  %232 = fsub double %1, %223
  %233 = fneg double %222
  %234 = fmul double %232, %233
  %235 = fmul double %232, %234
  %236 = tail call double @llvm.fmuladd.f64(double %231, double %230, double %235)
  %237 = tail call double @exp(double noundef %236) #3
  %238 = fmul double %228, %237
  %239 = fdiv double %219, %0
  %240 = fmul double %221, -2.000000e+00
  %241 = tail call double @llvm.fmuladd.f64(double %240, double %230, double %239)
  %242 = fmul double %241, %238
  %243 = fdiv double %220, %1
  %244 = fmul double %222, -2.000000e+00
  %245 = tail call double @llvm.fmuladd.f64(double %244, double %232, double %243)
  %246 = tail call double @llvm.fmuladd.f64(double %242, double %245, double %217)
  %247 = fsub double 1.000000e+00, %1
  %248 = fadd double %0, -1.000000e+00
  %249 = fmul double %248, %248
  %250 = fadd double %1, -1.000000e+00
  %251 = load double, ptr @coeffr4, align 16, !tbaa !5
  %252 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 8), align 8, !tbaa !5
  %253 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 16), align 16, !tbaa !5
  %254 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 24), align 8, !tbaa !5
  %255 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 32), align 16, !tbaa !5
  %256 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 40), align 8, !tbaa !5
  %257 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 48), align 16, !tbaa !5
  %258 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 56), align 8, !tbaa !5
  %259 = fdiv double 5.000000e-01, %254
  %260 = tail call double @pow(double noundef %249, double noundef %259) #3
  %261 = tail call double @llvm.fmuladd.f64(double %255, double %260, double %247)
  %262 = tail call double @pow(double noundef %249, double noundef %252) #3
  %263 = fmul double %256, %262
  %264 = tail call double @llvm.fmuladd.f64(double %261, double %261, double %263)
  %265 = fneg double %257
  %266 = fmul double %248, %265
  %267 = fneg double %258
  %268 = fmul double %250, %267
  %269 = fmul double %250, %268
  %270 = tail call double @llvm.fmuladd.f64(double %266, double %248, double %269)
  %271 = tail call double @exp(double noundef %270) #3
  %272 = tail call double @deltd(double noundef %0, double noundef %264, double noundef %255, double noundef %261, double noundef %254, double noundef %256, double noundef %252)
  %273 = tail call double @deltdd(double noundef %272, double noundef %0, double noundef %264, double noundef %255, double noundef %261, double noundef %254, double noundef %256, double noundef %252)
  %274 = tail call double @pow(double noundef %264, double noundef %253) #3
  %275 = tail call double @psit(double noundef %1, double noundef %271, double noundef %258)
  %276 = tail call double @psidt(double noundef %1, double noundef %0, double noundef %271, double noundef %257, double noundef %258)
  %277 = tail call double @llvm.fmuladd.f64(double %0, double %276, double %275)
  %278 = tail call double @deltabid(double noundef %253, double noundef %264, double noundef %272)
  %279 = fmul double %278, %0
  %280 = tail call double @psit(double noundef %1, double noundef %271, double noundef %258)
  %281 = fmul double %279, %280
  %282 = tail call double @llvm.fmuladd.f64(double %274, double %277, double %281)
  %283 = tail call double @deltabit(double noundef %253, double noundef %261, double noundef %264)
  %284 = tail call double @psid(double noundef %0, double noundef %271, double noundef %257)
  %285 = tail call double @llvm.fmuladd.f64(double %0, double %284, double %271)
  %286 = tail call double @llvm.fmuladd.f64(double %283, double %285, double %282)
  %287 = tail call double @deltabidt(double noundef %255, double noundef %253, double noundef %254, double noundef %264, double noundef %0, double noundef %261, double noundef %272)
  %288 = fmul double %287, %0
  %289 = tail call double @llvm.fmuladd.f64(double %288, double %271, double %286)
  %290 = tail call double @llvm.fmuladd.f64(double %251, double %289, double %246)
  %291 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 64), align 16, !tbaa !5
  %292 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 72), align 8, !tbaa !5
  %293 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 80), align 16, !tbaa !5
  %294 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 88), align 8, !tbaa !5
  %295 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 96), align 16, !tbaa !5
  %296 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 104), align 8, !tbaa !5
  %297 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 112), align 16, !tbaa !5
  %298 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 120), align 8, !tbaa !5
  %299 = fdiv double 5.000000e-01, %294
  %300 = tail call double @pow(double noundef %249, double noundef %299) #3
  %301 = tail call double @llvm.fmuladd.f64(double %295, double %300, double %247)
  %302 = tail call double @pow(double noundef %249, double noundef %292) #3
  %303 = fmul double %296, %302
  %304 = tail call double @llvm.fmuladd.f64(double %301, double %301, double %303)
  %305 = fneg double %297
  %306 = fmul double %248, %305
  %307 = fneg double %298
  %308 = fmul double %250, %307
  %309 = fmul double %250, %308
  %310 = tail call double @llvm.fmuladd.f64(double %306, double %248, double %309)
  %311 = tail call double @exp(double noundef %310) #3
  %312 = tail call double @deltd(double noundef %0, double noundef %304, double noundef %295, double noundef %301, double noundef %294, double noundef %296, double noundef %292)
  %313 = tail call double @deltdd(double noundef %312, double noundef %0, double noundef %304, double noundef %295, double noundef %301, double noundef %294, double noundef %296, double noundef %292)
  %314 = tail call double @pow(double noundef %304, double noundef %293) #3
  %315 = tail call double @psit(double noundef %1, double noundef %311, double noundef %298)
  %316 = tail call double @psidt(double noundef %1, double noundef %0, double noundef %311, double noundef %297, double noundef %298)
  %317 = tail call double @llvm.fmuladd.f64(double %0, double %316, double %315)
  %318 = tail call double @deltabid(double noundef %293, double noundef %304, double noundef %312)
  %319 = fmul double %318, %0
  %320 = tail call double @psit(double noundef %1, double noundef %311, double noundef %298)
  %321 = fmul double %319, %320
  %322 = tail call double @llvm.fmuladd.f64(double %314, double %317, double %321)
  %323 = tail call double @deltabit(double noundef %293, double noundef %301, double noundef %304)
  %324 = tail call double @psid(double noundef %0, double noundef %311, double noundef %297)
  %325 = tail call double @llvm.fmuladd.f64(double %0, double %324, double %311)
  %326 = tail call double @llvm.fmuladd.f64(double %323, double %325, double %322)
  %327 = tail call double @deltabidt(double noundef %295, double noundef %293, double noundef %294, double noundef %304, double noundef %0, double noundef %301, double noundef %312)
  %328 = fmul double %327, %0
  %329 = tail call double @llvm.fmuladd.f64(double %328, double %311, double %326)
  %330 = tail call double @llvm.fmuladd.f64(double %291, double %329, double %290)
  %331 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 128), align 16, !tbaa !5
  %332 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 136), align 8, !tbaa !5
  %333 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 144), align 16, !tbaa !5
  %334 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 152), align 8, !tbaa !5
  %335 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 160), align 16, !tbaa !5
  %336 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 168), align 8, !tbaa !5
  %337 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 176), align 16, !tbaa !5
  %338 = load double, ptr getelementptr inbounds (i8, ptr @coeffr4, i64 184), align 8, !tbaa !5
  %339 = fdiv double 5.000000e-01, %334
  %340 = tail call double @pow(double noundef %249, double noundef %339) #3
  %341 = tail call double @llvm.fmuladd.f64(double %335, double %340, double %247)
  %342 = tail call double @pow(double noundef %249, double noundef %332) #3
  %343 = fmul double %336, %342
  %344 = tail call double @llvm.fmuladd.f64(double %341, double %341, double %343)
  %345 = fneg double %337
  %346 = fmul double %248, %345
  %347 = fneg double %338
  %348 = fmul double %250, %347
  %349 = fmul double %250, %348
  %350 = tail call double @llvm.fmuladd.f64(double %346, double %248, double %349)
  %351 = tail call double @exp(double noundef %350) #3
  %352 = tail call double @deltd(double noundef %0, double noundef %344, double noundef %335, double noundef %341, double noundef %334, double noundef %336, double noundef %332)
  %353 = tail call double @deltdd(double noundef %352, double noundef %0, double noundef %344, double noundef %335, double noundef %341, double noundef %334, double noundef %336, double noundef %332)
  %354 = tail call double @pow(double noundef %344, double noundef %333) #3
  %355 = tail call double @psit(double noundef %1, double noundef %351, double noundef %338)
  %356 = tail call double @psidt(double noundef %1, double noundef %0, double noundef %351, double noundef %337, double noundef %338)
  %357 = tail call double @llvm.fmuladd.f64(double %0, double %356, double %355)
  %358 = tail call double @deltabid(double noundef %333, double noundef %344, double noundef %352)
  %359 = fmul double %358, %0
  %360 = tail call double @psit(double noundef %1, double noundef %351, double noundef %338)
  %361 = fmul double %359, %360
  %362 = tail call double @llvm.fmuladd.f64(double %354, double %357, double %361)
  %363 = tail call double @deltabit(double noundef %333, double noundef %341, double noundef %344)
  %364 = tail call double @psid(double noundef %0, double noundef %351, double noundef %337)
  %365 = tail call double @llvm.fmuladd.f64(double %0, double %364, double %351)
  %366 = tail call double @llvm.fmuladd.f64(double %363, double %365, double %362)
  %367 = tail call double @deltabidt(double noundef %335, double noundef %333, double noundef %334, double noundef %344, double noundef %0, double noundef %341, double noundef %352)
  %368 = fmul double %367, %0
  %369 = tail call double @llvm.fmuladd.f64(double %368, double %351, double %366)
  %370 = tail call double @llvm.fmuladd.f64(double %331, double %369, double %330)
  ret double %370
}

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