icn3d 3.8.1 → 3.9.0

package/icn3d.js

@@ -6862,10 +6862,11 @@ class Tube {
                     firstAtom = atom;
                 }
 
-                //if (index > 0 && (currentChain !== atom.chain || Math.abs(atom.coord.x - prevAtom.coord.x) > maxDist || Math.abs(atom.coord.y - prevAtom.coord.y) > maxDist || Math.abs(atom.coord.z - prevAtom.coord.z) > maxDist
-                //  || (currentResi + 1 !== atom.resi && (Math.abs(atom.coord.x - prevAtom.coord.x) > maxDist2 || Math.abs(atom.coord.y - prevAtom.coord.y) > maxDist2 || Math.abs(atom.coord.z - prevAtom.coord.z) > maxDist2) )
+                atom.structure + '_' + atom.chain + '_' + (parseInt(atom.resi) - 1).toString();
+
                 if (index > 0 && (currentChain !== atom.chain || Math.abs(atom.coord.x - prevAtom.coord.x) > maxDist || Math.abs(atom.coord.y - prevAtom.coord.y) > maxDist || Math.abs(atom.coord.z - prevAtom.coord.z) > maxDist
-                  || (parseInt(currentResi) + 1 < parseInt(atom.resi) && (Math.abs(atom.coord.x - prevAtom.coord.x) > maxDist2 || Math.abs(atom.coord.y - prevAtom.coord.y) > maxDist2 || Math.abs(atom.coord.z - prevAtom.coord.z) > maxDist2) )
+//                  || (parseInt(currentResi) + 1 < parseInt(atom.resi) && (Math.abs(atom.coord.x - prevAtom.coord.x) > maxDist2 || Math.abs(atom.coord.y - prevAtom.coord.y) > maxDist2 || Math.abs(atom.coord.z - prevAtom.coord.z) > maxDist2) && ic.firstAtomObjCls.getFirstCalphaAtomObj(ic.residues[resid]) && ic.firstAtomObjCls.getFirstCalphaAtomObj(ic.residues[resid]).ss == 'helix')
+                  || (parseInt(currentResi) + 1 < parseInt(atom.resi) && (Math.abs(atom.coord.x - prevAtom.coord.x) > maxDist2 || Math.abs(atom.coord.y - prevAtom.coord.y) > maxDist2 || Math.abs(atom.coord.z - prevAtom.coord.z) > maxDist2))
                   ) ) {
                     if(bHighlight !== 2) {
                         if(!isNaN(firstAtom.resi) && !isNaN(prevAtom.resi)) {
@@ -7459,7 +7460,15 @@ class Strand {
                 }
 
                 let maxDist = 6.0;
-                let bBrokenSs = (prevCoorCA && Math.abs(currentCA.x - prevCoorCA.x) > maxDist) || (prevCoorCA && Math.abs(currentCA.y - prevCoorCA.y) > maxDist) || (prevCoorCA && Math.abs(currentCA.z - prevCoorCA.z) > maxDist);
+                //let bBrokenSs = (prevCoorCA && Math.abs(currentCA.x - prevCoorCA.x) > maxDist) || (prevCoorCA && Math.abs(currentCA.y - prevCoorCA.y) > maxDist) || (prevCoorCA && Math.abs(currentCA.z - prevCoorCA.z) > maxDist);
+                let bBrokenSs = !atoms.hasOwnProperty(atom.serial) || (prevCoorCA && Math.abs(currentCA.x - prevCoorCA.x) > maxDist) || (prevCoorCA && Math.abs(currentCA.y - prevCoorCA.y) > maxDist) || (prevCoorCA && Math.abs(currentCA.z - prevCoorCA.z) > maxDist);
+
+                if(bBrokenSs && atom.ss === 'sheet') {
+                    bSheetSegment = true;
+                }
+                else if(bBrokenSs && atom.ss === 'helix') {
+                    bHelixSegment = true;
+                }
 
                 if ((atom.ssbegin || atom.ssend || (drawnResidueCount === totalResidueCount - 1) || bBrokenSs) && pnts[0].length > 0 && bSameChain) {
                     let atomName = 'CA';
@@ -11460,17 +11469,32 @@ class Surface {
                 //vertexColors.push(ic.atoms[verts[vb].atomid].color);
                 //vertexColors.push(ic.atoms[verts[vc].atomid].color);
 
-                colorArray[offset2++] = ic.atoms[verts[va].atomid].color.r;
-                colorArray[offset2++] = ic.atoms[verts[va].atomid].color.g;
-                colorArray[offset2++] = ic.atoms[verts[va].atomid].color.b;
+                if(type == 21 || type == 22 || type == 23) { // potential on surface
+                    colorArray[offset2++] = verts[va].color.r;
+                    colorArray[offset2++] = verts[va].color.g;
+                    colorArray[offset2++] = verts[va].color.b;
 
-                colorArray[offset2++] = ic.atoms[verts[vb].atomid].color.r;
-                colorArray[offset2++] = ic.atoms[verts[vb].atomid].color.g;
-                colorArray[offset2++] = ic.atoms[verts[vb].atomid].color.b;
+                    colorArray[offset2++] = verts[vb].color.r;
+                    colorArray[offset2++] = verts[vb].color.g;
+                    colorArray[offset2++] = verts[vb].color.b;
 
-                colorArray[offset2++] = ic.atoms[verts[vc].atomid].color.r;
-                colorArray[offset2++] = ic.atoms[verts[vc].atomid].color.g;
-                colorArray[offset2++] = ic.atoms[verts[vc].atomid].color.b;
+                    colorArray[offset2++] = verts[vc].color.r;
+                    colorArray[offset2++] = verts[vc].color.g;
+                    colorArray[offset2++] = verts[vc].color.b;
+                }
+                else {
+                    colorArray[offset2++] = ic.atoms[verts[va].atomid].color.r;
+                    colorArray[offset2++] = ic.atoms[verts[va].atomid].color.g;
+                    colorArray[offset2++] = ic.atoms[verts[va].atomid].color.b;
+    
+                    colorArray[offset2++] = ic.atoms[verts[vb].atomid].color.r;
+                    colorArray[offset2++] = ic.atoms[verts[vb].atomid].color.g;
+                    colorArray[offset2++] = ic.atoms[verts[vb].atomid].color.b;
+    
+                    colorArray[offset2++] = ic.atoms[verts[vc].atomid].color.r;
+                    colorArray[offset2++] = ic.atoms[verts[vc].atomid].color.g;
+                    colorArray[offset2++] = ic.atoms[verts[vc].atomid].color.b;
+                }
 
                 //var normals = [];
                 //normals.push(normalArrayIn[va]);
@@ -11722,7 +11746,7 @@ class ShareLink {
 
            if(!bPngHtml) {
                if(ic.bInputfile && !ic.bInputUrlfile) {
-                   var aaa = 1; //alert("Share Link does NOT work when the data is from custom files. Please save 'iCn3D PNG Image' in the File menu and open it in iCn3D.");
+                   var aaa = 1; //alert("Share Link does NOT work when the data are from custom files. Please save 'iCn3D PNG Image' in the File menu and open it in iCn3D.");
                    return;
                }
                if(bTooLong) {
@@ -14360,9 +14384,7 @@ class HBond {
                     }
 
                     //if(result > kMaxHBondEnergy) {
-                    if(atom.ss == 'helix' && atomHbond[j].ss == 'helix' && result > kMaxHBondEnergy) {
-                        continue;
-                    }
+                    if(atom.ss == 'helix' && atomHbond[j].ss == 'helix' && result > kMaxHBondEnergy) ;
                 }
               }
               else {
@@ -15511,7 +15533,8 @@ class GetGraph {
         let  strokewidth = '1';
         let  textcolor = '#000';
         let  fontsize = '6px'; // '6';
-        let  html = (bAfMap) ? "<g>" : "<g class='icn3d-node' resid='" + resid + "' >";
+        //let  html = (bAfMap) ? "<g>" : "<g class='icn3d-node' resid='" + resid + "' >";
+        let  html = "<g class='icn3d-node' resid='" + resid + "' >";
         html += "<title>" + node.id + "</title>";
         if(bVertical) {
             html += "<circle cx='" + y + "' cy='" + x + "' r='" + r + "' fill='" + color + "' stroke-width='" + strokewidth + "' stroke='" + strokecolor + "' resid='" + resid + "' />";
@@ -15554,10 +15577,10 @@ class GetGraph {
             else if(node.s == 'b') {
                 nodeArray2.push(node);
             }
-            //else if(node.s == 'ab') {
-            //    nodeArray1.push(node);
-            //    nodeArray2.push(node);
-            //}
+            else if(node.s == 'ab') {
+                nodeArray1.push(node);
+                nodeArray2.push(node);
+            }
         }
         // sort array
         nodeArray1.sort(function(a,b) {
@@ -15589,6 +15612,7 @@ class GetGraph {
               let  atom = ic.firstAtomObjCls.getFirstAtomObj(ic.residues[resid]);
               resid2color[resid] = atom.color.getHexString().toUpperCase();
           }
+
           let  target2resid = {};
           for(let i = 0, il = graphJson.nodes.length; i < il; ++i) {
               let  node = graphJson.nodes[i];
@@ -15614,6 +15638,7 @@ class GetGraph {
           }
           ic.graphStr = JSON.stringify(graphJson);
       }
+
       if(ic.bGraph) ic.drawGraphCls.drawGraph(ic.graphStr, ic.pre + 'dl_graph');
       if(ic.bLinegraph) ic.lineGraphCls.drawLineGraph(ic.graphStr);
       if(ic.bScatterplot) ic.lineGraphCls.drawLineGraph(ic.graphStr, true);
@@ -15945,7 +15970,7 @@ class LineGraph {
                     let resid2 = chainid2 + '_' + idArrayB[4];
 
                     let mapping1, mapping2;
-                    
+
                     if(ic.chainsMapping[chainid1] && ic.chainsMapping[chainid1][resid1]
                       && ic.chainsMapping[chainid2] && ic.chainsMapping[chainid2][resid2]) { 
                         mapping1 = (nodeA.s == "a") ? ic.chainsMapping[chainid1][resid1] : ic.chainsMapping[chainid2][resid2];
@@ -15961,7 +15986,7 @@ class LineGraph {
                     }
                 } 
             }
-            
+
             // set linkArraySplitCommon and nameHashSplitCommon
             // set linkArraySplitDiff and nameHashSplitDiff
             let separatorCommon = "=>", separatorDiff = "==>", postCommon = "-", postDiff = "--";
@@ -16077,17 +16102,19 @@ class LineGraph {
 
                 width = (maxWidth + 2) * (r + gap) + 2 * marginX;
             }
-            let id;
+            let  id, graphWidth;
             if(bScatterplot) {
                 ic.scatterplotWidth = 2 * width;
+                graphWidth = ic.scatterplotWidth;
                 id = me.scatterplotid;
             } else {
                 ic.linegraphWidth = 2 * width;
+                graphWidth = ic.linegraphWidth;
                 id = me.linegraphid;
             }
             html =(strucArray.length == 0) ? "No interactions found for each structure<br><br>" :
                 "2D integration graph for " + strucArray.length + " structure(s) <b>" + strucArray + "</b>. There are three sections: \"Interactions\", \"Common interactions\", and \"Different interactions\". Each section has " + strucArray.length + " graphs.<br><br>";
-            html += "<svg id='" + id + "' viewBox='0,0," + width + "," + heightAll + "'>";
+            html += "<svg id='" + id + "' viewBox='0,0," + width + "," + heightAll + "' width='" + graphWidth + "px'>";
 
             let  result, heightFinal = 0;            
  
@@ -16172,15 +16199,15 @@ class LineGraph {
         // draw common interaction
         let label, postfix;
         if(bCommonDiff == 0) {
-            label = "Interactions in structure ";
+            label = "Interactions in ";
             postfix = "";
         }
         else if(bCommonDiff == 1) {
-            label = "Common interactions in structure ";
+            label = "Common interactions in ";
             postfix = "_common";
         }
         else if(bCommonDiff == 2) {
-            label = "Different interactions in structure ";
+            label = "Different interactions in ";
             postfix = "_diff";
         }
 
@@ -16237,7 +16264,7 @@ class LineGraph {
         // draw label
         if(label) {
             height += textHeight;
-            html += "<text x='" + margin1 + "' y='" + height + "' style='font-size:8px; font-weight:bold'>" + label + "</text>";
+            html += "<text x='" + margin + "' y='" + height + "' style='font-size:8px; font-weight:bold'>" + label + "</text>";
         }
 
         let  h1 = 30 + height,
@@ -16388,8 +16415,13 @@ class LineGraph {
         
         if(bAfMap && ic.hex2skip[link.c]) ;
         else if(bAfMap && ic.hex2id[link.c]) {
-            let id = ic.hex2id[link.c];
-            html += "<use href='#" + id + "' x='" +(pos2.x - halfSize).toString() + "' y='" +(pos1.y - halfSize).toString() + "' />";
+            ic.hex2id[link.c];
+//            html += "<use href='#" + id + "' x='" +(pos2.x - halfSize).toString() + "' y='" +(pos1.y - halfSize).toString() + "' />";
+
+            //html += "<g class='icn3d-interaction' resid1='" + resid1 + "' resid2='" + resid2 + "' >";
+            //html += "<title>Interaction of residue " + node1.id + " with residue " + node2.id + "</title>";
+            html += "<rect class='icn3d-interaction' resid1='" + resid1 + "' resid2='" + resid2 + "' x='" +(pos2.x - halfSize).toString() + "' y='" +(pos1.y - halfSize).toString() + "' width='" + rectSize + "' height='" + rectSize + "' fill='" + strokecolor + "' stroke-width='" + linestrokewidth + "' stroke='" + strokecolor + "' />";
+            //html += "</g>";
         }
         else {
             html += "<g class='icn3d-interaction' resid1='" + resid1 + "' resid2='" + resid2 + "' >";
@@ -19353,19 +19385,19 @@ class Scap {
         //snp: 6M0J_E_484_K,6M0J_E_501_Y,6M0J_E_417_N
         let snpStr = '';
         let snpArray = snp.split(','); //stru_chain_resi_snp
-        let atomHash = {}, residArray = [];
+        let atomHash = {}, snpResidArray = [];
         for(let i = 0, il = snpArray.length; i < il; ++i) {
             let idArray = snpArray[i].split('_'); //stru_chain_resi_snp
 
             let resid = idArray[0] + '_' + idArray[1] + '_' + idArray[2];
             atomHash = me.hashUtilsCls.unionHash(atomHash, ic.residues[resid]);
-            residArray.push(resid);
+            snpResidArray.push(resid);
 
             snpStr += idArray[1] + '_' + idArray[2] + '_' + idArray[3];
             if(i != il -1) snpStr += ',';
         }
 
-        let selectSpec = ic.resid2specCls.residueids2spec(residArray);
+        let selectSpec = ic.resid2specCls.residueids2spec(snpResidArray);
         let select = "select " + selectSpec;
 
         let bGetPairs = false;
@@ -19374,7 +19406,7 @@ class Scap {
         let result = ic.showInterCls.pickCustomSphere_base(radius, atomHash, ic.atoms, false, false, undefined, select, bGetPairs);
 
 
-        residArray = Object.keys(result.residues);
+        let residArray = Object.keys(result.residues);
         ic.hAtoms = {};
         for(let index = 0, indexl = residArray.length; index < indexl; ++index) {
           let residueid = residArray[index];
@@ -19421,6 +19453,18 @@ class Scap {
               let bAddition = true;
               let hAtom1 = me.hashUtilsCls.cloneHash(ic.hAtoms);
 
+              // the wild type is the reference
+              for(let serial in hAtom1) {
+                  let atom = ic.atoms[serial];
+                  let chainid = atom.structure + '_' + atom.chain;
+                  let resid = chainid + '_' + atom.resi;
+
+                  if(!ic.chainsMapping.hasOwnProperty(chainid)) {
+                    ic.chainsMapping[chainid] = {};
+                  }
+                  ic.chainsMapping[chainid][resid] = me.utilsCls.residueName2Abbr(atom.resn) + atom.resi;
+              }
+
               ic.hAtoms = {};
               ic.loadPDBCls.loadPDB(pdbData, pdbid, false, false, bAddition);
               let hAtom2 = me.hashUtilsCls.cloneHash(ic.hAtoms);
@@ -19433,17 +19477,31 @@ class Scap {
 
               ic.opts['color'] = 'chain';
               ic.setColorCls.setColorByOptions(ic.opts, ic.dAtoms);
-
               for(let serial in hAtom2) {
                   let atom = ic.atoms[serial];
                   if(!atom.het) {
                       // use the same color as the wild type
-                      let resid = atom.structure.substr(0, 4) + '_' + atom.chain + '_' + atom.resi;
+                      let resid = atom.structure.substr(0, atom.structure.length - 1) + '_' + atom.chain + '_' + atom.resi;
 
                       let atomWT = ic.firstAtomObjCls.getFirstAtomObj(ic.residues[resid]);
                       ic.atoms[serial].color = atomWT.color;
                       ic.atomPrevColors[serial] = atomWT.color;
                   }
+
+                  let chainid = atom.structure + '_' + atom.chain;
+                  let resid = chainid + '_' + atom.resi;
+                  let residWT = atom.structure.substr(0, atom.structure.length - 1) + '_' + atom.chain + '_' + atom.resi;
+
+                  if(!ic.chainsMapping.hasOwnProperty(chainid)) {
+                    ic.chainsMapping[chainid] = {};
+                  }
+                  ic.chainsMapping[chainid][resid] = me.utilsCls.residueName2Abbr(atom.resn) + atom.resi;
+                  // use the wild type as reference
+
+                  if(snpResidArray.indexOf(residWT) != -1) {
+                      let atomWT = ic.firstAtomObjCls.getFirstAtomObj(ic.residues[residWT]);
+                      ic.chainsMapping[chainid][resid] = me.utilsCls.residueName2Abbr(atomWT.resn) + atomWT.resi;
+                  }
               }
 
               if(bPdb) {
@@ -26455,9 +26513,10 @@ class LoadScript {
               let  strArray = ic.commands[i].split("|||");
 
               if(Object.keys(ic.proteins).length > 0 && ic.mmdb_data === undefined &&(ic.bAjax3ddomain === undefined || !ic.bAjax3ddomain)) {
-                  $.when(thisClass.applyCommand3ddomain(strArray[0].trim())).then(function() {
-                      thisClass.execCommandsBase(i + 1, end, steps);
-                  });
+                  //$.when(thisClass.applyCommand3ddomain(strArray[0].trim())).then(function() {
+                    thisClass.applyCommand3ddomain(strArray[0].trim());    
+                    thisClass.execCommandsBase(i + 1, end, steps);
+                  //});
               }
               else {
                   if(Object.keys(ic.proteins).length > 0) {
@@ -26479,8 +26538,9 @@ class LoadScript {
                 &&(ic.bAjax3ddomain === undefined || !ic.bAjax3ddomain || ic.mmdb_data === undefined) ) {
                   $.when(thisClass.applyCommandClinvar(strArray[0].trim()))
                     .then(thisClass.applyCommandSnp(strArray[0].trim()))
-                    .then(thisClass.applyCommand3ddomain(strArray[0].trim()))
+                    //.then(thisClass.applyCommand3ddomain(strArray[0].trim()))
                     .then(function() {
+                      thisClass.applyCommand3ddomain(strArray[0].trim());
                       ic.annotationCls.setAnnoTabAll();
 
                       thisClass.execCommandsBase(i + 1, end, steps);
@@ -26499,8 +26559,9 @@ class LoadScript {
               else if(Object.keys(ic.proteins).length > 0 &&(ic.bAjaxClinvar === undefined || !ic.bAjaxClinvar)
                 &&(ic.bAjax3ddomain === undefined || !ic.bAjax3ddomain || ic.mmdb_data === undefined)) {
                   $.when(thisClass.applyCommandClinvar(strArray[0].trim()))
-                    .then(thisClass.applyCommand3ddomain(strArray[0].trim()))
+                    //.then(thisClass.applyCommand3ddomain(strArray[0].trim()))
                     .then(function() {
+                      thisClass.applyCommand3ddomain(strArray[0].trim());
                       ic.annotationCls.setAnnoTabAll();
 
                       thisClass.execCommandsBase(i + 1, end, steps);
@@ -26508,10 +26569,11 @@ class LoadScript {
               }
               else if(Object.keys(ic.proteins).length > 0 &&(ic.bAjax3ddomain === undefined || !ic.bAjax3ddomain || ic.mmdb_data === undefined)
                 &&(ic.bAjaxSnp === undefined || !ic.bAjaxSnp)) {
-                  $.when(thisClass.applyCommand3ddomain(strArray[0].trim()))
-                    .then(thisClass.applyCommandSnp(strArray[0].trim()))
+                  //$.when(thisClass.applyCommand3ddomain(strArray[0].trim()))
+                  $.when(thisClass.applyCommandSnp(strArray[0].trim()))
                     .then(function() {
-                      ic.annotationCls.setAnnoTabAll();
+                        thisClass.applyCommand3ddomain(strArray[0].trim());
+                        ic.annotationCls.setAnnoTabAll();
 
                       thisClass.execCommandsBase(i + 1, end, steps);
                   });
@@ -26533,12 +26595,14 @@ class LoadScript {
                   });
               }
               else if(Object.keys(ic.proteins).length > 0 &&(ic.bAjax3ddomain === undefined || !ic.bAjax3ddomain || ic.mmdb_data === undefined) ) {
-                  $.when(thisClass.applyCommand3ddomain(strArray[0].trim()))
-                    .then(function() {
+                  //$.when(thisClass.applyCommand3ddomain(strArray[0].trim()))
+                  // .then(function() {
+                      thisClass.applyCommand3ddomain(strArray[0].trim());
+
                       ic.annotationCls.setAnnoTabAll();
 
                       thisClass.execCommandsBase(i + 1, end, steps);
-                  });
+                  //});
               }
               else {
                   if(Object.keys(ic.proteins).length > 0) {
@@ -26805,12 +26869,13 @@ class LoadScript {
                         thisClass.applyCommand3ddomain(lastCommand);
                     }
                     else if(lastCommand.indexOf('set annotation all') == 0) {
-                        //$.when(thisClass.applyCommandAnnotationsAndCddSite(lastCommand))
-                        //    .then(thisClass.applyCommandSnpClinvar(lastCommand))
                         $.when(thisClass.applyCommandClinvar(lastCommand))
                             .then(thisClass.applyCommandSnp(lastCommand))
-                            .then(thisClass.applyCommand3ddomain(lastCommand));
-                        ic.annotationCls.setAnnoTabAll();
+                            //.then(thisClass.applyCommand3ddomain(lastCommand));
+                            .then(function() {
+                                thisClass.applyCommand3ddomain(lastCommand);
+                                ic.annotationCls.setAnnoTabAll();
+                            });
                     }
                     else if(lastCommand.indexOf('view interactions') == 0 && me.cfg.align !== undefined) {
                         thisClass.applyCommandViewinteraction(lastCommand);
@@ -27202,11 +27267,11 @@ class LoadScript {
       let  thisClass = this;
 
       // chain functions together
-      ic.deferred3ddomain = $.Deferred(function() {
+      //ic.deferred3ddomain = $.Deferred(function() {
           thisClass.applyCommand3ddomainBase(command);
-      }); // end of me.deferred = $.Deferred(function() {
+      //}); // end of me.deferred = $.Deferred(function() {
 
-      return ic.deferred3ddomain.promise();
+      //return ic.deferred3ddomain.promise();
     }
 
     applyCommandViewinteractionBase(command) { let  ic = this.icn3d, me = ic.icn3dui;
@@ -27545,6 +27610,1037 @@ class ResizeCanvas {
     }
 }
 
+/*
+ * @author Jiyao Wang <wangjiy@ncbi.nlm.nih.gov> / https://github.com/ncbi/icn3d
+ * Modified from Tom Madej's C++ code
+*/
+
+//import * as THREE from 'three';
+
+class Domain3d {
+    constructor(icn3d) {
+		this.icn3d = icn3d;
+
+        this.dcut = 8; // threshold for C-alpha interactions
+
+		this.MAX_SSE = 512;
+
+        //let this.ctc_cnt[this.MAX_SSE][this.MAX_SSE];		// contact count matrix
+        this.ctc_cnt = [];
+        for(let i = 0; i < this.MAX_SSE; ++i) {
+            this.ctc_cnt[i] = [];
+        }
+
+        //let this.elt_size[this.MAX_SSE];			// element sizes in residues
+        this.elt_size = [];
+        this.elt_size.length = this.MAX_SSE;
+
+        //let this.group_num[this.MAX_SSE];			// indicates required element groupings
+        this.group_num = [];
+        this.group_num.length = this.MAX_SSE;
+
+        this.split_ratio = 0.0;			//let // splitting ratio
+        this.min_size = 0;				// min required size of a domain
+        this.min_sse = 0;				// min number of SSEs required in a domain
+        this.max_csz = 0;				// max size of a cut, i.e. number of points
+        this.mean_cts = 0.0;				// mean number of contacts in a domain
+        this.c_delta = 0;				// cut set parameter
+        this.nc_fact = 0.0;				// size factor for internal contacts
+        //let this.elements[2*this.MAX_SSE];			// sets of this.elements to be split
+        this.elements = [];
+        this.elements.length = 2*this.MAX_SSE;
+
+        //let this.stack[this.MAX_SSE];			// this.stack of sets (subdomains) to split
+        this.stack = [];
+        this.stack.length = this.MAX_SSE;
+
+        this.top = 0;					// this.top of this.stack
+        //let this.curr_prt0[this.MAX_SSE];			// current part 0 this.elements
+        this.curr_prt0 = [];
+        this.curr_prt0.length = this.MAX_SSE;
+
+        //let this.curr_prt1[this.MAX_SSE];			// current part 1 this.elements
+        this.curr_prt1 = [];
+        this.curr_prt1.length = this.MAX_SSE;
+
+        this.curr_ne0 = 0;				// no. of this.elements in current part 0
+        this.curr_ne1 = 0;				// no. of this.elements in current part 1
+        this.curr_ratio = 0.0;			// current splitting ratio
+        this.curr_msize = 0;				// min of current part sizes
+        //let this.parts[2*this.MAX_SSE];			// final partition into domains
+        this.parts = [];
+        this.parts.length = 2*this.MAX_SSE;
+
+        this.np = 0;					// next free location in this.parts[]
+        this.n_doms = 0;				// number of domains
+        //let this.save_ratios[this.MAX_SSE];		// this.saved splitting ratios
+        this.save_ratios = [];
+        this.save_ratios.length = this.MAX_SSE;
+
+        this.saved = 0;				// number of this.saved ratios
+	}
+
+	// Partition the set of this.elements on this.top of the this.stack based on the input cut.
+	// If the partition is valid and the ratio is smaller than the current one, then
+	// save it as the best partition so far encountered.  Various criteria are
+	// employed for valid partitions, as described below.
+	//
+
+	//update_partition(int* cut, let k, let n) { let ic = this.icn3d, me = ic.icn3dui;
+	update_partition(cut, k, n) { let ic = this.icn3d; ic.icn3dui;
+		let i, il, j, t, nc0, nc1, ncx, ne, ne0, ne1, elts = [], prt = []; //int
+		let size0, size1, prt0 = [], prt1 = []; // int
+        prt0.length = this.MAX_SSE;
+        prt1.length = this.MAX_SSE;
+		let f, r0; //let
+
+		// this.elements from the this.top of the this.stack 
+		//elts = &this.elements[this.stack[this.top - 1]];
+		for(i = this.stack[this.top - 1], il = this.elements.length; i < il; ++i) {
+			elts.push(this.elements[i]);
+		}
+
+		// generate the partition based on the cut //
+		for (i = ne = ne0 = ne1 = 0, prt = prt0, t = -1; i < k; i++) {
+			// write the this.elements into prt //
+			for (j = t + 1; j <= cut[i]; j++)
+				prt[ne++] = elts[j];
+
+			t = cut[i];
+
+			// switch the partition //
+			if (prt == prt0) {
+				ne0 = ne;
+				prt = prt1;
+				ne = ne1;
+			}
+			else {
+				ne1 = ne;
+				prt = prt0;
+				ne = ne0;
+			}
+		}
+
+		// finish with the last part //
+		for (j = t + 1; j < n; j++)
+			prt[ne++] = elts[j];
+
+		if (prt == prt0)
+			ne0 = ne;
+		else
+			ne1 = ne;
+
+		// don't split into two teeny this.parts! //
+		if ((ne0 < this.min_sse) && (ne1 < this.min_sse))
+			return cut;
+
+		// check to see if the partition splits any required groups //
+		for (i = 0; i < ne0; i++) {
+			t = this.group_num[prt0[i]];
+
+			for (j = 0; j < ne1; j++) {
+				if (t == this.group_num[prt1[j]])
+					return cut;
+			}
+		}
+
+		// compute the sizes of the this.parts //
+		for (i = size0 = 0; i < ne0; i++)
+			size0 += this.elt_size[prt0[i]];
+
+		for (i = size1 = 0; i < ne1; i++)
+			size1 += this.elt_size[prt1[i]];
+
+		// count internal contacts for part 0 //
+		for (i = nc0 = 0; i < ne0; i++) {
+			for (j = i; j < ne0; j++)
+				nc0 += this.ctc_cnt[prt0[i]][prt0[j]];
+		}
+
+		// count internal contacts for part 1 //
+		for (i = nc1 = 0; i < ne1; i++) {
+			for (j = i; j < ne1; j++)
+				nc1 += this.ctc_cnt[prt1[i]][prt1[j]];
+		}
+
+		// check globularity condition //
+		if ((1.0 * nc0 / size0 < this.mean_cts) ||
+			(1.0 * nc1 / size1 < this.mean_cts))
+			return cut;
+
+		// to handle non-globular pieces make sure nc0, nc1, are large enough //
+		nc0 = Math.max(nc0, this.nc_fact*size0);
+		nc1 = Math.max(nc1, this.nc_fact*size1);
+
+		// count inter-part contacts //
+		for (i = ncx = 0; i < ne0; i++) {
+			t = prt0[i];
+
+			for (j = 0; j < ne1; j++)
+				ncx += this.ctc_cnt[t][prt1[j]];
+		}
+
+		// compute the splitting ratio //
+		f = Math.min(nc0, nc1);
+		r0 = 1.0 * ncx / (f + 1.0);
+
+		if ((r0 >= this.curr_ratio + 0.01) || (r0 > this.split_ratio))
+			return cut;
+
+		// If the difference in the ratios is insignificant then take the split
+		// that most evenly partitions the domain.
+
+		if ((r0 > this.curr_ratio - 0.01) && (Math.min(size0, size1) < this.curr_msize))
+				return cut;
+
+		// if we get to here then keep this split //
+		for (i = 0; i < ne0; i++)
+			this.curr_prt0[i] = prt0[i];
+
+		for (i = 0; i < ne1; i++)
+			this.curr_prt1[i] = prt1[i];
+
+		this.curr_ne0 = ne0;
+		this.curr_ne1 = ne1;
+		this.curr_ratio = r0;
+		this.curr_msize = Math.min(size0, size1);
+
+		return cut;
+
+	} // end update_partition //
+
+
+
+	// // Run through the possible cuts of size k for a set of this.elements of size n.
+	//  *
+	//  * To avoid small protrusions, no blocks of consecutive this.elements of length <= this.c_delta
+	//  * are allowed.  An example where this is desirable is as follows.  Let's say you
+	//  * have a protein with 2 subdomains, one of them an alpha-beta-alpha sandwich.  It
+	//  * could then happen that one of the helices in the sandwich domain might make more
+	//  * contacts with the other subdomain than with the sandwich.  The correct thing to
+	//  * do is to keep the helix with the rest of the sandwich, and the "this.c_delta rule"
+	//  * enforces this.
+	//  //
+
+	cut_size(k, n) { let ic = this.icn3d; ic.icn3dui;
+		let i, j, cok, cut0 = []; //int
+		cut0.length = this.MAX_SSE;
+
+		for (i = 0; i < k; i++)
+			cut0[i] = i;
+
+		// enumerate cuts of length k //
+		while (1) {
+			// check block sizes in the cut //
+			for (i = cok = 1; i < k; i++) {
+				if (cut0[i] - cut0[i - 1] <= this.c_delta) {
+					cok = 0;
+					break;
+				}
+			}
+			if (cok && (cut0[k - 1] < n - 1))
+				cut0 = this.update_partition(cut0, k, n);
+
+			// generate the next k-tuple of positions //
+			for (j = k - 1; (j >= 0) && (cut0[j] == n - k + j); j--);
+
+			if (j < 0) break;
+
+			cut0[j]++;
+
+			for (i = j + 1; i < k; i++)
+				cut0[i] = cut0[i - 1] + 1;
+		}
+
+	} // end cut_size //
+
+
+
+	// // Process the set of this.elements on this.top of the this.stack.  We generate cut sets in
+	//  * a limited size range, generally from 1 to 5.  For each cut the induced
+	//  * partition is considered and its splitting parameters computed.  The cut
+	//  * that yields the smallest splitting ratio is chosen as the correct one, if
+	//  * the ratio is low enough.  The subdomains are then placed on the this.stack for
+	//  * further consideration.
+	//  *
+	//  * Subdomains with < this.min_sse SSEs are not allowed to split further, however,
+	//  * it is possible to trim fewer than this.min_sse SSEs from a larger domain.  E.g.
+	//  * a chain with 7 SSEs can be split into a subdomain with 5 SSEs and another
+	//  * with 2 SSEs, but the one with 2 SSEs cannot be split further.
+	//  *
+	//  * Note that the invariant is, that this.stack[top] always points to the next free
+	//  * location in this.elements[].
+	//  //
+
+	process_set() { let ic = this.icn3d; ic.icn3dui;
+		let i, il, k, n, t, k0, elts = []; //int
+
+		// count the this.elements //
+		//elts = &this.elements[this.stack[this.top - 1]];
+		for(i = this.stack[this.top - 1], il = this.elements.length; i < il; ++i) {
+			elts.push(this.elements[i]);
+		}
+
+		//for (n = 0; *elts > -1; n++, elts++);
+		for (n = 0; n < elts.length && elts[n] > -1; n++);
+
+		// try various cut sizes //
+		k0 = Math.min(n - 1, this.max_csz);
+		this.curr_ne0 = this.curr_ne1 = 0;
+		this.curr_ratio = 100.0;
+
+		for (k = 1; k <= k0; k++)
+			this.cut_size(k, n);
+
+		// pop this.stack //
+		this.top--;
+
+		if (this.curr_ne0 == 0) {
+			// no split took place, save part //
+			t = this.stack[this.top];
+
+			//for (elts = &this.elements[t]; *elts > -1; elts++)
+			//	parts[np++] = *elts;
+
+			for (i = t; i < this.elements.length && this.elements[i] > -1; i++)
+				this.parts[this.np++] = this.elements[i];
+
+			this.parts[this.np++] = -1;
+			this.n_doms++;
+		}
+		else {
+			this.save_ratios[this.saved++] = this.curr_ratio;
+
+			if (this.curr_ne0 > this.min_sse) {
+				// push on part 0 //
+				t = this.stack[this.top];
+
+				for (i = 0; i < this.curr_ne0; i++)
+					this.elements[t++] = this.curr_prt0[i];
+
+				this.elements[t++] = -1;
+				this.stack[++this.top] = t;
+			}
+			else {
+				// save part 0 //
+				for (i = 0; i < this.curr_ne0; i++)
+					this.parts[this.np++] = this.curr_prt0[i];
+
+				this.parts[this.np++] = -1;
+				this.n_doms++;
+			}
+
+			if (this.curr_ne1 > this.min_sse) {
+				// push on part 1 //
+				t = this.stack[this.top];
+
+				for (i = 0; i < this.curr_ne1; i++)
+					this.elements[t++] = this.curr_prt1[i];
+
+				this.elements[t++] = -1;
+				this.stack[++this.top] = t;
+			}
+			else {
+				// save part 1 //
+				for (i = 0; i < this.curr_ne1; i++)
+					this.parts[this.np++] = this.curr_prt1[i];
+
+				this.parts[this.np++] = -1;
+				this.n_doms++;
+			}
+		}
+	} // end process_set //
+
+
+
+	// Main driver for chain splitting. //
+	//process_all(let n) { let ic = this.icn3d, me = ic.icn3dui;
+	process_all(n) { let ic = this.icn3d; ic.icn3dui;
+		let i; //int
+
+		// initialize the this.stack //
+		this.top = 1;
+		this.stack[0] = this.np = this.n_doms = 0;
+		this.saved = 0;
+
+		for (i = 0; i < n; i++)
+			this.elements[i] = i;
+
+		this.elements[n] = -1;
+
+		// recursively split the chain into domains //
+		while (this.top > 0) {
+			this.process_set();
+		}
+	} // end process_all //
+
+	// Output the domains.  For S we number the this.elements 1, 2, ..., n. //
+	//output(let n, int* prts) { let ic = this.icn3d, me = ic.icn3dui;
+	output(n) { let ic = this.icn3d; ic.icn3dui;
+		let i, k; //int
+		
+		let prts = [];
+
+		// zap the output array //
+		for (i = 0; i < 2*n; i++)
+			prts.push(0);
+
+		// now write out the subdomains //
+		for (i = k = 0; k < this.n_doms; i++) {
+			prts[i] = this.parts[i] + 1;
+
+			if (this.parts[i] < 0)
+				k++;
+		}
+
+		return prts;
+	} // end output //
+
+
+
+	// // S-interface to the chain-splitting program.
+	//  *
+	//  * Explanation of parameters:
+	//  *
+	//  *	ne - number of secondary structure this.elements (SSEs)
+	//  *	cts - contact count matrix
+	//  *	elt_sz - sizes of SSEs
+	//  *	grps - element group indicators
+	//  *	sratio - splitting ratio
+	//  *	msize - min size of a split domain
+	//  *	m_sse - min number of SSEs required in a split part
+	//  *	mcsz - max cut size, i.e. max number of split points
+	//  *	avg_cts - mean number of internal contacts for a domain
+	//  *	c_delt - cut set parameter
+	//  *	ncf0 - size factor for number of internal contacts
+	//  *	prts - output listing of domains
+	//  *	n_saved - number of this.saved splitting ratios
+	//  *	ratios - splitting ratios
+	//  *	ret - success/failure indicator
+	//  *	verb - flag to turn off/on splitting information
+	//  //
+
+	//new_split_chain(let ne, let sratio, let msize, let m_sse, let mcsz, let avg_cts,
+	//	let c_delt, let ncf0, int* prts, int* n_saved, let* ratios) { let ic = this.icn3d, me = ic.icn3dui;
+	new_split_chain(ne, sratio, msize, m_sse, mcsz, avg_cts,
+		c_delt, ncf0, prts, n_saved, ratios) { let ic = this.icn3d; ic.icn3dui;
+		let i; //int
+
+		this.split_ratio = sratio;
+		this.min_size = msize;
+		this.min_sse = m_sse;
+		this.max_csz = mcsz;
+		this.mean_cts = avg_cts;
+		this.c_delta = c_delt;
+		this.nc_fact = ncf0;
+		
+		this.process_all(ne);
+		//this.output(ne, prts);
+		this.parts = this.output(ne);
+		n_saved = this.saved;
+		for (i = 0; i < this.saved; i++)
+			ratios[i] = this.save_ratios[i];
+
+		return n_saved;
+
+	} // end new_split_chain //
+
+	//
+	// Actually, here is a better method that is also simple!
+	//
+	// If there are N atoms (residues) this algorithm should usually run in
+	// time O(N^4/3), and usually even much faster!  In very unusual cases
+	// it could take quadratic time.  The key idea is that atoms are not
+	// infinitely compressible, i.e. only a fixed number will fit in a given
+	// region of space.  So if the protein is roughly spherical, there will
+	// only be O(N^1/3) atoms close to any given diameter.  Therefore, a
+	// bound on the number of iterations of the inner loop is O(N^1/3).
+	//
+	// For an elongated protein that happens to have the x-axis normal to
+	// the long axis, then it is possible for the inner loop to take time
+	// O(N), in which case the whole takes O(N^2).  But this should rarely,
+	// if ever, occur in practice.  It would also be possible beforehand to
+	// choose the axis with the largest variance.
+	//
+
+	// typedef struct res_struct {
+	// 	let rnum;
+	// 	let x, y, z;
+	// } ResRec;
+
+	//list< pair< pair< int, let >, let > >
+	//c2b_AlphaContacts(let n0, let* x0, let* y0, let* z0,
+	//	const let incr = 4, const let dcut = 8.0) { let ic = this.icn3d, me = ic.icn3dui;
+	c2b_AlphaContacts(n0, x0, y0, z0, dcut) { let ic = this.icn3d; ic.icn3dui;
+		//if(!incr) incr = 4;
+		if(!dcut) dcut = this.dcut;
+
+		let list_cts = [], list_rr = [];
+
+		for (let i = 0; i < n0; i++) {
+			// don't include residues with missing coordinates
+			//if ((x0[i] == MissingCoord) || (y0[i] == MissingCoord) || (z0[i] == MissingCoord))
+			if (!x0[i]|| !y0[i] || !z0[i])
+				continue;
+
+			//ResRec rr0;
+			let rr0 = {};
+			rr0.rnum = i + 1;
+			rr0.x = x0[i];
+			rr0.y = y0[i];
+			rr0.z = z0[i];
+			list_rr.push(rr0);
+		}
+		
+		list_rr.sort(function(rr1, rr2) {
+				return rr1.x - rr2.x;
+			});
+					
+		//let rrit1, rrit2, rrbeg;
+		let i, j, len = list_rr.length;
+
+		//for (rrit1 = list_rr.begin(); rrit1 != list_rr.end(); rrit1++) {
+		for (i = 0; i < len; ++i) {	
+			//ResRec rr1 = *rrit1;
+			let rr1 = list_rr[i];
+			let x1 = rr1.x;
+			let y1 = rr1.y;
+			let z1 = rr1.z;
+			//rrbeg = rrit1;
+			//rrbeg++;
+
+			//for (rrit2 = rrbeg; rrit2 != list_rr.end(); rrit2++) {
+			for (j = i + 1; j < len; ++j) {	
+				//ResRec rr2 = *rrit2;
+				let rr2 = list_rr[j];
+				if ((rr1.rnum - rr2.rnum <= 3) && (rr2.rnum - rr1.rnum <= 3)) continue;
+				let x2 = rr2.x;
+				let y2 = rr2.y;
+				let z2 = rr2.z;
+
+				if (x2 > x1 + dcut)
+					break;
+
+				// x1 <= x2 <= x1 + dcut so compare
+				let sum = (x1 - x2)*(x1 - x2);
+				sum += (y1 - y2)*(y1 - y2);
+				sum += (z1 - z2)*(z1 - z2);
+				let d0 = Math.sqrt(sum);
+				if (d0 > dcut) continue;
+				//pair< pair< int, let >, let > lpair;
+				//pair< int, let > rpair;
+				let lpair = {}, rpair = {};
+
+				if (rr1.rnum < rr2.rnum) {
+					rpair.first = rr1.rnum;
+					rpair.second = rr2.rnum;
+				}
+				else {
+					rpair.first = rr2.rnum;
+					rpair.second = rr1.rnum;
+				}
+
+				lpair.first = rpair;
+				lpair.second = d0;
+				list_cts.push(lpair);
+			}
+		}
+
+		return list_cts;
+
+	} // end c2b_AlphaContacts
+
+
+
+	//
+	// Creates a table, actually a graph, of the contacts between SSEs.
+	//
+
+	//static map< pair< int, let >, let >
+	//c2b_ContactTable(vector<int>& v1, vector<int>& v2) { let ic = this.icn3d, me = ic.icn3dui;
+	c2b_ContactTable(v1, v2) { let ic = this.icn3d; ic.icn3dui;
+		let cmap = {};
+		let n0 = v1.length; //unsigned int
+
+		if (n0 != v2.length) {
+			// problem!
+
+			return cmap;
+		}
+
+		for (let i = 0; i < n0; i++) {
+			let e1 = v1[i];
+			let e2 = v2[i];
+			//pair<int, int> epr;
+			//let epr = {};
+			//epr.first = e1;
+			//epr.second = e2;
+			let epr = e1 + '_' + e2;
+
+			//if (cmap.count(epr) == 0) {
+			if (!cmap[epr]) {	
+				cmap[epr] = 1;
+			}
+			else
+				cmap[epr]++;
+		}
+
+		return cmap;
+
+	} // end c2b_ContactTable
+
+	
+	//https://www.geeksforgeeks.org/number-groups-formed-graph-friends/
+	countUtil(ss1, sheetNeighbor, existing_groups) {
+		this.visited[ss1] = true;
+		if(!this.groupnum2sheet[existing_groups]) this.groupnum2sheet[existing_groups] = [];
+		this.groupnum2sheet[existing_groups].push(parseInt(ss1));
+
+		for(let ss2 in sheetNeighbor[ss1]) {
+			if (!this.visited[ss2]) {
+				this.countUtil(ss2, sheetNeighbor, existing_groups);  
+			}
+		}  
+	}
+
+	//
+	// Residue ranges of the Vast domains, per protein chain.
+	//
+
+	//
+	// Subdomain definition rules are as follows; let m0 = minSSE:
+	//
+	//     1. A subdomain with <= m0 SSEs cannot be split.
+	//
+	//     2. A subdomain cannot be split into two this.parts, both with < m0 SSEs.
+	//
+	//     3. However, a subdomain can be trimmed, i.e. split into two this.parts,
+	//        one with < m0 SSEs.
+	//
+	//c2b_NewSplitChain(string asymId, let seqLen, let* x0, let* y0, let* z0) { let ic = this.icn3d, me = ic.icn3dui;
+	// x0, y0, z0: array of x,y,z coordinates of C-alpha atoms
+	c2b_NewSplitChain(chnid, dcut) { let ic = this.icn3d; ic.icn3dui;
+		let x0 = [], y0 = [], z0 = [];
+
+		//substruct: array of secondary structures, each of which has the keys: From (1-based), To (1-based), Sheet (0 or 1)
+		let substruct = [];
+		// sheets: array of sheets, each of which has the key: sheet_num (beta sandwich has two sheets, e.g., 0 and 1), adj_strand1 (not used), adj_strand2
+		let sheets = [];
+
+		let substructItem = {};
+		let resiOffset = 0;
+		for(let i = 0; i < ic.chainsSeq[chnid].length; ++i) {
+            let resi = ic.chainsSeq[chnid][i].resi;
+			if(i == 0) {
+				resiOffset = resi - 1;
+
+				for(let j = 0; j < resiOffset; ++j) {
+					x0.push(undefined);
+					y0.push(undefined);
+					z0.push(undefined);					
+				}
+			}
+
+			let resid = chnid + "_" + resi;
+			let atom = ic.firstAtomObjCls.getFirstCalphaAtomObj(ic.residues[resid]);
+
+			if(atom) {
+				x0.push(atom.coord.x);
+				y0.push(atom.coord.y);
+				z0.push(atom.coord.z);
+			}
+			else {
+				x0.push(undefined);
+				y0.push(undefined);
+				z0.push(undefined);
+			}
+
+			if(!atom) continue;           
+
+			if(atom.ssend) {
+				substructItem.To = parseInt(resi);
+				substructItem.Sheet = (atom.ss == 'sheet') ? true : false;
+				substruct.push(substructItem);
+				substructItem = {};				
+			}
+
+			// a residue could be both start and end. check ssend first, then check ssbegin 
+			if(atom.ssbegin) {
+				substructItem.From = parseInt(resi);
+			}
+        }
+
+		let nsse = substruct.length;
+
+		if (nsse <= 3)
+			// too small, can't split or trim
+			return;
+
+		if (nsse > this.MAX_SSE) {
+			// we have a problem...
+
+			return;
+		}
+
+		let seqLen = ic.chainsSeq[chnid].length + resiOffset;
+
+		// get a list of Calpha-Calpha contacts
+		///list< pair< pair< int, let >, let > >
+		let cts = this.c2b_AlphaContacts(seqLen, x0, y0, z0, dcut);
+
+		//
+		// Produce a "map" of the SSEs, i.e. vec_sse[i] = 0 means residue i + 1
+		// is in a loop, and vec_sse[i] = k means residue i + 1 belongs to SSE
+		// number k.
+		//
+		let vec_sse = []; //vector<int>
+
+		for (let i = 0; i < seqLen; i++)
+			vec_sse.push(0);
+			
+		let hasSheets = false;
+
+		//substruct: array of secondary structures, each of which has the keys: From, To, Sheet (0, 1)
+		for (let i = 0; i < substruct.length; i++) {
+			//SSE_Rec sserec = substruct[i];
+			let sserec = substruct[i];
+			let From = sserec.From;
+			let To = sserec.To;
+			this.elt_size[i] = To - From + 1;
+
+			// double-check indexing OK???
+			for (let j = From; j <= To; j++)
+				vec_sse[j - 1] = i + 1;
+
+			//if (sserec.Sheet > 0)
+			if (sserec.Sheet)
+				hasSheets = true;
+		}
+
+		// produce the SSE contact lists
+		let vec_cts1 = [], vec_cts2 = [], vec_cts1a = [], vec_cts2a = [];
+
+		//for (ctsit = cts.begin(); ctsit != cts.end(); ctsit++) {
+		for (let i = 0, il = cts.length; i < il; ++i) {
+			//pair< pair< int, let >, let > epr = *ctsit;
+			//pair< int, let > respair = epr.first;
+			let epr = cts[i];
+			let respair = epr.first;
+			let sse1 = vec_sse[respair.first - 1];
+			let sse2 = vec_sse[respair.second - 1];
+			// could be 0 or null
+			if ((sse1 <= 0) || (sse2 <= 0) || !sse1 || !sse2) continue;
+			vec_cts1.push(sse1);
+			vec_cts2.push(sse2);
+			if (sse1 == sse2) continue;
+			vec_cts1a.push(sse1);
+			vec_cts2a.push(sse2);
+		}
+
+		// this symmetrizes the contact data
+		for (let i = 0; i < vec_cts1a.length; i++) {
+			vec_cts1.push(vec_cts2a[i]);
+			vec_cts2.push(vec_cts1a[i]);
+		}
+
+		// add dummy contacts
+		for (let i = 0; i < nsse; i++) {
+			vec_cts1.push(i + 1);
+			vec_cts2.push(i + 1);
+		}
+		
+		// create contact counts from the contacts/interactions
+		//map< pair< int, let >, let > ctable = this.c2b_ContactTable(vec_cts1, vec_cts2);
+		let ctable = this.c2b_ContactTable(vec_cts1, vec_cts2);
+
+		// neighbor list of each sheet
+		let sheetNeighbor = {};
+		for(let pair in ctable) {
+			let ssPair = pair.split('_'); // 1-based
+			let ss1 = parseInt(ssPair[0]);
+			let ss2 = parseInt(ssPair[1]);
+
+			// both are sheets
+			if(substruct[ss1 - 1].Sheet && substruct[ss2 - 1].Sheet) {
+				if(!sheetNeighbor[ss1]) sheetNeighbor[ss1] = {};
+				if(!sheetNeighbor[ss2]) sheetNeighbor[ss2] = {};
+
+				sheetNeighbor[ss1][ss2] = 1;
+				sheetNeighbor[ss2][ss1] = 1;
+			}
+		}
+
+		//https://www.geeksforgeeks.org/number-groups-formed-graph-friends/
+		let existing_groups = 0;
+		let sheet2sheetnum = {};
+		this.groupnum2sheet = {};
+		this.visited = {};
+		for (let ss1 in sheetNeighbor) {
+			this.visited[ss1] = false;
+		}
+
+		// get this.groupnum2sheet
+		for (let ss1 in sheetNeighbor) {
+			// If not in any group.
+			if (this.visited[ss1] == false) {
+				existing_groups++;
+				 
+				this.countUtil(ss1, sheetNeighbor, existing_groups);
+			}
+		}
+
+		// get sheet2sheetnum
+		// each neighboring sheet willbe represented by the sheet with the smallest sse 
+		for(let groupnum in this.groupnum2sheet) {
+			let ssArray = this.groupnum2sheet[groupnum].sort();
+			for(let i = 0, il = ssArray.length; i < il; ++i) {
+				sheet2sheetnum[ssArray[i]] = ssArray[0];
+			}
+		}
+
+		for (let i = 0; i < nsse; i++) {
+			if(substruct[i].Sheet) {				
+				let sheetsItem = {};
+				if(sheet2sheetnum[i+1]) {
+					sheetsItem.sheet_num = sheet2sheetnum[i+1];
+					sheetsItem.adj_strand2 = 1; 
+					sheetsItem.sse = i + 1; 
+				}
+				else {
+					sheetsItem.sheet_num = 0;
+					sheetsItem.adj_strand2 = 0; 
+					sheetsItem.sse = i + 1; 
+				}
+
+				sheets.push(sheetsItem);
+			}
+		}
+
+		//
+		// Correct for dummy contacts; they're present to ensure that the
+		// table gives the right result in the possible case there is an
+		// element with no contacts.
+		//
+		for (let i = 0; i < nsse; i++) {
+			for (let j = 0; j < nsse; j++) {
+				//pair<int, int> epr;
+				//let epr = {};
+				//epr.first = i + 1;
+				//epr.second = j + 1;
+				let epr = (i+1).toString() + '_' + (j+1).toString();
+
+				//if (ctable.count(epr) == 0)
+				if (!ctable[epr])
+					this.ctc_cnt[i][j] = 0;
+				else {
+					let cnt = ctable[epr];
+					if (i == j) cnt--; // subtract dummy contact
+					this.ctc_cnt[i][j] = cnt;
+					this.ctc_cnt[j][i] = cnt;
+				}
+			}
+		}
+
+		let minStrand = 6;
+
+		if (hasSheets) {
+			//sheets: array of sheets, each of which has the key: sheet_num (number of strands), adj_strand1, adj_strand2
+
+			let cnt = 0;
+
+			for (let i = 0; i < sheets.length; i++) {
+				//BetaSheet_Rec bsrec = sheets[i];
+				let bsrec = sheets[i];
+
+				//if ((bsrec.sheet_num > 0) && (this.elt_size[i] >= minStrand) && (bsrec.adj_strand2 != 0))
+				if ((bsrec.sheet_num > 0) && (this.elt_size[bsrec.sse - 1] >= minStrand) && (bsrec.adj_strand2 != 0))
+					cnt++;
+			}
+
+			for (let i = 0; i < nsse; i++) {
+				//this.group_num[i] = (cnt == 0) ? i + 1 : 0;
+				this.group_num[i] = i + 1;
+			}
+
+			if (cnt> 0) {
+				for (let i = 0; i < sheets.length; i++) {
+					let bsrec = sheets[i];
+					this.group_num[bsrec.sse - 1] = bsrec.sheet_num;
+				}
+			}
+		}
+		else {
+			for (let i = 0; i < nsse; i++)
+				this.group_num[i] = i + 1;
+		}
+
+		let sratio = 0.25;
+		let minSize = 25;
+		let maxCsz = 4;
+		let avgCts = 0.0;
+		let ncFact = 0.0;
+		let cDelta = 3;
+		let minSSE = 3;
+
+		// call the domain splitter
+		this.parts = [];
+		this.parts.length = 2*this.MAX_SSE;
+		let ratios = [];
+		ratios.length = this.MAX_SSE;
+		let n_saved = 0;
+
+		for (let i = 0; i < nsse; i++) {
+			this.parts[2*i] = this.parts[2*i + 1] = 0;
+			ratios[i] = 0.0;
+		}
+
+		n_saved = this.new_split_chain(nsse, sratio, minSize, minSSE, maxCsz, avgCts, cDelta, ncFact, this.parts, n_saved, ratios);
+
+		// save domain data
+		//list< vector< let > > list_parts;
+		let list_parts = [];
+
+		if (n_saved > 0) {
+			// splits occurred...
+			let j = 0;
+			
+			for (let i = 0; i <= n_saved; i++) {
+				//vector<int> sselst;
+				let sselst = [];
+				//sselst.clear();
+
+				while (j < 2*nsse) {
+					let sse0 = this.parts[j++];
+
+					if (sse0 == 0) {
+						list_parts.push(sselst);
+						break;
+					}
+					else
+						sselst.push(sse0);
+				}
+			}
+		}
+		
+		list_parts.sort(function(v1, v2) {
+				return v1[0] - v2[0];
+			});
+
+		// determine residue ranges for each subdomain
+		let subdomains = [];
+
+		//for (lplet = list_parts.begin(); lplet != list_parts.end(); lpint++) {
+		for (let index = 0, indexl = list_parts.length; index < indexl; ++index) {
+			//vector<int> prts = *lpint;
+			let prts = list_parts[index];
+			//vector<int> resflags;
+			//resflags.clear();
+			let resflags = [];
+
+			// a domain must have at least 3 SSEs...
+			if (prts.length <= 2) continue;
+
+			for (let i = 0; i < seqLen; i++)
+				resflags.push(0);
+
+			for (let i = 0; i < prts.length; i++) {
+				let k = prts[i] - 1;
+
+				if ((k < 0) || (k >= substruct.length)) {
+					return;
+				}
+
+				//SSE_Rec sserec = substruct[k];
+				let sserec = substruct[k];
+				let From = sserec.From;
+				let To = sserec.To;
+
+				for (let j = From; j <= To; j++)
+					resflags[j - 1] = 1;
+
+				if ((k == 0) && (From > 1)) {
+					// residues with negative residue numbers will not be included
+					for (let j = 1; j < From; j++)
+						resflags[j - 1] = 1;
+				}
+
+				if ((k == substruct.length - 1) && (To < seqLen)) {
+					for (let j = To + 1; j <= seqLen; j++)
+						resflags[j - 1] = 1;
+				}
+
+				// left side
+				if (k > 0) {
+					//SSE_Rec sserec1 = substruct[k - 1];
+					let sserec1 = substruct[k - 1];
+					let To1 = sserec1.To;
+					//let ll = (int) floor(0.5*((let) (From - To1 - 1)));
+					let ll = parseInt(0.5 * (From - To1 - 1));
+
+					if (ll > 0) {
+						for (let j = From - ll; j <= From - 1; j++)
+							resflags[j - 1] = 1;
+					}
+				}
+
+				// right side
+				if (k < substruct.length - 1) {
+					//SSE_Rec sserec1 = substruct[k + 1];
+					let sserec1 = substruct[k + 1];
+					let From1 = sserec1.From;
+					//let ll = (int) ceil(0.5*((let) (From1 - To - 1)));
+					// let ft = From1 - To - 1;
+					// let ll = parseInt(ft/2);
+					// if (ft % 2 == 1) ll++;
+					let ll = parseInt(0.5 * (From1 - To - 1) + 0.5);
+
+					if (ll > 0) {
+						for (let j = To + 1; j <= To + ll; j++)
+							resflags[j - 1] = 1;
+					}
+				}
+			}
+
+			// extract the continguous segments
+			let inseg = false;
+			let startseg;
+			//vector<int> segments;
+			//segments.clear();
+			let segments = [];
+
+			for (let i = 0; i < seqLen; i++) {
+				let rf = resflags[i];
+
+				if (!inseg && (rf == 1)) {
+					// new segment starts here
+					startseg = i + 1;
+					inseg = true;
+					continue;
+				}
+
+				if (inseg && (rf == 0)) {
+					// segment ends
+					segments.push(startseg);
+					segments.push(i);
+					inseg = false;
+				}
+			}
+
+			// check for the last segment
+			if (inseg) {
+				segments.push(startseg);
+				segments.push(seqLen);
+			}
+
+			subdomains.push(segments);
+		}
+
+//console.log("subdomains: " + JSON.stringify(subdomains));
+
+		return subdomains;
+	} // end c2b_NewSplitChain
+}
+
 /**
  * @author Jiyao Wang <wangjiy@ncbi.nlm.nih.gov> / https://github.com/ncbi/icn3d
  */
@@ -29866,6 +30962,8 @@ class Picking {
               transformation.quaternion._y = parseFloat(ic.quaternion._y).toPrecision(5);
               transformation.quaternion._z = parseFloat(ic.quaternion._z).toPrecision(5);
               transformation.quaternion._w = parseFloat(ic.quaternion._w).toPrecision(5);
+
+              /*
               if(ic.bAddCommands) {
                   ic.commands.push('pickatom ' + atom.serial + '|||' + ic.transformCls.getTransformationStr(transformation));
                   ic.optsHistory.push(me.hashUtilsCls.cloneHash(ic.opts));
@@ -29877,6 +30975,9 @@ class Picking {
               if( $( "#" + ic.pre + "logtext" ).length )  {
                 $("#" + ic.pre + "logtext").val("> " + ic.logs.join("\n> ") + "\n> ").scrollTop($("#" + ic.pre + "logtext")[0].scrollHeight);
               }
+              */
+              me.htmlCls.clickMenuCls.setLogCmd('pickatom ' + atom.serial, true);
+
               // update the interaction flag
               ic.bSphereCalc = false;
               //me.htmlCls.clickMenuCls.setLogCmd('set calculate sphere false', true);
@@ -31253,11 +32354,16 @@ class ApplyCommand {
             ic.bGlycansCartoon = false;
         }
       }
+      else if(command.indexOf('save html') == 0) {
+        let  id = command.substr(command.lastIndexOf(' ') + 1);
+        me.htmlCls.eventsCls.saveHtml(id);
+      }
 
     // special, select ==========
 
       else if(command.indexOf('select displayed set') !== -1) {
-        ic.hAtoms = me.hashUtilsCls.cloneHash(ic.dAtoms);
+        //ic.hAtoms = me.hashUtilsCls.cloneHash(ic.dAtoms);
+        ic.hAtoms = me.hashUtilsCls.cloneHash(ic.viewSelectionAtoms);
         ic.hlUpdateCls.updateHlAll();
       }
       else if(command.indexOf('select prop') !== -1) {
@@ -32104,11 +33210,12 @@ class Resid2spec {
 
     atoms2spec(atomHash) {var ic = this.icn3d; ic.icn3dui;
         let  spec = "";
-
         let  i = 0;
         let  structureHash = {}, chainHash = {}, resiHash = {};
+
+        let atom;
         for(let serial in atomHash) {
-            let  atom = ic.atoms[serial];
+            atom = ic.atoms[serial];
             if(i > 0) {
                 spec += ' or ';
             }
@@ -34015,7 +35122,7 @@ class MmcifParser {
 
                                 // not all listed residues are considered missing, e.g., PDB ID 4OR2, only the firts four residues are considered missing
                                 if(!isNaN(resi) &&(prevMissingChain == '' ||(chain != prevMissingChain) ||(chain == prevMissingChain && resi > maxMissingResi)) ) {
-                                    chainMissingResidueArray[chainNum].push(resObject);
+                                    ic.chainMissingResidueArray[chainNum].push(resObject);
 
                                     maxMissingResi = resi;
                                     prevMissingChain = chain;
@@ -35167,7 +36274,7 @@ class ShowAnno {
                   // the missing residues at the end of the seq will be filled up in the API showNewTrack()
                   let nGap = 0;
                   ic.alnChainsSeq[chnid] = [];
-                  let offset =(ic.chainid2offset[chnid]) ? ic.chainid2offset[chnid] : 0;
+                  let offset =(ic.chainid2offset[chnid]) ? ic.chainid2offset[chnid] : 0;                
                   for(let i = 0, il = targetSeq.length; i < il; ++i) {
                       //text += ic.showSeqCls.insertGap(chnid, i, '-', true);
                       if(ic.targetGapHash.hasOwnProperty(i)) {
@@ -35540,7 +36647,8 @@ class AnnoDomain {
         let pdbArray = Object.keys(ic.structures);
         // show 3D domains
         let pdbid = pdbArray[index];
-        let url = me.htmlCls.baseUrl + "mmdb/mmdb_strview.cgi?v=2&program=icn3d&domain&molinfor&uid=" + pdbid;
+        me.htmlCls.baseUrl + "mmdb/mmdb_strview.cgi?v=2&program=icn3d&domain&molinfor&uid=" + pdbid;
+
         if(index == 0 && ic.mmdb_data !== undefined) {
             for(let chnid in ic.protein_chainid) {
                 if(chnid.indexOf(pdbid) !== -1) {
@@ -35556,6 +36664,43 @@ class AnnoDomain {
             }
         }
         else {
+            // calculate 3D domains on-the-fly
+            //ic.protein_chainid[chainArray[i]] 
+            let data = {};
+            data.domains = {};
+            for(let chainid in ic.chains) {
+                let structure = chainid.substr(0, chainid.indexOf('_'));
+                if(pdbid == structure && ic.protein_chainid.hasOwnProperty(chainid)) {
+                    data.domains[chainid] = {};
+                    data.domains[chainid].domains = [];
+                    
+                    let subdomains = ic.domain3dCls.c2b_NewSplitChain(chainid);
+                    for(let i = 0, il = subdomains.length; i < il; ++i) {
+                        // domain item: {"sdid":1722375,"intervals":[[1,104],[269,323]]}
+                        let domain = {};
+                        domain.intervals = [];
+
+                        for(let j = 0, jl = subdomains[i].length; j < jl; j += 2) {
+                            domain.intervals.push([subdomains[i][j], subdomains[i][j+1]]);
+                        }
+
+                        data.domains[chainid].domains.push(domain);
+                    }
+                }
+            }
+
+            ic.mmdb_dataArray[index] = data;
+            let bCalcDirect = true;
+            for(let chnid in ic.protein_chainid) {
+                if(chnid.indexOf(pdbid) !== -1) {
+                    thisClass.showDomainWithData(chnid, ic.mmdb_dataArray[index], bCalcDirect);
+                }
+            }
+
+            ic.bAjax3ddomain = true;
+            ic.bAjaxDoneArray[index] = true;
+
+/*            
             $.ajax({
               url: url,
               dataType: 'json',
@@ -35619,6 +36764,7 @@ class AnnoDomain {
                 return;
               }
             });
+*/            
         }
     }
 
@@ -35637,28 +36783,36 @@ class AnnoDomain {
             this.showDomainPerStructure(i);
         }
     }
-    showDomainWithData(chnid, data) { let ic = this.icn3d, me = ic.icn3dui;
+    showDomainWithData(chnid, data, bCalcDirect) { let ic = this.icn3d, me = ic.icn3dui;
         let html = '<div id="' + ic.pre + chnid + '_domainseq_sequence" class="icn3d-dl_sequence">';
         let html2 = html;
         let html3 = html;
         let domainArray, proteinname;
         let pos = chnid.indexOf('_');
         let chain = chnid.substr(pos + 1);
-        let molinfo = data.moleculeInfor;
-        let currMolid;
-        for(let molid in molinfo) {
-        if(molinfo[molid].chain === chain) {
-          currMolid = molid;
-          proteinname = molinfo[molid].name;
-          break;
-        }
-        }
-        if(currMolid !== undefined && data.domains[currMolid] !== undefined) {
-          domainArray = data.domains[currMolid].domains;
+
+        if(bCalcDirect) {
+            proteinname = chnid;
+            domainArray = data.domains[chnid].domains;
         }
-        if(domainArray === undefined) {
-          domainArray = [];
+        else {
+            let molinfo = data.moleculeInfor;
+            let currMolid;
+            for(let molid in molinfo) {
+                if(molinfo[molid].chain === chain) {
+                currMolid = molid;
+                proteinname = molinfo[molid].name;
+                break;
+                }
+            }
+            if(currMolid !== undefined && data.domains[currMolid] !== undefined) {
+                domainArray = data.domains[currMolid].domains;
+            }
+            if(domainArray === undefined) {
+                domainArray = [];
+            }
         }
+
         for(let index = 0, indexl = domainArray.length; index < indexl; ++index) {
             //var fulltitle = '3D domain ' +(index+1).toString() + ' of ' + proteinname + '(PDB ID: ' + data.pdbId + ')';
             let fulltitle = '3D domain ' +(index+1).toString() + ' of ' + proteinname;
@@ -35681,7 +36835,7 @@ class AnnoDomain {
                 }
 
                 // use the NCBI residue number, and convert to PDB residue number during selection
-                if(ic.bNCBI) {
+                if(ic.bNCBI || bCalcDirect) {
                     fromArray.push(domainFrom);
                     toArray.push(domainTo);
                 }
@@ -37780,6 +38934,18 @@ class Diagram2d {
             ic.bSelectResidue = false;
 */
         });
+
+        $(document).on("click", "#" + ic.pre + "dl_alignerrormap .icn3d-interaction", function(e) { thisClass.icn3d;
+            e.stopImmediatePropagation();
+
+            thisClass.clickInteraction(this);
+        });
+
+        $(document).on("click", "#" + ic.pre + "dl_alignerrormap .icn3d-node", function(e) { thisClass.icn3d;
+            e.stopImmediatePropagation();
+
+            thisClass.clickNode(this);
+        });
     }
 
     clickNode(node) {  let ic = this.icn3d, me = ic.icn3dui;
@@ -38459,8 +39625,11 @@ class Selection {
 
         for(let i in ic.chains) {
            ic.hAtoms = me.hashUtilsCls.unionHash(ic.hAtoms, ic.chains[i]);
-           ic.dAtoms = me.hashUtilsCls.unionHash(ic.dAtoms, ic.chains[i]);
+           //ic.dAtoms = me.hashUtilsCls.unionHash(ic.dAtoms, ic.chains[i]);
         }
+
+        ic.dAtoms = me.hashUtilsCls.cloneHash(ic.hAtoms);
+        ic.viewSelectionAtoms = me.hashUtilsCls.cloneHash(ic.hAtoms);
     }
 
     //Select a chain with the chain id "chainid" in the sequence dialog and save it as a custom selection with the name "commandname".
@@ -38695,6 +39864,7 @@ class Selection {
         if(Object.keys(ic.hAtoms).length == 0) this.selectAll_base();
 
         ic.dAtoms = me.hashUtilsCls.cloneHash(ic.hAtoms);
+        ic.viewSelectionAtoms = me.hashUtilsCls.cloneHash(ic.hAtoms);
 
         let  centerAtomsResults = ic.applyCenterCls.centerAtoms(me.hashUtilsCls.hash2Atoms(ic.dAtoms, ic.atoms));
         ic.maxD = centerAtomsResults.maxD;
@@ -40951,8 +42121,9 @@ class ApplySsbonds {
           let start, end;
 
           if(ic.bAlternate) {
-              start = ic.ALTERNATE_STRUCTURE;
-              end = ic.ALTERNATE_STRUCTURE + 1;
+              let nStructures = structureArray.length;
+              start = ic.ALTERNATE_STRUCTURE % nStructures;
+              end = ic.ALTERNATE_STRUCTURE % nStructures + 1;
           }
           else {
               start = 0;
@@ -44243,7 +45414,8 @@ class ClickMenu {
     //    clkMn2_selectdisplayed: function() {
         me.myEventCls.onIds("#" + me.pre + "mn2_selectdisplayed", "click", function(e) { let ic = me.icn3d;
            thisClass.setLogCmd("select displayed set", true);
-           ic.hAtoms = me.hashUtilsCls.cloneHash(ic.dAtoms);
+           //ic.hAtoms = me.hashUtilsCls.cloneHash(ic.dAtoms);
+           ic.hAtoms = me.hashUtilsCls.cloneHash(ic.viewSelectionAtoms);
            ic.hlUpdateCls.updateHlAll();
            //ic.drawCls.draw();
         });
@@ -45364,6 +46536,10 @@ class ClickMenu {
     //    },
     //    clkMn6_sidebyside: function() {
         me.myEventCls.onIds("#" + me.pre + "mn6_sidebyside", "click", function(e) { let ic = me.icn3d;
+           if(ic.bInputfile) {
+                var aaa = 1; //alert("Side-by-Side does NOT work when the input is from a local file.");
+                return;
+           }
            let url = ic.shareLinkCls.shareLinkUrl(undefined);
            //if(url.indexOf('http') !== 0) {
            //    var aaa = 1; //alert("The url is more than 4000 characters and may not work.");
@@ -45971,6 +47147,8 @@ class SetMenu {
         // show title at the top left corner
         html += me.htmlCls.divStr + "title' class='icn3d-commandTitle' style='font-size:1.2em; font-weight:normal; position:absolute; z-index:1; float:left; display:block; margin: 12px 0px 0px 40px; color:" + titleColor + "; width:" +(me.htmlCls.WIDTH - 40).toString() + "px'></div>";
         html += me.htmlCls.divStr + "viewer' style='position:relative; width:100%; height:100%; background-color: " + me.htmlCls.GREYD + ";'>";
+        // don't show legend in mobile
+        //html += me.htmlCls.divStr + "legend' class='icn3d-text icn3d-legend'></div>";
         html += me.htmlCls.divStr + "mnLogSection'>";
         html += "<div style='height: " + me.htmlCls.MENU_HEIGHT + "px;'></div>";
         html += "</div>";
@@ -48257,6 +49435,12 @@ class SetDialog {
         html += "<div style='width:500px'>";
         html += 'Please specify the mutations with a comma separated mutation list. Each mutation can be specified as "[PDB ID]_[Chain ID]_[Residue Number]_[One Letter Mutatnt Residue]". E.g., the mutation of N501Y in the E chain of PDB 6M0J can be specified as "6M0J_E_501_Y". <br/><br/>';
         html += "<div style='display:inline-block; width:110px'>Mutations: </div>" + me.htmlCls.inputTextStr + "id='" + me.pre + "mutationids' value='6M0J_E_484_K,6M0J_E_501_Y,6M0J_E_417_N' size=50><br/><br/>";
+
+        html += "<b>Data Source</b>: <select id='" + me.pre + "idsource'>";
+        html += "<option value='mmdbid' selected>PDB ID</option>";
+        html += "<option value='afid'>AlphaFold UniProt ID</option>";
+        html += "</select><br/><br/>";
+        
         html += me.htmlCls.buttonStr + "reload_mutation_3d' title='Show the mutations in 3D using the scap program'>3D with scap</button>";
         html += me.htmlCls.buttonStr + "reload_mutation_inter' style='margin-left:20px' title='Show the mutations in 3D and the change of interactions'>Interactions</button>";
         html += me.htmlCls.buttonStr + "reload_mutation_pdb' style='margin-left:20px' title='Show the mutations in 3D and export the PDB of the mutant within 10 angstrom'>PDB</button>";
@@ -48651,7 +49835,7 @@ class SetDialog {
 
         html += me.htmlCls.divStr + "dl_alignerrormap' style='background-color:white' class='" + dialogClass + "'>";
 
-        //html += me.htmlCls.divNowrapStr + "Hold Ctrl key to select multiple nodes." + me.htmlCls.space3 + "</div>";
+        html += me.htmlCls.divNowrapStr + "Hold Ctrl key to select multiple nodes." + me.htmlCls.space3 + "</div>";
       
         me.alignerrormapid = me.pre + 'alignerrormap';
         html += me.htmlCls.divNowrapStr + buttonStrTmp + me.alignerrormapid + '_svg">SVG</button>' + me.htmlCls.space2;
@@ -49122,6 +50306,18 @@ class Events {
        }
     }
 
+    saveHtml(id) { let me = this.icn3dui, ic = me.icn3d;
+        let html = '';
+        html += '<link rel="stylesheet" href="https:///structure.ncbi.nlm.nih.gov/icn3d/lib/jquery-ui-1.12.1.min.css">\n';
+        html += '<link rel="stylesheet" href="https:///structure.ncbi.nlm.nih.gov/icn3d/icn3d_full_ui.css">\n';
+        html += $("#" + id).html();
+        let idArray = id.split('_');
+        let idStr =(idArray.length > 2) ? idArray[2] : id;
+        let structureStr = Object.keys(ic.structures)[0];
+        if(Object.keys(ic.structures).length > 1) structureStr += '-' + Object.keys(ic.structures)[1];
+        ic.saveFileCls.saveFile(structureStr + '-' + idStr + '.html', 'html', encodeURIComponent(html));
+    }
+
     //Hold all functions related to click events.
     allEventFunctions() { let me = this.icn3dui, ic = me.icn3d;
         let thisClass = this;
@@ -49520,26 +50716,29 @@ class Events {
            e.preventDefault();
            if(!me.cfg.notebook) dialog.dialog( "close" );
            let mutationids = $("#" + me.pre + "mutationids").val();
-           let mmdbid = mutationids.substr(0, mutationids.indexOf('_'));
+           let idsource = $("#" + me.pre + "idsource").val();
+           let mmdbid = mutationids.substr(0, mutationids.indexOf('_'));           
            me.htmlCls.clickMenuCls.setLogCmd("3d of mutation " + mutationids, false);
            //window.open(me.htmlCls.baseUrl + 'icn3d/full.html?mmdbid=' + mmdbid + '&command=scap 3d ' + mutationids + '; select displayed set', '_blank');
-           window.open(hostUrl + '?mmdbid=' + mmdbid + '&command=scap 3d ' + mutationids + '; select displayed set', '_blank');
+           window.open(hostUrl + '?' + idsource + '=' + mmdbid + '&command=scap 3d ' + mutationids + '; select displayed set', '_blank');
         });
 
         me.myEventCls.onIds("#" + me.pre + "reload_mutation_pdb", "click", function(e) { me.icn3d;
            e.preventDefault();
            if(!me.cfg.notebook) dialog.dialog( "close" );
            let mutationids = $("#" + me.pre + "mutationids").val();
+           let idsource = $("#" + me.pre + "idsource").val();
            let mmdbid = mutationids.substr(0, mutationids.indexOf('_'));
            me.htmlCls.clickMenuCls.setLogCmd("pdb of mutation " + mutationids, false);
            //window.open(me.htmlCls.baseUrl + 'icn3d/full.html?mmdbid=' + mmdbid + '&command=scap pdb ' + mutationids + '; select displayed set', '_blank');
-           window.open(hostUrl + '?mmdbid=' + mmdbid + '&command=scap pdb ' + mutationids + '; select displayed set', '_blank');
+           window.open(hostUrl + '?' + idsource + '=' + mmdbid + '&command=scap pdb ' + mutationids + '; select displayed set', '_blank');
         });
 
         me.myEventCls.onIds("#" + me.pre + "reload_mutation_inter", "click", function(e) { let ic = me.icn3d;
            e.preventDefault();
            if(!me.cfg.notebook) dialog.dialog( "close" );
            let mutationids = $("#" + me.pre + "mutationids").val();
+           let idsource = $("#" + me.pre + "idsource").val();
 
            let mutationArray = mutationids.split(',');
            let residArray = [];
@@ -49561,7 +50760,7 @@ class Events {
            me.htmlCls.clickMenuCls.setLogCmd("interaction change of mutation " + mutationids, false);
            //window.open(me.htmlCls.baseUrl + 'icn3d/full.html?mmdbid=' + mmdbid + '&command=scap interaction ' + mutationids + '; select ' + selectSpec + ' | name test; line graph interaction pairs | selected non-selected | hbonds,salt bridge,interactions,halogen,pi-cation,pi-stacking | false | threshold 3.8 6 4 3.8 6 5.5; adjust dialog dl_linegraph; select displayed set', '_blank');
            //window.open(me.htmlCls.baseUrl + 'icn3d/full.html?mmdbid=' + mmdbid + '&command=scap interaction ' + mutationids, '_blank');
-           window.open(hostUrl + '?mmdbid=' + mmdbid + '&command=scap interaction ' + mutationids, '_blank');
+           window.open(hostUrl + '?' + idsource + '=' + mmdbid + '&command=scap interaction ' + mutationids, '_blank');
         });
 
     //    },
@@ -50780,18 +51979,12 @@ class Events {
         });
     //    },
     //    clickSaveDialog: function() {
-        $(document).on("click", ".icn3d-saveicon", function(e) { let ic = me.icn3d;
+        $(document).on("click", ".icn3d-saveicon", function(e) { me.icn3d;
            e.stopImmediatePropagation();
            let id = $(this).attr('pid');
-           let html = '';
-           html += '<link rel="stylesheet" href="https:///structure.ncbi.nlm.nih.gov/icn3d/lib/jquery-ui-1.12.1.min.css">\n';
-           html += '<link rel="stylesheet" href="https:///structure.ncbi.nlm.nih.gov/icn3d/icn3d_full_ui.css">\n';
-           html += $("#" + id).html();
-           let idArray = id.split('_');
-           let idStr =(idArray.length > 2) ? idArray[2] : id;
-           let structureStr = Object.keys(ic.structures)[0];
-           if(Object.keys(ic.structures).length > 1) structureStr += '-' + Object.keys(ic.structures)[1];
-           ic.saveFileCls.saveFile(structureStr + '-' + idStr + '.html', 'html', encodeURIComponent(html));
+
+           thisClass.saveHtml(id);
+           me.htmlCls.clickMenuCls.setLogCmd("save html " + id, true);
         });
     //    },
     //    clickHideDialog: function() {
@@ -52306,38 +53499,66 @@ class Alternate {
     // change the display atom when alternating
     //Show structures one by one.
     alternateStructures() { let ic = this.icn3d, me = ic.icn3dui;
-        let hAtomsCount = Object.keys(ic.hAtoms).length;
+        if(!ic.viewSelectionAtoms) {
+            ic.viewSelectionAtoms = me.hashUtilsCls.cloneHash(ic.dAtoms);
+        }
+
+        let viewSelectionAtomsCount = Object.keys(ic.viewSelectionAtoms).length;
         let allAtomsCount = Object.keys(ic.atoms).length;
 
         ic.dAtoms = {};
 
+        // alternate all displayed structures
+        let moleculeArray = Object.keys(ic.structures);
         // only alternate selected structures
-        //let moleculeArray = Object.keys(ic.structures);
-        let structureHash = {};
-        for(let i in ic.hAtoms) {
-            let structure = ic.atoms[i].structure;
-            structureHash[structure] = 1;
-        }
-        let moleculeArray = Object.keys(structureHash);
+        // let structureHash = {};
+        // for(let i in ic.hAtoms) {
+        //     let structure = ic.atoms[i].structure;
+        //     structureHash[structure] = 1;
+        // }
+        // let moleculeArray = Object.keys(structureHash);
 
         for(let i = 0, il = moleculeArray.length; i < il; ++i) {
             let structure = moleculeArray[i];
-            if(i > ic.ALTERNATE_STRUCTURE || (ic.ALTERNATE_STRUCTURE === il - 1 && i === 0) ) {
+            //if(i > ic.ALTERNATE_STRUCTURE || (ic.ALTERNATE_STRUCTURE === il - 1 && i === 0) ) {
+            let bChoose;
+            if(ic.bShift) {
+                // default ic.ALTERNATE_STRUCTURE = -1
+                if(ic.ALTERNATE_STRUCTURE < 0) ic.ALTERNATE_STRUCTURE = 1;
+
+                bChoose = (i == ic.ALTERNATE_STRUCTURE % il - 1) 
+                  || (ic.ALTERNATE_STRUCTURE % il === 0 && i === il - 1);
+            } 
+            else {
+                bChoose = (i == ic.ALTERNATE_STRUCTURE % il + 1) 
+                  || (ic.ALTERNATE_STRUCTURE % il === il - 1 && i === 0);
+            }
+
+            if(bChoose) {
                 for(let k in ic.structures[structure]) {
                     let chain = ic.structures[structure][k];
                     ic.dAtoms = me.hashUtilsCls.unionHash(ic.dAtoms, ic.chains[chain]);
                 }
 
-                ic.ALTERNATE_STRUCTURE = i;
+                //ic.ALTERNATE_STRUCTURE = i;
+                if(ic.bShift) {
+                    --ic.ALTERNATE_STRUCTURE;
+                }
+                else {
+                    ++ic.ALTERNATE_STRUCTURE;
+                }
+
+                if(ic.ALTERNATE_STRUCTURE < 0) ic.ALTERNATE_STRUCTURE += il;
 
                 $("#" + ic.pre + "title").html(structure);
 
                 break;
             }
-        }
+        } 
 
-        if(hAtomsCount < allAtomsCount) {
-            ic.dAtoms = me.hashUtilsCls.intHash(ic.dAtoms, ic.hAtoms);
+        if(viewSelectionAtomsCount < allAtomsCount) {
+            //ic.dAtoms = me.hashUtilsCls.intHash(ic.dAtoms, ic.hAtoms);
+            ic.dAtoms = me.hashUtilsCls.intHash(ic.dAtoms, ic.viewSelectionAtoms);
 
             ic.bShowHighlight = false;
             ic.opts['rotationcenter'] = 'highlight center';
@@ -53067,17 +54288,10 @@ class ContactMap {
         let bContactMap = true;
 
         if(bAfMap) { // cleaned the code by using "use" in SVG, but didn't improve rendering
-            let  factor = 1;
-            let  r = 3 * factor;
-            let  rectSize = 2 * r;
 
             ic.hex2id = {};
             let threshold = 29.0 / max;
             ic.hex2skip = {}; // do not display any error larger than 29 angstrom
-
-            html += "<defs>";
-
-            let linestrokewidth = 1;
             let nRef = 1000;
             for(let i = 0; i < nRef; ++i) {
                 let ratio = 1.0 * i / nRef;
@@ -53087,7 +54301,6 @@ class ContactMap {
                 let gHex = (g.length == 1) ? '0' + g : g;
                 let bHex = rHex;
                 let color = rHex + gHex + bHex;
-                let strokecolor = "#" + color;
 
                 let idRect = me.pre + "afmap_" + i;
 
@@ -53097,11 +54310,11 @@ class ContactMap {
                 }
                 
                 //html += "<g id='" + id + "'>";
-                html += "<rect id='" + idRect + "' x='0' y='0' width='" + rectSize + "' height='" + rectSize + "' fill='" 
-                    + strokecolor + "' stroke-width='" + linestrokewidth + "' stroke='" + strokecolor + "' />";
+//                html += "<rect id='" + idRect + "' x='0' y='0' width='" + rectSize + "' height='" + rectSize + "' fill='" 
+//                    + strokecolor + "' stroke-width='" + linestrokewidth + "' stroke='" + strokecolor + "' />";
                 //html += "</g>"
             }
-            html += "</defs>";
+//            html += "</defs>";
         }
 
         html += ic.lineGraphCls.drawScatterplot_base(nodeArray1, nodeArray2, linkArray, name2node, 0, bContactMap, undefined, undefined, bAfMap);
@@ -54878,7 +56091,7 @@ class Control {
 
             else if(e.keyCode === 65 ) { // A, alternate
                if(Object.keys(ic.structures).length > 1) {
-                   ic.alternateCls.alternateWrapper();
+                 ic.alternateCls.alternateWrapper();
                }
             }
 
@@ -55425,6 +56638,7 @@ class iCn3D {
     this.annoSnpClinVarCls = new AnnoSnpClinVar(this);
     this.annoSsbondCls = new AnnoSsbond(this);
     this.annoTransMemCls = new AnnoTransMem(this);
+    this.domain3dCls = new Domain3d(this);
 
     this.addTrackCls = new AddTrack(this);
     this.annotationCls = new Annotation(this);
@@ -55677,7 +56891,7 @@ class iCn3DUI {
     //even when multiple iCn3D viewers are shown together.
     this.pre = this.cfg.divid + "_";
 
-    this.REVISION = '3.8.1';
+    this.REVISION = '3.9.0';
 
     // In nodejs, iCn3D defines "window = {navigator: {}}"
     this.bNode = (Object.keys(window).length < 2) ? true : false;
@@ -56192,6 +57406,7 @@ exports.Delphi = Delphi;
 exports.DensityCifParser = DensityCifParser;
 exports.Diagram2d = Diagram2d;
 exports.Dialog = Dialog;
+exports.Domain3d = Domain3d;
 exports.Draw = Draw;
 exports.DrawGraph = DrawGraph;
 exports.Dsn6Parser = Dsn6Parser;