int psede_transf_ones_apply(double *target, int size, int stride, int howmany,
                            int dist, void *params) {
  int i, j;
  for (i = 0; i < howmany; ++i) {
    for (j = 0; j < size; ++j) {
      target[i * dist + j * stride] = 1.0;
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<psede_transf_ones_apply>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1184 <psede_transf_ones_apply+0x84>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1171 <psede_transf_ones_apply+0x71>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xea2(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1133 <psede_transf_ones_apply+0x33>
jmp    1176 <psede_transf_ones_apply+0x76>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1120 <psede_transf_ones_apply+0x20>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<psede_transf_ones_apply>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %esi,-0x64(%rsp)
mov    %rdi,-0x50(%rsp)
test   %ecx,%ecx
jle    136f <psede_transf_ones_apply+0x26f>
mov    %r8d,%eax
movslq %edx,%rdx
movslq %r8d,%rsi
mov    %ecx,%eax
mov    %rax,-0x10(%rsp)
mov    -0x64(%rsp),%eax
mov    %eax,%ecx
mov    %ecx,%edi
and    $0xfffffffc,%edi
mov    %rdi,-0x28(%rsp)
add    $0xfffffffffffffffc,%rdi
mov    %rdi,-0x30(%rsp)
shr    $0x2,%rdi
add    $0x1,%rdi
cmp    $0x4,%eax
setae  %al
cmp    $0x1,%edx
sete   %bl
and    %al,%bl
mov    %bl,-0x65(%rsp)
mov    %edi,%eax
and    $0x3,%eax
mov    %rax,-0x48(%rsp)
and    $0xfffffffffffffffc,%rdi
mov    %rdi,-0x40(%rsp)
mov    %rcx,-0x58(%rsp)
mov    %ecx,%r15d
and    $0x3,%r15d
mov    %rsi,-0x8(%rsp)
lea    0x0(,%rsi,8),%rax
mov    %rax,-0x18(%rsp)
mov    %rdx,%r8
shl    $0x7,%r8
mov    %rdx,%rbp
shl    $0x5,%rbp
mov    -0x50(%rsp),%r10
lea    0x10(%r10),%rax
mov    %rax,-0x38(%rsp)
mov    %r15,-0x20(%rsp)
neg    %r15
movaps 0xe46(%rip),%xmm0
mov    %rdx,-0x60(%rsp)
lea    0x0(,%rdx,8),%r11
xor    %edi,%edi
movabs $0x3ff0000000000000,%r14
xor    %r9d,%r9d
jmp    11f9 <psede_transf_ones_apply+0xf9>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r9
add    -0x18(%rsp),%r10
add    -0x8(%rsp),%rdi
cmp    -0x10(%rsp),%r9
je     136f <psede_transf_ones_apply+0x26f>
cmpl   $0x0,-0x64(%rsp)
jle    11e0 <psede_transf_ones_apply+0xe0>
cmpb   $0x0,-0x65(%rsp)
je     1220 <psede_transf_ones_apply+0x120>
cmpq   $0xc,-0x30(%rsp)
jae    1228 <psede_transf_ones_apply+0x128>
xor    %ecx,%ecx
jmp    1279 <psede_transf_ones_apply+0x179>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %r12d,%r12d
jmp    12c3 <psede_transf_ones_apply+0x1c3>
mov    -0x40(%rsp),%rbx
mov    %r10,%rsi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups %xmm0,(%rsi)
movups %xmm0,0x10(%rsi)
lea    (%rsi,%rbp,1),%rax
movups %xmm0,(%rsi,%rbp,1)
movups %xmm0,0x10(%rsi,%rbp,1)
lea    (%rax,%rbp,1),%rdx
movups %xmm0,0x0(%rbp,%rax,1)
movups %xmm0,0x10(%rbp,%rax,1)
movups %xmm0,0x0(%rbp,%rdx,1)
movups %xmm0,0x10(%rbp,%rdx,1)
add    $0x10,%rcx
add    %r8,%rsi
add    $0xfffffffffffffffc,%rbx
jne    1240 <psede_transf_ones_apply+0x140>
cmpq   $0x0,-0x48(%rsp)
je     12b0 <psede_transf_ones_apply+0x1b0>
imul   -0x60(%rsp),%rcx
add    %rdi,%rcx
mov    -0x38(%rsp),%rax
lea    (%rax,%rcx,8),%rax
mov    -0x48(%rsp),%rcx
nopl   0x0(%rax,%rax,1)
movups %xmm0,-0x10(%rax)
movups %xmm0,(%rax)
add    %rbp,%rax
add    $0xffffffffffffffff,%rcx
jne    12a0 <psede_transf_ones_apply+0x1a0>
mov    -0x28(%rsp),%rax
mov    %rax,%r12
cmp    -0x58(%rsp),%rax
je     11e0 <psede_transf_ones_apply+0xe0>
mov    %r12,%rcx
not    %rcx
add    -0x58(%rsp),%rcx
cmpq   $0x0,-0x20(%rsp)
je     1302 <psede_transf_ones_apply+0x202>
mov    -0x60(%rsp),%rax
imul   %r12,%rax
add    %rdi,%rax
mov    -0x50(%rsp),%rdx
lea    (%rdx,%rax,8),%rdx
xor    %eax,%eax
nopl   (%rax)
mov    %r14,(%rdx)
add    $0xffffffffffffffff,%rax
add    %r11,%rdx
cmp    %rax,%r15
jne    12f0 <psede_transf_ones_apply+0x1f0>
sub    %rax,%r12
cmp    $0x3,%rcx
jb     11e0 <psede_transf_ones_apply+0xe0>
mov    -0x58(%rsp),%rcx
sub    %r12,%rcx
lea    0x3(%r12),%rax
mov    -0x60(%rsp),%rdx
imul   %rdx,%rax
lea    (%r10,%rax,8),%r13
lea    0x2(%r12),%rax
imul   %rdx,%rax
lea    (%r10,%rax,8),%rsi
mov    %rdx,%rax
imul   %r12,%rax
lea    (%r10,%rax,8),%rbx
add    $0x1,%r12
imul   %rdx,%r12
lea    (%r10,%r12,8),%rdx
xor    %eax,%eax
nopl   0x0(%rax)
mov    %r14,(%rbx,%rax,1)
mov    %r14,(%rdx,%rax,1)
mov    %r14,(%rsi,%rax,1)
mov    %r14,0x0(%r13,%rax,1)
add    %rbp,%rax
add    $0xfffffffffffffffc,%rcx
jne    1350 <psede_transf_ones_apply+0x250>
jmp    11e0 <psede_transf_ones_apply+0xe0>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @psede_transf_ones_apply(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, i8* %5) {
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i8*
%13 = alloca i32
%14 = alloca i32
store double* %0, double** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store i8* %5, i8** %12
store i32 0, i32* %13
br label %15
15: ; preds = %39, %6
%16 = load i32, i32* %13
%17 = load i32, i32* %10
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %42
19: ; preds = %15
store i32 0, i32* %14
br label %20
20: ; preds = %35, %19
%21 = load i32, i32* %14
%22 = load i32, i32* %8
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %38
24: ; preds = %20
%25 = load double*, double** %7
%26 = load i32, i32* %13
%27 = load i32, i32* %11
%28 = mul nsw i32 %26, %27
%29 = load i32, i32* %14
%30 = load i32, i32* %9
%31 = mul nsw i32 %29, %30
%32 = add nsw i32 %28, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %25, i64 %33
store double 1.000000e+00, double* %34
br label %35
35: ; preds = %24
%36 = load i32, i32* %14
%37 = add nsw i32 %36, 1
store i32 %37, i32* %14
br label %20
38: ; preds = %20
br label %39
39: ; preds = %38
%40 = load i32, i32* %13
%41 = add nsw i32 %40, 1
store i32 %41, i32* %13
br label %15
42: ; preds = %15
ret i32 0
}

LLVM IR (-O2):
define i32 @psede_transf_ones_apply(double* nocapture writeonly %0, i32 %1, i32 %2, i32 %3, i32 %4, i8* nocapture readnone %5) {
%7 = icmp sgt i32 %3, 0
br i1 %7, label %8, label %120
8: ; preds = %6
%9 = icmp sgt i32 %1, 0
%10 = sext i32 %2 to i64
%11 = sext i32 %4 to i64
%12 = zext i32 %3 to i64
%13 = zext i32 %1 to i64
%14 = and i64 %13, 4294967292
%15 = add nsw i64 %14, -4
%16 = lshr exact i64 %15, 2
%17 = add nuw nsw i64 %16, 1
%18 = icmp ugt i32 %1, 3
%19 = icmp eq i32 %2, 1
%20 = and i1 %18, %19
%21 = and i64 %13, 4294967292
%22 = and i64 %17, 3
%23 = icmp ult i64 %15, 12
%24 = and i64 %17, 9223372036854775804
%25 = icmp eq i64 %22, 0
%26 = icmp eq i64 %21, %13
%27 = and i64 %13, 3
%28 = icmp eq i64 %27, 0
br label %29
29: ; preds = %8, %117
%30 = phi i64 [ 0, %8 ], [ %118, %117 ]
br i1 %9, label %31, label %117
31: ; preds = %29
%32 = mul nsw i64 %30, %11
br i1 %20, label %33, label %82
33: ; preds = %31
br i1 %23, label %67, label %34
34: ; preds = %33, %34
%35 = phi i64 [ %64, %34 ], [ 0, %33 ]
%36 = phi i64 [ %65, %34 ], [ 0, %33 ]
%37 = mul nsw i64 %35, %10
%38 = add nsw i64 %37, %32
%39 = getelementptr inbounds double, double* %0, i64 %38
%40 = bitcast double* %39 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %40
%41 = getelementptr inbounds double, double* %39, i64 2
%42 = bitcast double* %41 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %42
%43 = or i64 %35, 4
%44 = mul nsw i64 %43, %10
%45 = add nsw i64 %44, %32
%46 = getelementptr inbounds double, double* %0, i64 %45
%47 = bitcast double* %46 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %47
%48 = getelementptr inbounds double, double* %46, i64 2
%49 = bitcast double* %48 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %49
%50 = or i64 %35, 8
%51 = mul nsw i64 %50, %10
%52 = add nsw i64 %51, %32
%53 = getelementptr inbounds double, double* %0, i64 %52
%54 = bitcast double* %53 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %54
%55 = getelementptr inbounds double, double* %53, i64 2
%56 = bitcast double* %55 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %56
%57 = or i64 %35, 12
%58 = mul nsw i64 %57, %10
%59 = add nsw i64 %58, %32
%60 = getelementptr inbounds double, double* %0, i64 %59
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %61
%62 = getelementptr inbounds double, double* %60, i64 2
%63 = bitcast double* %62 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %63
%64 = add nuw i64 %35, 16
%65 = add nuw i64 %36, 4
%66 = icmp eq i64 %65, %24
br i1 %66, label %67, label %34
67: ; preds = %34, %33
%68 = phi i64 [ 0, %33 ], [ %64, %34 ]
br i1 %25, label %81, label %69
69: ; preds = %67, %69
%70 = phi i64 [ %78, %69 ], [ %68, %67 ]
%71 = phi i64 [ %79, %69 ], [ 0, %67 ]
%72 = mul nsw i64 %70, %10
%73 = add nsw i64 %72, %32
%74 = getelementptr inbounds double, double* %0, i64 %73
%75 = bitcast double* %74 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %75
%76 = getelementptr inbounds double, double* %74, i64 2
%77 = bitcast double* %76 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %77
%78 = add nuw i64 %70, 4
%79 = add i64 %71, 1
%80 = icmp eq i64 %79, %22
br i1 %80, label %81, label %69
81: ; preds = %69, %67
br i1 %26, label %117, label %82
82: ; preds = %31, %81
%83 = phi i64 [ 0, %31 ], [ %21, %81 ]
%84 = xor i64 %83, -1
%85 = add nsw i64 %84, %13
br i1 %28, label %95, label %86
86: ; preds = %82, %86
%87 = phi i64 [ %92, %86 ], [ %83, %82 ]
%88 = phi i64 [ %93, %86 ], [ 0, %82 ]
%89 = mul nsw i64 %87, %10
%90 = add nsw i64 %89, %32
%91 = getelementptr inbounds double, double* %0, i64 %90
store double 1.000000e+00, double* %91
%92 = add nuw nsw i64 %87, 1
%93 = add i64 %88, 1
%94 = icmp eq i64 %93, %27
br i1 %94, label %95, label %86
95: ; preds = %86, %82
%96 = phi i64 [ %83, %82 ], [ %92, %86 ]
%97 = icmp ult i64 %85, 3
br i1 %97, label %117, label %98
98: ; preds = %95, %98
%99 = phi i64 [ %115, %98 ], [ %96, %95 ]
%100 = mul nsw i64 %99, %10
%101 = add nsw i64 %100, %32
%102 = getelementptr inbounds double, double* %0, i64 %101
store double 1.000000e+00, double* %102
%103 = add nuw nsw i64 %99, 1
%104 = mul nsw i64 %103, %10
%105 = add nsw i64 %104, %32
%106 = getelementptr inbounds double, double* %0, i64 %105
store double 1.000000e+00, double* %106
%107 = add nuw nsw i64 %99, 2
%108 = mul nsw i64 %107, %10
%109 = add nsw i64 %108, %32
%110 = getelementptr inbounds double, double* %0, i64 %109
store double 1.000000e+00, double* %110
%111 = add nuw nsw i64 %99, 3
%112 = mul nsw i64 %111, %10
%113 = add nsw i64 %112, %32
%114 = getelementptr inbounds double, double* %0, i64 %113
store double 1.000000e+00, double* %114
%115 = add nuw nsw i64 %99, 4
%116 = icmp eq i64 %115, %13
br i1 %116, label %117, label %98
117: ; preds = %95, %98, %81, %29
%118 = add nuw nsw i64 %30, 1
%119 = icmp eq i64 %118, %12
br i1 %119, label %120, label %29
120: ; preds = %117, %6
ret i32 0
}