Exemplo n.º 1
0
 @Override
 public double[] applyAll(VariableArray valAry, Map<Object, Object> cache) {
   int len = valAry.length();
   double[] r = valAry.get("r");
   double[] s = valAry.get("s");
   double[] t = valAry.get("t");
   double[] rlt = new double[len];
   double cr = vt[funIndex][0];
   double cs = vt[funIndex][1];
   double ct = vt[funIndex][2];
   for (int i = 0; i < len; i++)
     rlt[i] = (cr * r[i] + 1.0) * (cs * s[i] + 1.0) * (ct * t[i] + 1.0) / 8.0;
   return rlt;
 }
Exemplo n.º 2
0
 @Override
 public double[] applyAll(VariableArray valAry, Map<Object, Object> cache) {
   int len = valAry.length();
   double[] r = valAry.get("r");
   double[] s = valAry.get("s");
   double[] t = valAry.get("t");
   double[] rlt = new double[len];
   double cr = vt[funIndex][0];
   double cs = vt[funIndex][1];
   double ct = vt[funIndex][2];
   if ("r".equals(var))
     for (int i = 0; i < len; i++) rlt[i] = cr * (cs * s[i] + 1.0) * (ct * t[i] + 1.0) / 8.0;
   else if ("s".equals(var))
     for (int i = 0; i < len; i++) rlt[i] = (cr * r[i] + 1.0) * cs * (ct * t[i] + 1.0) / 8.0;
   else if ("t".equals(var))
     for (int i = 0; i < len; i++) rlt[i] = (cr * r[i] + 1.0) * (cs * s[i] + 1.0) * ct / 8.0;
   else throw new FutureyeException();
   return rlt;
 }
Exemplo n.º 3
0
    @Override
    public double[] applyAll(VariableArray valAry, Map<Object, Object> cache) {
      double[] detJ = null;
      double[][][] J = null;
      int len = valAry.length();
      double[] rlt = new double[len];
      if (cache != null) {
        detJ = (double[]) cache.get(1);
        J = (double[][][]) cache.get(2);
      }
      if (detJ == null || J == null) {
        J = new double[3][3][len];
        J[0][0] = x_r.applyAll(valAry, cache);
        J[0][1] = x_s.applyAll(valAry, cache);
        J[0][2] = x_t.applyAll(valAry, cache);
        J[1][0] = y_r.applyAll(valAry, cache);
        J[1][1] = y_s.applyAll(valAry, cache);
        J[1][2] = y_t.applyAll(valAry, cache);
        J[2][0] = z_r.applyAll(valAry, cache);
        J[2][1] = z_s.applyAll(valAry, cache);
        J[2][2] = z_t.applyAll(valAry, cache);
        // @see ./doc/invA33.png
        detJ = Utils.determinant(J);
        if (cache != null) {
          cache.put(1, detJ);
          cache.put(2, J);
        }
      }

      if ("r".equals(rst)) {
        if ("x".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[1][1][i] * J[2][2][i] - J[1][2][i] * J[2][1][i]) / detJ[i];
        else if ("y".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[0][2][i] * J[2][1][i] - J[0][1][i] * J[2][2][i]) / detJ[i];
        else if ("z".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[0][1][i] * J[1][2][i] - J[0][2][i] * J[1][1][i]) / detJ[i];
      } else if ("s".equals(rst)) {
        if ("x".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[1][2][i] * J[2][0][i] - J[1][0][i] * J[2][2][i]) / detJ[i];
        else if ("y".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[0][0][i] * J[2][2][i] - J[0][2][i] * J[2][0][i]) / detJ[i];
        else if ("z".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[0][2][i] * J[1][0][i] - J[0][0][i] * J[1][2][i]) / detJ[i];
      } else if ("t".equals(rst)) {
        if ("x".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[1][0][i] * J[2][1][i] - J[1][1][i] * J[2][0][i]) / detJ[i];
        else if ("y".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[0][1][i] * J[2][0][i] - J[0][0][i] * J[2][1][i]) / detJ[i];
        else if ("z".equals(xyz))
          for (int i = 0; i < len; i++)
            rlt[i] = (J[0][0][i] * J[1][1][i] - J[0][1][i] * J[1][0][i]) / detJ[i];
      } else {
        throw new FutureyeException();
      }
      return rlt;
    }