void laplacian(const double *v1, double *v2, const int dim_m, const int dim_n) {
  int offset = dim_m + 1;
  int my_rank = 0;
  int num_threads = 1;
#pragma omp parallel private(my_rank, num_threads)
  {
    int n_block = num_threads;
    int n_block_size = (int)ceil((float)(dim_n - 2) / n_block);
    int n_offset = my_rank * n_block_size;
    n_block_size = (((dim_n - 2 - my_rank * n_block_size) < (n_block_size))
                        ? (dim_n - 2 - my_rank * n_block_size)
                        : (n_block_size));
    laplacian_v6(v1 + offset + n_offset * dim_m, v2 + offset + n_offset * dim_m,
                 dim_m - 2, n_block_size, dim_m);
    asm volatile("mfence" ::: "memory");
  }
}


x86 assembly (clang, -O0):
<laplacian>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%eax
sub    $0x2,%eax
cvtsi2ss %eax,%xmm0
cvtsi2ssl -0x28(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
imul   -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%eax
sub    $0x2,%eax
mov    -0x20(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
sub    %ecx,%eax
cmp    -0x2c(%rbp),%eax
jge    11af <laplacian+0x8f>
mov    -0x18(%rbp),%eax
sub    $0x2,%eax
mov    -0x20(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
jmp    11b5 <laplacian+0x95>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x30(%rbp),%eax
imul   -0x14(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x10(%rbp),%rsi
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x30(%rbp),%eax
imul   -0x14(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x14(%rbp),%edx
sub    $0x2,%edx
mov    -0x2c(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    $0x0,%al
call   1040 <laplacian_v6@plt>
mfence
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<laplacian>:
push   %rax
add    $0xfffffffe,%ecx
cvtsi2ss %ecx,%xmm0
cvttss2si %xmm0,%eax
cmp    %eax,%ecx
cmovge %eax,%ecx
movslq %edx,%r8
lea    (%rdi,%r8,8),%rdi
add    $0x8,%rdi
lea    (%rsi,%r8,8),%rsi
add    $0x8,%rsi
lea    -0x2(%r8),%edx
xor    %eax,%eax
call   1030 <laplacian_v6@plt>
mfence
pop    %rax
ret

LLVM IR (-O0):
define void @laplacian(double* %0, double* %1, i32 %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%15 = load i32, i32* %7
%16 = add nsw i32 %15, 1
store i32 %16, i32* %9
store i32 0, i32* %10
store i32 1, i32* %11
%17 = load i32, i32* %11
store i32 %17, i32* %12
%18 = load i32, i32* %8
%19 = sub nsw i32 %18, 2
%20 = sitofp i32 %19 to float
%21 = load i32, i32* %12
%22 = sitofp i32 %21 to float
%23 = fdiv float %20, %22
%24 = fpext float %23 to double
%25 = call double @llvm.ceil.f64(double %24)
%26 = fptosi double %25 to i32
store i32 %26, i32* %13
%27 = load i32, i32* %10
%28 = load i32, i32* %13
%29 = mul nsw i32 %27, %28
store i32 %29, i32* %14
%30 = load i32, i32* %8
%31 = sub nsw i32 %30, 2
%32 = load i32, i32* %10
%33 = load i32, i32* %13
%34 = mul nsw i32 %32, %33
%35 = sub nsw i32 %31, %34
%36 = load i32, i32* %13
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %45
38: ; preds = %4
%39 = load i32, i32* %8
%40 = sub nsw i32 %39, 2
%41 = load i32, i32* %10
%42 = load i32, i32* %13
%43 = mul nsw i32 %41, %42
%44 = sub nsw i32 %40, %43
br label %47
45: ; preds = %4
%46 = load i32, i32* %13
br label %47
47: ; preds = %45, %38
%48 = phi i32 [ %44, %38 ], [ %46, %45 ]
store i32 %48, i32* %13
%49 = load double*, double** %5
%50 = load i32, i32* %9
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load i32, i32* %14
%54 = load i32, i32* %7
%55 = mul nsw i32 %53, %54
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %52, i64 %56
%58 = load double*, double** %6
%59 = load i32, i32* %9
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load i32, i32* %14
%63 = load i32, i32* %7
%64 = mul nsw i32 %62, %63
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %61, i64 %65
%67 = load i32, i32* %7
%68 = sub nsw i32 %67, 2
%69 = load i32, i32* %13
%70 = load i32, i32* %7
%71 = call i32 (double*, double*, i32, i32, i32, ...) bitcast (i32 (...)* @laplacian_v6 to i32 (double*, double*, i32, i32, i32, ...)*)(double* %57, double* %66, i32 %68, i32 %69, i32 %70)
ret void
}

LLVM IR (-O2):
define void @laplacian(double* %0, double* %1, i32 %2, i32 %3) {
%5 = add nsw i32 %3, -2
%6 = sitofp i32 %5 to float
%7 = fptosi float %6 to i32
%8 = icmp slt i32 %5, %7
%9 = select i1 %8, i32 %5, i32 %7
%10 = add nsw i32 %2, 1
%11 = sext i32 %10 to i64
%12 = getelementptr inbounds double, double* %0, i64 %11
%13 = getelementptr inbounds double, double* %1, i64 %11
%14 = add nsw i32 %2, -2
%15 = tail call i32 (double*, double*, i32, i32, i32, ...) bitcast (i32 (...)* @laplacian_v6 to i32 (double*, double*, i32, i32, i32, ...)*)(double* %12, double* %13, i32 %14, i32 %9, i32 %2)
tail call void asm sideeffect "mfence", "~{memory},~{dirflag},~{fpsr},~{flags}"()
ret void
}