Wolfgang,
It queries the boundary conditions as vector of length n_components, just as needed. I do not think so, maybe i be in mistake so if please clarify for me, notice that i do not mean a vector field like velocity vector, but multi-field that each one can be a multi-component vector 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. PS: notice that i do not enquiry about difference of "n_solution_vectors" and "n_components" but logically they should be different, e.g., stimator of density based compressible flow solver. Cheers RT
_______________________________________________
