missing operator== for block diagonal LDLT

Author: peddie

include/albatross/src/linalg/block_diagonal.hpp

@@ -68,6 +68,8 @@ struct BlockDiagonalLDLT {
   Eigen::Index rows() const;
 
   Eigen::Index cols() const;
+
+  bool operator==(const BlockDiagonalLDLT &other) const;
 };
 
 struct BlockDiagonal {
@@ -248,6 +250,11 @@ inline const BlockDiagonalLDLT &BlockDiagonalLDLT::adjoint() const {
   return *this;
 }
 
+inline bool
+BlockDiagonalLDLT::operator==(const BlockDiagonalLDLT &other) const {
+  return blocks == other.blocks;
+}
+
 /*
  * Block Diagonal
  */