[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[ff3d-users] Does function "fuzzySearch()" work correctly
From: |
Robert Li |
Subject: |
[ff3d-users] Does function "fuzzySearch()" work correctly |
Date: |
Thu, 19 Feb 2004 09:45:41 -0800 (PST) |
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
- [ff3d-users] Does function "fuzzySearch()" work correctly,
Robert Li <=