func AddScaledTo(dst, y []float64, alpha float64, s []float64) []float64 {
	if len(dst) != len(s) || len(dst) != len(y) {
		panic("floats: lengths of slices do not match")
	}
	f64.AxpyUnitaryTo(dst, alpha, s, y)
	return dst
}
AddScaledTo performs dst = y + alpha * s, where alpha is a scalar,
and dst, y and s are all slices.
It panics if the lengths of dst, y, and s are not equal.

At the return of the function, dst[i] = y[i] + alpha * s[i]
package floats
import (
	"errors"
	"math"
	"sort"
	"strconv"
	"gonum.org/v1/gonum/internal/asm/f64"
)
func Add(dst, s []float64) {
	if len(dst) != len(s) {
		panic("floats: length of the slices do not match")
	}
	f64.AxpyUnitaryTo(dst, 1, s, dst)
}
func AddTo(dst, s, t []float64) []float64 {
	if len(s) != len(t) {
		panic("floats: length of adders do not match")
	}
	if len(dst) != len(s) {
		panic("floats: length of destination does not match length of adder")
	}
	f64.AxpyUnitaryTo(dst, 1, s, t)
	return dst
}
func AddConst(c float64, dst []float64) {
	for i := range dst {
		dst[i] += c
	}
}
func AddScaled(dst []float64, alpha float64, s []float64) {
	if len(dst) != len(s) {
		panic("floats: length of destination and source to not match")
	}
	f64.AxpyUnitaryTo(dst, alpha, s, dst)
}