# ScalarField division

 Register Blogs Members List Search Today's Posts Mark Forums Read

 August 26, 2007, 17:57 I do the following operation b #1 Senior Member   Maka Mohu Join Date: Mar 2009 Posts: 305 Rep Power: 10 I do the following operation but it gives wrong result (1.22005e-312): const scalarField& ry = patch().deltaCoeffs(); // magnitude of patch-normal distance from face to cell centers: const scalarField& yn= 1.0/ry; but if I make a loop and use: const scalarField& yn= 1.0/ry[facei]; it gives correct result. How can I do mathematical operations that contain scalarFields and scalars without looping? Thanks. Best regards, Maka

 August 26, 2007, 18:16 Classic mistake: reference to #2 Senior Member   Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,783 Rep Power: 22 Classic mistake: reference to a temporary: if you do: scalarField yn= 1.0/ry; all will be well. Operation 1/ry will return a COPY, i.e. a new field. You took a reference to it, which is a memory error - that's the little & bit after the scalarField. I bet you already know this - if not, time for a careful look at a good C++ book. Got it? Hrv __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

 August 27, 2007, 05:10 Thanks Hrv for the explanation #3 Senior Member   Maka Mohu Join Date: Mar 2009 Posts: 305 Rep Power: 10 Thanks Hrv for the explanation. Best regards, Maka.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post xiao OpenFOAM Running, Solving & CFD 9 July 22, 2010 03:23 maka OpenFOAM Bugs 9 February 19, 2009 10:43 vitke OpenFOAM Running, Solving & CFD 5 September 1, 2008 05:35 wese OpenFOAM Running, Solving & CFD 0 March 14, 2008 08:50 mbeaudoin OpenFOAM Meshing & Mesh Conversion 10 February 22, 2007 08:43

All times are GMT -4. The time now is 17:40.