[ff3d-users] Does function "fuzzySearch()" work correctly

[Robert Li]

Dear Stephane Del Pino,
  when I read a mesh from input file and traced the
source code. I found that cell number is different in
assembleMatrix()and MeshOfTetrahedra::find().
  Fox example, it is 0 in assembleMatrix(), while it
becomes 234 in find().
  I think they should have same cell number, because
find()is called in cell loop of
assembleMatrix()indircetly.
  Thanks!
Robert

>>>output>>>>>>>>>>>>>>>>>>>>>>>>
Using standard Finite Element Method discretization:
Solving ...
reserve
Discretizing the PDE Problem ... Assembling matrix
__currentCellNumber=0
119  204  203  47
guessedCellNumber=234


void assembleMatrix()
  {
    Timer t;

    ffout(2) << "Assembling matrix\n";
    t.start();
  
    switch ((__A).type()) {
    case BaseMatrix::doubleHashedMatrix: {

      DoubleHashedMatrix& A =
dynamic_cast<DoubleHashedMatrix&>(__A);

      for(typename MeshType::const_iterator
icell(__mesh);
          !(icell.end()); ++icell) {
        const ElementGeometry& K = *icell;
                int nn=__mesh.numberOfCells();////????


        size_t
indices[FiniteElement::numberOfDegreesOfFreedom];
        for (size_t l=0;
l<FiniteElement::numberOfDegreesOfFreedom; ++l) {
          indices[l] = __mesh.vertexNumber(K(l));

          printf("%ld  ", indices[l]);
        }
        printf("\n");

        ConformTransformation T(K);
        JacobianTransformation J(T);
        generatesElementaryMatrix(__eSet,J);

        TinyVector<3> massCenter;
        T.value(0.5,0.5,0.5, massCenter);

        for(typename
DiscretizedOperators<ElementaryMatrixType>
              ::iterator iOperator =
__discretizedOperators.begin();
            iOperator != __discretizedOperators.end();
++iOperator) {
          const real_t coef =
((*iOperator).second)(massCenter);

...


2)
MeshOfTetrahedra::const_iterator
MeshOfTetrahedra::find(const double& x,
                       const double& y,
                       const double& z) const
{
  TinyVector<3,real_t> X(x,y,z);
  Octree<size_t>::iterator i =
(*__octree).fuzzySearch(X);

  size_t guessedCellNumber = (*i).value();////????
printf("guessedCellNumber=%ld\n", guessedCellNumber);
...


__________________________________
Do you Yahoo!?
Yahoo! Mail SpamGuard - Read only the mail you want.
http://antispam.yahoo.com/tools