-- FreeFem++ v4.6 (Thu Apr  2 15:47:38 CEST 2020 - git v4.6)
 Load: lg_fem lg_mesh lg_mesh3 eigenvalue 
    1 : //  Discussion:
    2 : //
    3 : //    Solve the Poisson equation over the L-shaped region.
    4 : //
    5 : //    Export the data, using ffmatlib, for graphic processing by MATLAB or Octave.
    6 : //
    7 : //  Licensing:
    8 : //
    9 : //    Copyright (C) 2018 Chloros2 <chloros2@gmx.de>
   10 : //
   11 : //    This program is free software: you can redistribute it and/or modify it
   12 : //    under the terms of the GNU General Public License as published by
   13 : //    the Free Software Foundation, either version 3 of the License, or
   14 : //    (at your option) any later version.
   15 : //
   16 : //    This program is distributed in the hopeC that it will be useful, but
   17 : //    WITHOUT ANY WARRANTY; without even the implied warranty of
   18 : //    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   19 : //    GNU General Public License for more details.
   20 : //
   21 : //    You should have received a copy of the GNU General Public License
   22 : //    along with this program.  If not, see
   23 : //    <https://www.gnu.org/licenses/>.
   24 : //
   25 : //  Modified:
   26 : //
   27 : //    2018-05-19
   28 : //
   29 : // Author: 
   30 : //
   31 : //   Chloros2 <chloros2@gmx.de>
   32 : //
   33 : include "ffmatlib.idp"/*ffmatlib.idp Export Macros (Plot Solutions for Matlab and Octave)
    2 :  *
    3 :  * Author: Chloros2 <chloros2@gmx.de>
    4 :  * Created: 2018-12-21
    5 :  *
    6 :  *   This file is part of the ffmatlib which is hosted at
    7 :  *   https://github.com/samplemaker/freefem_matlab_octave_plot
    8 :  *
    9 :  * Copyright (C) 2018 Chloros2 <chloros2@gmx.de>
   10 :  *
   11 :  * This program is free software: you can redistribute it and/or modify it
   12 :  * under the terms of the GNU General Public License as published by
   13 :  * the Free Software Foundation, either version 3 of the License, or
   14 :  * (at your option) any later version.
   15 :  *
   16 :  * This program is distributed in the hope that it will be useful, but
   17 :  * WITHOUT ANY WARRANTY; without even the implied warranty of
   18 :  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   19 :  * GNU General Public License for more details.
   20 :  *
   21 :  * You should have received a copy of the GNU General Public License
   22 :  * along with this program.  If not, see
   23 :  * <https://www.gnu.org/licenses/>.
   24 :  *
   25 :  */
   26 : 
   27 : //Save finite element connectivity (j-th degree of freedom for i-th element)
   28 : macro ffSaveVh(Th, Vh, filename){
   29 #   {
   30 #   ofstream file(filename);
   31 #   file.precision(10);
   32 #   for (int i=0; i<Th.nt; i++){
   33 #     for (int j=0; j<Vh.ndofK; j++){
   34 #       file << Vh(i,j) << "\n";
   35 #     }
   36 #   }
   37 #   }
   38 # }  )  //EOM
   39 : 
   40 : //Save FE Space data
   41 : macro ffSaveData(u, filename){
   42 #   {
   43 #   ofstream file(filename);
   44 #   file.precision(10);
   45 #   for (int j=0; j<u[].n; j++){
   46 #     file << u[][j] << endl;
   47 #   }
   48 #   }
   49 # }  )  //EOM
   50 : 
   51 : macro ffSaveData2(u1, u2, filename){
   52 #   {
   53 #   ofstream file(filename);
   54 #   int datalen=u1[].n;
   55 #   if (u2[].n!=datalen){
   56 #     cout << "error: arguments must have same size" << endl;
   57 #     exit(1);
   58 #   }
   59 #   file.precision(10);
   60 #   for (int j=0; j<datalen; j++){
   61 #     file << u1[][j] << " " << u2[][j] << endl;
   62 #   }
   63 #   }
   64 # }  )  //EOM
   65 : 
   66 : macro ffSaveData3(u1, u2, u3, filename){
   67 #   {
   68 #   ofstream file(filename);
   69 #   int datalen=u1[].n;
   70 #   if ((u2[].n!=datalen) | (u3[].n!=datalen)){
   71 #     cout << "error: arguments must have same size" << endl;
   72 #     exit(1);
   73 #   }
   74 #   file.precision(10);
   75 #   for (int j=0; j<datalen; j++){
   76 #     file << u1[][j] << " " << u2[][j] << " " << u3[][j] << endl;
   77 #   }
   78 #   }
   79 # }  )  //EOM
   80 : 
   81 : macro ffSaveData4(u1, u2, u3, u4, filename){
   82 #   {
   83 #   ofstream file(filename);
   84 #   int datalen=u1[].n;
   85 #   if ((u2[].n!=datalen) | (u3[].n!=datalen) | (u4[].n!=datalen)){
   86 #     cout << "error: arguments must have same size" << endl;
   87 #     exit(1);
   88 #   }
   89 #   file.precision(10);
   90 #   for (int j=0; j<datalen; j++){
   91 #     file << u1[][j] << " " << u2[][j] << " " << u3[][j] << " " << u4[][j] << endl;
   92 #   }
   93 #   }
   94 # }  )  //EOM
   95 : 
   34 : 
   35 : cout << endl;
   36 : cout << "poisson_ell:" << endl;
   37 : cout << "  FreeFem++ version." << endl;
   38 : cout << "  Solve the Poisson equation over the L-shaped re
  ... : gion." << endl;
   39 : cout << "  Using the ffmatlib interface, export the data f
  ... : or graphic" << endl;
   40 : cout << "  processing by MATLAB/Octave." << endl;
   41 : 
   42 : 
   43 : verbosity=0;
   44 : border a(t=0,1){x=t;y=0;label=1;};
   45 : border b(t=0,0.5){x=1;y=t;label=1;};
   46 : border c(t=0,0.5){x=1-t;y=0.5;label=1;};
   47 : border d(t=0.5,1){x=0.5;y=t;label=1;};
   48 : border e(t=0.5,1){x=1-t;y=1;label=1;};
   49 : border f(t=0,1){x=0;y=1-t;label=1;};
   50 : 
   51 : int nn=40;
   52 : mesh Th = buildmesh ( a(nn) + b(.5*nn) + c(.5*nn) +d(.5*nn) + e(.5*nn) + f(nn));
   53 : 
   54 : fespace Vh(Th,P2);
   55 : //fespace Vh(Th,P1);
   56 : 
   57 : Vh u,v;
   58 : 
   59 : func rho=1.0;
   60 : 
   61 : solve Poisson(u,v,solver=LU)=int2d(Th)(dx(u)*dx(v) + dy(u)*dy(v))
   62 :                            - int2d(Th)(rho*v)+on(a,b,c,d,e,f,u=0);
   63 : 
   64 : plot(Th,u);
   65 : 
   66 : cout << endl;
   67 : cout << "NbBoundaryElements:	" <<  Th.nbe << endl;
   68 : cout << "NbTriangles:		" <<  Th.nt << endl;
   69 : cout << "NbVertices:		" <<  Th.nv << endl;
   70 : cout << "nDoF:			" << Vh.ndof << endl;
   71 : cout << endl;
   72 : 
   73 : real lambda = 1.0;
   74 : Vh qx, qy;
   75 : qx = -lambda*dx(u);
   76 : qy = -lambda*dy(u);
   77 : //
   78 : //  Write out data for later processing by MATLAB/Octave.
   79 : //
   80 : savemesh ( Th, "poisson_ell.msh" );
   81 : ffSaveVh ( Th, Vh, "poisson_ell_vh.txt" )
   29 @   
   30 @    
   31 @   
   32 @       
   33 @         
   34 @           
   35 @     
   36 @   
   37 @   
   38 @  {
   29 @   {
   30 @   ofstream file( "poisson_ell_vh.txt");
   31 @   file.precision(10);
   32 @   for (int i=0; i< Th.nt; i++){
   33 @     for (int j=0; j< Vh.ndofK; j++){
   34 @       file <<  Vh(i,j) << "\n";
   35 @     }
   36 @   }
   37 @   }
   38 @ } ;
   82 : ffSaveData3 ( u, qx, qy, "poisson_ell_data.txt" )
   67 @   
   68 @    
   69 @    
   70 @      
   71 @         
   72 @     
   73 @   
   74 @   
   75 @       
   76 @                 
   77 @   
   78 @   
   79 @  {
   67 @   {
   68 @   ofstream file( "poisson_ell_data.txt");
   69 @   int datalen= u[].n;
   70 @   if (( qx[].n!=datalen) | ( qy[].n!=datalen)){
   71 @     cout << "error: arguments must have same size" << endl;
   72 @     exit(1);
   73 @   }
   74 @   file.precision(10);
   75 @   for (int j=0; j<datalen; j++){
   76 @     file <<  u[][j] << " " <<  qx[][j] << " " <<  qy[][j] << endl;
   77 @   }
   78 @   }
   79 @ } ;
   83 : //
   84 : //  Terminate.
   85 : //
   86 : cout << "\n";
   87 : cout << "poisson_ell:\n";
   88 : cout << "  Normal end of execution.\n";
   89 : 
   90 : exit ( 0 );
   91 : 
   92 :  sizestack + 1024 =2952  ( 1928 )


poisson_ell:
  FreeFem++ version.
  Solve the Poisson equation over the L-shaped region.
  Using the ffmatlib interface, export the data for graphic
  processing by MATLAB/Octave.

NbBoundaryElements:	160
NbTriangles:		2696
NbVertices:		1429
nDoF:			5553


poisson_ell:
  Normal end of execution.
  current line = 90
exit(0)
 err code 0 ,  mpirank 0