diff --git a/mrdna/segmentmodel.py b/mrdna/segmentmodel.py index 4093096586b29cddfc3c31783fc5bb4033ee05ea..3c7c05c7ff2ee5cf18d118033c2b1fa76f4e1e63 100644 --- a/mrdna/segmentmodel.py +++ b/mrdna/segmentmodel.py @@ -1010,17 +1010,17 @@ class SingleStrandedSegment(Segment): seg2 = other.container end1 = self.end3 end2 = other - # assert( other.type_ == "end5" ) - if (other.type_ is not "end5"): - print("WARNING: code does not prevent connecting 3prime to 3prime, etc") + assert(other.type_ != "end3") + # if (other.type_ is not "end5"): + # print("WARNING: code does not prevent connecting 3prime to 3prime, etc") else: seg1 = other.container seg2 = self end1 = other end2 = self.start - # assert( other.type_ == "end3" ) - if (other.type_ is not "end3"): - print("WARNING: code does not prevent connecting 3prime to 3prime, etc") + assert(other.type_ != "end5") + # if (other.type_ is not "end3"): + # print("WARNING: code does not prevent connecting 3prime to 3prime, etc") ## Remove other connections involving these points if end1.connection is not None: @@ -1840,7 +1840,8 @@ class SegmentModel(ArbdModel): for b1,b2 in intra_beads: # assert( not np.isclose( np.linalg.norm(b1.collapsedPosition() - b2.collapsedPosition()), 0 ) ) if np.linalg.norm(b1.collapsedPosition() - b2.collapsedPosition()) < 1: - print("WARNING: some beads are very close") + # print("WARNING: some beads are very close") + ... parent = self._getParent(b1,b2)