; LLVM IR from chrr/libECP :: src (transformations)

; Unoptimized IR (clang -disable-llvm-passes):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

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

; Function Attrs: nounwind allocsize(0,1)
declare noalias ptr @calloc(i64 noundef, i64 noundef) #1

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

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

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

29:                                               ; preds = %55, %4
  %30 = load i32, ptr %9, align 4, !tbaa !9
  %31 = load i32, ptr %6, align 4, !tbaa !9
  %32 = icmp sle i32 %30, %31
  br i1 %32, label %33, label %58

33:                                               ; preds = %29
  %34 = load i32, ptr %9, align 4, !tbaa !9
  %35 = mul nsw i32 3, %34
  %36 = add nsw i32 %35, 2
  %37 = sitofp i32 %36 to double
  %38 = load ptr, ptr %8, align 8, !tbaa !5
  %39 = load i32, ptr %9, align 4, !tbaa !9
  %40 = add nsw i32 %39, 3
  %41 = sext i32 %40 to i64
  %42 = getelementptr inbounds double, ptr %38, i64 %41
  %43 = load double, ptr %42, align 8, !tbaa !11
  %44 = fmul double %37, %43
  %45 = load ptr, ptr %8, align 8, !tbaa !5
  %46 = load i32, ptr %9, align 4, !tbaa !9
  %47 = sext i32 %46 to i64
  %48 = getelementptr inbounds double, ptr %45, i64 %47
  %49 = load double, ptr %48, align 8, !tbaa !11
  %50 = fmul double 1.200000e+01, %49
  %51 = fdiv double %44, %50
  %52 = fptosi double %51 to i32
  %53 = load i32, ptr %26, align 4, !tbaa !9
  %54 = add nsw i32 %53, %52
  store i32 %54, ptr %26, align 4, !tbaa !9
  br label %55

55:                                               ; preds = %33
  %56 = load i32, ptr %9, align 4, !tbaa !9
  %57 = add nsw i32 %56, 1
  store i32 %57, ptr %9, align 4, !tbaa !9
  br label %29, !llvm.loop !13

58:                                               ; preds = %29
  %59 = load i32, ptr %26, align 4, !tbaa !9
  %60 = sext i32 %59 to i64
  %61 = call noalias ptr @calloc(i64 noundef %60, i64 noundef 8) #6
  store ptr %61, ptr %27, align 8, !tbaa !5
  store ptr %61, ptr %28, align 8, !tbaa !5
  store i32 0, ptr %9, align 4, !tbaa !9
  br label %62

62:                                               ; preds = %322, %58
  %63 = load i32, ptr %9, align 4, !tbaa !9
  %64 = load i32, ptr %6, align 4, !tbaa !9
  %65 = icmp sle i32 %63, %64
  br i1 %65, label %66, label %325

66:                                               ; preds = %62
  store i32 0, ptr %11, align 4, !tbaa !9
  br label %67

67:                                               ; preds = %318, %66
  %68 = load i32, ptr %11, align 4, !tbaa !9
  %69 = load i32, ptr %9, align 4, !tbaa !9
  %70 = add nsw i32 %69, 1
  %71 = load i32, ptr %9, align 4, !tbaa !9
  %72 = add nsw i32 %71, 2
  %73 = mul nsw i32 %70, %72
  %74 = sdiv i32 %73, 2
  %75 = icmp slt i32 %68, %74
  br i1 %75, label %76, label %321

76:                                               ; preds = %67
  %77 = load i32, ptr %9, align 4, !tbaa !9
  %78 = sub nsw i32 %77, 1
  %79 = add nsw i32 %78, 1
  %80 = load i32, ptr %9, align 4, !tbaa !9
  %81 = sub nsw i32 %80, 1
  %82 = add nsw i32 %81, 2
  %83 = mul nsw i32 %79, %82
  %84 = load i32, ptr %9, align 4, !tbaa !9
  %85 = sub nsw i32 %84, 1
  %86 = add nsw i32 %85, 3
  %87 = mul nsw i32 %83, %86
  %88 = sdiv i32 %87, 6
  %89 = load i32, ptr %11, align 4, !tbaa !9
  %90 = add nsw i32 %88, %89
  %91 = mul nsw i32 %90, 3
  store i32 %91, ptr %13, align 4, !tbaa !9
  %92 = load ptr, ptr %7, align 8, !tbaa !5
  %93 = load i32, ptr %13, align 4, !tbaa !9
  %94 = sext i32 %93 to i64
  %95 = getelementptr inbounds i32, ptr %92, i64 %94
  %96 = load i32, ptr %95, align 4, !tbaa !9
  store i32 %96, ptr %17, align 4, !tbaa !9
  %97 = load ptr, ptr %7, align 8, !tbaa !5
  %98 = load i32, ptr %13, align 4, !tbaa !9
  %99 = add nsw i32 %98, 1
  %100 = sext i32 %99 to i64
  %101 = getelementptr inbounds i32, ptr %97, i64 %100
  %102 = load i32, ptr %101, align 4, !tbaa !9
  store i32 %102, ptr %18, align 4, !tbaa !9
  %103 = load ptr, ptr %7, align 8, !tbaa !5
  %104 = load i32, ptr %13, align 4, !tbaa !9
  %105 = add nsw i32 %104, 2
  %106 = sext i32 %105 to i64
  %107 = getelementptr inbounds i32, ptr %103, i64 %106
  %108 = load i32, ptr %107, align 4, !tbaa !9
  store i32 %108, ptr %19, align 4, !tbaa !9
  %109 = load ptr, ptr %8, align 8, !tbaa !5
  %110 = load i32, ptr %17, align 4, !tbaa !9
  %111 = sext i32 %110 to i64
  %112 = getelementptr inbounds double, ptr %109, i64 %111
  %113 = load double, ptr %112, align 8, !tbaa !11
  %114 = load ptr, ptr %8, align 8, !tbaa !5
  %115 = load i32, ptr %18, align 4, !tbaa !9
  %116 = sext i32 %115 to i64
  %117 = getelementptr inbounds double, ptr %114, i64 %116
  %118 = load double, ptr %117, align 8, !tbaa !11
  %119 = fmul double %113, %118
  %120 = load ptr, ptr %8, align 8, !tbaa !5
  %121 = load i32, ptr %19, align 4, !tbaa !9
  %122 = sext i32 %121 to i64
  %123 = getelementptr inbounds double, ptr %120, i64 %122
  %124 = load double, ptr %123, align 8, !tbaa !11
  %125 = fmul double %119, %124
  %126 = load ptr, ptr %8, align 8, !tbaa !5
  %127 = load i32, ptr %17, align 4, !tbaa !9
  %128 = mul nsw i32 2, %127
  %129 = sext i32 %128 to i64
  %130 = getelementptr inbounds double, ptr %126, i64 %129
  %131 = load double, ptr %130, align 8, !tbaa !11
  %132 = load ptr, ptr %8, align 8, !tbaa !5
  %133 = load i32, ptr %18, align 4, !tbaa !9
  %134 = mul nsw i32 2, %133
  %135 = sext i32 %134 to i64
  %136 = getelementptr inbounds double, ptr %132, i64 %135
  %137 = load double, ptr %136, align 8, !tbaa !11
  %138 = fmul double %131, %137
  %139 = load ptr, ptr %8, align 8, !tbaa !5
  %140 = load i32, ptr %19, align 4, !tbaa !9
  %141 = mul nsw i32 2, %140
  %142 = sext i32 %141 to i64
  %143 = getelementptr inbounds double, ptr %139, i64 %142
  %144 = load double, ptr %143, align 8, !tbaa !11
  %145 = fmul double %138, %144
  %146 = fdiv double %125, %145
  %147 = call double @sqrt(double noundef %146) #5
  store double %147, ptr %24, align 8, !tbaa !11
  store i32 0, ptr %12, align 4, !tbaa !9
  br label %148

148:                                              ; preds = %314, %76
  %149 = load i32, ptr %12, align 4, !tbaa !9
  %150 = load i32, ptr %9, align 4, !tbaa !9
  %151 = add nsw i32 %150, 1
  %152 = load i32, ptr %9, align 4, !tbaa !9
  %153 = add nsw i32 %152, 2
  %154 = mul nsw i32 %151, %153
  %155 = sdiv i32 %154, 2
  %156 = icmp slt i32 %149, %155
  br i1 %156, label %157, label %317

157:                                              ; preds = %148
  %158 = load i32, ptr %9, align 4, !tbaa !9
  %159 = sub nsw i32 %158, 1
  %160 = add nsw i32 %159, 1
  %161 = load i32, ptr %9, align 4, !tbaa !9
  %162 = sub nsw i32 %161, 1
  %163 = add nsw i32 %162, 2
  %164 = mul nsw i32 %160, %163
  %165 = load i32, ptr %9, align 4, !tbaa !9
  %166 = sub nsw i32 %165, 1
  %167 = add nsw i32 %166, 3
  %168 = mul nsw i32 %164, %167
  %169 = sdiv i32 %168, 6
  %170 = load i32, ptr %12, align 4, !tbaa !9
  %171 = add nsw i32 %169, %170
  %172 = mul nsw i32 %171, 3
  store i32 %172, ptr %13, align 4, !tbaa !9
  %173 = load ptr, ptr %7, align 8, !tbaa !5
  %174 = load i32, ptr %13, align 4, !tbaa !9
  %175 = sext i32 %174 to i64
  %176 = getelementptr inbounds i32, ptr %173, i64 %175
  %177 = load i32, ptr %176, align 4, !tbaa !9
  store i32 %177, ptr %20, align 4, !tbaa !9
  %178 = load ptr, ptr %7, align 8, !tbaa !5
  %179 = load i32, ptr %13, align 4, !tbaa !9
  %180 = add nsw i32 %179, 1
  %181 = sext i32 %180 to i64
  %182 = getelementptr inbounds i32, ptr %178, i64 %181
  %183 = load i32, ptr %182, align 4, !tbaa !9
  store i32 %183, ptr %21, align 4, !tbaa !9
  %184 = load ptr, ptr %7, align 8, !tbaa !5
  %185 = load i32, ptr %13, align 4, !tbaa !9
  %186 = add nsw i32 %185, 2
  %187 = sext i32 %186 to i64
  %188 = getelementptr inbounds i32, ptr %184, i64 %187
  %189 = load i32, ptr %188, align 4, !tbaa !9
  store i32 %189, ptr %22, align 4, !tbaa !9
  %190 = load i32, ptr %17, align 4, !tbaa !9
  %191 = load i32, ptr %20, align 4, !tbaa !9
  %192 = add nsw i32 %190, %191
  store i32 %192, ptr %14, align 4, !tbaa !9
  %193 = load i32, ptr %18, align 4, !tbaa !9
  %194 = load i32, ptr %21, align 4, !tbaa !9
  %195 = add nsw i32 %193, %194
  store i32 %195, ptr %15, align 4, !tbaa !9
  %196 = load i32, ptr %19, align 4, !tbaa !9
  %197 = load i32, ptr %22, align 4, !tbaa !9
  %198 = add nsw i32 %196, %197
  store i32 %198, ptr %16, align 4, !tbaa !9
  %199 = load i32, ptr %14, align 4, !tbaa !9
  %200 = srem i32 %199, 2
  %201 = icmp eq i32 %200, 0
  br i1 %201, label %202, label %309

202:                                              ; preds = %157
  %203 = load i32, ptr %15, align 4, !tbaa !9
  %204 = srem i32 %203, 2
  %205 = icmp eq i32 %204, 0
  br i1 %205, label %206, label %309

206:                                              ; preds = %202
  %207 = load i32, ptr %16, align 4, !tbaa !9
  %208 = srem i32 %207, 2
  %209 = icmp eq i32 %208, 0
  br i1 %209, label %210, label %309

210:                                              ; preds = %206
  %211 = load ptr, ptr %8, align 8, !tbaa !5
  %212 = load i32, ptr %20, align 4, !tbaa !9
  %213 = sext i32 %212 to i64
  %214 = getelementptr inbounds double, ptr %211, i64 %213
  %215 = load double, ptr %214, align 8, !tbaa !11
  %216 = load ptr, ptr %8, align 8, !tbaa !5
  %217 = load i32, ptr %21, align 4, !tbaa !9
  %218 = sext i32 %217 to i64
  %219 = getelementptr inbounds double, ptr %216, i64 %218
  %220 = load double, ptr %219, align 8, !tbaa !11
  %221 = fmul double %215, %220
  %222 = load ptr, ptr %8, align 8, !tbaa !5
  %223 = load i32, ptr %22, align 4, !tbaa !9
  %224 = sext i32 %223 to i64
  %225 = getelementptr inbounds double, ptr %222, i64 %224
  %226 = load double, ptr %225, align 8, !tbaa !11
  %227 = fmul double %221, %226
  %228 = load ptr, ptr %8, align 8, !tbaa !5
  %229 = load i32, ptr %20, align 4, !tbaa !9
  %230 = mul nsw i32 2, %229
  %231 = sext i32 %230 to i64
  %232 = getelementptr inbounds double, ptr %228, i64 %231
  %233 = load double, ptr %232, align 8, !tbaa !11
  %234 = load ptr, ptr %8, align 8, !tbaa !5
  %235 = load i32, ptr %21, align 4, !tbaa !9
  %236 = mul nsw i32 2, %235
  %237 = sext i32 %236 to i64
  %238 = getelementptr inbounds double, ptr %234, i64 %237
  %239 = load double, ptr %238, align 8, !tbaa !11
  %240 = fmul double %233, %239
  %241 = load ptr, ptr %8, align 8, !tbaa !5
  %242 = load i32, ptr %22, align 4, !tbaa !9
  %243 = mul nsw i32 2, %242
  %244 = sext i32 %243 to i64
  %245 = getelementptr inbounds double, ptr %241, i64 %244
  %246 = load double, ptr %245, align 8, !tbaa !11
  %247 = fmul double %240, %246
  %248 = fdiv double %227, %247
  %249 = call double @sqrt(double noundef %248) #5
  store double %249, ptr %25, align 8, !tbaa !11
  %250 = load ptr, ptr %8, align 8, !tbaa !5
  %251 = load i32, ptr %14, align 4, !tbaa !9
  %252 = sext i32 %251 to i64
  %253 = getelementptr inbounds double, ptr %250, i64 %252
  %254 = load double, ptr %253, align 8, !tbaa !11
  %255 = load ptr, ptr %8, align 8, !tbaa !5
  %256 = load i32, ptr %15, align 4, !tbaa !9
  %257 = sext i32 %256 to i64
  %258 = getelementptr inbounds double, ptr %255, i64 %257
  %259 = load double, ptr %258, align 8, !tbaa !11
  %260 = fmul double %254, %259
  %261 = load ptr, ptr %8, align 8, !tbaa !5
  %262 = load i32, ptr %16, align 4, !tbaa !9
  %263 = sext i32 %262 to i64
  %264 = getelementptr inbounds double, ptr %261, i64 %263
  %265 = load double, ptr %264, align 8, !tbaa !11
  %266 = fmul double %260, %265
  %267 = load double, ptr %24, align 8, !tbaa !11
  %268 = fmul double %266, %267
  %269 = load double, ptr %25, align 8, !tbaa !11
  %270 = fmul double %268, %269
  %271 = load ptr, ptr %8, align 8, !tbaa !5
  %272 = load i32, ptr %14, align 4, !tbaa !9
  %273 = sdiv i32 %272, 2
  %274 = sext i32 %273 to i64
  %275 = getelementptr inbounds double, ptr %271, i64 %274
  %276 = load double, ptr %275, align 8, !tbaa !11
  %277 = load ptr, ptr %8, align 8, !tbaa !5
  %278 = load i32, ptr %15, align 4, !tbaa !9
  %279 = sdiv i32 %278, 2
  %280 = sext i32 %279 to i64
  %281 = getelementptr inbounds double, ptr %277, i64 %280
  %282 = load double, ptr %281, align 8, !tbaa !11
  %283 = fmul double %276, %282
  %284 = load ptr, ptr %8, align 8, !tbaa !5
  %285 = load i32, ptr %16, align 4, !tbaa !9
  %286 = sdiv i32 %285, 2
  %287 = sext i32 %286 to i64
  %288 = getelementptr inbounds double, ptr %284, i64 %287
  %289 = load double, ptr %288, align 8, !tbaa !11
  %290 = fmul double %283, %289
  %291 = fdiv double %270, %290
  store double %291, ptr %23, align 8, !tbaa !11
  %292 = load i32, ptr %9, align 4, !tbaa !9
  %293 = sub nsw i32 0, %292
  store i32 %293, ptr %10, align 4, !tbaa !9
  br label %294

294:                                              ; preds = %305, %210
  %295 = load i32, ptr %10, align 4, !tbaa !9
  %296 = load i32, ptr %9, align 4, !tbaa !9
  %297 = icmp sle i32 %295, %296
  br i1 %297, label %298, label %308

298:                                              ; preds = %294
  %299 = load double, ptr %23, align 8, !tbaa !11
  %300 = load ptr, ptr %5, align 8, !tbaa !5
  %301 = load double, ptr %300, align 8, !tbaa !11
  %302 = load ptr, ptr %28, align 8, !tbaa !5
  %303 = load double, ptr %302, align 8, !tbaa !11
  %304 = call double @llvm.fmuladd.f64(double %299, double %301, double %303)
  store double %304, ptr %302, align 8, !tbaa !11
  br label %305

305:                                              ; preds = %298
  %306 = load i32, ptr %10, align 4, !tbaa !9
  %307 = add nsw i32 %306, 1
  store i32 %307, ptr %10, align 4, !tbaa !9
  br label %294, !llvm.loop !15

308:                                              ; preds = %294
  br label %309

309:                                              ; preds = %308, %206, %202, %157
  %310 = load ptr, ptr %28, align 8, !tbaa !5
  %311 = getelementptr inbounds double, ptr %310, i32 1
  store ptr %311, ptr %28, align 8, !tbaa !5
  %312 = load ptr, ptr %5, align 8, !tbaa !5
  %313 = getelementptr inbounds double, ptr %312, i32 1
  store ptr %313, ptr %5, align 8, !tbaa !5
  br label %314

314:                                              ; preds = %309
  %315 = load i32, ptr %12, align 4, !tbaa !9
  %316 = add nsw i32 %315, 1
  store i32 %316, ptr %12, align 4, !tbaa !9
  br label %148, !llvm.loop !16

317:                                              ; preds = %148
  br label %318

318:                                              ; preds = %317
  %319 = load i32, ptr %11, align 4, !tbaa !9
  %320 = add nsw i32 %319, 1
  store i32 %320, ptr %11, align 4, !tbaa !9
  br label %67, !llvm.loop !17

321:                                              ; preds = %67
  br label %322

322:                                              ; preds = %321
  %323 = load i32, ptr %9, align 4, !tbaa !9
  %324 = add nsw i32 %323, 1
  store i32 %324, ptr %9, align 4, !tbaa !9
  br label %62, !llvm.loop !18

325:                                              ; preds = %62
  %326 = load ptr, ptr %27, align 8, !tbaa !5
  call void @llvm.lifetime.end.p0(i64 8, ptr %28) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %27) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %26) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %25) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %24) #5
  call void @llvm.lifetime.end.p0(i64 8, ptr %23) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %22) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %21) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %20) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %19) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %18) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %17) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %16) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %15) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %14) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %13) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %12) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %11) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %10) #5
  call void @llvm.lifetime.end.p0(i64 4, ptr %9) #5
  ret ptr %326
}

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

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

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


; Optimized IR (opt -passes='default<O3>'):
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-i128:128-f80:128-n8:16:32:64-S128"
target triple = "x86_64-unknown-linux-gnu"

; Function Attrs: mustprogress nofree nounwind willreturn allockind("alloc,zeroed") allocsize(0,1) memory(inaccessiblemem: readwrite)
declare noalias noundef ptr @calloc(i64 noundef, i64 noundef) local_unnamed_addr #0

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

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

; Function Attrs: nofree nounwind uwtable
define dso_local noalias noundef ptr @TM_sph2cart(ptr nocapture noundef readonly %0, i32 noundef %1, ptr nocapture noundef readonly %2, ptr nocapture noundef readonly %3) local_unnamed_addr #3 {
  %.not114 = icmp slt i32 %1, 0
  br i1 %.not114, label %._crit_edge.thread, label %.lr.ph.preheader

._crit_edge.thread:                               ; preds = %4
  %5 = tail call noalias ptr @calloc(i64 noundef 0, i64 noundef 8) #5
  br label %._crit_edge137

.lr.ph.preheader:                                 ; preds = %4
  %6 = add nuw i32 %1, 1
  %wide.trip.count = zext i32 %6 to i64
  %min.iters.check = icmp ult i32 %1, 3
  br i1 %min.iters.check, label %.lr.ph.preheader168, label %vector.ph

vector.ph:                                        ; preds = %.lr.ph.preheader
  %n.vec = and i64 %wide.trip.count, 4294967292
  br label %vector.body

vector.body:                                      ; preds = %vector.body, %vector.ph
  %index = phi i64 [ 0, %vector.ph ], [ %index.next, %vector.body ]
  %vec.phi = phi <2 x i32> [ zeroinitializer, %vector.ph ], [ %25, %vector.body ]
  %vec.phi161 = phi <2 x i32> [ zeroinitializer, %vector.ph ], [ %26, %vector.body ]
  %vec.ind = phi <2 x i32> [ <i32 0, i32 1>, %vector.ph ], [ %vec.ind.next, %vector.body ]
  %7 = mul <2 x i32> %vec.ind, <i32 3, i32 3>
  %8 = mul <2 x i32> %vec.ind, <i32 3, i32 3>
  %9 = add <2 x i32> %7, <i32 2, i32 2>
  %10 = add <2 x i32> %8, <i32 8, i32 8>
  %11 = uitofp <2 x i32> %9 to <2 x double>
  %12 = uitofp <2 x i32> %10 to <2 x double>
  %13 = getelementptr inbounds double, ptr %3, i64 %index
  %14 = getelementptr inbounds i8, ptr %13, i64 24
  %15 = getelementptr inbounds i8, ptr %13, i64 40
  %wide.load = load <2 x double>, ptr %14, align 8, !tbaa !5
  %wide.load163 = load <2 x double>, ptr %15, align 8, !tbaa !5
  %16 = fmul <2 x double> %wide.load, %11
  %17 = fmul <2 x double> %wide.load163, %12
  %18 = getelementptr inbounds i8, ptr %13, i64 16
  %wide.load164 = load <2 x double>, ptr %13, align 8, !tbaa !5
  %wide.load165 = load <2 x double>, ptr %18, align 8, !tbaa !5
  %19 = fmul <2 x double> %wide.load164, <double 1.200000e+01, double 1.200000e+01>
  %20 = fmul <2 x double> %wide.load165, <double 1.200000e+01, double 1.200000e+01>
  %21 = fdiv <2 x double> %16, %19
  %22 = fdiv <2 x double> %17, %20
  %23 = fptosi <2 x double> %21 to <2 x i32>
  %24 = fptosi <2 x double> %22 to <2 x i32>
  %25 = add <2 x i32> %vec.phi, %23
  %26 = add <2 x i32> %vec.phi161, %24
  %index.next = add nuw i64 %index, 4
  %vec.ind.next = add <2 x i32> %vec.ind, <i32 4, i32 4>
  %27 = icmp eq i64 %index.next, %n.vec
  br i1 %27, label %middle.block, label %vector.body, !llvm.loop !9

middle.block:                                     ; preds = %vector.body
  %bin.rdx = add <2 x i32> %26, %25
  %28 = tail call i32 @llvm.vector.reduce.add.v2i32(<2 x i32> %bin.rdx)
  %cmp.n = icmp eq i64 %n.vec, %wide.trip.count
  br i1 %cmp.n, label %._crit_edge, label %.lr.ph.preheader168

.lr.ph.preheader168:                              ; preds = %middle.block, %.lr.ph.preheader
  %indvars.iv.ph = phi i64 [ 0, %.lr.ph.preheader ], [ %n.vec, %middle.block ]
  %.0100116.ph = phi i32 [ 0, %.lr.ph.preheader ], [ %28, %middle.block ]
  br label %.lr.ph

.lr.ph:                                           ; preds = %.lr.ph.preheader168, %.lr.ph
  %indvars.iv = phi i64 [ %indvars.iv.next, %.lr.ph ], [ %indvars.iv.ph, %.lr.ph.preheader168 ]
  %.0100116 = phi i32 [ %41, %.lr.ph ], [ %.0100116.ph, %.lr.ph.preheader168 ]
  %29 = trunc i64 %indvars.iv to i32
  %30 = mul i32 %29, 3
  %31 = add i32 %30, 2
  %32 = uitofp nneg i32 %31 to double
  %33 = getelementptr inbounds double, ptr %3, i64 %indvars.iv
  %34 = getelementptr inbounds i8, ptr %33, i64 24
  %35 = load double, ptr %34, align 8, !tbaa !5
  %36 = fmul double %35, %32
  %37 = load double, ptr %33, align 8, !tbaa !5
  %38 = fmul double %37, 1.200000e+01
  %39 = fdiv double %36, %38
  %40 = fptosi double %39 to i32
  %41 = add nsw i32 %.0100116, %40
  %indvars.iv.next = add nuw nsw i64 %indvars.iv, 1
  %exitcond.not = icmp eq i64 %indvars.iv.next, %wide.trip.count
  br i1 %exitcond.not, label %._crit_edge, label %.lr.ph, !llvm.loop !13

._crit_edge:                                      ; preds = %.lr.ph, %middle.block
  %.lcssa160 = phi i32 [ %28, %middle.block ], [ %41, %.lr.ph ]
  %42 = sext i32 %.lcssa160 to i64
  %43 = tail call noalias ptr @calloc(i64 noundef %42, i64 noundef 8) #5
  br i1 %.not114, label %._crit_edge137, label %.lr.ph129

.loopexit113:                                     ; preds = %180
  %indvars.iv.next143 = add nuw nsw i32 %indvars.iv142, 1
  %exitcond151.not = icmp eq i32 %.1102134, %1
  br i1 %exitcond151.not, label %._crit_edge137, label %.lr.ph129, !llvm.loop !14

.lr.ph129:                                        ; preds = %._crit_edge, %.loopexit113
  %indvars.iv142 = phi i32 [ %indvars.iv.next143, %.loopexit113 ], [ 1, %._crit_edge ]
  %.0136 = phi ptr [ %177, %.loopexit113 ], [ %43, %._crit_edge ]
  %.097135 = phi ptr [ %178, %.loopexit113 ], [ %0, %._crit_edge ]
  %.1102134 = phi i32 [ %45, %.loopexit113 ], [ 0, %._crit_edge ]
  %44 = shl nuw i32 %.1102134, 1
  %45 = add nuw nsw i32 %.1102134, 1
  %46 = add nuw nsw i32 %.1102134, 2
  %47 = mul nuw nsw i32 %45, %46
  %48 = lshr i32 %47, 1
  %factor.op.mul = mul i32 %46, %.1102134
  %.reass = mul i32 %factor.op.mul, %45
  %49 = udiv i32 %.reass, 6
  %50 = sub nsw i32 0, %.1102134
  %51 = zext nneg i32 %49 to i64
  %52 = zext nneg i32 %48 to i64
  %53 = and i32 %44, 6
  %54 = icmp ult i32 %.1102134, 4
  br label %55

55:                                               ; preds = %.lr.ph129, %180
  %indvars.iv148 = phi i64 [ 0, %.lr.ph129 ], [ %indvars.iv.next149, %180 ]
  %.1128 = phi ptr [ %.0136, %.lr.ph129 ], [ %177, %180 ]
  %.198127 = phi ptr [ %.097135, %.lr.ph129 ], [ %178, %180 ]
  %56 = add nuw nsw i64 %indvars.iv148, %51
  %.idx = mul i64 %56, 12
  %57 = getelementptr inbounds i8, ptr %2, i64 %.idx
  %58 = load i32, ptr %57, align 4, !tbaa !15
  %59 = getelementptr inbounds i8, ptr %57, i64 4
  %60 = load i32, ptr %59, align 4, !tbaa !15
  %61 = getelementptr inbounds i8, ptr %57, i64 8
  %62 = load i32, ptr %61, align 4, !tbaa !15
  %63 = sext i32 %58 to i64
  %64 = getelementptr inbounds double, ptr %3, i64 %63
  %65 = load double, ptr %64, align 8, !tbaa !5
  %66 = sext i32 %60 to i64
  %67 = getelementptr inbounds double, ptr %3, i64 %66
  %68 = load double, ptr %67, align 8, !tbaa !5
  %69 = fmul double %65, %68
  %70 = sext i32 %62 to i64
  %71 = getelementptr inbounds double, ptr %3, i64 %70
  %72 = load double, ptr %71, align 8, !tbaa !5
  %73 = fmul double %69, %72
  %74 = shl nsw i32 %58, 1
  %75 = sext i32 %74 to i64
  %76 = getelementptr inbounds double, ptr %3, i64 %75
  %77 = load double, ptr %76, align 8, !tbaa !5
  %78 = shl nsw i32 %60, 1
  %79 = sext i32 %78 to i64
  %80 = getelementptr inbounds double, ptr %3, i64 %79
  %81 = load double, ptr %80, align 8, !tbaa !5
  %82 = fmul double %77, %81
  %83 = shl nsw i32 %62, 1
  %84 = sext i32 %83 to i64
  %85 = getelementptr inbounds double, ptr %3, i64 %84
  %86 = load double, ptr %85, align 8, !tbaa !5
  %87 = fmul double %82, %86
  %88 = fdiv double %73, %87
  %89 = tail call double @sqrt(double noundef %88) #6
  br label %90

90:                                               ; preds = %55, %.loopexit
  %indvars.iv145 = phi i64 [ 0, %55 ], [ %indvars.iv.next146, %.loopexit ]
  %.2124 = phi ptr [ %.1128, %55 ], [ %177, %.loopexit ]
  %.299123 = phi ptr [ %.198127, %55 ], [ %178, %.loopexit ]
  %91 = add nuw nsw i64 %indvars.iv145, %51
  %.idx153 = mul i64 %91, 12
  %92 = getelementptr inbounds i8, ptr %2, i64 %.idx153
  %93 = load i32, ptr %92, align 4, !tbaa !15
  %94 = getelementptr inbounds i8, ptr %92, i64 4
  %95 = load i32, ptr %94, align 4, !tbaa !15
  %96 = getelementptr inbounds i8, ptr %92, i64 8
  %97 = load i32, ptr %96, align 4, !tbaa !15
  %98 = add nsw i32 %93, %58
  %99 = add nsw i32 %95, %60
  %100 = add nsw i32 %97, %62
  %101 = and i32 %98, 1
  %102 = icmp eq i32 %101, 0
  %103 = and i32 %99, 1
  %104 = icmp eq i32 %103, 0
  %or.cond = select i1 %102, i1 %104, i1 false
  %105 = and i32 %100, 1
  %106 = icmp eq i32 %105, 0
  %or.cond112 = select i1 %or.cond, i1 %106, i1 false
  br i1 %or.cond112, label %.lr.ph120, label %.loopexit

.lr.ph120:                                        ; preds = %90
  %107 = sext i32 %93 to i64
  %108 = getelementptr inbounds double, ptr %3, i64 %107
  %109 = load double, ptr %108, align 8, !tbaa !5
  %110 = sext i32 %95 to i64
  %111 = getelementptr inbounds double, ptr %3, i64 %110
  %112 = load double, ptr %111, align 8, !tbaa !5
  %113 = fmul double %109, %112
  %114 = sext i32 %97 to i64
  %115 = getelementptr inbounds double, ptr %3, i64 %114
  %116 = load double, ptr %115, align 8, !tbaa !5
  %117 = fmul double %113, %116
  %118 = shl nsw i32 %93, 1
  %119 = sext i32 %118 to i64
  %120 = getelementptr inbounds double, ptr %3, i64 %119
  %121 = load double, ptr %120, align 8, !tbaa !5
  %122 = shl nsw i32 %95, 1
  %123 = sext i32 %122 to i64
  %124 = getelementptr inbounds double, ptr %3, i64 %123
  %125 = load double, ptr %124, align 8, !tbaa !5
  %126 = fmul double %121, %125
  %127 = shl nsw i32 %97, 1
  %128 = sext i32 %127 to i64
  %129 = getelementptr inbounds double, ptr %3, i64 %128
  %130 = load double, ptr %129, align 8, !tbaa !5
  %131 = fmul double %126, %130
  %132 = fdiv double %117, %131
  %133 = tail call double @sqrt(double noundef %132) #6
  %134 = sext i32 %98 to i64
  %135 = getelementptr inbounds double, ptr %3, i64 %134
  %136 = load double, ptr %135, align 8, !tbaa !5
  %137 = sext i32 %99 to i64
  %138 = getelementptr inbounds double, ptr %3, i64 %137
  %139 = load double, ptr %138, align 8, !tbaa !5
  %140 = fmul double %136, %139
  %141 = sext i32 %100 to i64
  %142 = getelementptr inbounds double, ptr %3, i64 %141
  %143 = load double, ptr %142, align 8, !tbaa !5
  %144 = fmul double %140, %143
  %145 = fmul double %89, %144
  %146 = fmul double %133, %145
  %147 = ashr exact i32 %98, 1
  %148 = sext i32 %147 to i64
  %149 = getelementptr inbounds double, ptr %3, i64 %148
  %150 = load double, ptr %149, align 8, !tbaa !5
  %151 = ashr exact i32 %99, 1
  %152 = sext i32 %151 to i64
  %153 = getelementptr inbounds double, ptr %3, i64 %152
  %154 = load double, ptr %153, align 8, !tbaa !5
  %155 = fmul double %150, %154
  %156 = ashr exact i32 %100, 1
  %157 = sext i32 %156 to i64
  %158 = getelementptr inbounds double, ptr %3, i64 %157
  %159 = load double, ptr %158, align 8, !tbaa !5
  %160 = fmul double %155, %159
  %161 = fdiv double %146, %160
  %162 = load double, ptr %.299123, align 8, !tbaa !5
  %.2.promoted = load double, ptr %.2124, align 8, !tbaa !5
  br label %163

163:                                              ; preds = %163, %.lr.ph120
  %164 = phi double [ %.2.promoted, %.lr.ph120 ], [ %165, %163 ]
  %.0103118.prol = phi i32 [ %50, %.lr.ph120 ], [ %166, %163 ]
  %prol.iter = phi i32 [ 0, %.lr.ph120 ], [ %prol.iter.next, %163 ]
  %165 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %164)
  %166 = add i32 %.0103118.prol, 1
  %prol.iter.next = add i32 %prol.iter, 1
  %prol.iter.cmp.not = icmp eq i32 %prol.iter, %53
  br i1 %prol.iter.cmp.not, label %.prol.loopexit, label %163, !llvm.loop !17

.prol.loopexit:                                   ; preds = %163
  br i1 %54, label %..loopexit_crit_edge, label %.lr.ph120.new

.lr.ph120.new:                                    ; preds = %.prol.loopexit, %.lr.ph120.new
  %167 = phi double [ %175, %.lr.ph120.new ], [ %165, %.prol.loopexit ]
  %.0103118 = phi i32 [ %176, %.lr.ph120.new ], [ %166, %.prol.loopexit ]
  %168 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %167)
  %169 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %168)
  %170 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %169)
  %171 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %170)
  %172 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %171)
  %173 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %172)
  %174 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %173)
  %175 = tail call double @llvm.fmuladd.f64(double %161, double %162, double %174)
  %176 = add i32 %.0103118, 8
  %exitcond144.not.7 = icmp eq i32 %176, %indvars.iv142
  br i1 %exitcond144.not.7, label %..loopexit_crit_edge, label %.lr.ph120.new, !llvm.loop !19

..loopexit_crit_edge:                             ; preds = %.lr.ph120.new, %.prol.loopexit
  %.lcssa = phi double [ %165, %.prol.loopexit ], [ %175, %.lr.ph120.new ]
  store double %.lcssa, ptr %.2124, align 8, !tbaa !5
  br label %.loopexit

.loopexit:                                        ; preds = %..loopexit_crit_edge, %90
  %177 = getelementptr inbounds i8, ptr %.2124, i64 8
  %178 = getelementptr inbounds i8, ptr %.299123, i64 8
  %indvars.iv.next146 = add nuw nsw i64 %indvars.iv145, 1
  %179 = icmp ult i64 %indvars.iv.next146, %52
  br i1 %179, label %90, label %180, !llvm.loop !20

180:                                              ; preds = %.loopexit
  %indvars.iv.next149 = add nuw nsw i64 %indvars.iv148, 1
  %181 = icmp ult i64 %indvars.iv.next149, %52
  br i1 %181, label %55, label %.loopexit113, !llvm.loop !21

._crit_edge137:                                   ; preds = %.loopexit113, %._crit_edge.thread, %._crit_edge
  %182 = phi ptr [ %5, %._crit_edge.thread ], [ %43, %._crit_edge ], [ %43, %.loopexit113 ]
  ret ptr %182
}

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare i32 @llvm.vector.reduce.add.v2i32(<2 x i32>) #4

attributes #0 = { mustprogress nofree nounwind willreturn allockind("alloc,zeroed") allocsize(0,1) memory(inaccessiblemem: readwrite) "alloc-family"="malloc" "no-trapping-math"="true" "stack-protector-buffer-size"="8" "target-cpu"="x86-64" "target-features"="+cmov,+cx8,+fxsr,+mmx,+sse,+sse2,+x87" "tune-cpu"="generic" }
attributes #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 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #3 = { nofree 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 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #5 = { nounwind allocsize(0,1) }
attributes #6 = { 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, !11, !12}
!10 = !{!"llvm.loop.mustprogress"}
!11 = !{!"llvm.loop.isvectorized", i32 1}
!12 = !{!"llvm.loop.unroll.runtime.disable"}
!13 = distinct !{!13, !10, !12, !11}
!14 = distinct !{!14, !10}
!15 = !{!16, !16, i64 0}
!16 = !{!"int", !7, i64 0}
!17 = distinct !{!17, !18}
!18 = !{!"llvm.loop.unroll.disable"}
!19 = distinct !{!19, !10}
!20 = distinct !{!20, !10}
!21 = distinct !{!21, !10}


; Effective O3 passes: EarlyCSEPass,GVNPass,GlobalOptPass,IndVarSimplifyPass,InferFunctionAttrsPass,InstCombinePass,InstSimplifyPass,JumpThreadingPass,LCSSAPass,LICMPass,LoopInstSimplifyPass,LoopLoadEliminationPass,LoopRotatePass,LoopSimplifyPass,LoopUnrollPass,LoopVectorizePass,PostOrderFunctionAttrsPass,ReassociatePass,SCCPPass,SROAPass,SimplifyCFGPass,TailCallElimPass