{{{
Index: src/jacobian.cpp
===================================================================
--- src/jacobian.cpp (revision 29354)
+++ src/jacobian.cpp (working copy)
@@ -30,7 +30,8 @@
}

Jacobian::Jacobian(const Jacobian& arg):
- size(arg.columns())
+ size(arg.columns()),
+ nr_blocks(arg.nr_blocks)
{
twists = new Twist[size*nr_blocks];
for(unsigned int i=0;i @@ -45,8 +46,17 @@
twists[i]=arg.twists[i];
return *this;
}
-

+ bool Jacobian::operator == (const Jacobian& rhs)
+ {
+ bool rc = true; // prove otherwise
+ for(unsigned int i=0;i + {
+ rc &= (twists[i] == rhs.twists[i]);
+ }
+ return rc;
+ }
+
Jacobian::~Jacobian()
{
delete [] twists;
Index: src/jacobian.hpp
===================================================================
--- src/jacobian.hpp (revision 29354)
+++ src/jacobian.hpp (working copy)
@@ -41,6 +41,7 @@
Jacobian(const Jacobian& arg);

Jacobian& operator=(const Jacobian& arg);
+ bool operator == (const Jacobian& rhs);

~Jacobian();
}}}