routine calcul matrice

tempdata:
    .db 0,0,0,0,0,0,0,0,0,0
        .db 0,0,0,0,0,0,0,0,0,0
GRotatesur_X:
    .db 0,0,0
gCalculMatrix:
    ld ixl,128
    ld hl,GRotatesur_X
        ld b,(hl)
        inc hl
        ld c,(hl)
        inc hl
        ld d,(hl)
        ;50

        ld hl,SinTable/256
        ld l,c
        ld a,(hl)
        ld (gMatrix+2),a    ;optimisation
    ;34
        ld a,b
        add a,c
        ld e,a
        ld l,a
        ld a,(hl)
        ld (tempdata+1),a       ;sin(x+y)
        ld a,ixl
        sub c
        ld l,a
        ld a,(hl)
        ld (tempdata),a ;cos(x+y)
        ld a,e
        add a,d
        ld l,a
        ld a,(hl)
        ld (tempdata+13),a  ;sin(x+y+z)
        ld a, ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+12),a  ;cos(x+y+z)
        ld a,e
        sub d
        ld l,a
        ld a,(hl)
        ld (tempdata+15),a  ;sin(x+y-z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+14),a  ;cos(x+y-z)
    ;196


        ld a,b
        sub c
        ld e,a
        ld l,a
        ld a,(hl)
        ld (tempdata+3),a   ;sin(x-y)
        ld a,ixl
        sub e
        ld l,a
        ld a,(hl)
        ld (tempdata+2),a   ;cos(x-y)
        ld a,e
        add d
        ld l,a
        ld a,(hl)
        ld (tempdata+17),a  ;sin(x-y+z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+16),a  ;cos(x-y+z)
        ld a,e
        sub d
        ld l,a
        ld a,(hl)
        ld (tempdata+19),a  ;sin(x-y-z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+18),a  ;cos(x-y-z)
    ;196

        ld a,b
        add a,d
        ld l,a
        ld a,(hl)
        ld (tempdata+9),a   ;sin(x+z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+8),a   ;cos(x+z)
        ld a,b
        sub d
        ld l,a
        ld a,(hl)
        ld (tempdata+11),a  ;sin(x-z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+10),a  ;cos(x-z)
    ;128
        ld a,c
        add a,d
        ld l,a
        ld a,(hl)
        ld (tempdata+5),a   ;sin(y+z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+4),a   ;cos(y+z)
        ld a,c
        sub d
        ld l,a
        ld a,(hl)
        ld (tempdata+7),a   ;sin(y-z)
        ld a,ixl
        sub l
        ld l,a
        ld a,(hl)
        ld (tempdata+6),a   ;cos(y-z)
    ;128

    ;total: 732

    ;calculs
    ld ix, tempdata
    ;14
    ;matrice 0 et 1
    ld hl,(tempdata+4)
    ld de,(tempdata+6)
    ld a,l
    add a,e
    rra
    ld l,a
    ld a,h
    add a,d
    rra
    ld h,a
    ld (gMatrix),hl
    ;84
    ;matrice 4 et 5
    ld hl,(tempdata+12)
    ld de,(tempdata+18)
    ld bc,(tempdata+10)
    ld a,l
    sub (ix+14)
    add a,(ix+16)
    sub e
    ld l,a
    ld a,b
    add a,(ix+9)
    add a,a
    ld b,a
    ld a,l
    sub b
    scf
    ccf
    rra
    rra
    ld l,a
    ld a,(tempdata+8)
    add a,c
    add a,a
    sub (ix+17)
    add a, d
    sub h
    add a,(ix+15)
    scf
    ccf
    rra
    rra
    ld h,a
    ld (gMatrix+4),hl
    ;252
    ;matrice 6 et 10
    ld hl,(tempdata)
    ld de,(tempdata+2)
    xor a
    sub d
    sub h
    rra
    ld (gMatrix+6),a
    ld a,e
    add a,l
    rra
    ld (gMatrix+10),a
    ;90
    ;matrice 8 et 9
    ld hl,(tempdata+12)
    ld de,(tempdata+18)
        ld a,(tempdata+10)
    sub (ix+8)
    add a,a
    sub h
    add a, (ix+17)
    sub(ix+15)
    add a,d
    scf
    ccf
    rra
    rra
    ld c,a
    ld a,b
    add a,(ix+14)
    sub l
    sub e
    add a,(ix+16)
    scf
    ccf
    rra
    rra
    ld h,a
    ld l,c
    ld (gMatrix+8),a
    ;211
    ret
    ;1436