 CollOfScalar DeviceGrid::norm_of_cells(const thrust::device_vector<int>& cells,
 				       const bool full) const {
     if (full) {
-	CollOfScalar out(number_of_cells_,0);
+	CollOfScalar out(number_of_cells_);
 	cudaStatus_ = cudaMemcpy( out.data(), cell_volumes_, 
 				  sizeof(double)*number_of_cells_,
 				  cudaMemcpyDeviceToDevice);

	return out;
    }
    else {
	CollOfScalar out(faces.size(),0);
	kernelSetup s = out.setup();
	const int* faces_ptr = thrust::raw_pointer_cast( &faces[0] );
	normKernel<<<s.grid, s.block>>>( out.data(), faces_ptr, faces.size(),

	return out;
    }
    else {
	CollOfScalar out(cells.size());
	kernelSetup s = out.setup();
	const int* cells_ptr = thrust::raw_pointer_cast( &cells[0] );
	normKernel<<<s.grid, s.block>>>( out.data(), cells_ptr, cells.size(),
