https://github.com/geodynamics/citcoms

sort by:
Revision Author Date Message Commit Date
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:26 UTC
3ee0608 Fixed typo. 13 October 2009, 16:37:49 UTC
0271108 z_layer needs 4 elements 08 October 2009, 18:01:00 UTC
4e0b2f4 added only_check_vel_convergence to pyre version. 07 October 2009, 04:39:33 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:28 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:34 UTC
2622434 Refactoring z_layer input validation. 05 October 2009, 19:05:33 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:04 UTC
a3b2b41 Experimental traction surface BC from grds. 02 October 2009, 23:16:39 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:25 UTC
e407de4 clean up unused viscosity variables 28 September 2009, 20:05:38 UTC
bb0c1ea Fixed a bug in pseudo free surface, reported by Robert Moucha. 28 September 2009, 20:05:03 UTC
f173a47 remove duplicating prototype declaration 28 September 2009, 20:04:38 UTC
5f0f67f Added pyre binding for write_q_files 03 September 2009, 17:24:47 UTC
3024d9f Added a warning message about tracer + coupling 03 September 2009, 17:24:19 UTC
37b376d Forgot one valid check for Uzawa iterations. 25 August 2009, 19:54:47 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:36 UTC
57faf89 test implementation of steady state strain-rate weakending via PDEPV psrw=on parameter. 21 August 2009, 20:32:56 UTC
4b24572 Add Murnaghan's integrated linear EoS as reference_state=2 11 August 2009, 17:48:33 UTC
9157604 Add gather_r() to gather all radial coordinates 11 August 2009, 17:47:36 UTC
a0d609e Adding augmented Lagrangian component of stiffness matrix in cbf topo 11 August 2009, 17:35:23 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:56 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:08 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:42 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:36 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:52 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:38 UTC
fb1101a This version of precise strain-rates may now be worth testing out. 24 July 2009, 01:03:47 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:41 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:26 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:11 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:38 UTC
9d02f10 New scaffolding with annotations. 16 July 2009, 03:13:33 UTC
25fc342 Bug fixes... still doesn't work right. 16 July 2009, 00:49:42 UTC
60515e6 Wrote code to allocate & initialize device memory, and copy the result from the device. 15 July 2009, 21:54:24 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:53 UTC
362c2e3 Parallel n_assemble_del2_u() takes shape. 15 July 2009, 00:14:53 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:02 UTC
5825a00 include density (rho) scaling in convective heat flux (necessary for compressible flow) 10 July 2009, 21:08:14 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:10 UTC
0ee4a4c Only one processor complains about removing rigid rotations while prescribing plate motions. 09 July 2009, 17:27:41 UTC
67d7472 Make sure "-arch sm_13" is given -- for double-precision floating-point. Added missing #ifndef USE_CUDA. 06 July 2009, 22:13:05 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:00 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:31 UTC
678e5fc The 'temp1' array is unnecessary since nproc == 1. 30 June 2009, 22:33:19 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:40 UTC
87a1bca Proposed CUDA kernel for multigrid solver. 29 June 2009, 20:18:25 UTC
d01a1c7 forgot to add these files 25 June 2009, 16:51:10 UTC
f038904 Update ChangeLog to r15180 15 June 2009, 18:21:41 UTC
9a9cbd3 Added COND_GGRD flags; removed redundant LIBS 10 June 2009, 20:42:43 UTC
c5915dc Checking malloc.h earlier to initialize even when not using pyre 10 June 2009, 20:41:31 UTC
055dd50 checking netcdf header and library 10 June 2009, 17:49:35 UTC
aa8a0ac Append GMTHOME/lib to LD path 10 June 2009, 04:11:35 UTC
55d0112 Adding ggrd compilation support 09 June 2009, 23:47:35 UTC
03dafa6 One more feature and mention where to find the features in the cookbooks. 08 June 2009, 21:03:32 UTC
234f0ef Added a safe guard 08 June 2009, 18:11:24 UTC
ba6a738 Reordering parameter 04 June 2009, 21:04:33 UTC
bc9a15c Update SVTInlet to inherit from BaseSVTInlet 03 June 2009, 21:08:01 UTC
56e2c06 adding files to makefile 03 June 2009, 21:07:04 UTC
b3dfd01 update to current coupler 02 June 2009, 22:58:45 UTC
7fd3100 reducing the accuracy of solvers, otherwise, they will converge to the wrong result 02 June 2009, 22:58:24 UTC
748c305 Added parameters 'amending_outflow' and 'exchange_pressure' to help the convergence of esolver When 'amending_outflow' is set to true, the imposed velocity BC will be amended slightly to be divergence-free. (The divergence, e.g. outflow, is caused by the combination of solver inaccuracy and interpolation inaccuracy.) When 'exchange_pressure' is set to true, the initial pressure (at element level) of the embedded solver is taken from the pressure solution of the containing solver. 02 June 2009, 22:56:46 UTC
179463a Boundary for global mesh, not sure whether it works with interpolator or not 02 June 2009, 22:09:25 UTC
2ef0c3d init'd E->output.seismic and E->output.coord_bin 02 June 2009, 21:58:40 UTC
101123f refining cookbook10 02 June 2009, 21:58:23 UTC
e1af994 CUDA-ized e_assemble_del2_u(), and ran it successfully -- and slowly! -- under the device emulator. Added a new 'configure' option: "--with-cuda". 28 May 2009, 01:50:40 UTC
500e3f5 Parallelized strip_bcs_from_residual(). 27 May 2009, 02:32:03 UTC
511d3f3 Sketched data-parallel version of e_assemble_del2_u(). 27 May 2009, 01:28:55 UTC
060b73c Added blob_bc_persist parameter (default = off) to allow assigning blob-like temperature anomalies to BCs. 26 May 2009, 15:21:42 UTC
5a6accf Fixed typo in blob temperature assignment. Changed conversion of coordinates, now assignment works. 26 May 2009, 15:03:36 UTC
accc9b4 add changes in r14999 21 May 2009, 23:16:27 UTC
f99c960 ignoring the degree-0 and degree-1 components of geoid 15 May 2009, 22:32:48 UTC
1670e55 add more chaanges in v3.1.0 10 May 2009, 23:18:05 UTC
a622619 Change version number for v3.1.0 alpha release 10 May 2009, 16:37:31 UTC
658d80a Summary of changes in this release 10 May 2009, 16:35:07 UTC
573f5d0 Adding cookbook for synthetic seismogram, mostly copied from cookbook7 01 May 2009, 22:19:35 UTC
7a399ef Depth is measured from r=1.0, not outer radius, because in coupled model, the outer radius can be < 1.0. 01 May 2009, 22:18:00 UTC
c9c08bb Printing a big warning when plate velocity is imposed but remove_angular_momentum or remove_rigid_rotation is on. 01 May 2009, 22:17:27 UTC
24392fd Removing angular momentum in pseudo free surface case 01 May 2009, 22:15:05 UTC
5ab9dda ignore anything after ';' 01 May 2009, 22:14:29 UTC
b215c86 explicitly add Python.h to requirement 01 May 2009, 22:13:26 UTC
05da830 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 28 April 2009, 22:53:38 UTC
2064543 Moved malloc/free calls outside of kernel. (The amount of memory allocated appears to be invariant.) 27 April 2009, 20:15:10 UTC
fa871b7 Consolidated calls to malloc/free. 27 April 2009, 18:39:15 UTC
2d995dd Added flag to allow using the more precise method of computing strain-rates throughout the globe, and not just close to the poles. precise_strain_rate=on will select this, default is "odd", i.e. poles only as before. 24 April 2009, 18:37:49 UTC
0eb5c37 Don't call through function pointers -- otherwise, 'nvcc' crashes! 23 April 2009, 02:32:01 UTC
185800d Proposed CUDA kernel for conjugate gradient solver. 23 April 2009, 02:22:19 UTC
a6bc514 Avoid core dump in record() when E->fp hasn't been set yet. 22 April 2009, 00:59:20 UTC
e45eab4 The new Python.h header check reorders this script somewhat. CIT_CHECK_PYTHON_HEADER must follow AC_PROG_CC, and precede CIT_PYTHON_EGG_SETUP. 21 April 2009, 22:29:26 UTC
c46afd8 Issue45: Check for Python.h. 21 April 2009, 19:33:59 UTC
756b97b Added viscosity options 9 and 10 for TDPEPV which are identical to 3 and 8 but temperature is not limited to [0;1]. This is cumbersome, but perhaps needed for backward compatibility. 17 April 2009, 03:37:28 UTC
9ea8a44 Move parameter "mantle_temp" from CitcomS.solver.param to CitcomS.solver.ic. Renamed E->control.lith_age_mantle_temp to E->control.mantle_temp. 16 April 2009, 18:24:28 UTC
0591133 update cookbook8 with the manual 16 April 2009, 18:23:06 UTC
487eeea Add -DUSE_HDF5 to CPPFLAGS in module/Exchanger 16 April 2009, 18:22:22 UTC
075b734 change the viscosity section to be the same as that in cookbook3. 16 April 2009, 18:21:50 UTC
2c35d25 Removed domain extents, since they are determined by coor_file. 16 April 2009, 18:20:08 UTC
a7389af Removed instructions to 'make install', 'make clean' and 'make distclean'. Added dependencies. 16 April 2009, 18:07:53 UTC
a2a39ef Fixed a bug in cgrad solver introduced in r13270. 09 April 2009, 20:53:58 UTC
7dbc76c Left debugging output in, now cleaned up 09 April 2009, 20:39:24 UTC
d5b89ca Fixed file output issue which caused vtk log file I/O to bomb during Citcom finalize. 09 April 2009, 20:37:59 UTC
back to top