I got varbvsmixoptimize working finally in R.

I also discovered a major bug in my R code where

x <- 1 + 2 + 3
     + 4 + 5

does *not* set x to the sum of numbers 1 through 5.