e4438939-2b38-4b8f-a77c-1cb58c8d1d58 91 4 5 false false 1 decimal

&[DATE] &[TIME] - &[FILENAME]

B #

Intersection of two cylinders

Draghilev method (n = 3 X0 tmin tmax N D3 k 13 range : u k el x k el : b 1 - (f u Jacobian 1233 submatrix * : 13 mat A f u Jacobian 1212 submatrix : Δ 3 el A det - : 12 mat Δ 1 el b A 1222 submatrix augment det - : Δ 2 el A 1211 submatrix b augment det - : D k el Δ k el Δ norme / : t x D D : 13 mat X0 0 tmax N t x D rkfixed 51 line :

b m t f k 1 b rows range L1 k el m k row b k row stack π 36 / t * Ω2 γ * * eval : S k el L1 k el 1 col L1 k el 2 col augment 1010^1010^1010^1010^1010^16 mat transpose eval : k 1 ≡ F S k el eval : F F S k el stack eval : if 31 line for F 21 line :

f x 1 el (2^x 3 el (2^+ 4 - x 1 el (1 - (2^x 2 el (2^+ 1 - 21 mat : X01 21013 -^013 mat transpose : X0 00213 mat transpose : 13 mat

x0 t2 N F x0 0 t2 N D3 1 N 24 submatrix :

B X01 4.1 - 41 F :

B1 X01 8.1 - 81 F :

B2 X0 8.1 - 81 F :

B3 X01 8.1 81 F :

B4 X01 4.1 41 F :

γ 0.5757 - 0.1732 0.7991 0.5677 0.6187 - 0.543 0.5885 0.7663 0.2579 33 mat :

Матрица преобразования для поворота вокруг оси y angle Ω2 angle cos 0 angle sin 010 angle sin - 0 angle cos 33 mat :

t ints1 B1 π 36 / t * Ω2 γ * * 1 col B1 π 36 / t * Ω2 γ * * 2 col augment :

t ints2 B3 π 36 / t * Ω2 * γ * 1 col B3 π 36 / t * Ω2 * γ * 2 col augment :

t Lints t ints1 t ints2 stack eval :

M B 1 col B rows 1 matrix 3.2 + B 3 col augment :

M2 B2 1 col B2 rows 1 matrix 3.2 - B2 3 col augment :

M4 B3 1 col B3 rows 1 matrix 3.2 + B3 3 col augment :

M1 B1 1 col B1 2 col B1 rows 1 matrix 3.2 + augment :

M3 B3 1 col B3 2 col B3 rows 1 matrix 3.2 - augment :

t A1 t Lints 11 mat B2 M2 t f B M t f B4 M4 t f stack eval mat2sys.1 eval :

t A2 t Lints 11 mat B1 M1 t f B3 M3 t f B M t f B2 M2 t f B4 M4 t f stack eval mat2sys.1 eval :

step1 072 range :

t A1 t A2