> e.g., in "error_estimator.cc" line 1322 we consider single component case
> but can be multi-field, vector "g" stor one scaler (not vector) per QP
> (quadrature ponit), when we take bc, we just take one bc for all field, and
> as you can see in line 1331, we loop over "n_solution_vectors" and exert
>
> per_thread_data.phi[n][point][0] -= g[point]
>
> actually it should be something like:
>
> per_thread_data.phi[n][point][0] -= g[n][point];
>
> i.e., do not contrast between various "n"
>
> in the same manner we have this scenario for multi-component vector fields
> around line 1333.
Ah, I see. Yes, we use the same boundary conditions for all solution vectors.
I wouldn't call it a bug (the code doesn't do anything different from what
the documentation says) but we could indeed improve the documentation and say
explicitly that this is the case.
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [EMAIL PROTECTED]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
