code: public Matrix4d transfToOrthonormal(Matrix4d m) {
  		Vector3d trans = new Vector3d(m.m03,m.m13,m.m23);
  		transfToOrthonormal(trans);
  
  		Matrix3d rot = new Matrix3d();
  		m.getRotationScale(rot);
  		// see Giacovazzo section 2.E, eq. 2.E.1
  		// Rprime = MT-1 * R * MT
  		rot.mul(this.getMTranspose());
  		rot.mul(this.getMTransposeInv(),rot);
  
  		return new Matrix4d(rot,trans,1.0);
  	}
docs: Transform given Matrix4d in crystal basis to the orthonormal basis using
  the PDB axes convention (NCODE=1)
  @param m
  @return