int dicom_hdr(short int *src, int nhdr, int *ncol, int *nrow, int *npix) {
  int i;
  *ncol = src[223 + (nhdr - 496) / 2];
  *nrow = src[218 + (nhdr - 496) / 2];
  *npix = (*ncol) * (*nrow);
  return (2 * (*npix) + nhdr);
}


x86 assembly (clang, -O0):
<dicom_hdr>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1f0,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
add    $0xdf,%ecx
movslq %ecx,%rcx
movswl (%rax,%rcx,2),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x8(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1f0,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x38(%rbp),%rax
add    $0xda,%ecx
movslq %ecx,%rcx
movswl (%rax,%rcx,2),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
mov    (%rax),%ecx
mov    -0x20(%rbp),%rax
imul   (%rax),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
shl    $0x1,%eax
add    -0xc(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<dicom_hdr>:
mov    %esi,%eax
add    $0xfffffe10,%eax
shr    $0x1f,%eax
add    %esi,%eax
add    $0xfffffe10,%eax
sar    %eax
movslq %eax,%r9
movswl 0x1be(%rdi,%r9,2),%eax
mov    %eax,(%rdx)
movswl 0x1b4(%rdi,%r9,2),%eax
mov    %eax,(%rcx)
imul   (%rdx),%eax
mov    %eax,(%r8)
lea    (%rsi,%rax,2),%eax
ret

LLVM IR (-O0):
define i32 @dicom_hdr(i16* %0, i32 %1, i32* %2, i32* %3, i32* %4) {
%6 = alloca i16*
%7 = alloca i32
%8 = alloca i32*
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca i32
store i16* %0, i16** %6
store i32 %1, i32* %7
store i32* %2, i32** %8
store i32* %3, i32** %9
store i32* %4, i32** %10
%12 = load i16*, i16** %6
%13 = load i32, i32* %7
%14 = sub nsw i32 %13, 496
%15 = sdiv i32 %14, 2
%16 = add nsw i32 223, %15
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds i16, i16* %12, i64 %17
%19 = load i16, i16* %18
%20 = sext i16 %19 to i32
%21 = load i32*, i32** %8
store i32 %20, i32* %21
%22 = load i16*, i16** %6
%23 = load i32, i32* %7
%24 = sub nsw i32 %23, 496
%25 = sdiv i32 %24, 2
%26 = add nsw i32 218, %25
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds i16, i16* %22, i64 %27
%29 = load i16, i16* %28
%30 = sext i16 %29 to i32
%31 = load i32*, i32** %9
store i32 %30, i32* %31
%32 = load i32*, i32** %8
%33 = load i32, i32* %32
%34 = load i32*, i32** %9
%35 = load i32, i32* %34
%36 = mul nsw i32 %33, %35
%37 = load i32*, i32** %10
store i32 %36, i32* %37
%38 = load i32*, i32** %10
%39 = load i32, i32* %38
%40 = mul nsw i32 2, %39
%41 = load i32, i32* %7
%42 = add nsw i32 %40, %41
ret i32 %42
}

LLVM IR (-O2):
define i32 @dicom_hdr(i16* nocapture %0, i32 %1, i32* nocapture %2, i32* nocapture writeonly %3, i32* nocapture writeonly %4) {
%6 = add nsw i32 %1, -496
%7 = sdiv i32 %6, 2
%8 = add nsw i32 %7, 223
%9 = sext i32 %8 to i64
%10 = getelementptr inbounds i16, i16* %0, i64 %9
%11 = load i16, i16* %10
%12 = sext i16 %11 to i32
store i32 %12, i32* %2
%13 = add nsw i32 %7, 218
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds i16, i16* %0, i64 %14
%16 = load i16, i16* %15
%17 = sext i16 %16 to i32
store i32 %17, i32* %3
%18 = load i32, i32* %2
%19 = mul nsw i32 %18, %17
store i32 %19, i32* %4
%20 = shl nsw i32 %19, 1
%21 = add nsw i32 %20, %1
ret i32 %21
}