- for (iy = 0; iy < ny; iy++)
- vOut[ix][iy] = sqrt (vIn[ix][iy]);
- histString = "Sqrt transformation";
+ for (iy = 0; iy < ny; iy++)\r
+ if (vIn[ix][iy] < 0)
+ vOut[ix][iy] = sqrt (-vIn[ix][iy]);\r
+ else\r
+ vOut[ix][iy] = sqrt (vIn[ix][iy]);\r
+ histString = "Square root transformation";