451 void gLevelsetPoints::RbfOp(int p, int index, const fullMatrix<double> &cntrs,  in RbfOp()  argument
455   fullMatrix<double> rbfMatB = generateRbfMat(p, index, cntrs, nodes);  in RbfOp()
459     rbfInvA = generateRbfMat(0, index, cntrs, cntrs);  in RbfOp()
466   D.resize(nodes.size1(), cntrs.size1());  in RbfOp()
471                                  const fullMatrix<double> &cntrs,  in evalRbfDer()  argument
479   RbfOp(p, index, cntrs, nodes, D, isLocal);  in evalRbfDer()
483 void gLevelsetPoints::setup_level_set(const fullMatrix<double> &cntrs,  in setup_level_set()  argument
487   int numNodes = cntrs.size1();  in setup_level_set()
501     cntrsPlus(i, 0) = cntrs(i, 0);  in setup_level_set()
502     cntrsPlus(i, 1) = cntrs(i, 1);  in setup_level_set()
503     cntrsPlus(i, 2) = cntrs(i, 2);  in setup_level_set()
506         sqrt((cntrs(i, 0) - cntrs(j, 0)) * (cntrs(i, 0) - cntrs(j, 0)) +  in setup_level_set()
507              (cntrs(i, 1) - cntrs(j, 1)) * (cntrs(i, 1) - cntrs(j, 1)) +  in setup_level_set()
508              (cntrs(i, 2) - cntrs(j, 2)) * (cntrs(i, 2) - cntrs(j, 2)));  in setup_level_set()
514   cntrsPlus(numNodes, 0) = cntrs(0, 0) + dist_max;  in setup_level_set()
515   cntrsPlus(numNodes, 1) = cntrs(0, 1) + dist_max;  in setup_level_set()
516   cntrsPlus(numNodes, 2) = cntrs(0, 2) + dist_max;  in setup_level_set()
519   evalRbfDer(1, 1, cntrsPlus, cntrs, ONES, sx, true);  in setup_level_set()
520   evalRbfDer(2, 1, cntrsPlus, cntrs, ONES, sy, true);  in setup_level_set()
521   evalRbfDer(3, 1, cntrsPlus, cntrs, ONES, sz, true);  in setup_level_set()
535       level_set_nodes(i, j) = cntrs(i, j);  in setup_level_set()
536       level_set_nodes(i + numNodes, j) = cntrs(i, j) - delta * norms(i, j);  in setup_level_set()
537       level_set_nodes(i + 2 * numNodes, j) = cntrs(i, j) + delta * norms(i, j);  in setup_level_set()