1NOTES on building phasta.: 2the usual command line is 3 4ARCHOS needs to be specified when ever the build target is different from the native 5ARCHOS of the build machine. 6 7gmake [ VERSION=xyz ] [ COMPRESSIBLE=1 ] [ VERS=opt ] 8PARALLEL=[ mpich2, mpigm,mpivmi,sunmpi,sgimpi,decmpi,ibmmpi] 9 10VERSION needs to be specified if the source base is in the dirtory structure 11phSolver/xyz/phSolver 12 13PARALLEL can take any of the values specified above based on the target architechture 14and mpi build one would like to use. 15 16if PARALLEL=mpich2 is specified, we also need to specify a MPICOMM type ch_shmem or 17ch_p4. the comm type also has a -icc extension for the intel compiler based builds on 18linix. 19 20at the moment phasta expectes the libraries phastaIO, shapeFunction and metis to have a 21-$(PARALLEL) extension on them...simple symbolic links should suffice. 22 23 249/17/2007 KJ 25After severe neglect we are back to version control. We start with a rollback to 1.9_7b and then bring AMG into the main trunk. 26EXPORT environment variable AMG to include (AMG=1) or exclude (AMG=0) AMG specific code. Note if you include you have to either compile your own PARPACK libraries or link to those compiled by someone else (sunc2 or kjansen). They are not put into SVN because they the only changes that have been made are those contained within the download instructions. Note there is also now a directory under phSolver called AMG. If you are looking to build AMG code you will need to execute the Makefile in that directory between the common build and the incompressible build. 27 28 29 305/3/2006 KJ 31Fixed several places that compressible code was broken in the 3 previous commits. 32Modified the Makefile so that it gets the MKL libraries properly under their current placement. This is a temporary fix that hopefully will move into the buildUtil environment soon. 33 34 354/20/2006 EB 36two files were added: symline.c and shptri.c 37 381/16/2006 EB TAGGED as v1_9_8 39phasta works with 2D meshes for scalar solve (advection diffusion) with specified initial velocity. 40 41July 28 2005 TAGGED AS v1_9_7b 42Stanford deformable wall and generalized outflow boundary conditions merged onto top of the trunk. Ditched the use of common.h by .c functions (instead used extern of phasta_iotype. Also adopted Nathan's changes to partition.cc. 43 44CVS: ---------------------------------------------------------------------- 45CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 46CVS: 47CVS: Committing in . 48CVS: 49CVS: Modified Files: 50CVS: CHANGES common/asbnabi.f common/bctint.f common/common.h 51CVS: common/common_c.h common/fillsparse.f common/gendat.f 52CVS: common/genini.f common/input.config common/input.f 53CVS: common/input_fform.cc common/local_mass.f 54CVS: common/new_interface.c common/partition.cc common/phasta.cc 55CVS: common/pointer.f common/proces.f common/pvsqbi.f common/qpbc.f 56CVS: common/readnblk.f common/rwvelb.f common/solvecon.f 57CVS: common/timeseries.f incompressible/Makefile 58CVS: incompressible/advLES.f incompressible/asbflx.f 59CVS: incompressible/asbmfg.f incompressible/asigmr.f 60CVS: incompressible/bc3lhs.f incompressible/bc3res.f 61CVS: incompressible/bflux.f incompressible/e3b.f 62CVS: incompressible/e3bvar.f incompressible/e3ivar.f 63CVS: incompressible/elmgmr.f incompressible/forces.f 64CVS: incompressible/itrPC.f incompressible/itrdrv.f 65CVS: incompressible/lesSparse.f incompressible/solfar.f 66CVS: incompressible/stats.f incompressible/usr.c 67CVS: ---------------------------------------------------------------------- 68 69======= 70 71June 27 2005 JM/OS 72v1_9_7 -tagged commit 73introduced ybar: time averaged fields: 74ybar(:,1) - ybar(:,3) is average velocity components 75ybar(:,4) is average pressure 76ybar(:,5) is average speed 77turned on upon setting: 78"Print ybar: True" 79dumped into restart files with tag "ybar" 80required to perform avergage speed hessian based aniso adaptivity 81 82 83 84May 18 2005 JM 85small but essential change in genbkb.f 86reads now 'double' instead of 'integer' as the nbc values ! 87 88 89February 3 2005 90v1_9_6 -tagged commit 91communication fixed 92 93wall shear stress computation implemented 94wall shear stress print out into restart (via tag 'wall shear stresses' ) is 95turned on by setting 'Print Wall Fluxes' to 'True' 96 97ifdef cplusplus added to common_c.h to allow new_interface.c to use it 98normals fixed so that bflux performs properly on curved surfaces (f3lhs.f) 99GETINTPNTS function specialization for intel platforms removed (getIntPnts.c) 100iotype for arrays can now be binary or ascii (set by user in solver.inp). 101error data arrays are written into restart now 102syntax cleanup on lessKey function for C++ standards (partition.cc) 103debugger behaviour made more uniform across platforms (phasta.cc) 104phasta became a function in common that is called by main in either incompressible or compressible directory (phasta.cc) 105fixed a rho scaling problem (e3lhs.f and e3stab.f) 106forces_s properly handled for parallel use (forces.f) 107 108CVS: Modified Files: 109CVS: phSolver//common/f3lhs.f phSolver//incompressible/bflux.f 110 111 112 113v1_9_5 September 24th 2003 114Makefiles changed to so that executables are always relinked. 115newshape.cc indented for readability 116 117September 12th 2003 118v1_9_4 119 120make file changes to avoid automatic override of ARCHOS 121removed unnecessary defines. 122removed all references to MPI to buildUtil, all parallel builds must 123specify a value of PARALLEL and potentially an MPICOMM(ch_shmem/ch_p4) 124 125PARALLEL can take values ( mpich2, mpigm, ibmmpi, mpivmi, summpi, 126 sgimpi, decmpi ) 127mpivmi is for ia64_linux( mpich built on VMI ) 128 129MPICOMM needs to be specified only for PARALLEL=mpich2 130 131MPINATIVE is nolonger a valid build option, one needs to use 132PARALLEL={arch}mpi 133 134September 11th 2003 135v1_9_3 136 137Unstructured SPEBC is add to phasta 138CVS: Modified Files: 139CVS: common/common.h common/common_c.h common/gendat.f 140CVS: common/genini.f common/genscale.f common/getvel.f 141CVS: common/input.config common/input_fform.cc common/proces.f 142CVS: common/readnblk.f common/renum.f common/spebc.f 143CVS: compressible/bc3lhs.f compressible/bc3res.f 144CVS: compressible/itrPC.f compressible/itrdrv.f 145CVS: Added Files: 146CVS: common/elem-search.f common/eqn_plane.f common/get_coeff.f 147CVS: ---------------------------------------------------------------------- 148 149July 18th 2003 150v1_9_0 151 152Branch for windows returns to the tree. This version also uses SCOREC 153makefiles and as such has a makefile in each of the three 154subdirectories (common makes a library that is used by incompressible 155or compressible). Many files have changed but almost all in a superficial way 156 157 1581)removal of stdcall for intel which was a digital fortran need (now using Intel compiler for the fortran). 159 1602) SCOREC makefiles wanted cc instead of cpp 161 1623) mpistub.f90 is gone so now you have to get mpich for windows 163 164There are bug fixes for surface tension relative to the earlier version as 165well and restart.latest.<nproc> should again work.. 166 167 168CVS: Modified Files: 169CVS: Makefile common/auxmpi.h common/commu.f common/ctypes.f 170CVS: common/gltri.c common/new_interface.c common/proces.f 171CVS: common/qpbc.f common/readnblk.f common/shp6w.c common/shphex.c 172CVS: common/shppyr.c common/shptet.c common/symhex.c 173CVS: common/sympyr.c common/symquad.c common/symquadw.c 174CVS: common/symtet.c common/symtri.c common/symtripyr.c 175CVS: common/symwdg.c common/tmrc.c compressible/get_h.c 176CVS: compressible/gltet.c incompressible/getSol.c 177CVS: incompressible/lestools.c incompressible/usr.c 178CVS: incompressible/usr.h incompressible/write_hessian.c 179CVS: Added Files: 180CVS: common/Cinput.cc common/Makefile common/ValType.cc 181CVS: common/flush.c common/input_fform.cc common/missing_win32.c 182CVS: common/newshape.cc common/partition.cc common/phasta.cc 183CVS: compressible/Makefile incompressible/Makefile 184CVS: compressible/stats.f 185CVS: Removed Files: 186CVS: common/Cinput.cpp common/ValType.cpp common/input_fform.cpp 187CVS: common/mpistub.f90 common/newshape.cpp common/partition.cpp 188CVS: common/phasta.cpp incompressible/mpistub.f90 189 190June 13th 2003 191v1_8_0 192 193File: Makefile Status: Modified 194File: auxmpi.h Status: Modified 195File: commu.f Status: Modified 196File: ctypes.f Status: Modified 197File: genbkb.f Status: Modified 198File: genblk.f Status: Modified 199File: input.config Status: Modified 200File: input_fform.cpp Status: Modified 201File: phasta.cpp Status: Modified 202File: readnblk.f Status: Modified 203File: bflux.f Status: Modified 204 205Needs NSpre 3.6.2 or higher. 206 2071) NSpre_Parallel now integrated with phSolver 208 209 1/ now phSolver can start from a 1 proc output and partition the problem 210 2/ it still retains the capability to run on already partitioned problems 211 3/ mesh dual info now in phasta I/O format 212 4/ dual.dat ceases to exist, dual is now written to geombc.dat.1 213 2142) SONFATH is now specified in solver.inp 215 2163) I/O format specified in solver.inp 217 2184) A new MPI data type added to avoid commu error for error indicators. 219 kdof=14 corresponds to n = 10 220 221------------------------------------------------------------------------------- 222 223May 22, 2003 224v1_7_2 225File: CHANGES Status: Modified 226File: Makefile Status: Modified 227File: itrdrv.f Status: Modified 228File: les.h Status: Modified 229File: usr.c Status: Modified 230File: usr.h Status: Modified 231 2321) implemented restarting of les Solves, the aperm vector is now saved to 233 the restart file and again read in before starting a solve. 234 2352) the myfles routines now call the C routines inside insted of fortran 236 routines. 237 238------------------------------------------------------------------------ 239 240May 22, 2003 241v1_7_1 242CVS: 243CVS: Modified Files: 244CVS: CHANGES common/common.h common/common_c.h common/input.config 245CVS: common/input_fform.cpp common/qpbc.f common/readnblk.f 246CVS: incompressible/bflux.f incompressible/e3b.f 247CVS: incompressible/e3ivar.f incompressible/e3lhs.f 248CVS: incompressible/e3q.f incompressible/e3stab.f 249CVS: incompressible/elmgmr.f incompressible/itrdrv.f 250CVS: incompressible/solfar.f 251CVS: ---------------------------------------------------------------------- 252 2531) (IC) bflux now uses phastaIO and flux array is written into restart file 2542) (IC) bflux correctly accounts for periodicity 2553) (IC) bugs in discontinuous field reconstruction (second derivative of velocity and normal) were resolved which should fix surface tension bug 2564) (IC) fixed lumped mass scalar bug 2575) ascii I/O is now supported with iotype variable set in solver.inp (almost) 2586) somewhat improved messages in readnblk.f with regard to non-fatal reading errors. 2597) (IC) fixed a bug in the lhs for the scalar solver. It was not accounting for the phi_{,t} term in the stabilization tangent. 2608) (IC) scalar can be run as Galerkin. After 7) we have verified that scalar advection diffusion converges to requested solver tolerance (as a linear problem should) 261 262 263March 31, 2003 264Fixed a bug in compressible/itrdrv.f: the second argument to the final call of 265 Bflux was ac. To get restart behavior consistent with continuous running, we 266 now pass acold instead. 267 268March 10, 2003 269Fixed a bug in common/settauw.f: iLES.ge.0 was being used as a condition upon 270 which we use LES quantities for near-wall modeling. Sadly, these quantities 271 are only present when iLES.gt.0, so this is the condition we adopt instead. 272 273 274 275v1_7_0 276 277* new phastaIO libarary introduced 278* linux debugging support included 279 280Modified Files: 281-------------- 282 Makefile 283 common/new_interface.c 284 common/phasta.cpp 285 common/readnblk.f 286 incompressible/itrdrv.f 287 288Removed Files: 289------------- 290 common/genbcb.f 291 292 293-------------------------------------------------------------------------- 294CVS: ---------------------------------------------------------------------- 295CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 296CVS: 297CVS: Committing in . 298CVS: 299CVS: Modified Files: 300CVS: CHANGES common/bctint.f common/blkdta.f common/common.h 301CVS: common/common_c.h common/input.config common/input_fform.cpp 302CVS: common/proces.f compressible/e3massl.f compressible/e3source.f 303CVS: compressible/itrdrv.f common/new_interface.c 304CVS: ---------------------------------------------------------------------- 305 306 3071) both Analytical and IC sponge enabled 308 changes to input.config include 309# Body Force Option: Cooling Analytic # ibody=4 => matflag(5,n) 310# Body Force Option: Cooling Initial Condition # ibody=5 => matflag(5,n) 311 3122) Prandtl number (and thereby conductivity) is set in input.config (or solver.inp) 313 Prandtl Number: 0.72 # Compressible code sets diffusivity with this 314 315 3163) Ability to set sponge on one equation at a time 317 Sponge for Continuity Equation : False 318 Sponge for x Momentum Equation : False 319 Sponge for y Momentum Equation : False 320 Sponge for z Momentum Equation : False 321 Sponge for Energy Equation : False 322} 3234) Upper limit set on sponge coefficient 324 325 Maximum Value of Sponge Parameter: 1.0 326 3275) Changes to io interface ANIL DESCRIBE ( Anil: no changes, memory leaks fixed , code 328 indented properly ) 329 330 link and unlink used instead of ln -s for making restart.latest.* 331 332June 21 2002 v1_6_3 333CVS: ---------------------------------------------------------------------- 334CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 335CVS: 336CVS: Committing in . 337CVS: 338CVS: Modified Files: 339CVS: CHANGES Makefile Test_C/solver.inp common/bardmc.f 340CVS: common/common.h common/common_c.h common/genbkb.f 341CVS: common/genblk.f common/genscale.f common/genshp.f 342CVS: common/input.config common/input_fform.cpp common/newshape.cpp 343CVS: common/proces.f common/readnblk.f compressible/asigmr.f 344CVS: compressible/bc3res.f compressible/e3.f compressible/e3conv.f 345CVS: compressible/e3juel.f compressible/e3ls.f 346CVS: compressible/e3wmlt.f compressible/elmgmr.f 347CVS: compressible/elmmfg.f compressible/itrdrv.f 348CVS: compressible/localt.f compressible/solgmr.f 349CVS: compressible/solmfg.f incompressible/asbflx.f 350CVS: incompressible/asigmr.f incompressible/e3.f 351CVS: incompressible/e3b.f incompressible/e3dc.f 352CVS: incompressible/e3ivar.f incompressible/e3res.f 353CVS: incompressible/elmgmr.f incompressible/itrdrv.f 354CVS: ---------------------------------------------------------------------- 355 356Cleaned up a number of bugs and made it DEC compliant (mostly harmless "bugs") 357 358Modified the way scalar parameters are handled: 359 a) solver tolerance 360 b) number of solves per LHS formation 361 c) limiting 362 363In the incompressible code, temperature is also afforded the same type of descriptors 364 365The easiest way to understand the above two changes is to compare input.config of this version to an earliear version. 366 367To add the capability to change the number of solves per LHS formation to scalars required rewiring the use of iprec and ilhs back into the scalar routines. This was "cleaned up" throughout the code with iprec and ilhs flags being set at the itrdrv.f level with the exception of matrix free where it changes. 368 369Returned the option to use pressure outflow boundary condition for the compressible code in the "entropy variable W" form. Default is to simply zero res(:,1) but if you set 370 Entropy Form of Pressure Constraint on Weight Space: 1 # 1 turns it on 371you will get what the name implies (see bc3res.f). 372 373Fixed a bug that prevented you from solving the second scalar without solveing the first. We were trying too hard to save memory when allocating the LHS matrix. 374 375Some of the changes above have "taken over" arrays previously used for iteration sequences. As this capability has long been depricated, in this version forward it is officially pronounced DEAD. 376 377May 25 2002 tagged as v1_6_3Almost 378CVS: ---------------------------------------------------------------------- 379CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 380CVS: 381CVS: Committing in . 382CVS: 383CVS: Modified Files: 384CVS: CHANGES common/input_fform.cpp compressible/getdiff.f 385CVS: compressible/itrbc.f 386 387Fixed substantial bug in itrbcSclr.f that was keeping axisymmetric RANS from converging 388 389Fixed solver tolerance error problem that was present when the incompressible code solved temperature equation. 390 391Fixed e3tau_nd.f bug with DC 392 393May 8 2002 tagged as v1_6_2 394CVS: ---------------------------------------------------------------------- 395CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 396CVS: 397CVS: Committing in . 398CVS: 399CVS: Modified Files: 400CVS: Test_C/solver.inp common/asbwmod.f common/aveprep.f 401CVS: common/bardmc.f common/bctint.f common/common.h 402CVS: common/common_c.h common/filtprep.f common/genblk.f 403CVS: common/gendat.f common/genini.f common/gensvb.f common/getdmc.f 404CVS: common/getstrl.f common/input.config common/input.f 405CVS: common/input_fform.cpp common/pointer.f common/proces.f 406CVS: common/readnblk.f common/rwvelb.f common/settauw.f 407CVS: compressible/asiflx.f compressible/asigmr.f 408CVS: compressible/asimfg.f compressible/asiq.f 409CVS: compressible/asires.f compressible/bflux.f compressible/e3.f 410CVS: compressible/e3b.f compressible/e3bvar.f compressible/e3ivar.f 411CVS: compressible/e3q.f compressible/e3ql.f compressible/e3qvar.f 412CVS: compressible/e3source.f compressible/e3tau.f 413CVS: compressible/e3visc.f compressible/elmgmr.f 414CVS: compressible/elmmfg.f compressible/getdiff.f 415CVS: compressible/itrdrv.f compressible/itrres.f 416CVS: compressible/solgmr.f incompressible/advLES.f 417CVS: incompressible/asbflx.f incompressible/asbmfg.f 418CVS: incompressible/asigmr.f incompressible/asiq.f 419CVS: incompressible/bflux.f incompressible/e3.f 420CVS: incompressible/e3b.f incompressible/e3bvar.f 421CVS: incompressible/e3ivar.f incompressible/e3q.f 422CVS: incompressible/e3ql.f incompressible/e3res.f 423CVS: incompressible/e3source.f incompressible/e3stab.f 424CVS: incompressible/elmgmr.f incompressible/filters.f 425CVS: incompressible/itrdrv.f incompressible/rstatic.f 426CVS: Added Files: 427CVS: incompressible/getdiff.f 428CVS: Removed Files: 429CVS: compressible/scatnu.f incompressible/e3visc.f 430CVS: incompressible/scatnu.f 431CVS: ---------------------------------------------------------------------- 432 433The above routines: 434 4351) added auto-sponge wherein a flow can be forced to remain close to 436 its initial condition by a cooling body force where the target flow 437 state is set to be the initial condition 438 4392) First cut at a DES (Detatched Eddy Simulation) model 440 added. Basically the code solves both RANS (Spalart-Allmaras) and any 441 of the LES models and uses a user specified edge length as a switch to 442 govern which local eddy viscosity to use. This necessitated: 443 444 iturb replaced by iRANS and iLES which are both set when DES is the 445 selected model. 446 eLES the edge length set by "DES Edge Length" in solver.inp 447 blend of viscosities in getdiff 448 4493) First steps towards a common getDiff (e3visc.f removed from incompressible 450 and renamed getDiff 451 4524) xmudml is now xmudmi (since it is at quadrature points not local nodes) 453 and is now correctly dimensioned throughout code (along with strl) 454 4555) scatnu made its way to common 456 4576) gensvb had a bug that distributed piecewise constant boundary flux 458 data (what is provided by NSpre) to all nshlb shapefunction 459 coefficients (which means that quadratic and higher modes had this 460 value). Now it is only distributed to vertices. 461 462CVS: ---------------------------------------------------------------------- 463CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 464CVS: 465CVS: Committing in . 466CVS: 467CVS: Modified Files: 468CVS: Makefile Test_C/solver.inp common/common.h common/common_c.h 469CVS: common/input.config common/input_fform.cpp common/localy.f 470CVS: common/proces.f compressible/e3bdg_nd.f compressible/e3ivar.f 471CVS: compressible/itrdrv.f incompressible/asbflx.f 472CVS: incompressible/asbmfg.f incompressible/asigmr.f 473CVS: incompressible/asiq.f incompressible/e3.f incompressible/e3b.f 474CVS: incompressible/e3bvar.f incompressible/e3ivar.f 475CVS: incompressible/e3q.f incompressible/e3res.f 476CVS: incompressible/e3source.f incompressible/e3stab.f 477CVS: incompressible/e3visc.f incompressible/itrPC.f 478CVS: incompressible/itrdrv.f 479CVS: Added Files: 480CVS: common/elm3keps.f common/turbke.f 481CVS: ---------------------------------------------------------------------- 482 483The above routines changed to incorporate the k-epsilon model into the incompressible code. Other minor fixes to a misplaced tdecay in input_fform.cpp. 484Note also that scalar solver tolerances must now be individually set (default for all is 1.0e-3) in both the compressible and incompressible code (see input.config). Incompressible scalars can also now be solved using the conservation restoring velocity. In this way if one can work with the advective form and still have conservative scalars (as is the case with flow). 485 486 487 488 4894/23/02 490CVS: ---------------------------------------------------------------------- 491CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 492CVS: 493CVS: Committing in . 494CVS: 495CVS: Modified Files: 496CVS: phSolver/CHANGES phSolver/Makefile phSolver/Test_C/solver.inp 497CVS: phSolver/common/bctint.f phSolver/common/common.h 498CVS: phSolver/common/common_c.h phSolver/common/input.config 499CVS: phSolver/common/input_fform.cpp 500CVS: phSolver/common/new_interface.c phSolver/common/rwvelb.f 501CVS: phSolver/compressible/asigmr.f phSolver/compressible/asires.f 502CVS: phSolver/compressible/au1mfg.f phSolver/compressible/au2mfg.f 503CVS: phSolver/compressible/e3.f phSolver/compressible/e3ls.f 504CVS: phSolver/compressible/e3source.f phSolver/compressible/e3tau.f 505CVS: phSolver/compressible/e3visc.f phSolver/compressible/itrdrv.f 506CVS: phSolver/compressible/itrfdi.f phSolver/compressible/itrres.f 507CVS: phSolver/compressible/solmfg.f 508CVS: phSolver/incompressible/asigmr.f 509CVS: phSolver/incompressible/e3ivar.f 510CVS: phSolver/incompressible/itrdrv.f 511CVS: Added Files: 512CVS: phSolver/common/timedata.f phSolver/compressible/e3bdg_nd.f 513CVS: phSolver/compressible/e3eig1.f phSolver/compressible/e3eig1.f~ 514CVS: Removed Files: 515CVS: phSolver/compressible/e3cool.f 516CVS: ---------------------------------------------------------------------- 517 518Re-instated matrix tau and extended it to higher order cubics 519Added timeseries through arbitrary points including higher order 520eliminated e3cool.f (this technique is in e3source now) 521 522Changed restart naming convention. Now a soft link is created and updated for the last restart (restart.lastest.[1..nproc]). This can be nice for script chained jobs at ncsa sinc tarred solution files can go one way while the restart needed to restart the calculation is put elsewhere (without having to know its name). To disable this feature simply remove restart.latest.* and old behavior of looking for numstart.dat and reading the restart that it points to is recovered. 523 524To facilitate transitions from RANS to LES a flag was created to turn off source terms to scalar (turbulence model) equations. 525A second flag was created to force the decay of the scalar fields according to 526 527s_{n+1}=df*s{n} 528 529where df is the decay factor. 530 531Source term for a harmonic oscillator has been added. 532 533Compressible code now also uses minIters and maxIters like ACUSIM 534solver. WARNING this may change regression test values unless you set 535minIters to be 1 as it was before. WARNING II note that the word ACUSIM has been removed from the key phrases 536 537 Minimum Number of Iterations per Nonlinear Iteration: 5 # minIters 538 Maximum Number of Iterations per Nonlinear Iteration: 50 # maxIter 539 540Rational for change is that you typically spent quite bit of time forming the LHS, if you have a loose tolerance that is satified quickly you might as well solve it a bit better after investing the CPU time in forming a LHS. This argument is diminished if you are able to reform LHS infrequently. Currently, only minIters is used in the compressible code but consolidation effort underway to address duplicity of this with nGMRES and Kspace. 541 5424/22/02 543CVS: ---------------------------------------------------------------------- 544CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 545CVS: 546CVS: Committing in . 547CVS: 548CVS: Modified Files: 549CVS: CHANGES common/asithf.f common/common.h common/common_c.h 550CVS: common/filtprep.f common/genshp.f common/getdmc.f 551CVS: common/getstrl.f common/input.config common/input_fform.cpp 552CVS: common/proces.f compressible/asbflx.f compressible/asbmfg.f 553CVS: compressible/asbres.f compressible/asigmr.f compressible/e3.f 554CVS: compressible/e3dc.f compressible/e3ivar.f compressible/e3q.f 555CVS: compressible/e3source.f compressible/getdiff.f 556CVS: compressible/itrbc.f incompressible/advLES.f 557CVS: incompressible/asbmfg.f incompressible/asigmr.f 558CVS: incompressible/asiq.f incompressible/e3.f incompressible/e3b.f 559CVS: incompressible/e3dc.f incompressible/e3ivar.f 560CVS: incompressible/e3lhs.f incompressible/e3q.f 561CVS: incompressible/e3ql.f incompressible/e3qvar.f 562CVS: incompressible/e3res.f incompressible/e3source.f 563CVS: incompressible/e3stab.f incompressible/e3visc.f 564CVS: incompressible/elmgmr.f incompressible/itrdrv.f 565CVS: incompressible/solfar.f 566CVS: Added Files: 567CVS: common/solvecon.f 568CVS: Removed Files: 569CVS: compressible/solvecon.f 570CVS: ---------------------------------------------------------------------- 571 572Significant memory savings for dynamic model 573 introduced maxnint to common (it tells the maximum numberof integration 574 points FOR A SPECIFIC RUN which allows savings of memory) 575split turbvar block to avoid strange behavior 576ivart was still being used, now iconvflow is used correctly in its place 577Major cleanup of two phase flow work including addition of surface tension to incompressible flow 578 579 580bugfixes 581shapefunctions mis-dimensioned in e3q 582MAXSH=32 NOT 125 583several asb* and e3q routines had benign misdimensioning problems (benign because array was not used). 584bugs in DC lhs 585 5864/18/02 587CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 588CVS: 589CVS: Committing in . 590CVS: 591CVS: Modified Files: 592CVS: CHANGES Makefile common/asadj.f common/bardmc.f 593CVS: common/bctint.f common/fillsparse.f common/genadj.f 594CVS: common/genbc.f common/gendat.f common/genini.f 595CVS: common/getIntPnts.c common/getdmc.f common/proces.f 596CVS: common/settauw.f common/symtripyr.c common/symwdg.c 597CVS: compressible/bc3lhs.f compressible/e3.f compressible/elmgmr.f 598CVS: compressible/get_h.c compressible/gltet.c 599CVS: compressible/itrdrv.f compressible/solvecon.f 600CVS: incompressible/advLES.f incompressible/bc3lhs.f 601CVS: incompressible/e3res.f incompressible/elmgmr.f 602CVS: incompressible/forces.f incompressible/itrdrv.f 603CVS: Added Files: 604CVS: common/qpbc.f common/rerun_check.f 605CVS: Removed Files: 606CVS: compressible/qpbc.f incompressible/qpbc.f 607CVS: ---------------------------------------------------------------------- 608Bug fixes 6091) locblkb bug in settauw.f 6102) cubic quadrature rule fixed 6113) qpbc moved to common 6124) removed the inline absolute values on pointers to ien (this confuses people at the very least and may have been a performance hit in fortran90) 613 614KAPL REPORTED CHANGES 615 K1) bardmc.f and getdmc.f nsons(nshg)->nsons(nfath) 616 K2) x and shape functions passed into BCint.f 617 K3) genbc.f now uses nenl (properly) instead of nen 618 K4) genini gets shape functions (for K2) 619 K5) dec-> decalp for ifdefing in c routines 620 K6) 1 point quadrature added for wedges 621 K7) extra arrays removed from solvecon.f 622 K8) rotating frame added back to incompressible code 623 K9) surface forces added to forces.f 624K10) itrdrv now checks at each step to see if job should be stopped 625 626 62711/15/01 628 629Fixed a bug that was preventing us from using the temporally averaged 630dynamic model without wall functions. rwvelb.f now writes 631bar.<stepno>.<procno> instead of velb.<stepno>.<procno>. In that file 632you may find the y field averaged over the past runs, the dynamic 633model numerater and denominator coefficient or both. 634 635This commit also brought in changes to the sponge as needed for the 636jet simulation (from 1 diameter to 1.5 diameters downstream). Changes 637only in e3source.f 638 639 640CVS: ---------------------------------------------------------------------- 641CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 642CVS: 643CVS: Committing in common 644CVS: 645CVS: Modified Files: 646CVS: Cinput.cpp GaussLegendreSimplex.c asbnabi.f asithf.f common.h 647CVS: common_c.h dtn.f filtprep.f genbkb.f gendat.f genibc.f 648CVS: getvel.f input.config input_fform.cpp new_interface.c 649CVS: phasta.cpp proces.f settauw.f shp6w.c shphex.c shptet.c 650CVS: Added Files: 651CVS: bardmc.f cmass.f getdmc.f getstrl.f vec_func.h 652CVS: Removed Files: 653CVS: les.h 654CVS: ---------------------------------------------------------------------- 655 656 657CVS: ---------------------------------------------------------------------- 658CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 659CVS: 660CVS: Committing in compressible 661CVS: 662CVS: Modified Files: 663CVS: bc3per.f e3.f e3b.f e3bvar.f e3cool.f e3dc.f e3ivar.f e3mtrx.f 664CVS: e3source.f e3tau.f elmgmr.f elmmfg.f get_h.c getthm.f itrbc.f 665CVS: itrdrv.f qpbc.f rstat.f solgmr.f 666CVS: Removed Files: 667CVS: bardmc.f getdmc.f getstrl.f vec_func.h 668CVS: ---------------------------------------------------------------------- 669 670 671CVS: ---------------------------------------------------------------------- 672CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 673CVS: 674CVS: Committing in incompressible 675CVS: 676CVS: Modified Files: 677CVS: bflux.f e3b.f elmgmr.f itrdrv.f qpbc.f scatnu.f usr.c 678CVS: Removed Files: 679CVS: Input.h auxmpi.h bardmc.f getdmc.f getstrl.f advLES.f filters.f 680CVS: ---------------------------------------------------------------------- 681 682 683not yet tagged changes 684 685commu can now handle multiple communications per processor boundary 686(bi-directional). 687 688parallel scalars now work (for sgi only as it was a bug in the solver 689library and other architectures have not been updated yet). 690appropriate changes were also made to file usr.c, routine myflesnew 691 692 693STEP sequence size removed from input.config since it is calculated now 694 695CVS: ---------------------------------------------------------------------- 696CVS: Enter Log. Lines beginning with `CVS:' are removed automatically 697CVS: 698CVS: Committing in common 699CVS: 700CVS: Modified Files: 701CVS: commu.f input.config 702CVS: ---------------------------------------------------------------------- 703CVS: Committing in incompressible 704CVS: 705CVS: Modified Files: 706CVS: les.h usr.c 707CVS: ---------------------------------------------------------------------- 708 709For previous changes see phSolverI v1_6_0 and phSolverC v1_5_3 710