; スケルトンからアウトラインへの変換を行なう |
; スケルトンからアウトラインへの変換を行なう |
; |
; |
|
|
|
; |
|
; pointsのn番目を取り,floatに変換する |
|
; |
(defun point-n (n points) |
(defun point-n (n points) |
(let ((point (nth n points))) |
(let ((point (nth n points))) |
`(,(float (car point)),(float (cadr point)) .,(cddr point)))) |
`(,(float (car point)),(float (cadr point)) .,(cddr point)))) |
|
|
|
; |
|
; points全体をfloatに変換する |
|
; |
(defun floatlist (list) |
(defun floatlist (list) |
(mapcar list |
(mapcar list |
(function (lambda (x) `(,(float (car x)),(float (cadr x)).,(cddr x)))))) |
(function (lambda (x) `(,(float (car x)),(float (cadr x)).,(cddr x)))))) |
|
; |
|
; 見てのとおり |
|
; |
(defun appendrev (a b) (append a (reverse b))) |
(defun appendrev (a b) (append a (reverse b))) |
|
; |
|
; |
|
; |
(defun skeleton2list (l tag) |
(defun skeleton2list (l tag) |
(setq l (rm-limit l)) |
(setq l (rm-limit l)) |
(let ((func (get-def 'typehook tag))) |
(let ((func (get-def 'typehook tag))) |
(*$ 0.5 sqrt51 costheta) |
(*$ 0.5 sqrt51 costheta) |
(*$ 0.5 sqrt35 cosfai)))))) |
(*$ 0.5 sqrt35 cosfai)))))) |
|
|
|
; |
|
; Gauss-Seidel 法により三重対角行列の解を求めているが, |
|
; 優対角行列でない場合は問題があり |
|
; LU分解の方が良い? |
|
; |
(defun gs (n array x b) |
(defun gs (n array x b) |
(do ((i 0 (1+ i))) |
(do ((i 0 (1+ i))) |
((> i 10)) |
((> i 10)) |