sort by:
Revision Author Date Message Commit Date
0c84e99 Read in sdepv_misfit regardless of rheology. 14 September 2010, 01:08 UTC
bf59b6b Added option to start with isotropic viscosity during first anisotropic iteration. 13 September 2010, 14:07 UTC
a112ea8 Forgot part of the rotation matrix, minor fixes else, still experimental. 13 September 2010, 12:01 UTC
ab35881 Made sure current version compiled nicely without additional flags. Anisotropic viscosity is set up more flexibly. 13 September 2010, 07:13 UTC
b5d6c30 Made sure that viscosity prefactors get assigned for material dependent viscosity only (TDEPV = off) 07 September 2010, 07:05 UTC
a345eef Added random number generator initialization to ensure exact reproducibility of all runs, including those with tracers. Earlier reported solver issues might be related to our cluster, at least all revisions seems to give the same answer now. 07 September 2010, 02:22 UTC
1cf026d Checking in current version because of need for debugging. Revision 17172, compared to (presumably) 17153 give strange multigrid convergence behavior. Checking why that is now. 07 September 2010, 01:40 UTC
1efa65c Trial implementation of anisotropic viscosity. 06 September 2010, 05:20 UTC
4ecc25c Added support of multigrid solver in Exchanger. Lifted the solver assertion in Coupler. Cookbook9 can be converted to multigrid solver once appropriate parameters (*solver.mesher.levels=2 and *solver.vsolver.Solver=multigrid) are added. However, the embedded solver will converge incorrectly. Setting esolver.vsolver.accuracy=6e-3 seems to fix the problem. 30 August 2010, 20:41 UTC
bd28b8f And a typo fix. 25 May 2010, 02:47 UTC
f71b317 toplayerbc > 0 will now assign BCs for all noes with r >= toplayerbc_r, which defaults to 0.984303876942 i.e. 100 km depth. (Still experimental.) 24 May 2010, 18:00 UTC
ed7e8ae Improvements for internal BCs. 12 March 2010, 00:41 UTC
b732e74 This logic should make more sense for internal BCs, but still testing. 10 March 2010, 19:55 UTC
8a2427e Experimental implementation of a single node layer for internal boundary condition assignment, assigned if toplayerbc < 0, to nodes at layer noz+toplayerbc, but only for the top level multigrid, for now. 10 March 2010, 19:38 UTC
70a834f Fix a misleading comment 26 February 2010, 21:04 UTC
4ff7faa Fix an error when running autoreconf -i 02 February 2010, 23:41 UTC
f5834c5 Changed internal variable name for use_se_smoothing, as per Eh's suggestion. 01 February 2010, 19:43 UTC
2a3b029 Added old CitcomS element -> nodes -> element viscosity smoothing with new use_ne_visc_smooth parameter, default is off. 29 January 2010, 14:48 UTC
8bb44b5 Added pyre input for the new parameters added in r16138. 21 January 2010, 21:23 UTC
3ad8c8c Replaced fixed accuracy setting. Bug reported by Rob Moucha. 21 January 2010, 21:22 UTC
c4154bc Added opendx macro to visualize opt (e.g. composition) fields. This fixed Issue99. 21 January 2010, 20:11 UTC
0dc0383 Preliminary support for the absolute tracer method. Only tested with one flavor, reading in from a file, and the initial volume of the tracers is HARD-CODED! In addition, the code automaticaly truncates the composition to one, although this should be a switch for the user to specify. 17 January 2010, 01:19 UTC
110bd1c Implementation of two, perhaps entirely unnecessary features, in search of why certain prescribed velocity models do not converge like they should. - allow for each inner loop solution to have a higher accuracy than the outer, Uzawa loop. This is set by new parameter "inner_accuracy_scale", which pre-multiplies the inner accuracy (imp) setting (defaults to unity) - remove any rigid body rotation at each Uzawa iteration step. Parmaeter "inner_remove_rigid_rotation", by default off. For this, I also added an experimental routine assign_v_to_vector() in order to reasign the V velocities to the U solution vector after removing the net rotations from V (how is supposed to work for the regular operational mode?) 16 January 2010, 02:47 UTC
33e0d0b Replaced TRUE and FALSE. Added function prototypes. 09 December 2009, 02:25 UTC
ef97d2c Added two boolean parameters for controlling convergence criterion: check_continuity_convergence (default=on): include div/v in the criterion. check_pressure_convergence (default=on): include dp/p in the criterion. Retired the parameter: only_check_vel_convergence. 08 December 2009, 21:44 UTC
7efded4 removed unused variables in Drive_solvers.c 08 December 2009, 21:43 UTC
e829bf5 Merged solve_constrained_flow_iterative_pseudo_surf() with void solve_constrained_flow_iterative(); assemble_forces_pseudo_surf() wi Removed used python bindings in module/stokes_solver.c 08 December 2009, 21:43 UTC
0c6082f Truly apply internal boundary conditions on horizontal surfaces for E->mesh.toplayerbc > 0 (default, of course, 0) 08 December 2009, 17:19 UTC
025e452 Modified handling of internal BCs, as per Eh's suggestion, now handled in a separate function. 04 December 2009, 23:16 UTC
21c984f Removed metioning the version number in README. Added the url of the manual. 04 December 2009, 21:57 UTC
e20fdfc Suggested unified treatment of horizontal_bc (moved to BC_util) for review by Eh Improved handling of grd-read velocity and traction boundary conditions. 04 December 2009, 17:44 UTC
45144ae Temporary fix for regular operations, for Eh's review. Ggrd assignment still needs tending to. 04 December 2009, 02:33 UTC
1f61de1 Experimental implementation of internal stress or velocity boundary conditions. As of right now, toplayerbc > 0 (default = 0) will assign the surface boundary condition to all nodes within the first toplayerbc layers. If BCs are assigned from file, this will only work for the grd method, not the regular file I/O. 03 December 2009, 23:51 UTC
adf8b19 Exit criterion for incompressibility iteration allowed for velocity and pressure convergence at finite compressibility > accuracy. Changed criterion to loop while (convergence < 2) OR (compressibility < imp) Please review Also modified the sloppy iteration branch, where pressure isn't checked. 21 November 2009, 00:15 UTC
5f38ecb Clean up #include 20 November 2009, 21:58 UTC
aca8438 Fixed two minor errors 20 November 2009, 21:58 UTC
0ba4cb4 Created a new array E-> for the radii of all layers. This makes gather_r() obsolete. 19 November 2009, 23:11 UTC
7ae00fb refactor viscosity definition 19 November 2009, 23:09 UTC
8e9f046 'localhost' can abbr. as '-' 19 November 2009, 23:08 UTC
09cea7d Changed defaults for z_410, z_lmantle, z_lith to more Earth like parameters. Does this break any tests? 14 October 2009, 18:26 UTC
3ee0608 Fixed typo. 13 October 2009, 16:37 UTC
0271108 z_layer needs 4 elements 08 October 2009, 18:01 UTC
4e0b2f4 added only_check_vel_convergence to pyre version. 07 October 2009, 04:39 UTC
4b0faab Modified handling of ggrd_mat_control, > 0 will select all layers above, < 0 will select individual layer == - ggrd_mat_control 06 October 2009, 18:48 UTC
ee214e5 Add pyre binding of z_layer, r_grid_layers, rr, and nr. replace MAX_MAT with CITCOM_MAX_VISC_LAYER. 05 October 2009, 19:06 UTC
2622434 Refactoring z_layer input validation. 05 October 2009, 19:05 UTC
599ab80 Dealing with tracers that belong to the current processor but cannot find its element. When itracer_warning is off, the code will ignore the lost tracers. When itracer_warning is on (default), the code will abort. 05 October 2009, 19:05 UTC
a3b2b41 Experimental traction surface BC from grds. 02 October 2009, 23:16 UTC
11d52ac visc_layer_control: allowing each radial layer has its own viscosity parameters When "visc_layer_control" is true, the file "visc_layer_file" is read in. The format of visc_layer_file: param_nameA Aval_1 Aval_2 ... Aval_nel param_nameB Bval_1 Bval_2 ... Bval_nel . . . The first line is assigned to the basal layer (CMB), and the last line to the surface layer. Note that each parameter has nodez-1 values. 28 September 2009, 20:08 UTC
e407de4 clean up unused viscosity variables 28 September 2009, 20:05 UTC
bb0c1ea Fixed a bug in pseudo free surface, reported by Robert Moucha. 28 September 2009, 20:05 UTC
f173a47 remove duplicating prototype declaration 28 September 2009, 20:04 UTC
5f0f67f Added pyre binding for write_q_files 03 September 2009, 17:24 UTC
3024d9f Added a warning message about tracer + coupling 03 September 2009, 17:24 UTC
37b376d Forgot one valid check for Uzawa iterations. 25 August 2009, 19:54 UTC
ba42e48 Ensure that Uzawa solver does not converge if low level, MG solver does not yield a valid solution. 25 August 2009, 17:31 UTC
57faf89 test implementation of steady state strain-rate weakending via PDEPV psrw=on parameter. 21 August 2009, 20:32 UTC
4b24572 Add Murnaghan's integrated linear EoS as reference_state=2 11 August 2009, 17:48 UTC
9157604 Add gather_r() to gather all radial coordinates 11 August 2009, 17:47 UTC
a0d609e Adding augmented Lagrangian component of stiffness matrix in cbf topo 11 August 2009, 17:35 UTC
e5a69cd The function collect_terms() turns out to be quite expensive (over 70% according to gprof). Perform it once per level, caching the result. 04 August 2009, 01:45 UTC
6957a07 Fixed the remaining bug in n_assemble_del2_u(). Glancing at "global_defs.h", I had assumed that VBX, VBY, and VBZ were defined in the natrual way: #define VBX 0x2 #define VBY 0x4 #define VBZ 0x8 I.e., "0x1 << doff"... but today I looked at it more closely; Y & Z are reversed!!! #define VBX 0x2 #define VBZ 0x4 #define VBY 0x8 With this fix, the CUDA multigrid solver seems to work correctly. 29 July 2009, 20:16 UTC
6a06281 Ensure debugging output for ggrd_vtop mixed BC if more than one CPU with depth is used. Previous check in of precise_strain_rate=on mode seems to be OK, stress, geoid, and dynamic topo output (for cbf topo off) have been tested. 29 July 2009, 04:05 UTC
b9278a2 Experiments show that using a block size of 1 leads to poor performance. Using the recommended minimum of 64, gauss_seidel_0 and gauss_seidel_1 are about 3x to 4x faster. 28 July 2009, 01:00 UTC
b500001 Fixed stupid reduce bug that affected several routines. (This is what I get for coping & pasting example code without thinking.) Now the only remaining broken routine is n_assemble_del2_u(). 27 July 2009, 23:44 UTC
9a97294 Wrote two new versions of gauss_seidel_2() -- neither of which is very parallel. Both produce similar wrong answers. 25 July 2009, 02:48 UTC
fb1101a This version of precise strain-rates may now be worth testing out. 24 July 2009, 01:03 UTC
1c4fd5d Corrected call to get_ba in Topo_gravity for precise_strain_rate, still testing. Rest should be cosmetic. 23 July 2009, 18:40 UTC
b2a7cb3 This junk is my effort to understand gauss_seidel_2(). Because of the way E->temp[] is used, it is quite difficult (impossible?) to parallelize. I even wrote code to generate a "dot" graph of the expression tree. Of course, this tree is too big to plot in its entirety. With neq == 2187, the tree has 87479 nodes and is 355 levels deep; each level can have anywhere from a handful to several hundred nodes, with over a thousand nodes in each of the bottommost layers. 23 July 2009, 03:50 UTC
02a3b02 Tentative check in of stress computation modification such that Eh can take a look of what we're working on. 22 July 2009, 20:17 UTC
5bfdbd3 Created drop-in, host-side replacements for each of the multigrid kernel functions. After swapping them in/out, it seems n_assemble_del2_u(), gauss_seidel_2(), and gauss_seidel_3() are all buggy :-( Both gauss_seidel_0() and gauss_seidel_1() appear to be OK... but that isn't worth writing home about, because they are both trivial. 17 July 2009, 21:22 UTC
9d02f10 New scaffolding with annotations. 16 July 2009, 03:13 UTC
25fc342 Bug fixes... still doesn't work right. 16 July 2009, 00:49 UTC
60515e6 Wrote code to allocate & initialize device memory, and copy the result from the device. 15 July 2009, 21:54 UTC
0bef9e8 Parallelized the last fragment of the gauss_seidel() main loop as gauss_seidel_3(). It turns out to be very similar to the second half of n_assemble_del2_u(). Tore down scaffolding. 15 July 2009, 01:58 UTC
362c2e3 Parallel n_assemble_del2_u() takes shape. 15 July 2009, 00:14 UTC
3e27612 Pardon my dust... tally_n_assemble_del2_u() is instrumention and printfs I wrote in order to understand n_assemble_del2_u(). As with e_assemble_del2_u(), it appears I need one block for each Au[i] element, with a reduction in shared memory. 10 July 2009, 22:15 UTC
5825a00 include density (rho) scaling in convective heat flux (necessary for compressible flow) 10 July 2009, 21:08 UTC
e0d4a44 Second material element file only read in and processed when (pos_age) criteria met (i.e. consistent with vel b.c. approach). 09 July 2009, 20:40 UTC
0ee4a4c Only one processor complains about removing rigid rotations while prescribing plate motions. 09 July 2009, 17:27 UTC
67d7472 Make sure "-arch sm_13" is given -- for double-precision floating-point. Added missing #ifndef USE_CUDA. 06 July 2009, 22:13 UTC
735ae7c Sketched data-parallel version of gauss_seidel(). I haven't figured out how to parallelize the following: Au[C[i]] += ... Ad[C[j]] += ... 02 July 2009, 18:41 UTC
43a1aa4 Fixed a pair of mistakes I made: under cgrad, don't replace solve_del2_u(); under multigrid, 'level' is not a constant. 01 July 2009, 19:09 UTC
678e5fc The 'temp1' array is unnecessary since nproc == 1. 30 June 2009, 22:33 UTC
9743b5e For my own sanity, simplified code assuming the following: E->sphere.caps_per_proc == 1 E->parallel.nproc == 1 E->mesh.levmax == 0 E->mesh.nsd == 3 29 June 2009, 22:10 UTC
87a1bca Proposed CUDA kernel for multigrid solver. 29 June 2009, 20:18 UTC
d01a1c7 forgot to add these files 25 June 2009, 16:51 UTC
f038904 Update ChangeLog to r15180 15 June 2009, 18:21 UTC
9a9cbd3 Added COND_GGRD flags; removed redundant LIBS 10 June 2009, 20:42 UTC
c5915dc Checking malloc.h earlier to initialize even when not using pyre 10 June 2009, 20:41 UTC
055dd50 checking netcdf header and library 10 June 2009, 17:49 UTC
aa8a0ac Append GMTHOME/lib to LD path 10 June 2009, 04:11 UTC
55d0112 Adding ggrd compilation support 09 June 2009, 23:47 UTC
03dafa6 One more feature and mention where to find the features in the cookbooks. 08 June 2009, 21:03 UTC
234f0ef Added a safe guard 08 June 2009, 18:11 UTC
ba6a738 Reordering parameter 04 June 2009, 21:04 UTC
bc9a15c Update SVTInlet to inherit from BaseSVTInlet 03 June 2009, 21:08 UTC
56e2c06 adding files to makefile 03 June 2009, 21:07 UTC
b3dfd01 update to current coupler 02 June 2009, 22:58 UTC
7fd3100 reducing the accuracy of solvers, otherwise, they will converge to the wrong result 02 June 2009, 22:58 UTC
back to top