compute_layer(input_activation[][],input_indices[][], weights[][][][], output_activation[][], output_index[][]) {
	for each ix,iy	// image coordinates
		float[] tmp_output;
		clear(tmp_output);
		for each fx,fy,fz // filter coordinates  .. WITHOUT SPARSE INPUT, THIS IS 4D TRAVERSAL, NOT 3D
			tmp_output[fz]+= input[ix][iy] * weights[fx][fy][ input_indicies[ix][iy] ][fz];

		output_index[ix][iy]=index_of_max_value_of(tmp_output);
		output_activation[ix][iy] = max_value_of(tmp_output)
}