interpolate_complex.edp

/* Author: Chloros2 <chloros2@gmx.de>
 * Created: 2018-05-19
 *
 * Copyright (C) 2018 Chloros2 <chloros2@gmx.de>
 *
 * This program is free software: you can redistribute it and/or modify it
 * under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hopeC that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see
 * <https://www.gnu.org/licenses/>.
 *
 */

include "ffmatlib.idp"

mesh Th = square(10, 10, [2*x*pi,2*y*pi]);

fespace Vh(Th, P1);
Vh<complex>  uC=(x+1i*y)*(x+1i*y);

savemesh(Th,"complex_test.msh");
ffSaveVh(Th,Vh,"complex_vh.txt");
ffSaveData(uC,"complex_data.txt");