(Closed) project(a*b, v) get weird result also


129
views
0
5 months ago by
Hi:

from dolfin import *

mesh=UnitSquareMesh(5,5)

V=FunctionSpace(mesh, "Lagrange",2)

a_=Function(V)
b_=Function(V)

a_.vector()[:]=1
a_.vector()[25]=-1
a_.vector().array()
array([ 1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1., -1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.])

b_.vector()[:]=a_.vector().array()>=0
b_.vector().array()
array([ 1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  0.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,  1.,
        1.,  1.,  1.,  1.])

c_=project(a_*b_,V)
c_.vector().array()
array([ 1.00031668,  0.99999569,  1.00047384,  1.0022144 ,  1.00093042,
        0.99992455,  0.99884712,  1.00406464,  1.00255549,  1.00448146,
        1.00218528,  0.99884577,  0.99992672,  1.00002007,  1.03493756,
        0.98977631,  1.05153577,  1.05715593,  1.01121546,  0.98977235,
        1.00405916,  1.00031645,  1.0002511 ,  0.99998577,  1.00013188,
       -1.09687929,  1.05716495,  1.05717414,  1.05716445,  1.00087187,
        1.03490742,  0.99989058,  1.00449171,  1.00393634,  1.05150069,
        0.98980301,  1.00012144,  1.00003169,  1.00013686,  0.99998191,
        1.03492624,  1.05150318,  0.98980991,  1.00445898,  1.0011249 ,
        0.99977202,  0.99907742,  1.0012775 ,  1.03488204,  0.9897977 ,
        0.98979908,  1.00451004,  1.0109353 ,  1.00392932,  0.99907808,
        0.99997348,  1.00010122,  1.0000652 ,  1.00402368,  1.00121613,
        0.99982798,  1.00029666,  0.99997029,  0.99977047,  1.00013193,
        1.01100618,  1.0039191 ,  0.9990719 ,  0.99913846,  1.00109692,
        1.00082801,  1.00003396,  0.99989165,  1.00001037,  0.99988216,
        1.00018995,  1.00007131,  0.99998961,  1.0000099 ,  1.00212773,
        1.00131193,  0.99953989,  1.00006429,  1.00000665,  1.00006605,
        0.99999772,  1.00114253,  1.00001591,  0.99989392,  1.0000245 ,
        1.00016197,  1.00006992,  0.999977  ,  1.00001227,  1.00000298,
        1.00000412,  1.00056137,  1.00010307,  0.99994833,  1.00003603,
        1.00009259,  0.99999777,  0.99999755,  1.00000733,  1.00002342,
        0.99999958,  1.00000855,  1.00000105,  1.00000174,  1.0000673 ,
        0.99999921,  0.99999961,  1.0000141 ,  1.00000509,  0.99999899,
        1.00000418,  1.0000056 ,  0.9999997 ,  1.00000033,  1.00000069,
        1.00000143])​


This is similar to the post https://www.allanswered.com/post/nqzmq/project-absfunction-get-weird-result/

Any advise would be appreciated!

Thanks

Community: FEniCS Project

This post is closed as a duplicate. Similar post can be found at:

written 5 months ago by Michal Habera  
Please login to add an answer/comment or follow this question.
The thread is closed. No new answer/comment may be added.

Similar posts:
Search »