FEniCS Project

How to change vector values of a Function

Asked by Xujun Zhao on 2013-03-12

a beginner's problem for C++ implementation:

I created two functions in the same function space, say
Function f1(V), f2(V)

We know there are two main components for a function, vector and functionspace. Here f1 and f2's vector should have the same length. Now I want to change the vector values of f2, for example, f2.vector() = 0.5*(f1.vector() + f2.vector()). what operations do I need to do?

Thank you for your answers.
XZ

Question information

Language:: English Edit question

Status:: Answered

For:: FEniCS Project Edit question

Assignee:: No assignee Edit question

Last query:: 2013-03-15

Last reply:: 2013-03-16

Link existing bug

Revision history for this message

Xujun Zhao (xzhao99) said on 2013-03-15:

OK, I check the source codes and found the solution. Although I am not sure if this is optimal, it does work! Hope it helps to other ones who want to use the same functions.

using vector->get_local() to pick out each components of function vectors
do any operations: + - * / to get an intermediate variable w
using vector->set() to change the vector value into w
finally vector-apply() to finalize the assemble (I am not sure this is necessary, because I also tried not to use "apply" and found that the vector values were changed. Anyone can give me an explanation? )