{"action":"created","issue":{"url":"https://api.github.com/repos/geodynamics/aspect/issues/871","repository_url":"https://api.github.com/repos/geodynamics/aspect","labels_url":"https://api.github.com/repos/geodynamics/aspect/issues/871/labels{/name}","comments_url":"https://api.github.com/repos/geodynamics/aspect/issues/871/comments","events_url":"https://api.github.com/repos/geodynamics/aspect/issues/871/events","html_url":"https://github.com/geodynamics/aspect/issues/871","id":154553715,"number":871,"title":"\"geometric average\" material averaging crashes in 3D","user":{"login":"jperryhouts","id":6118266,"avatar_url":"https://avatars.githubusercontent.com/u/6118266?v=3","gravatar_id":"","url":"https://api.github.com/users/jperryhouts","html_url":"https://github.com/jperryhouts","followers_url":"https://api.github.com/users/jperryhouts/followers","following_url":"https://api.github.com/users/jperryhouts/following{/other_user}","gists_url":"https://api.github.com/users/jperryhouts/gists{/gist_id}","starred_url":"https://api.github.com/users/jperryhouts/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/jperryhouts/subscriptions","organizations_url":"https://api.github.com/users/jperryhouts/orgs","repos_url":"https://api.github.com/users/jperryhouts/repos","events_url":"https://api.github.com/users/jperryhouts/events{/privacy}","received_events_url":"https://api.github.com/users/jperryhouts/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":7,"created_at":"2016-05-12T18:46:17Z","updated_at":"2016-05-13T11:35:42Z","closed_at":null,"body":"Maybe just on my machine, but maybe something we should look into before 1.4 release? Can be reproduced with this input:\r\n\r\n(works if you change to Dimension=2 or Material averaging = anything else)\r\n```\r\n set Dimension                               = 3\r\n set Output directory                        = build/output\r\n \r\n subsection Geometry model\r\n   set Model name = chunk\r\n   subsection Chunk\r\n     set Chunk inner radius      = 1\r\n     set Chunk outer radius      = 2\r\n     set Chunk minimum latitude  = 37\r\n     set Chunk maximum latitude  = 49\r\n     set Chunk minimum longitude = -132\r\n     set Chunk maximum longitude = -116\r\n   end\r\n end\r\n \r\n subsection Model settings\r\n   set Fixed temperature boundary indicators     = inner, outer\r\n   set Tangential velocity boundary indicators   = outer\r\n   set Zero velocity boundary indicators         = inner\r\n end\r\n \r\n subsection Material model\r\n   set Material averaging    = geometric average\r\n   set Model name            = simple\r\n end\r\n \r\n subsection Initial conditions\r\n   set Model name = adiabatic\r\n \r\n   subsection Adiabatic\r\n     set Age bottom boundary layer = 1e6\r\n     set Age top boundary layer = 5e6\r\n     set Amplitude = 100\r\n     set Position = center\r\n     set Radius = 100e3\r\n     set Subadiabaticity = 0\r\n   end\r\n end\r\n \r\n subsection Boundary temperature model\r\n   set Model name    = initial temperature\r\n end\r\n \r\n subsection Gravity model\r\n   set Model name    = radial earth-like\r\n end\r\n \r\n subsection Mesh refinement\r\n   set Initial global refinement             = 3\r\n   set Initial adaptive refinement           = 0\r\n end\r\n \r\n subsection Termination criteria\r\n   set Termination criteria = end step\r\n   set End step = 0\r\n end\r\n```"},"comment":{"url":"https://api.github.com/repos/geodynamics/aspect/issues/comments/219019510","html_url":"https://github.com/geodynamics/aspect/issues/871#issuecomment-219019510","issue_url":"https://api.github.com/repos/geodynamics/aspect/issues/871","id":219019510,"user":{"login":"tjhei","id":1531285,"avatar_url":"https://avatars.githubusercontent.com/u/1531285?v=3","gravatar_id":"","url":"https://api.github.com/users/tjhei","html_url":"https://github.com/tjhei","followers_url":"https://api.github.com/users/tjhei/followers","following_url":"https://api.github.com/users/tjhei/following{/other_user}","gists_url":"https://api.github.com/users/tjhei/gists{/gist_id}","starred_url":"https://api.github.com/users/tjhei/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/tjhei/subscriptions","organizations_url":"https://api.github.com/users/tjhei/orgs","repos_url":"https://api.github.com/users/tjhei/repos","events_url":"https://api.github.com/users/tjhei/events{/privacy}","received_events_url":"https://api.github.com/users/tjhei/received_events","type":"User","site_admin":false},"created_at":"2016-05-13T11:35:42Z","updated_at":"2016-05-13T11:35:42Z","body":"[0]PETSC ERROR: Caught signal number 8 FPE: Floating Point Exception,probably divide by zero\r\n```\r\n*** Timestep 0:  t=0 years\r\n   Solving temperature system... 0 iterations.\r\n\r\nProgram received signal SIGFPE, Arithmetic exception.\r\naspect::Assemblers::CompleteEquations<3>::local_assemble_stokes_system_incompressible (this=0x154fff0, pressure_scaling=1.4275000000000002e+27, rebuild_stokes_matrix=true, scratch=..., data=...) at /ssd/aspect-local/source/simulator/assembly.cc:1061\r\n1061\t                    data.local_matrix(i,j) += ( (use_tensor ?\r\n(gdb) p eta\r\n$2 = inf\r\n```\r\nit turns out the product of the viscosities overflows to infinity because\r\n```\r\np values_out\r\n$14 = std::vector of length 27, capacity 27 = {5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24,\r\n 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+2\r\n4, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e\r\n+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.000000000000000\r\n5e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000\r\n005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24, 5.0000000000000005e+24}\r\n```\r\n:-)"}}