Beispiel #1
0
  private double getDerivative(Sextad sextad, Sigma sigma) {
    Node a = sigma.getA();
    Node b = sigma.getB();

    Node n1 = sextad.getI();
    Node n2 = sextad.getJ();
    Node n3 = sextad.getK();
    Node n4 = sextad.getL();
    Node n5 = sextad.getM();
    Node n6 = sextad.getN();

    if (a == n1) {
      if (b == n4) {
        return r(n2, n5) * r(n3, n6) - r(n2, n6) * r(n3, n5);

      } else if (b == n5) {
        return -r(n2, n4) * r(n3, n6) + r(n3, n4) * r(n2, n6);

      } else if (b == n6) {
        return r(n2, n4) * r(n3, n5) - r(n3, n4) * r(n2, n5);
      }

    } else if (a == n2) {
      if (b == n4) {
        return r(n3, n5) * r(n1, n6) - r(n1, n5) * r(n3, n6);

      } else if (b == n5) {
        return r(n1, n4) * r(n3, n6) - r(n3, n4) * r(n1, n6);

      } else if (b == n6) {
        return -r(n1, n4) * r(n3, n5) + r(n3, n4) * r(n1, n5);
      }

    } else if (a == n3) {
      if (b == n4) {
        return r(n1, n5) * r(n2, n6) - r(n2, n5) * r(n1, n6);

      } else if (b == n5) {
        return -r(n1, n4) * r(n2, n6) + r(n2, n4) * r(n1, n6);

      } else if (b == n6) {
        return r(n1, n4) * r(n2, n5) - r(n2, n4) * r(n1, n5);
      }
    }

    // symmetry
    n6 = sextad.getI();
    n5 = sextad.getJ();
    n4 = sextad.getK();
    n3 = sextad.getL();
    n2 = sextad.getM();
    n1 = sextad.getN();

    if (a == n1) {
      if (b == n4) {
        return r(n2, n5) * r(n3, n6) - r(n2, n6) * r(n3, n5);

      } else if (b == n5) {
        return -r(n2, n4) * r(n3, n6) + r(n3, n4) * r(n2, n6);

      } else if (b == n6) {
        return r(n2, n4) * r(n3, n5) - r(n3, n4) * r(n2, n5);
      }

    } else if (a == n2) {
      if (b == n4) {
        return r(n3, n5) * r(n1, n6) - r(n1, n5) * r(n3, n6);

      } else if (b == n5) {
        return r(n1, n4) * r(n3, n6) - r(n3, n4) * r(n1, n6);

      } else if (b == n6) {
        return -r(n1, n4) * r(n3, n5) + r(n3, n4) * r(n1, n5);
      }

    } else if (a == n3) {
      if (b == n4) {
        return r(n1, n5) * r(n2, n6) - r(n2, n5) * r(n1, n6);

      } else if (b == n5) {
        return -r(n1, n4) * r(n2, n6) + r(n2, n4) * r(n1, n6);

      } else if (b == n6) {
        return r(n1, n4) * r(n2, n5) - r(n2, n4) * r(n1, n5);
      }
    }

    return 0.0;
  }