Mercurial > fem-fenics-eugenio

comparison src/assemble.cc @ 168:67944f307560

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove the verbose debug output
author gedeone-octave <marcovass89@hotmail.it>
date Sun, 06 Oct 2013 22:31:48 +0100
parents 3895c8d4c066
children 8ea37cfc7a14
comparison
equal deleted inserted replaced
167:3895c8d4c066 168:67944f307560
72 const dolfin::Form & a = frm.get_form (); 72 const dolfin::Form & a = frm.get_form ();
73 a.check (); 73 a.check ();
74 74
75 if (a.rank () == 2) 75 if (a.rank () == 2)
76 { 76 {
77 std::cout << "Assembling Matrix from the bilinear form..."
78 << std::endl;
79 dolfin::parameters["linear_algebra_backend"] = "uBLAS"; 77 dolfin::parameters["linear_algebra_backend"] = "uBLAS";
80 dolfin::Matrix A; 78 dolfin::Matrix A;
81 dolfin::assemble (A, a); 79 dolfin::assemble (A, a);
82 80
83 for (std::size_t i = 1; i < nargin; ++i) 81 for (std::size_t i = 1; i < nargin; ++i)
137 retval(0) = sm; 135 retval(0) = sm;
138 } 136 }
139 137
140 else if (a.rank () == 1) 138 else if (a.rank () == 1)
141 { 139 {
142 std::cout << "Assembling Vector from the linear form..."
143 << std::endl;
144 dolfin::Vector A; 140 dolfin::Vector A;
145 dolfin::assemble (A, a); 141 dolfin::assemble (A, a);
146 142
147 for (std::size_t i = 1; i < nargin; ++i) 143 for (std::size_t i = 1; i < nargin; ++i)
148 { 144 {
173 retval(0) = myb; 169 retval(0) = myb;
174 } 170 }
175 171
176 else if (a.rank () == 0) 172 else if (a.rank () == 0)
177 { 173 {
178 std::cout << "Assembling double from the functional form..."
179 << std::endl;
180 double b = dolfin::assemble (a); 174 double b = dolfin::assemble (a);
181 retval(0) = octave_value (b); 175 retval(0) = octave_value (b);
182 } 176 }
183 177
184 else 178 else
187 } 181 }
188 } 182 }
189 } 183 }
190 else if (nargout == 2) 184 else if (nargout == 2)
191 { 185 {
192 std::cout << "Assemble: apply boundary condition to a vector for a nonlinear problem..."
193 << std::endl;
194
195 if (nargin < 2) 186 if (nargin < 2)
196 print_usage (); 187 print_usage ();
197 else 188 else
198 { 189 {
199 if (args(0).type_id () == form::static_type_id ()) 190 if (args(0).type_id () == form::static_type_id ())
206 const dolfin::Form & a = frm.get_form (); 197 const dolfin::Form & a = frm.get_form ();
207 a.check (); 198 a.check ();
208 199
209 if (a.rank () == 1) 200 if (a.rank () == 1)
210 { 201 {
211 std::cout << "Assembling Vector from the linear form..."
212 << std::endl;
213 dolfin::Vector A; 202 dolfin::Vector A;
214 dolfin::assemble (A, a); 203 dolfin::assemble (A, a);
215 204
216 dolfin::Vector x (myx.length ()); 205 dolfin::Vector x (myx.length ());
217 for (std::size_t i = 0; i < myx.length (); ++i) 206 for (std::size_t i = 0; i < myx.length (); ++i)