dolfin team mailing list archive

[ndl@xxxxxxxxxxxxxx]

Hi.
I'm having some "dimension" troubles when applying some DirichletBC's.

I've something like this:
-------------------------------
MyFunctionSpace V(mesh); //Th=P1+P1
SubSpace V0(V,0); //S
SubSpace V1(V,1); //P
...
DirichletBoundary dirichletboundary;
DBC bc(t);
DirichletBC dbc(V1,bc,dirichletboundary);
.....
MyBilinearForm a(V1,V1);
Function p(V1);
MyLinearForm L(V1);
VariationalProblem pde(a,L,dbc);
pde.solve(p);
----------------------------------------------
The V dimension is 24136 and V1 is half:12068.
PETSc for instance gives me:
[0]PETSC ERROR: Argument out of range!
[0]PETSC ERROR: row 12068 out of range!
...........
No problems applying DirichletBC to the V systems.
After digging the mailing list I found a reference to a bug:
"2008-12-09 [DOLFIN-dev] [HG DOLFIN] Bug fix for dimension check in
DirichletBC."
Could this be related?
Or am I doing something wrong when defining the sub-problem?

Thanks in advance.

-- 
Nuno David Lopes

e-mail:ndl@xxxxxxxxxxxxxx        (FCUL/CMAF)
           nlopes@xxxxxxxxxxxxxxx    (ISEL)
http://ptmat.ptmat.fc.ul.pt/%7Endl/

Sat Mar 14 01:26:04 WET 2009