Hello, When computing an integral in parallel, I know one needs to scale things and do a global reduction. There are examples in previous posts (for example here ) with the integral of 2 vectors, say: Vh<complex> scaledVecu = Vecu; for[i, v : A.D] scaledVecu[] *= v; complex wrong = int2d(Th)(scaled…

This is in fact independent of the number of components. Keep in mind that in FreeFEM, whenever you do u[], it selects all degrees of freedom. So whether u is a function defined on a fespace Vh(Th, [Pk, Qk, Tk, Zk, Wk]) or fespace Vh(Th, Pk), the syntax from your snippet will indeed scale all values…

Thank you for your reply. Sorry, my question was not clear: I was not talking of the number of components (for instance [ux] vs [ux,uy] vs [ux,uy,uz] etc.), but the number of terms that are multiplied in the integrand: int2d(Th)(u*v) vs int2d(Th)(u*v*w).

Integral in parallel

General Discussion

prj April 7, 2022, 7:07am 4

I see. Only one of the term has to be scaled.

Poisson Equation (Parallel Implementation)

Topic		Replies	Views
Compute matrix scalar product without creating new variable General Discussion	0	337	March 4, 2021
mpiAllReduce(uReduce[], uG[], mpiCommWorld, mpiSUM) General Discussion	3	345	June 11, 2021
Integral over nested meshes General Discussion	13	512	October 17, 2022
Memory issue with Parallel computing General Discussion	17	815	July 1, 2021
Computing 3D integral with mixed/vector element variables is very slow General Discussion	0	414	April 20, 2019

Integral in parallel

Related topics