Removed leftover print statements from previous debugging efforts

chemplot/ramachandran.go

@@ -29,14 +29,13 @@ import (
 	"fmt"
 	"image/color"
 	"math"
-	"strings"
-	"image/color"
+//	"strings"
 	"github.com/gonum/plot"
 	"github.com/gonum/plot/plotter"
 	"github.com/gonum/plot/vg"
 	"github.com/gonum/plot/vg/draw"
-	"github.com/rmera/gochem"
-	"github.com/rmera/gochem/v3"
+//	"github.com/rmera/gochem"
+//	"github.com/rmera/gochem/v3"
 )
 
 const (
@@ -121,7 +120,7 @@ func RamaPlotParts(data [][][]float64, tag [][]int, title, plotname string) erro
 			}
 			//set the colors
 			r, g, b := colors(key, len(data))
-			fmt.Println("DATA POINT", key, "color", r, g, b)
+		//	fmt.Println("DATA POINT", key, "color", r, g, b)
 			s.GlyphStyle.Color = color.RGBA{R: r, B: b, G: g, A: 255}
 			//The tagging procedure is a bit complex.
 			p.Add(s)

clifford.go

@@ -35,7 +35,7 @@ import (
 	"runtime"
 )
 
-import "fmt" //debug
+//import "fmt" //debug
 
 type paravector struct {
 	Real  float64
@@ -278,7 +278,7 @@ func RotateP(Target, Res, axis, Rv, Rvrev, tmp1, tmp2, tmp3, tmp4, itmp1, itmp2,
 			//this "print" causes a data race but it shouldn't matter, as it's only for debugging purposes.
 			//removing the print removes the race warning, so there isn't apparently any data race going on.
 			//fmt.Println("WTF",r,c,i,tmp2,"\n") /////////
-			fmt.Println("") ////////////
+			//fmt.Println("") ////////////
 			t2 := tmp2.VecView(i)
 			t4 := tmp4.VecView(i)
 			pv := paravectorFromVector(t2, t4)

files.go

@@ -612,9 +612,9 @@ func xyzReadSnap(xyz *bufio.Reader, ReadTopol bool) (*v3.Matrix, []*Atom, error)
 	}
 	//This could be done faster if done in the same loop where the coords are read
 	//Instead of having another loop just for them.
-	for k, i := range errs {
+	for _, i := range errs {
 		if i != nil {
-			fmt.Println("line", line, k)
+		//	fmt.Println("line", line, k)
 			return nil, nil, CError{i.Error(), []string{"strconv.ParseFloat", "xyzReadSnap"}}
 		}
 	}

handy.go

@@ -313,25 +313,25 @@ func MakeWater(a1, a2 *v3.Matrix, distance, angle float64, oxygen bool) *v3.Matr
 	dist := v3.Zeros(1)
 	dist.Sub(a1, a2)
 	a1a2dist := dist.Norm(0)
-	fmt.Println("ala2dist", a1a2dist, distance) ////////////////7777
+//	fmt.Println("ala2dist", a1a2dist, distance) ////////////////7777
 	w.Scale(distance+a1a2dist, w)
 	w.Add(w, a1)
 	for i := 0; i <= 1; i++ {
 		o := water.VecView(0)
 		w = water.VecView(i + 1)
 		w.Copy(o)
-		fmt.Println("w1", w) ////////
+//		fmt.Println("w1", w) ////////
 		w.Sub(w, a2)
-		fmt.Println("w12", w) ///////////////
+//		fmt.Println("w12", w) ///////////////
 		w.Unit(w)
-		fmt.Println("w4", w)
+//		fmt.Println("w4", w)
 		w.Scale(WaterOHDist+distance, w)
-		fmt.Println("w3", w, WaterOHDist, distance)
+//		fmt.Println("w3", w, WaterOHDist, distance)
 		o.Sub(o, a2)
 		t, _ := v3.NewMatrix([]float64{0, 0, 1})
 		upp := v3.Zeros(1)
 		upp.Cross(w, t)
-		fmt.Println("upp", upp, w, t)
+//		fmt.Println("upp", upp, w, t)
 		upp.Add(upp, o)
 		upp.Add(upp, a2)
 		//water.SetMatrix(3,0,upp)

qm/turbomole.go

@@ -445,7 +445,7 @@ func (O *TMHandle) Run(wait bool) (err error) {
 	os.Chdir(O.inputname)
 	defer os.Chdir("..")
 	filename := strings.Fields(O.command)
-	fmt.Println("nohup " + O.command + " > " + filename[0] + ".out")
+	//fmt.Println("nohup " + O.command + " > " + filename[0] + ".out")
 	command := exec.Command("sh", "-c", "nohup "+O.command+" >"+filename[0]+".out")
 	if wait == true {
 		err = command.Run()

v3/gocoords.go

@@ -99,7 +99,7 @@ func (F *Matrix) DelVec(A *Matrix, i int) {
 	tempF1.Copy(tempA1)
 	//now the other part
 	//	if i != ar-1 {
-	fmt.Println("options", ar, i, ar-i-1)
+	//fmt.Println("options", ar, i, ar-i-1)
 	if i < ar-1 {
 		tempA2 := A.View(i+1, 0, ar-i-1, ac) //The magic happens here
 		tempF2 := F.View(i, 0, ar-i-1, fc)

v3/gonum_purego.go

@@ -206,7 +206,7 @@ func EigenWrap(in *VecMatrix, epsilon float64) (*VecMatrix, []float64, error) {
 				return eig.evecs, evals[:], notOrthogonal
 			}
 		}
-		fmt.Println("VECTORS", eig.evecs) //////////////////////////
+	//	fmt.Println("VECTORS", eig.evecs) //////////////////////////
 		if math.Abs(vectori.Norm(0)-1) > epsilon {
 			//Of course I could just normalize the vectors instead of complaining.
 			//err= fmt.Errorf("Vectors not normalized %s",err.Error())
@@ -314,7 +314,7 @@ func (F *Dense) Copy(A Matrix) {
 	ar, ac := A.Dims()
 	fr, fc := F.Dims()
 	if ac > fc || ar > fr {
-		fmt.Println(ar, ac, fr, fc)
+	//	fmt.Println(ar, ac, fr, fc)
 		panic(gnErrShape)
 	}