astroimsim-geometry 0.1.19

Geometry package for astroimsim
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57

use crate::grid2d::PlotPoint;
use plotpy::Plot;
use astroimsim_geometry::coordinate_system::Coordinates::RELATIVE;
use astroimsim_geometry::coordinate_system::{CoordinateSystem, Coordinates};
use astroimsim_geometry::grid2d::GRID2D;
use astroimsim_geometry::grid2d::PlotPoint::Given;
use astroimsim_geometry::points::Point;
use crate::coordinate_system::Coordinates::ABSOLUTE;
use crate::projection_test::test_project;

pub mod coordinate_system;
pub mod grid1d;
pub mod grid2d;
pub mod points;
mod projection_test;

fn main() {

    let mut plot = Plot::new();
    let coords = CoordinateSystem::new((7.0,-6.2),(-0.4,4.0),(0.0,0.0), "test_1");
    let grid = GRID2D::new_empty((9,9),(1.0,1.0),(0.0,0.0), 0.0001,RELATIVE(coords.clone()));
   // println!("rid number 50 {:?}", grid.locate(50).to_absolute().values());
    let point_2 = Point::new(1.5,2.2,Coordinates::RELATIVE(coords.clone()));
   // println!("{:?}", point_2);
    let point = Point::new(2.6, -4.2,Coordinates::ABSOLUTE);
   // println!("point from absolute{:?}",CoordinateSystem::point_from_absolute(&coords.clone(),point));
   // println!("{:?}", point_2.to_absolute());
    //println!("{:?}", point_2.convert(&RELATIVE(coords.clone())));

   // println!("{:?}",grid.fit_grid(&point_2));


    grid.plot_points(&mut plot, Given(point_2));
    coords.plot(&mut plot,"green");
    plot.legend();
    plot.save("test");
    plot.show("slefj").unwrap();



    /*
let grid = GRID2D::new_from_width((9000,9000),(1.0,1.0), (0.0,0.0), 0.01,RELATIVE(coords.clone()));
println!("{:?}",grid.locate(grid.num_points-1));


let grid = GRID2D::new_from_width((6,6),(1.0,1.0), (0.0,0.0), 0.01,RELATIVE(coords.clone()));
println!("{:?}",grid.locate(grid.num_points-1))

     */

}

pub fn norm(vect:(f64,f64))->(f64,f64){
    let norm = (vect.0.powi(2)+ vect.1.powi(2)).sqrt();
    (vect.0/norm,vect.1/norm)

}