thrust

Just replace double …. with double& ….. 3 solved overload operator [ ] to enable element assignment