// Copyright (C) 2015-2017 Mikael Berthe <mikael@lilotux.net>. All rights reserved.
// Use of this source code is governed by the MIT license,
// which can be found in the LICENSE file.
//
// Thanks to atbrask's blog post <http://www.atbrask.dk/?p=98> for the
// protocol details.

// gobm65 is a Beurer BM65 Blood Pressure Monitor CLI reader.

package main

// Installation:
//
// % go get hg.lilotux.net/golang/mikael/gobm65
//
// Examples:
//
// Get help:
// % gobm65 --help
//
// Get records and display the average:
// % gobm65 --average
// ... display more statistics:
// % gobm65 --stats
//
// Display the latest 3 records with the average:
// % gobm65 -l 3 --average
// Display all records since a specific date:
// % gobm65 --since "2016-06-01"
// Display all records of the last 7 days:
// % gobm65 --since "$(date "+%F" -d "7 days ago")"
//
// Display statistics for morning records:
// % gobm65 --from-time 06:00 --to-time 12:00 --stats
// One can invert times to get night data:
// % gobm65 --from-time 21:00 --to-time 09:00
//
// Display the last/first 10 records in JSON:
// % gobm65 -l 10 --format json
//
// Save the records to a JSON file:
// % gobm65 -o data_u2.json
//
// Read a JSON file and display average of the last 3 records:
// % gobm65 -i data_u2.json -l 3 --average
// % gobm65 -i data_u2.json -l 3 --stats
// Read a JSON file, merge with device records, and save to another file:
// % gobm65 -i data_u2.json --merge -o data_u2-new.json

import (
	"encoding/json"
	"fmt"
	"io"
	"io/ioutil"
	"log"
	"math"
	"sort"
	"time"

	flag "github.com/docker/docker/pkg/mflag"
	"github.com/tarm/serial"
)

type measurement struct {
	Header    int
	Systolic  int
	Diastolic int
	Pulse     int
	Month     int
	Day       int
	Hour      int
	Minute    int
	Year      int
}

type simpleTime struct {
	hour, minute int
}

func getData(s io.ReadWriteCloser, buf []byte, size int) (int, error) {
	t := 0
	b := buf
	for t < size {
		n, err := s.Read(b[t:])
		if err != nil {
			log.Fatal(err) // XXX
			return t, err
		}
		//log.Printf("(%d bytes) %q\n", n, b[t:t+1])
		t = t + n
	}
	return t, nil
}

func fetchData(dev string) (items []measurement, err error) {
	c := &serial.Config{Name: dev, Baud: 4800}

	var s *serial.Port
	s, err = serial.OpenPort(c)
	if err != nil {
		return items, err
	}

	// =================== Handshake =====================
	q := []byte("\xaa")
	//log.Printf("Query: %q\n", q)
	log.Println("Starting handshake...")
	n, err := s.Write(q)
	if err != nil {
		return items, err
	}

	buf := make([]byte, 128)
	n, err = getData(s, buf, 1)
	if err != nil {
		return items, err
	}
	if n == 1 && buf[0] == '\x55' {
		log.Println("Handshake successful.")
	} else {
		log.Printf("(%d bytes) %q\n", n, buf[:n])
		s.Close()
		return items, fmt.Errorf("handshake failed")
	}

	// =================== Desc =====================
	q = []byte("\xa4")
	//log.Printf("Query: %q\n", q)
	log.Println("Requesting device description...")
	n, err = s.Write(q)
	if err != nil {
		return items, err
	}

	n, err = getData(s, buf, 32)
	log.Printf("DESC> %q\n", buf[:n])

	// =================== Count =====================
	q = []byte("\xa2")
	//log.Printf("Query: %q\n", q)
	log.Println("Requesting data counter...")
	n, err = s.Write(q)
	if err != nil {
		return items, err
	}

	n, err = getData(s, buf, 1)
	if err != nil {
		return items, err
	}
	var nRecords int
	if n == 1 {
		log.Printf("%d item(s) available.", buf[0])
		nRecords = int(buf[0])
	} else {
		log.Printf("(%d bytes) %q\n", n, buf[:n])
		return items, fmt.Errorf("no measurement found")
	}

	// =================== Records =====================
	for i := 0; i < nRecords; i++ {
		q = []byte{'\xa3', uint8(i + 1)}
		//log.Printf("Query: %q\n", q)
		//log.Printf("Requesting measurement %d...", i+1)
		n, err = s.Write(q)
		if err != nil {
			return items, err
		}

		n, err = getData(s, buf, 9)
		//log.Printf("DESC> %q\n", buf[:n])

		var data measurement
		data.Header = int(buf[0])
		data.Systolic = int(buf[1]) + 25
		data.Diastolic = int(buf[2]) + 25
		data.Pulse = int(buf[3])
		data.Month = int(buf[4])
		data.Day = int(buf[5])
		data.Hour = int(buf[6])
		data.Minute = int(buf[7])
		data.Year = int(buf[8]) + 2000
		items = append(items, data)
	}

	s.Close()
	return items, nil
}

func loadFromJSONFile(filename string) (items []measurement, err error) {
	data, err := ioutil.ReadFile(filename)
	if err != nil {
		return items, err
	}

	err = json.Unmarshal(data, &items)
	return items, err
}

func mergeItems(newItems, oldItems []measurement) []measurement {
	var result []measurement

	// Sort method: isLater returns true if mi's date is later or
	// equal to mj's date.
	isLater := func(mi, mj measurement) bool {
		switch {
		case mi.Year < mj.Year:
			return false
		case mi.Year > mj.Year:
			return true
		case mi.Month < mj.Month:
			return false
		case mi.Month > mj.Month:
			return true
		case mi.Day < mj.Day:
			return false
		case mi.Day > mj.Day:
			return true
		case mi.Hour < mj.Hour:
			return false
		case mi.Hour > mj.Hour:
			return true
		case mi.Minute < mj.Minute:
			return false
		default:
			return true
		}
	}

	sort.Slice(oldItems, func(i, j int) bool {
		return isLater(oldItems[i], oldItems[j])
	})
	sort.Slice(newItems, func(i, j int) bool {
		return isLater(newItems[i], newItems[j])
	})

	appendIfMissing := func(l []measurement, m measurement) []measurement {
		var i int
		for i = range l {
			if !isLater(l[i], m) {
				break
			}
			if l[i] == m { // Duplicate
				return l
			}
		}
		if i == len(l) {
			return append(l[:i], m)
		}
		r := append(l[:i], append(l[i:], m)...)
		return r
	}

	for _, item := range newItems {
		result = appendIfMissing(result, item)
	}
	for _, item := range oldItems {
		result = appendIfMissing(result, item)
	}
	return result
}

func parseDate(dateStr string) (date time.Time, err error) {
	if dateStr == "" {
		return
	}

	var yy, mm, dd, h, m, s int
	n, e := fmt.Sscanf(dateStr, "%d-%d-%d %d:%d:%d", &yy, &mm, &dd, &h, &m, &s)
	if e != nil && n < 3 {
		err = e
		return
	}
	if n < 6 {
		log.Printf("Date parsed with only %d fields\n", n)
	}
	date = time.Date(yy, time.Month(mm), dd, h, m, s, 0, time.Local)
	return
}

func parseTime(timeStr string) (t simpleTime, err error) {
	_, err = fmt.Sscanf(timeStr, "%d:%d", &t.hour, &t.minute)
	return
}

func average(items []measurement) (measurement, error) {
	var avgMeasure measurement
	var avgCount int

	for _, data := range items {
		avgMeasure.Systolic += data.Systolic
		avgMeasure.Diastolic += data.Diastolic
		avgMeasure.Pulse += data.Pulse
		avgCount++
	}

	roundDivision := func(a, b int) int {
		return int(0.5 + float64(a)/float64(b))
	}

	if avgCount == 0 {
		return avgMeasure, fmt.Errorf("cannot compute average: empty set")
	}

	avgMeasure.Systolic = roundDivision(avgMeasure.Systolic, avgCount)
	avgMeasure.Diastolic = roundDivision(avgMeasure.Diastolic, avgCount)
	avgMeasure.Pulse = roundDivision(avgMeasure.Pulse, avgCount)

	return avgMeasure, nil
}

func intMedian(numbers []int) int {
	middle := len(numbers) / 2
	med := numbers[middle]
	if len(numbers)%2 == 0 {
		med = (med + numbers[middle-1]) / 2
	}
	return med
}

func median(items []measurement) (measurement, error) {
	var med measurement
	if len(items) == 0 {
		return med, fmt.Errorf("cannot compute average: empty set")
	}

	var sys, dia, pul []int
	for _, data := range items {
		sys = append(sys, data.Systolic)
		dia = append(dia, data.Diastolic)
		pul = append(pul, data.Pulse)
	}

	sort.Ints(sys)
	sort.Ints(dia)
	sort.Ints(pul)

	med.Systolic = intMedian(sys)
	med.Diastolic = intMedian(dia)
	med.Pulse = intMedian(pul)

	return med, nil
}

func stdDeviation(items []measurement) (measurement, error) {
	var sDev measurement

	if len(items) <= 1 {
		return sDev, fmt.Errorf("cannot compute deviation: set too small")
	}

	var sumSys, sumDia, sumPul float64
	avg, err := average(items)
	if err != nil {
		return sDev, err
	}

	for _, data := range items {
		sumSys += math.Pow(float64(data.Systolic-avg.Systolic), 2)
		sumDia += math.Pow(float64(data.Diastolic-avg.Diastolic), 2)
		sumPul += math.Pow(float64(data.Pulse-avg.Pulse), 2)
	}

	sDev.Systolic = int(math.Sqrt(sumSys / float64(len(items)-1)))
	sDev.Diastolic = int(math.Sqrt(sumDia / float64(len(items)-1)))
	sDev.Pulse = int(math.Sqrt(sumPul / float64(len(items)-1)))

	return sDev, nil
}

func main() {
	inFile := flag.String([]string{"-input-file", "i"}, "", "Input JSON file")
	outFile := flag.String([]string{"-output-file", "o"}, "", "Output JSON file")
	limit := flag.Uint([]string{"-limit", "l"}, 0, "Limit number of items to N first")
	since := flag.String([]string{"-since"}, "",
		"Filter records from date (YYYY-mm-dd HH:MM:SS)")
	format := flag.String([]string{"-format", "f"}, "", "Output format (csv, json)")
	avg := flag.Bool([]string{"-average", "a"}, false, "Compute average")
	stats := flag.Bool([]string{"-stats"}, false, "Compute statistics")
	merge := flag.Bool([]string{"-merge", "m"}, false,
		"Try to merge input JSON file with fetched data")
	device := flag.String([]string{"-device", "d"}, "/dev/ttyUSB0", "Serial device")
	fromTime := flag.String([]string{"-from-time"}, "", "Select records after time (HH:MM)")
	toTime := flag.String([]string{"-to-time"}, "", "Select records bofore time (HH:MM)")

	var startTime, endTime simpleTime

	flag.Parse()

	switch *format {
	case "":
		if *outFile == "" {
			*format = "csv"
		}
		break
	case "json", "csv":
		break
	default:
		log.Fatal("Unknown output format.  Possible choices are csv, json.")
	}

	if *fromTime != "" {
		if t, err := parseTime(*fromTime); err != nil {
			log.Fatal("Cannot parse 'from' time: ", err)
		} else {
			startTime = t
		}
	}
	if *toTime != "" {
		if t, err := parseTime(*toTime); err != nil {
			log.Fatal("Cannot parse 'to' time: ", err)
		} else {
			endTime = t
		}
	}

	startDate, err := parseDate(*since)
	if err != nil {
		log.Fatal("Could not parse date: ", err)
	}

	var items []measurement

	// Read data

	if *inFile == "" {
		// Read from device
		if items, err = fetchData(*device); err != nil {
			log.Fatal(err)
		}
	} else {
		// Read from file
		var fileItems []measurement
		if fileItems, err = loadFromJSONFile(*inFile); err != nil {
			log.Fatal(err)
		}
		if *merge {
			if items, err = fetchData(*device); err != nil {
				log.Fatal(err)
			}
			items = mergeItems(items, fileItems)
		} else {
			items = fileItems
		}
	}

	// Apply filters

	if !startDate.IsZero() {
		log.Printf("Filtering out records before %v...\n", startDate)
		for i := range items {
			iDate := time.Date(items[i].Year, time.Month(items[i].Month),
				items[i].Day, items[i].Hour, items[i].Minute, 0, 0,
				time.Local)
			if iDate.Sub(startDate) < 0 {
				items = items[0:i]
				break
			}
		}
	}

	if *fromTime != "" || *toTime != "" {
		log.Println("Filtering hours...")

		compare := func(m measurement, t simpleTime) int {
			if m.Hour*60+m.Minute < t.hour*60+t.minute {
				return -1
			}
			if m.Hour*60+m.Minute > t.hour*60+t.minute {
				return 1
			}
			return 0
		}

		inv := false
		if *fromTime != "" && *toTime != "" &&
			startTime.hour*60+startTime.minute > endTime.hour*60+endTime.minute {
			inv = true
		}

		var newItems []measurement
		for _, data := range items {
			if inv {
				if compare(data, startTime) == -1 && compare(data, endTime) == 1 {
					continue
				}
				newItems = append(newItems, data)
				continue
			}
			if *fromTime != "" && compare(data, startTime) == -1 {
				continue
			}
			if *toTime != "" && compare(data, endTime) == 1 {
				continue
			}
			newItems = append(newItems, data)
		}
		items = newItems
	}

	if *limit > 0 && len(items) > int(*limit) {
		items = items[0:*limit]
	}

	// Done with filtering

	if *format == "csv" {
		for i, data := range items {
			fmt.Printf("%d;%x;%d-%02d-%02d %02d:%02d;%d;%d;%d\n",
				i+1, data.Header,
				data.Year, data.Month, data.Day,
				data.Hour, data.Minute,
				data.Systolic, data.Diastolic, data.Pulse)
		}
	}

	if *stats {
		*avg = true
	}

	if *avg && len(items) > 0 {
		avgMeasure, err := average(items)
		if err != nil {
			log.Println("Error:", err)
		} else {
			fmt.Printf("Average: %d;%d;%d\n", avgMeasure.Systolic,
				avgMeasure.Diastolic, avgMeasure.Pulse)
		}
	}

	if *stats && len(items) > 1 {
		d, err := stdDeviation(items)
		if err != nil {
			log.Println("Error:", err)
		} else {
			fmt.Printf("Standard deviation: %d;%d;%d\n",
				d.Systolic, d.Diastolic, d.Pulse)
		}
	}
	if *stats && len(items) > 0 {
		m, err := median(items)
		if err != nil {
			log.Println("Error:", err)
		} else {
			fmt.Printf("Median values: %d;%d;%d\n",
				m.Systolic, m.Diastolic, m.Pulse)
		}
	}

	if *format == "json" || *outFile != "" {
		rawJSON, err := json.MarshalIndent(items, "", "  ")
		if err != nil {
			log.Fatal("Error:", err)
		}

		if *format == "json" {
			fmt.Println(string(rawJSON))
		}
		if *outFile != "" {
			err = ioutil.WriteFile(*outFile, rawJSON, 0600)
			if err != nil {
				log.Println("Could not write output file:", err)
			}
		}
	}
}