2 I I / 61

Size: px

Start display at page:

Download "2 I I / 61"

しょうすけおおはし
5 years ago
Views:

1 2 I I / 61

2 I / 61

3 I / 61

4 7/13 2 7/20 I 7/27 II I / 61

5 π-computer gnuplot MobaXterm Wiki PC X11 DISPLAY I / 61

6 Mac 1. X11 XQuartz.app 2. ssh -Y my 3. π-computer gnuplot I / 61

7 gnuplot> gnuplot> plot sin(x) I / 61

8 cp -r /tmp/170713/ your_directory ls -l 4 data, sierpinski, src01d, src02d cd your_directory/src01d I / 61

9 1D data1d generator.f95 module I / 61

10 const m module const m i m p l i c i t none! < f o r t r a n c o n s t a n t s > i n t e g e r, parameter : : SP = s e l e c t e d r e a l k i n d ( 6 ) i n t e g e r, parameter : : DP = s e l e c t e d r e a l k i n d ( 1 5 )! <s i m u l a t i o n c o n s t a n t s > i n t e g e r, parameter : : NN = 101 end module const m I / 61

11 f95 data1d_generator.f95 &&./a.out > data.1d make data.1d make I / 61

12 data.1d I / 61

13 gnuplot # 1 x, y 1 plot x 2 y I / 61

14 gnuplot gnuplot> plot data.1d w lp PDF w lp with linespoints linepoint w linespoints I / 61

15 I / 61

16 gnuplot gnuplot gnuplot script src01d graph1d.gp # # A simple sample gnuplot script: graph1d.gp # # # set size square # same side lengths for x and y # set size 0.65, 1 # other way to set regular rectan set xlabel x # x-axis set ylabel func # y-axis set xrange[-10:10] # x coords range set yrange[-1:1] # y coords range set title graph 1d plot data.1d w lp pause -1 I / 61

17 gnuplot quit shell gnuplot graph1d.gp make graph I / 61

18 Makefile make target source target: source [ ] command_to_make_target_from_source target source source target I / 61

19 Makefile make target: source01 [ ] command_to_make_target_from_source01 source01: source02 [ ] command_to_make_source01_from_source02 source02: source03 [ ]... I / 61

20 I / 61

21 recsin(x, n) = sin(sin(sin( (sin(x))) ) anim data generator.f95 I / 61

22 anim data generator.f95 make anim.data 50 I / 61

23 gnuplot gnuplot pause t t pause t t I / 61

24 set xlabel x set ylabel func set xrange[-10:10] set yrange[-1:1] set title anim 1d plot data.1d.000 w lp pause 5 plot data.1d.001 w lp pause 0.5 plot data.1d.002 w lp pause 0.5 plot data.1d.003 w lp.. pause -1 # x-axis # y-axis # x coords range # y coords range # Enter I / 61

25 100 plot 100 : anim plotscript generator.f95 make anim.gp anim.gp I / 61

26 make anim I / 61

27 I / 61

28 2 1 gnuplot 2 plot splot gnuplot > splot [ f(x,y)] gnuplot > splot x**2 + y**2 I / 61

29 splot f(x,y) I / 61

30 splot gnuplot splot exp(-x**2-y**2) set isosamples 50 # replot set xrange[-2:2] # x replot set yrange[-2:2] # y replot I / 61

31 show view # set view 60,0 # x 60 replot set view 60,30 # x 60 z 30 I / 61

32 I / 61

33 set hidden3d replot # I / 61

34 set contour base # replot set cntrparam levels 30 # replot I / 61

35 unset surface replot # I / 61

36 set pm3d at b replot # draw with colored contour I / 61

37 2D I / 61

38 2D x00 y00 x01 y00 x02 y00 x09 y00 x00 y01 x01 y01 x09 y01 x09 y09 I / 61

39 cd [your directory]/src02d I / 61

40 heat4.f95 I / 61

41 2D x,y s u b r o u t i n e p r i n t p r o f i l e 2 d ( p, j j, f ) t y p e ( r a n k s t ), i n t e n t ( i n ) : : p t y p e ( s p a n t ), i n t e n t ( i n ) : : j j r e a l (DP), d i m e n s i o n ( 0 : NGRID+1, & r e a l (DP), j j%s t t 1: j j%end+1), i n t e n t ( i n ) : : f d i m e n s i o n ( 0 : NGRID+1,0:NGRID+1) & : : f g l o b a l! 2d p r o f to be saved i n t e g e r : : c o u n t e r = 0! has s a v e a t t r i b. t y p e ( s p a n t ) : : j j 2! used f o r f g l o b a l c h a r a c t e r ( l e n =4) : : s e r i a l n u m! put on f i l e name c h a r a c t e r ( l e n = ), parameter : : base =.. / data /temp. 2 d. i n t e g e r : : i, j I / 61

42 2D j j 2 = a d j u s t j s t a r t a n d j e n d ( p, j j ) w r i t e ( s e r i a l n u m, ( i 4. 4 ) ) c o u n t e r f g l o b a l ( :, : ) = s e t p r o f 2 d ( j j, j j 2, f ) i f ( p%myrank==0 ) then open ( 1 0, f i l e =base // s e r i a l n u m ) do j = 0, NGRID+1 do i = 0, NGRID+1 w r i t e ( 1 0, ) i, j, f g l o b a l ( i, j ) end do w r i t e ( 1 0, )! g n u p l o t r e q u i r e s a b l a n k l i n e h e r e. end do c l o s e ( 1 0 ) end i f c o u n t e r = c o u n t e r + 1 end s u b r o u t i n e p r i n t p r o f i l e 2 d I / 61

43 1. mpifrtpx heat4.f95 2. pjsub heat4.sh../data temp.2d.0000 make data OK I / 61

44 E E E E-002 I / 61../data temp.2d E E E E E

45 heat4 lines.gp 2 gnuplot # # A sample gnuplot script: heat4_lines.gp # # [ line contours ] # # set size square # same side lengths for x and y set size 0.65, 1 # same side lengths for x and y set xlabel "i" # x-axis set ylabel "j" # y-axis set xrange[0:50] # i-grid min & max set yrange[0:50] # j-grid min & max set nosurface # do not show surface plot unset ztics # do not show z-tics set contour base # enables contour lines set cntrparam levels 10 # draw 10 contours set view 0,0 # view from the due north set title "Temperature" I / 61

46 gnuplot UNIX gnuplot heat4 lines.gp I / 61

47 I / 61

48 gnuplot heat4 colors.gp I / 61

49 heat4 colors.gp # # A sample gnuplot script: heat4_plot_contour_colors.gp # # [ color contours ] # # set size square # same side lengths for x and y set size 0.65, 1 # same side lengths for x and y set xlabel "i" # x-axis set ylabel "j" # y-axis set xrange[0:50] # i-grid min & max set yrange[0:50] # j-grid min & max set palette defined (0 "blue", 0.15 "red", 0.3 "yellow") set nosurface # do not show surface plot unset ztics # do not show z-tics set pm3d at b # draw with colored contour set view 0,0 # view from the due north set title "Temperature " splot "../data/temp.2d.0000" using 1:2:3 pause -1 I / 61

50 I / 61

51 bird s eye view I / 61

52 gnuplot plot4_plot_birdseyeview.gp # # a sample gnuplot script: plot4_plot_birdseyeview.gp # # [ Bird"s Eye View ] # # set size square # same side lengths for x and y set size 0.65, 1 set xlabel "i" # x-axis set ylabel "j" # y-axis set xrange[0:50] # i-grid min & max set yrange[0:50] # j-grid min & max set contour base # enables contour lines set cntrparam levels 10 # draw 10 contours # set palette defined (0 "blue", 0.15 "red", 0.3 "yellow") # set pm3d # draw with colored contour set title "Temperature " splot "../data/temp.2d.0000" using 1:2:3 w l pause -1 I / 61

53 I / 61

54 splot view heat4 plot rotating birdseyeview generator.f95 gfortran heat4 plot rotating birdseyeview generator.f95./a.out > anim.gp anim.gp gnuplot anim.gp 5 make anim OK make clean make data../data make anim I / 61

55 I / 61

56 cd [your directory]/sierpinski OpenMP 2 sierp.f95 I / 61

57 sierp.f95 sierp2.f95 omp section 3 module private public I / 61

58 sierp2.f95 enshu.sh frtpx -Kopenmp sierp2.f95 pjsub enshu.sh make data Makefile : data data.099 data.??? make graph I / 61

59 I / 61

60 sierp2.f95 1/3 (x n+1, y n+1 ) = (x n /2 + 1, y n /2) 1/3 (x n+1, y n+1 ) = (x n /2 1, y n /2) 1/3 (x n+1, y n+1 ) = (x n /2, y n + 1) 1/10 (x n+1, y n+1 ) = (x n /2 + 1, y n /2) 1/10 (x n+1, y n+1 ) = (x n /2 1, y n /2) 8/10 (x n+1, y n+1 ) = (x n /2, y n + 1) sierp3.f95 gnuplot I / 61

61 π-computer diff sierp2.f95 sierp3.f95 mail kage 7/19 23:59 I / 61

2 A I / 58

2 A I / 58 2 A 2018.07.12 I 2 2018.07.12 1 / 58 I 2 2018.07.12 2 / 58 π-computer gnuplot 5/31 1 π-computer -X ssh π-computer gnuplot I 2 2018.07.12 3 / 58 gnuplot> gnuplot> plot sin(x) I 2 2018.07.12 4 / 58 cp -r