// Path: src/com/rcwang/seal/util/StringFactory.java
// public static class SID implements Iterable<String> {
//     
//   String[] array;
//     
//   public SID(List<String> strs) {
//     array = strs.toArray(new String[strs.size()]);
//   }
// 
//   public SID(SID id) {
//     this(id.array);
//   }
// 
//   public SID(String[] array) {
//     this.array = array.clone();
//   }
//     
//   public void append(SID sid) {
//     append(sid.array);
//   }
//     
//   public void append(String s) {
//     append(StringFactory.toID(s));
//   }
//     
//   public void append(String[] array) {
//     String[] array2 = new String[this.array.length+array.length];
//     System.arraycopy(this.array, 0, array2, 0, this.array.length);
//     System.arraycopy(array, 0, array2, this.array.length, array.length);
//     this.array = array2;
//   }
//     
//   @Override
//   public boolean equals(Object obj) {
//     if (this == obj) return true;
//     if (obj == null) return false;
//     if (getClass() != obj.getClass()) return false;
//     final SID other = (SID) obj;
//     if (!Arrays.equals(array, other.array))
//       return false;
//     return true;
//   }
//     
//   @Override
//   public int hashCode() {
//     return 31 + Arrays.hashCode(array);
//   }
//     
//   public Iterator<String> iterator() {
//     return Arrays.asList(array).iterator();
//   }
//     
//   public int length() {
//     int length = 0;
//     for (String s : this)
//       length += s.length();
//     return length;
//   }
// 
//   @Override
//   public String toString() {
//     return StringFactory.toName(this);
//   }
//     
//   public EntityLiteral toLiteral() { return new EntityLiteral(StringFactory.toName(this)); }
//     
//   public SID toLower() {
//   	return StringFactory.toID(new EntityLiteral(StringFactory.toName(this).toLowerCase()));
//   }
// }


import java.util.HashSet;
import java.util.Set;
import java.util.Map.Entry;
import org.apache.log4j.Logger;
import com.rcwang.seal.util.StringFactory.SID;

/**************************************************************************
 * Developed by Language Technologies Institute, Carnegie Mellon University
 * Written by Richard Wang (rcwang#cs,cmu,edu)
 **************************************************************************/
package com.rcwang.seal.util;




public class MatrixUtil {
  
  public static Logger log = Logger.getLogger(MatrixUtil.class);
  
  public static SparseMatrix dotProduct(SparseMatrix sm1, SparseMatrix sm2) {
    SparseMatrix sm = new SparseMatrix();
    for (SparseVector col2 : sm2.getColumns())
      sm.addColumn(dotProduct(sm1, col2));
    return sm;
  }
  
  public static SparseVector dotProduct(SparseMatrix sm, SparseVector sv1) {
    SparseVector sv = new SparseVector(sv1.id);
    for (SparseVector row : sm.getRows())
      sv.put(row.id, dotProduct(row, sv1));
    return sv;
  }
  
  public static double dotProduct(SparseVector sv1, SparseVector sv2) {
    if (sv2.size() < sv1.size()) {
      SparseVector sv = sv2;
      sv2 = sv1; sv1 = sv;
    }
    double sum = 0;

// Path: src/com/rcwang/seal/util/StringFactory.java
// public static class SID implements Iterable<String> {
//     
//   String[] array;
//     
//   public SID(List<String> strs) {
//     array = strs.toArray(new String[strs.size()]);
//   }
// 
//   public SID(SID id) {
//     this(id.array);
//   }
// 
//   public SID(String[] array) {
//     this.array = array.clone();
//   }
//     
//   public void append(SID sid) {
//     append(sid.array);
//   }
//     
//   public void append(String s) {
//     append(StringFactory.toID(s));
//   }
//     
//   public void append(String[] array) {
//     String[] array2 = new String[this.array.length+array.length];
//     System.arraycopy(this.array, 0, array2, 0, this.array.length);
//     System.arraycopy(array, 0, array2, this.array.length, array.length);
//     this.array = array2;
//   }
//     
//   @Override
//   public boolean equals(Object obj) {
//     if (this == obj) return true;
//     if (obj == null) return false;
//     if (getClass() != obj.getClass()) return false;
//     final SID other = (SID) obj;
//     if (!Arrays.equals(array, other.array))
//       return false;
//     return true;
//   }
//     
//   @Override
//   public int hashCode() {
//     return 31 + Arrays.hashCode(array);
//   }
//     
//   public Iterator<String> iterator() {
//     return Arrays.asList(array).iterator();
//   }
//     
//   public int length() {
//     int length = 0;
//     for (String s : this)
//       length += s.length();
//     return length;
//   }
// 
//   @Override
//   public String toString() {
//     return StringFactory.toName(this);
//   }
//     
//   public EntityLiteral toLiteral() { return new EntityLiteral(StringFactory.toName(this)); }
//     
//   public SID toLower() {
//   	return StringFactory.toID(new EntityLiteral(StringFactory.toName(this).toLowerCase()));
//   }
// }


// Path: src/com/rcwang/seal/util/MatrixUtil.java
import java.util.HashSet;
import java.util.Set;
import java.util.Map.Entry;
import org.apache.log4j.Logger;
import com.rcwang.seal.util.StringFactory.SID;

/**************************************************************************
 * Developed by Language Technologies Institute, Carnegie Mellon University
 * Written by Richard Wang (rcwang#cs,cmu,edu)
 **************************************************************************/
package com.rcwang.seal.util;




public class MatrixUtil {
  
  public static Logger log = Logger.getLogger(MatrixUtil.class);
  
  public static SparseMatrix dotProduct(SparseMatrix sm1, SparseMatrix sm2) {
    SparseMatrix sm = new SparseMatrix();
    for (SparseVector col2 : sm2.getColumns())
      sm.addColumn(dotProduct(sm1, col2));
    return sm;
  }
  
  public static SparseVector dotProduct(SparseMatrix sm, SparseVector sv1) {
    SparseVector sv = new SparseVector(sv1.id);
    for (SparseVector row : sm.getRows())
      sv.put(row.id, dotProduct(row, sv1));
    return sv;
  }
  
  public static double dotProduct(SparseVector sv1, SparseVector sv2) {
    if (sv2.size() < sv1.size()) {
      SparseVector sv = sv2;
      sv2 = sv1; sv1 = sv;
    }
    double sum = 0;

    for (Entry<SID, Cell> entry : sv1) {
