program main c*********************************************************************72 c cc ellipsoid_grid_test() tests ellipsoid_grid(). c c Licensing: c c This code is distributed under the MIT license. c c Modified: c c 11 November 2011 c c Author: c c John Burkardt c implicit none call timestamp ( ) write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'ellipsoid_grid_test():' write ( *, '(a)' ) ' FORTRAN77 version' write ( *, '(a)' ) ' Test ellipsoid_grid().' call ellipsoid_grid_test01 ( ) c c Terminate. c write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'ellipsoid_grid_test():' write ( *, '(a)' ) ' Normal end of execution.' write ( *, '(a)' ) ' ' call timestamp ( ) return end subroutine ellipsoid_grid_test01 ( ) c*********************************************************************72 c cc ELLIPSOID_GRID_TEST01 tests ELLIPSOID_GRID. c c Licensing: c c This code is distributed under the MIT license. c c Modified: c c 11 November 2011 c c Author: c c John Burkardt c implicit none double precision c(3) character * ( 80 ) filename integer n integer ng double precision r(3) double precision xyz(3,1200) write ( *, '(a)' ) ' ' write ( *, '(a)' ) 'TEST01:' write ( *, '(a)' ) & ' ELLIPSOID_GRID can define a grid of points' write ( *, '(a)' ) ' with N+1 points on the minor half axis,' write ( *, '(a)' ) ' based on any ellipsoid.' n = 4 r(1) = 2.0D+00 r(2) = 1.0D+00 r(3) = 1.5D+00 c(1) = 1.0D+00 c(2) = 2.0D+00 c(3) = 1.5D+00 write ( *, '(a)' ) ' ' write ( *, '(a)' ) ' We use N = ', n write ( *, '(a,g14.6,a,g14.6,a,g14.6,a)' ) & ' Radius R = (', r(1), ',', r(2), ',', r(3), ')' write ( *, '(a,g14.6,a,g14.6,a,g14.6,a)' ) & ' Center C = (', c(1), ',', c(2), ',', c(3), ')' call ellipsoid_grid_count ( n, r, c, ng ) write ( *, '(a)' ) ' ' write ( *, '(a,i8)' ) ' Number of grid points will be ', ng call ellipsoid_grid ( n, r, c, ng, xyz ) call r83vec_print_part ( ng, xyz, 20, & ' Part of the grid point array:' ) filename = 'ellipsoid_grid_test01.xyz' call r8mat_write ( filename, 3, ng, xyz ) write ( *, '(a)' ) ' ' write ( *, '(a)' ) & ' Data written to the file "' // trim ( filename ) // '".' return end