You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
101 lines
2.1 KiB
Go
101 lines
2.1 KiB
Go
package main
|
|
|
|
import (
|
|
"database/sql"
|
|
_ "github.com/mattn/go-sqlite3"
|
|
"math"
|
|
"os"
|
|
)
|
|
|
|
func setupDb() {
|
|
db, err := sql.Open("sqlite3", "./exo.db")
|
|
if err != nil {
|
|
bail(E_No_DB, "unable to open database: %v", err)
|
|
}
|
|
defer db.Close()
|
|
|
|
stmnt := `create table if not exists planets (
|
|
id integer not null primary key autoincrement,
|
|
name text,
|
|
x integer,
|
|
y integer,
|
|
z integer,
|
|
planets integer
|
|
);`
|
|
_, err = db.Exec(stmnt)
|
|
if err != nil {
|
|
log_error("couldn't create table: %v", err)
|
|
return
|
|
}
|
|
|
|
n, err := countPlanets(db)
|
|
log_info("num planets: %v error: %v", n, err)
|
|
if n == 0 {
|
|
fi, err := os.Open(dataPath)
|
|
if err != nil {
|
|
bail(E_No_Data, "unable to open data path: %v", err)
|
|
}
|
|
c := make(chan exoSystem)
|
|
go speckStream(fi, c)
|
|
for planet := range c {
|
|
planet.Store(db)
|
|
}
|
|
}
|
|
|
|
idx := indexPlanets(db)
|
|
log_info("%v", idx)
|
|
fillEdges(db, idx)
|
|
|
|
stmnt = `create table if not exists edges (
|
|
id_1 integer,
|
|
id_2 integer,
|
|
distance real
|
|
);`
|
|
_, err = db.Exec(stmnt)
|
|
if err != nil {
|
|
log_error("couldn't create distance table: %v", err)
|
|
return
|
|
}
|
|
|
|
}
|
|
|
|
func sq(x float64) float64 {
|
|
return x * x
|
|
}
|
|
|
|
func dist3d(x1, y1, z1, x2, y2, z2 float64) float64 {
|
|
return math.Sqrt(sq(x1-x2) + sq(y1-y2) + sq(z1-z2))
|
|
}
|
|
|
|
func planetDistance(p1, p2 exoSystem) float64 {
|
|
return dist3d(p1.x, p1.y, p1.z, p2.x, p2.y, p2.z)
|
|
}
|
|
|
|
func indexPlanets(db *sql.DB) map[int]exoSystem {
|
|
rows, err := db.Query(`select * from planets`)
|
|
if err != nil {
|
|
log_error("unable to select all planets: %v", err)
|
|
return nil
|
|
}
|
|
defer rows.Close()
|
|
planets := make(map[int]exoSystem, 551)
|
|
for rows.Next() {
|
|
var id int
|
|
p := exoSystem{}
|
|
if err := rows.Scan(&id, &p.name, &p.x, &p.y, &p.z, &p.planets); err != nil {
|
|
log_info("unable to scan planet row: %v", err)
|
|
continue
|
|
}
|
|
planets[id] = p
|
|
}
|
|
return planets
|
|
}
|
|
|
|
func fillEdges(db *sql.DB, planets map[int]exoSystem) {
|
|
for i := 0; i < len(planets); i++ {
|
|
for j := i + 1; j < len(planets); j++ {
|
|
log_info("distance from %s to %s: %v", planets[i].name, planets[j].name, planetDistance(planets[i], planets[j]))
|
|
}
|
|
}
|
|
}
|