float getPhase(Point3f pt) {
		ptPsi.set(pt);
		getCalcPoint(ptPsi);
		switch (colorPhase) {
		case 0:
			return (hydrogenAtomPsiAt(ptPsi, psi_n, psi_l, psi_m) > 0 ? 1 : -1);
		case -1:
		case 1:
			return (pt.x > 0 ? 1 : -1);
		case 2:
			return (pt.y > 0 ? 1 : -1);
		case 3:
			return (pt.z > 0 ? 1 : -1);
		case 4:
			return (pt.x * pt.y > 0 ? 1 : -1);
		case 5:
			return (pt.y * pt.z > 0 ? 1 : -1);
		case 6:
			return (pt.x * pt.z > 0 ? 1 : -1);
		case 7:
			return (pt.x * pt.x - pt.y * pt.y > 0 ? 1 : -1);
		case 8:
			return (pt.z * pt.z * 2f - pt.x * pt.x - pt.y * pt.y > 0 ? 1 : -1);
		}
		return 1;
	}