#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ mat444Multiply(double **, double **, double **);
void mat44Isometry(double rot[3][3], double trans[3], double isom[4][4]) {
  double trans44[4][4] = {1, 0, 0, trans[0], 0, 1, 0, trans[1],
                          0, 0, 1, trans[2], 0, 0, 0, 1};
  double rot44[4][4] = {
      rot[0][0], rot[0][1], rot[0][2], 0, rot[1][0], rot[1][1], rot[1][2], 0,
      rot[2][0], rot[2][1], rot[2][2], 0, 0,         0,         0,         1};
  mat444Multiply(trans44, rot44, isom);
}


x86 assembly (clang, -O0):
<mat44Isometry>:
push   %rbp
mov    %rsp,%rbp
sub    $0x120,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
lea    -0xa0(%rbp),%rdi
movsd  0xeba(%rip),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
add    $0x20,%rdi
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  0xe78(%rip),%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
add    $0x20,%rdi
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xe41(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
add    $0x20,%rdi
xor    %esi,%esi
mov    $0x20,%edx
call   1040 <memset@plt>
mov    0x2e47(%rip),%rax
mov    %rax,-0x40(%rbp)
mov    0x2e44(%rip),%rax
mov    %rax,-0x38(%rbp)
mov    0x2e41(%rip),%rax
mov    %rax,-0x30(%rbp)
mov    0x2e3e(%rip),%rax
mov    %rax,-0x28(%rbp)
lea    -0x120(%rbp),%rdi
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
add    $0x20,%rdi
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
add    $0x20,%rdi
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
add    $0x20,%rdi
xor    %esi,%esi
mov    $0x20,%edx
call   1040 <memset@plt>
mov    0x2d63(%rip),%rax
mov    %rax,-0xc0(%rbp)
mov    0x2d5d(%rip),%rax
mov    %rax,-0xb8(%rbp)
mov    0x2d57(%rip),%rax
mov    %rax,-0xb0(%rbp)
mov    0x2d51(%rip),%rax
mov    %rax,-0xa8(%rbp)
lea    -0xa0(%rbp),%rdi
lea    -0x120(%rbp),%rsi
mov    -0x18(%rbp),%rdx
call   1030 <mat444Multiply@plt>
add    $0x120,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mat44Isometry>:
sub    $0x108,%rsp
movabs $0x3ff0000000000000,%rax
mov    %rax,0x80(%rsp)
xorps  %xmm0,%xmm0
movups %xmm0,0x88(%rsp)
movsd  (%rsi),%xmm1
movsd  %xmm1,0x98(%rsp)
xorps  %xmm1,%xmm1
movhps 0xeb5(%rip),%xmm1
movaps %xmm1,0xa0(%rsp)
movq   $0x0,0xb0(%rsp)
movsd  0x8(%rsi),%xmm1
movsd  %xmm1,0xb8(%rsp)
movaps %xmm0,0xc0(%rsp)
mov    %rax,0xd0(%rsp)
movsd  0x10(%rsi),%xmm0
movsd  %xmm0,0xd8(%rsp)
movups 0xe76(%rip),%xmm0
movaps %xmm0,0xe0(%rsp)
movups 0xe77(%rip),%xmm1
movaps %xmm1,0xf0(%rsp)
movups (%rdi),%xmm2
movaps %xmm2,(%rsp)
movsd  0x10(%rdi),%xmm2
movsd  %xmm2,0x10(%rsp)
movq   $0x0,0x18(%rsp)
movups 0x18(%rdi),%xmm2
movaps %xmm2,0x20(%rsp)
movsd  0x28(%rdi),%xmm2
movsd  %xmm2,0x30(%rsp)
movq   $0x0,0x38(%rsp)
movups 0x30(%rdi),%xmm2
movaps %xmm2,0x40(%rsp)
movsd  0x40(%rdi),%xmm2
movsd  %xmm2,0x50(%rsp)
movq   $0x0,0x58(%rsp)
movaps %xmm0,0x60(%rsp)
movaps %xmm1,0x70(%rsp)
lea    0x80(%rsp),%rdi
mov    %rsp,%rsi
call   1030 <mat444Multiply@plt>
add    $0x108,%rsp
ret

LLVM IR (-O0):
@constinit = private global [4 x double] [double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.000000e+00]
@constinit.1 = private global [4 x double] [double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.000000e+00]
define void @mat44Isometry([3 x double]* %0, double* %1, [4 x double]* %2) {
%4 = alloca [3 x double]*
%5 = alloca double*
%6 = alloca [4 x double]*
%7 = alloca [4 x [4 x double]]
%8 = alloca [4 x [4 x double]]
store [3 x double]* %0, [3 x double]** %4
store double* %1, double** %5
store [4 x double]* %2, [4 x double]** %6
%9 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %7, i64 0, i64 0
%10 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 0, i64 0
store double 1.000000e+00, double* %10
%11 = getelementptr inbounds double, double* %10, i64 1
store double 0.000000e+00, double* %11
%12 = getelementptr inbounds double, double* %11, i64 1
store double 0.000000e+00, double* %12
%13 = getelementptr inbounds double, double* %12, i64 1
%14 = load double*, double** %5
%15 = getelementptr inbounds double, double* %14, i64 0
%16 = load double, double* %15
store double %16, double* %13
%17 = getelementptr inbounds [4 x double], [4 x double]* %9, i64 1
%18 = getelementptr inbounds [4 x double], [4 x double]* %17, i64 0, i64 0
store double 0.000000e+00, double* %18
%19 = getelementptr inbounds double, double* %18, i64 1
store double 1.000000e+00, double* %19
%20 = getelementptr inbounds double, double* %19, i64 1
store double 0.000000e+00, double* %20
%21 = getelementptr inbounds double, double* %20, i64 1
%22 = load double*, double** %5
%23 = getelementptr inbounds double, double* %22, i64 1
%24 = load double, double* %23
store double %24, double* %21
%25 = getelementptr inbounds [4 x double], [4 x double]* %17, i64 1
%26 = getelementptr inbounds [4 x double], [4 x double]* %25, i64 0, i64 0
store double 0.000000e+00, double* %26
%27 = getelementptr inbounds double, double* %26, i64 1
store double 0.000000e+00, double* %27
%28 = getelementptr inbounds double, double* %27, i64 1
store double 1.000000e+00, double* %28
%29 = getelementptr inbounds double, double* %28, i64 1
%30 = load double*, double** %5
%31 = getelementptr inbounds double, double* %30, i64 2
%32 = load double, double* %31
store double %32, double* %29
%33 = getelementptr inbounds [4 x double], [4 x double]* %25, i64 1
%34 = bitcast [4 x double]* %33 to i8*
%35 = getelementptr inbounds [4 x double], [4 x double]* %33, i64 0, i64 0
%36 = bitcast [4 x double]* %33 to i8*
%37 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %8, i64 0, i64 0
%38 = getelementptr inbounds [4 x double], [4 x double]* %37, i64 0, i64 0
%39 = load [3 x double]*, [3 x double]** %4
%40 = getelementptr inbounds [3 x double], [3 x double]* %39, i64 0
%41 = getelementptr inbounds [3 x double], [3 x double]* %40, i64 0, i64 0
%42 = load double, double* %41
store double %42, double* %38
%43 = getelementptr inbounds double, double* %38, i64 1
%44 = load [3 x double]*, [3 x double]** %4
%45 = getelementptr inbounds [3 x double], [3 x double]* %44, i64 0
%46 = getelementptr inbounds [3 x double], [3 x double]* %45, i64 0, i64 1
%47 = load double, double* %46
store double %47, double* %43
%48 = getelementptr inbounds double, double* %43, i64 1
%49 = load [3 x double]*, [3 x double]** %4
%50 = getelementptr inbounds [3 x double], [3 x double]* %49, i64 0
%51 = getelementptr inbounds [3 x double], [3 x double]* %50, i64 0, i64 2
%52 = load double, double* %51
store double %52, double* %48
%53 = getelementptr inbounds double, double* %48, i64 1
store double 0.000000e+00, double* %53
%54 = getelementptr inbounds [4 x double], [4 x double]* %37, i64 1
%55 = getelementptr inbounds [4 x double], [4 x double]* %54, i64 0, i64 0
%56 = load [3 x double]*, [3 x double]** %4
%57 = getelementptr inbounds [3 x double], [3 x double]* %56, i64 1
%58 = getelementptr inbounds [3 x double], [3 x double]* %57, i64 0, i64 0
%59 = load double, double* %58
store double %59, double* %55
%60 = getelementptr inbounds double, double* %55, i64 1
%61 = load [3 x double]*, [3 x double]** %4
%62 = getelementptr inbounds [3 x double], [3 x double]* %61, i64 1
%63 = getelementptr inbounds [3 x double], [3 x double]* %62, i64 0, i64 1
%64 = load double, double* %63
store double %64, double* %60
%65 = getelementptr inbounds double, double* %60, i64 1
%66 = load [3 x double]*, [3 x double]** %4
%67 = getelementptr inbounds [3 x double], [3 x double]* %66, i64 1
%68 = getelementptr inbounds [3 x double], [3 x double]* %67, i64 0, i64 2
%69 = load double, double* %68
store double %69, double* %65
%70 = getelementptr inbounds double, double* %65, i64 1
store double 0.000000e+00, double* %70
%71 = getelementptr inbounds [4 x double], [4 x double]* %54, i64 1
%72 = getelementptr inbounds [4 x double], [4 x double]* %71, i64 0, i64 0
%73 = load [3 x double]*, [3 x double]** %4
%74 = getelementptr inbounds [3 x double], [3 x double]* %73, i64 2
%75 = getelementptr inbounds [3 x double], [3 x double]* %74, i64 0, i64 0
%76 = load double, double* %75
store double %76, double* %72
%77 = getelementptr inbounds double, double* %72, i64 1
%78 = load [3 x double]*, [3 x double]** %4
%79 = getelementptr inbounds [3 x double], [3 x double]* %78, i64 2
%80 = getelementptr inbounds [3 x double], [3 x double]* %79, i64 0, i64 1
%81 = load double, double* %80
store double %81, double* %77
%82 = getelementptr inbounds double, double* %77, i64 1
%83 = load [3 x double]*, [3 x double]** %4
%84 = getelementptr inbounds [3 x double], [3 x double]* %83, i64 2
%85 = getelementptr inbounds [3 x double], [3 x double]* %84, i64 0, i64 2
%86 = load double, double* %85
store double %86, double* %82
%87 = getelementptr inbounds double, double* %82, i64 1
store double 0.000000e+00, double* %87
%88 = getelementptr inbounds [4 x double], [4 x double]* %71, i64 1
%89 = bitcast [4 x double]* %88 to i8*
%90 = getelementptr inbounds [4 x double], [4 x double]* %88, i64 0, i64 0
%91 = bitcast [4 x double]* %88 to i8*
%92 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %7, i64 0, i64 0
%93 = bitcast [4 x double]* %92 to double**
%94 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %8, i64 0, i64 0
%95 = bitcast [4 x double]* %94 to double**
%96 = load [4 x double]*, [4 x double]** %6
%97 = bitcast [4 x double]* %96 to double**
%98 = call i32 @mat444Multiply(double** %93, double** %95, double** %97)
ret void
}

LLVM IR (-O2):
@constinit.1 = [4 x double] [double 0.000000e+00, double 0.000000e+00, double 0.000000e+00, double 1.000000e+00]
define void @mat44Isometry([3 x double]* nocapture %0, double* nocapture %1, [4 x double]* %2) {
%4 = alloca [4 x [4 x double]]
%5 = alloca [4 x [4 x double]]
%6 = bitcast [4 x [4 x double]]* %4 to i8*
%7 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 0, i64 0
store double 1.000000e+00, double* %7
%8 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 0, i64 1
%9 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 0, i64 3
%10 = bitcast double* %8 to i8*
%11 = load double, double* %1
store double %11, double* %9
%12 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 1, i64 0
%13 = bitcast double* %12 to <2 x double>*
store <2 x double> <double 0.000000e+00, double 1.000000e+00>, <2 x double>* %13
%14 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 1, i64 2
store double 0.000000e+00, double* %14
%15 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 1, i64 3
%16 = getelementptr inbounds double, double* %1, i64 1
%17 = load double, double* %16
store double %17, double* %15
%18 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 2, i64 0
%19 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 2, i64 2
%20 = bitcast double* %18 to i8*
store double 1.000000e+00, double* %19
%21 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 2, i64 3
%22 = getelementptr inbounds double, double* %1, i64 2
%23 = load double, double* %22
store double %23, double* %21
%24 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %4, i64 0, i64 3
%25 = bitcast [4 x double]* %24 to i8*
%26 = bitcast [4 x [4 x double]]* %5 to i8*
%27 = bitcast [3 x double]* %0 to <2 x double>*
%28 = load <2 x double>, <2 x double>* %27
%29 = bitcast [4 x [4 x double]]* %5 to <2 x double>*
store <2 x double> %28, <2 x double>* %29
%30 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 0, i64 2
%31 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 0, i64 2
%32 = load double, double* %31
store double %32, double* %30
%33 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 0, i64 3
store double 0.000000e+00, double* %33
%34 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 1, i64 0
%35 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 1, i64 0
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = bitcast double* %34 to <2 x double>*
store <2 x double> %37, <2 x double>* %38
%39 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 1, i64 2
%40 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 1, i64 2
%41 = load double, double* %40
store double %41, double* %39
%42 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 1, i64 3
store double 0.000000e+00, double* %42
%43 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 2, i64 0
%44 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 2, i64 0
%45 = bitcast double* %44 to <2 x double>*
%46 = load <2 x double>, <2 x double>* %45
%47 = bitcast double* %43 to <2 x double>*
store <2 x double> %46, <2 x double>* %47
%48 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 2, i64 2
%49 = getelementptr inbounds [3 x double], [3 x double]* %0, i64 2, i64 2
%50 = load double, double* %49
store double %50, double* %48
%51 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 2, i64 3
%52 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %5, i64 0, i64 3
%53 = bitcast [4 x double]* %52 to i8*
%54 = bitcast double* %51 to i64*
store i64 0, i64* %54
%55 = bitcast [4 x [4 x double]]* %4 to double**
%56 = bitcast [4 x [4 x double]]* %5 to double**
%57 = bitcast [4 x double]* %2 to double**
%58 = call i32 @mat444Multiply(double** %55, double** %56, double** %57)
ret void
}