ブランコに乗った水槽

#include "udf.h"

#define RADIUS  4.0  /* [m] */
#define PERIOD  8.0  /* [sec] */
#define OMG     (2.0*M_PI/PERIOD)   /* [rad/sec] */

#define ANGLEMAX  30.0   /* [deg] */
#define THETAMAX  (ANGLEMAX*M_PI/180.0)  /* [rad] */

#define GRAV   9.81  /* [m/sec^2] */

/************************************************************/
DEFINE_CG_MOTION(pendulum,dt,vel,omega,time,dtime)
/************************************************************/
{
   real theta,omgt,thetadot;

   omgt = OMG*time;
   theta = THETAMAX * sin(omgt);

   thetadot = THETAMAX * OMG * cos(omgt);

   vel[0] = RADIUS * cos(theta) * thetadot;
   vel[1] = RADIUS * sin(theta) * thetadot;
   vel[2] = 0;

   omega[0] = 0;
   omega[1] = 0;
   omega[2] = thetadot;
}

/************************************************************/
DEFINE_SOURCE(gravy, c, t, dS, eqn)
/************************************************************/
{
   real theta,source;

   source = - C_R(c,t) * GRAV;

   dS[eqn] = 0;
   return source;
}

(command1 argument1 argument2 ...)

> (define x 10) ; Defines a global variable
> (set! x 20) ;Set! for redefinitions of global variables
> (display x) ; Prints the value of the variable on the prompt
20
> (let ((x 2) (y 3)) (* x y))
6

> (define a1 #t) ; Define the variable as: true
> (define a2 #f) ; Define the variable as: false
> (define c1 #\A) ; Define the variable as the uppercase letter : A
> (define d1 "teststring") ; Define the variable as the string : teststring
> (define e1 ’testsymbol ) ; Define the variable as the symbol : testsymbol
> (symbol-bound? ’e1 (the-environment)) ;Check , if variable e1 is bound
#t
> (symbol-assigned? ’e1 (the-environment)) ;Check , if variable e1 is assigned (contains a value)
#t

> (rpgetvar ’time-step) ;時間ステップ番号を得る
0
> (rpsetvar ’time-step 0) ;時間ステップ番号を強制的に0に設定する

> (cxgetvar 'contour/filled?) ; 塗りつぶしコンター図の指定
#t
> (cxsetvar 'contour/filled? #f)
#f

> (define x 10) ; Define variable x as 10
> (define y 20) ; Define variable y as 20
> (define z1 (+ x y)) ; Calculates x+y= 30
> (- x y) ; Calculate x-y= -10
> (* x y) ; Calculate x*y= 200
> (/ x y) ; Calculate x/y= 0.5
> (/ x y 4) ; Calculate x/y/4= 0.125
> (exp 2) ; Calculate e?2=7.389
> (expt 2 5) ; Calculate 2?5=32
> (log 7.389) ; Calculate log(7.389) =1.9999
> (sin 0) ; Calculate sin(0 rad) = 0, same for cos and tan
> (asin 1) ; Calculate arcsin(1) = pi/2, same for acos and atan

> (define l1 ’(a b c)) ; Define l1 as list with elements a, b and c
> (car l1) ; Return the first element of a list
> (cdr l1) ; Return all elements of a list , except the first
> (length l1) ; Return the length (3) of a list
> (list -ref l1 2) ; Return the second element of a list

> (= a b) ;Check , if a and b are equal , use for numbers
> (< a b) ;Check , if a is less than b, use for numbers
> (<= a b) ;Check , if a is less than or equal to b, use for numbers
> (positive? a) ;Check , if a is a positive number , (also : zero? and negative?)
> (eqv? a b) ;Check , if values of a and b are equal , use for non-numeric values
> (eqal? a b) ;Check , if values of a and b are equal , use for lists , vectors etc.
> (not a) ; Invert logical value of a
> (and ((= 3 3) (eqv? a b))) ;#t, if 3 is equal to 3 AND values of a and b
> (or ((= 3 3) (eqv? a b))) ;#t, if 3 is equal to 3 OR values of a and b
> (if (< 2 3) ’yes ’no) ;If the test expression is #t, yes is returned otherwise no
> (if (< 2 3) ; Begin enables more than one expression per body
> (begin (+ 10 20) (+ 30 40)) ; True body
> (begin (+ 1 2) (+ 3 4))) ; False body

(do ((x1 x1start (+ x1 x1delta))) ((> x1 x1end)) ... loop-body ... )

> (do ((x1 1 (+ x1 1))) ; カウンターx1の初期値、増分値を定義
((> x1 5)) ; 終了条件を定義
(display (* x1 10))) ; 実行文を記述
1020304050

> (define sVariable1 "bananas")
> (define iVariable1 4)
> (define iVariable2 13)
> (define fpVariable1 0.001)
> (format "Monkeys like ?a." sVariable1) ; Insert sVariable1 in the string Monkeys like bananas.
> (format "Monkeys like ?a." "bananas") ; Insert "bananas" in the string Monkeys like bananas.
> (format "Monkeys like ?s." sVariable1) ; String specifier inserts the string with " into the string.
Monkeys like "bananas".
> (format "Elephants can reach a height of ?d meters/?d feet." iVariable1 iVariable2) ; Usage of more than one format specifier
Elephants can reach a height of 4 meters/13 feet.
> (format " Elephants can reach a height of ?02d meters/?04d feet." iVariable1 iVariable2 ) ; Fill zeros on the left side
Elephants can reach a height of 04 meters/0013 feet.
> (format "A millisecond is specified as ?f seconds." fpVariable1 ) ;Use a floating point variable with no extra specifications
A millisecond is specified as 0 seconds.
> (format "A millisecond is specified as ?4.3f seconds." fpVariable1 ) ;Specify explicitly 4 digits , 3 after the decimal
A millisecond is specified as 0.001 seconds.

(ti-menu-load-string (format #f "/views/camera/position  ~a ~a 1" xc yc) )

 ;
 ;  ファイル localview.scm
 ;
(define maxstep 4000)
(define ident "vof")
(define head (string-append ident "-1-"))
(define tail ".h5")

; 計算結果読み込みのための無名関数を定義
(define rcd-step (lambda (step) (read-case-data (string-append head step tail))))


(define radius 4)
(define period 8)
(define anglemax 30)
(define pi 3.14159265358979)
(define twopi (* 2 pi))
(define thetamax (/ (* anglemax pi) 180))
(define omega    (/ twopi period))

; 初期基準点(position)座標 (cameraで使用)
(define xini 0.0)
(define yini 0.57)

; テキストボックス表示、文字サイズを拡大
(ti-menu-load-string (format #f "/disp/set/window/text/visible y") )
(ti-menu-load-string (format #f "/disp/set/window/text/date n") )
(ti-menu-load-string (format #f "/disp/set/window/text/application n") )
(cxsetvar 'frame/title/font (format #f "size=0.3wru"))



;;; -------------------------------------- 時間ステップループ
; (0ステップから4000ステップまで20ステップ刻み)
(do ((step 0 (+ step 20))) ((> step maxstep))
(define mystep (format #f "~05d" step))

; 計算結果ファイルの読み込み
(rcd-step mystep)


; 時刻データの取得
(define time (rpgetvar 'flow-time))

; theta, sin(theta), cos(theta)の計算
(define theta (* thetamax (sin (* omega time))))
(define sn (sin theta))
(define cs (cos theta))

; 角度theta回転後の(xini,yini)座標位置の計算
; xc = (R-yi)*sin + xi*cos
; yc = R - ((R-yi)*cos - xi*sin)
(define xc (+
             (* (- radius yini) sn)
             (* xini cs)
           ))
(define yc (- radius
               (-
                 (* (- radius yini) cs)
                 (* xini sn)
               )
            ))

; camera表示範囲、上向き方向の指定
(define fldx 1.5)
(define fldy 1.5)
(define upvx (- sn))
(define upvy cs)

; TUIコマンドを schemeコマンドに変換して実行
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; vof
(ti-menu-load-string (format #f "/disp/obj/disp vof") )

(ti-menu-load-string (format #f "/views/camera/position  ~a ~a 1" xc yc) )
(ti-menu-load-string (format #f "/views/camera/target    ~a ~a 0" xc yc) )
(ti-menu-load-string (format #f "/views/camera/up-vector ~a ~a 0" upvx upvy) )
(ti-menu-load-string (format #f "/views/camera/field     ~a ~a  " fldx fldy) )
(ti-menu-load-string (format #f "/disp/save-pic/lvof-~05d.png" step) )

(ti-menu-load-string (format #f "/views/restore-view view-0" ) )
(ti-menu-load-string (format #f "/disp/save-pic/gvof-~05d.png" step) )

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; vect
(ti-menu-load-string (format #f "/disp/obj/edit vof-local-vect scale auto-scale n scale-f 0.5 q q q") )
(ti-menu-load-string (format #f "/disp/obj/disp sc-local-vect") )

(ti-menu-load-string (format #f "/views/camera/position  ~a ~a 1" xc yc) )
(ti-menu-load-string (format #f "/views/camera/target    ~a ~a 0" xc yc) )
(ti-menu-load-string (format #f "/views/camera/up-vector ~a ~a 0" upvx upvy) )
(ti-menu-load-string (format #f "/views/camera/field     ~a ~a  " fldx fldy) )
(ti-menu-load-string (format #f "/disp/save-pic/lvect-~05d.png" step) )

(ti-menu-load-string (format #f "/disp/obj/edit vof-vect scale auto-scale n scale-f 0.2 q q q") )
(ti-menu-load-string (format #f "/disp/obj/disp sc-vect") )
(ti-menu-load-string (format #f "/views/restore-view view-0" ) )
(ti-menu-load-string (format #f "/disp/save-pic/gvect-~05d.png" step) )

)

; コンソールウインドウから実行する場合
(load "localview.scm")