@lofcz/tikz 0.0.0-stage → 0.1.0

package/runtime/dist/bundles/latex-extra/files/tex/generic/pgfplots/pgfplotsticks.code.tex

@@ -0,0 +1,2589 @@
+%--------------------------------------------
+%
+% Package pgfplots
+%
+% Provides a user-friendly interface to create function plots (normal
+% plots, semi-logplots and double-logplots).
+%
+% It is based on Till Tantau's PGF package.
+%
+% Copyright 2007/2008 by Christian Feuersänger.
+%
+% This program is free software: you can redistribute it and/or modify
+% it under the terms of the GNU General Public License as published by
+% the Free Software Foundation, either version 3 of the License, or
+% (at your option) any later version.
+%
+% This program is distributed in the hope that it will be useful,
+% but WITHOUT ANY WARRANTY; without even the implied warranty of
+% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+% GNU General Public License for more details.
+%
+% You should have received a copy of the GNU General Public License
+% along with this program.  If not, see <http://www.gnu.org/licenses/>.
+%
+%--------------------------------------------
+
+\newif\ifpgfplots@log@tick@isminor@tick@pos
+% Checks whether the tick position given as #1.#2=log10(T) belongs to
+% T=i*10^j with an integer i>1.
+%
+% If T=i*10^j,  \ifpgfplots@log@tick@isminor@tick@pos will be set to true and
+% \pgfmathresult will contain T.
+%
+% Otherwise, \ifpgfplots@log@tick@isminor@tick@pos will be set to false and
+% pgfmathresult to #1.#2
+%
+% Arguments:
+% #1.#2  the value log10(T)
+%
+% Implementation:
+% if T = i*10^j,  log10(T) = log10(i) + j.
+% That means if log10(T) in \Z,  we have T = 10^j. If not, we need to
+% check wether i is an integer. Please note that log10(i) < 1.
+%
+% Further note: log(T) < 0 <=>  j<0.
+% In case j<0, we have
+%   #1.#2 = j + log(i)
+%         = - ( -j - log(i) )
+%         = - ( -j - 1  + (1-log(i)) )
+%         = #1 '.' #2 [ up to the '0.'
+% that means #1 = j-1  and #2 = 1-log(i).
+\def\pgfplots@is@log@tick@a@minor@tick@pos#1.#2\relax{%
+	\pgfmathapproxequalto@{#1.#2}{#1.0}%
+	\ifpgfmathcomparison
+		% in MOST cases, this here will be true:
+		\pgfplots@log@tick@isminor@tick@posfalse
+		\def\pgfmathresult{#1.#2}%
+	\else
+		% I guess this won't happen too often. In fact, it's a very
+		% special case.
+		\begingroup
+		\c@pgf@counta=#1\relax
+		\ifdim#1.#2pt<0pt %
+			\advance\c@pgf@counta by-1
+			\pgfmathsubtract@{1}{0.#2}%
+			\expandafter\pgfplots@is@log@tick@a@minor@tick@pos@IDENTIFY@LOGi\pgfmathresult\relax
+			\ifpgfplots@log@tick@isminor@tick@pos
+				\aftergroup\pgfplots@log@tick@isminor@tick@postrue
+				\edef\pgfmathresult{\pgfmathresult e\the\c@pgf@counta}%
+			\else
+			\aftergroup\pgfplots@log@tick@isminor@tick@posfalse
+				\def\pgfmathresult{#1.#2}%
+			\fi
+		\else
+			\pgfplots@is@log@tick@a@minor@tick@pos@IDENTIFY@LOGi0.#2\relax
+			\ifpgfplots@log@tick@isminor@tick@pos
+				\aftergroup\pgfplots@log@tick@isminor@tick@postrue
+				\edef\pgfmathresult{\pgfmathresult e\the\c@pgf@counta}%
+			\else
+				\aftergroup\pgfplots@log@tick@isminor@tick@posfalse
+				\def\pgfmathresult{#1.#2}%
+			\fi
+		\fi
+		\pgfmath@smuggleone\pgfmathresult
+		\endgroup
+	\fi
+}
+
+% expects a positive number.
+\def\pgfplots@is@log@tick@a@minor@tick@pos@IDENTIFY@LOGi0.#1\relax{%
+	\pgfplots@log@tick@isminor@tick@postrue
+	\pgfmathapproxequalto@{0.#1}{0.3010299956639}%
+	\ifpgfmathcomparison
+		\def\pgfmathresult{2}%
+	\else
+		\pgfmathapproxequalto@{0.#1}{0.4771212547196}%
+		\ifpgfmathcomparison
+			\def\pgfmathresult{3}%
+		\else
+			\pgfmathapproxequalto@{0.#1}{0.6020599913279}%
+			\ifpgfmathcomparison
+				\def\pgfmathresult{4}%
+			\else
+				\pgfmathapproxequalto@{0.#1}{0.698970004}%
+				\ifpgfmathcomparison
+					\def\pgfmathresult{5}%
+				\else
+					\pgfmathapproxequalto@{0.#1}{0.7781512503}%
+					\ifpgfmathcomparison
+						\def\pgfmathresult{6}%
+					\else
+						\pgfmathapproxequalto@{0.#1}{0.8450980400}%
+						\ifpgfmathcomparison
+							\def\pgfmathresult{7}%
+						\else
+							\pgfmathapproxequalto@{0.#1}{0.9030899869}%
+							\ifpgfmathcomparison
+								\def\pgfmathresult{8}%
+							\else
+								\pgfmathapproxequalto@{0.#1}{0.954242509439}%
+								\ifpgfmathcomparison
+									\def\pgfmathresult{9}%
+								\else
+									\pgfplots@log@tick@isminor@tick@posfalse
+								\fi
+							\fi
+						\fi
+					\fi
+				\fi
+			\fi
+		\fi
+	\fi
+}
+
+% Checks whether we need to create a separate 'tick scale label',
+% a node with ' * 10^3' on the side of the axis:
+%
+% PRECONDITION:
+%    Axis limits for #1 are given. I need their values before any data
+%    scale transformation has been applied.
+%    If
+%    	\pgfplots@#1min@unscaled@as@float
+%    and
+%    	\pgfplots@#1max@unscaled@as@float
+%    exist; I will use these macros.
+%    Otherwise, I will use \pgfplots@#1min and \pgfplots@#1max;
+%    assuming that no data scale transformation is active.
+%    FIXME : does that need further attention?
+\def\pgfplots@init@scaled@tick@for#1{%
+	\global\def\pgfplots@glob@TMPa{0}%
+	\expandafter\pgfplotslistcheckempty\csname pgfplots@prepared@tick@positions@major@#1\endcsname
+	\ifpgfplotslistempty
+		% we have no tick labels. Omit the tick scale label as well!
+	\else
+	\begingroup
+	\ifcase\csname pgfplots@scaled@ticks@#1@choice\endcsname\relax
+	% CASE 0 : scaled #1 ticks=false: do nothing here.
+	\or
+		% CASE 1 : scaled #1 ticks=true:
+		%--------------------------------
+		% the \pgfplots@xmin@unscaled@as@float  is set just before the data
+		% scale transformation is initialised.
+		%
+		% The variables are empty if there is no datascale transformation.
+		\expandafter\let\expandafter\pgfplots@cur@min@unscaled\csname pgfplots@#1min@unscaled@as@float\endcsname
+		\expandafter\let\expandafter\pgfplots@cur@max@unscaled\csname pgfplots@#1max@unscaled@as@float\endcsname
+		%
+		\ifx\pgfplots@cur@min@unscaled\pgfutil@empty
+			\edef\pgfplots@loc@TMPa{\csname pgfplots@#1min\endcsname}%
+			\expandafter\pgfmathfloatparsenumber\expandafter{\pgfplots@loc@TMPa}%
+			\let\pgfplots@cur@min@unscaled=\pgfmathresult
+			\edef\pgfplots@loc@TMPa{\csname pgfplots@#1max\endcsname}%
+			\expandafter\pgfmathfloatparsenumber\expandafter{\pgfplots@loc@TMPa}%
+			\let\pgfplots@cur@max@unscaled=\pgfmathresult
+		\fi
+		%
+		\expandafter\pgfmathfloat@decompose@E\pgfplots@cur@min@unscaled\relax\pgfmathfloat@a@E
+		\expandafter\pgfmathfloat@decompose@E\pgfplots@cur@max@unscaled\relax\pgfmathfloat@b@E
+		\pgfplots@init@scaled@tick@normalize@exponents
+		\ifnum\pgfmathfloat@b@E<\pgfmathfloat@a@E
+			\pgfmathfloat@b@E=\pgfmathfloat@a@E
+		\fi
+		\xdef\pgfplots@glob@TMPa{\pgfplots@scale@ticks@above@exponent}%
+		\ifnum\pgfplots@glob@TMPa<\pgfmathfloat@b@E
+			% ok, scale it:
+			\multiply\pgfmathfloat@b@E by-1
+			\xdef\pgfplots@glob@TMPa{\the\pgfmathfloat@b@E}%
+		\else
+			\xdef\pgfplots@glob@TMPa{\pgfplots@scale@ticks@below@exponent}%
+			\ifnum\pgfplots@glob@TMPa>\pgfmathfloat@b@E
+				% ok, scale it:
+				\multiply\pgfmathfloat@b@E by-1
+				\xdef\pgfplots@glob@TMPa{\the\pgfmathfloat@b@E}%
+			\else
+				% no scaling necessary:
+				\xdef\pgfplots@glob@TMPa{0}%
+			\fi
+		\fi
+	\or
+		% CASE 2 : scaled #1 ticks=base 10:
+		%--------------------------------
+		\c@pgf@counta=\csname pgfplots@scaled@ticks@#1@arg\endcsname\relax
+		%\multiply\c@pgf@counta by-1
+		\xdef\pgfplots@glob@TMPa{\the\c@pgf@counta}%
+	\or
+		% CASE 3 : scaled #1 ticks=real:
+		%--------------------------------
+		\pgfmathfloatparsenumber{\csname pgfplots@scaled@ticks@#1@arg\endcsname}%
+		\global\let\pgfplots@glob@TMPa=\pgfmathresult
+	\or
+		% CASE 4 : scaled #1 ticks=manual:
+		\expandafter\global\expandafter\let\expandafter\pgfplots@glob@TMPa\csname pgfplots@scaled@ticks@#1@arg\endcsname
+	\fi
+	\endgroup
+	\fi
+	\expandafter\let\csname pgfplots@tick@scale@#1\endcsname=\pgfplots@glob@TMPa%
+}
+
+% Handles the case that one of the limits is 0 (or unbounded, although
+% that might not be a use-case at all).
+%
+% Note that 0 = 0*10^A (naturally). In this case, A can be undefined,
+% and we want to use B's exponent (only) for the decision here.
+%
+% INPUT:
+%  - \pgfplots@cur@min@unscaled,
+%  - \pgfplots@cur@max@unscaled,
+%  - \pgfmathfloat@a@E
+%  - \pgfmathfloat@b@E
+%
+% Output:
+%  normalized values of \pgfmathfloat@a@E and \pgfmathfloat@b@E
+\def\pgfplots@init@scaled@tick@normalize@exponents{%
+	\pgfmathfloatgetflags\pgfplots@cur@min@unscaled\pgfmathfloat@a@S
+	\pgfmathfloatgetflags\pgfplots@cur@max@unscaled\pgfmathfloat@b@S
+	\ifcase\pgfmathfloat@a@S%
+		% min = 0.
+		\ifcase\pgfmathfloat@b@S
+			% max =0
+			% normalize to 0 * 10^0 !
+			\pgfmathfloat@a@E=0 %
+			\pgfmathfloat@b@E=0 %
+		\or
+			% max > 0
+			% since 0 = 0 * 10^A for any A, tick scaling is based on
+			% B only.
+			\pgfmathfloat@a@E=\pgfmathfloat@b@E
+		\or
+			% max < 0
+			\pgfmathfloat@a@E=\pgfmathfloat@b@E
+		\else
+			% max is unbounded. normalize exponent to something
+			% useful:
+			\pgfmathfloat@a@E=0 %
+			\pgfmathfloat@b@E=0 %
+		\fi
+	\or
+		% min>0
+		\ifcase\pgfmathfloat@b@S
+			% max =0
+			% since 0 = 0 * 10^B for any B, tick scaling is based on
+			% A only.
+			\pgfmathfloat@b@E=\pgfmathfloat@a@E
+		\or
+			% max > 0
+		\or
+			% max < 0
+		\else
+			% max is unbounded. normalize exponent to something
+			% useful:
+			\pgfmathfloat@b@E=\pgfmathfloat@a@E
+		\fi
+	\or
+		% min<0
+		\ifcase\pgfmathfloat@b@S
+			% max =0
+			% since 0 = 0 * 10^B for any B, tick scaling is based on
+			% A only.
+			\pgfmathfloat@b@E=\pgfmathfloat@a@E
+		\or
+			% max > 0
+		\or
+			% max < 0
+		\else
+			% max is unbounded. normalize exponent to something
+			% useful:
+			\pgfmathfloat@b@E=\pgfmathfloat@a@E
+		\fi
+	\else
+		% min is unbounded:
+		% normalize somehow.
+		\pgfmathfloat@a@E=0 %
+		\pgfmathfloat@b@E=0 %
+	\fi
+}
+
+% x-axis tick labels for #1th tick
+% #1: the axis (x,y or z)
+% #2: the value
+% #3,#4+#5: arguments for \pgfplotspointonorientedsurfaceabwithbshift
+% #5: ticknumber
+\def\pgfplots@show@ticklabel#1#2(#3,#4+#5)#6{%
+	\begingroup%
+	\edef\pgfmathresult{#2}%
+	\pgfkeyslet{/pgfplots/sloped/at position}\pgfmathresult%
+	\csname ifpgfplots@#1ticklabel@interval\endcsname
+		% Special case for the INTERVAL feature:
+		% we have to do additional work here.
+		\pgfmathparse{#3}%
+		\edef\pgfplots@show@ticklabel@coord@x@new{\pgfmathresult}%
+		\pgfmathparse{#4}%
+		\edef\pgfplots@show@ticklabel@coord@y@new{\pgfmathresult}%
+		%
+		\pgfplots@show@ticklabel@{#1}{#2}%
+		\let\nexttick=\tick
+		\ifx\pgfplots@show@ticklabel@LASTTICK\pgfutil@empty
+			% its the first call. Simply remember arguments and wait
+			% for interval boundary before proceeding.
+		\else
+			% acquire options of first interval boundary:
+			\pgfplots@show@ticklabel@LASTTICK
+			% compute new node position:
+			\pgfmathparse{0.5*(\pgfplots@show@ticklabel@coord@x + \pgfplots@show@ticklabel@coord@x@new)}%
+			\let\pgfplots@show@ticklabel@coord@x=\pgfmathresult%
+			\pgfmathparse{0.5*(\pgfplots@show@ticklabel@coord@y + \pgfplots@show@ticklabel@coord@y@new)}%
+			\let\pgfplots@show@ticklabel@coord@y=\pgfmathresult%
+			\let\ticknum=\pgfplots@show@ticklabel@num\relax%
+			\let\tick=\pgfplots@show@ticklabel@tick%
+			\pgfplots@show@ticklabel@@{#1}
+				{\pgfplotspointonorientedsurfaceabwithbshift{\pgfplots@show@ticklabel@coord@x}{\pgfplots@show@ticklabel@coord@y}{#5}}%
+		\fi
+		\xdef\pgfplots@show@ticklabel@LASTTICK{%
+			\noexpand\def\noexpand\pgfplots@show@ticklabel@tick{\nexttick}%
+			\noexpand\def\noexpand\pgfplots@show@ticklabel@coord@x{\pgfplots@show@ticklabel@coord@x@new}%
+			\noexpand\def\noexpand\pgfplots@show@ticklabel@coord@y{\pgfplots@show@ticklabel@coord@y@new}%
+			\noexpand\edef\noexpand\pgfplots@show@ticklabel@num{#6}%
+		}%
+	\else
+		\let\ticknum=#6\relax%
+		\pgfplots@show@ticklabel@{#1}{#2}%
+		\pgfplots@show@ticklabel@@{#1}{\pgfplotspointonorientedsurfaceabwithbshift{#3}{#4}{#5}}%
+	\fi
+	\endgroup
+}
+
+% #1: the axis (x,y or z)
+% #2: the location where to place it.
+\def\pgfplots@show@ticklabel@@#1#2{%
+	% Typeset the label!
+	\pgfinterruptboundingbox
+		% What makes this complicated is the 'ticklabel cs' feature.
+		% What we need is to compute the MAXIMUM LENGTH over each tick
+		% label IN DIRECTION OF THE OUTER NORMAL.
+		%
+		% This needs to
+		% 1. enable bounding box computation even in case of
+		% 'overlay',
+		% 2. projection of the bounding box in direction of the outer
+		% normal,
+		% 3. update of the bounding box if 'overlay' is not active.
+		\begingroup
+		%
+		% prepare step (1.):
+		\pgfkeysalso{%
+			/tikz/every node/.append code={%
+				\ifpgf@relevantforpicturesize
+					\gdef\pgfplots@show@ticklabel@@update@BB{1}%
+				\else
+					\gdef\pgfplots@show@ticklabel@@update@BB{0}%
+				\fi
+				\pgf@relevantforpicturesizetrue
+			}%
+		}%
+		%
+		% Compute and remember the position '#2':
+		\pgf@process{#2}%
+		\edef\pgfplots@ticklabel@at@x{\the\pgf@x}%
+		\edef\pgfplots@ticklabel@at@y{\the\pgf@y}%
+		%
+		% ok, generate the label!
+		\node at (\pgfplots@ticklabel@at@x,\pgfplots@ticklabel@at@y) {\csname pgfplots@#1ticklabel\endcsname};%
+		%
+		% compute the label's dimensions, step (2.):
+		\pgfplots@ticklabel@maxtickdimen@updateforcurrentpath
+			{#1}
+			{\pgf@x=\pgfplots@ticklabel@at@x\space\pgf@y=\pgfplots@ticklabel@at@y\space}%%
+		%
+		% prepare step (3.): update of bounding box:
+		\if\pgfplots@show@ticklabel@@update@BB1%
+			\xdef\pgfplots@glob@TMPa{%
+				\pgf@xa=\the\pgf@picminx\space
+				\pgf@xb=\the\pgf@picminy\space
+				\pgf@ya=\the\pgf@picmaxx\space
+				\pgf@yb=\the\pgf@picmaxy\space
+				\noexpand\pgf@protocolsizes{\pgf@xa}{\pgf@xb}%
+				\noexpand\pgf@protocolsizes{\pgf@ya}{\pgf@yb}%
+				\noexpand\pgf@resetpathsizes
+			}%
+		\else
+			\global\let\pgfplots@glob@TMPa=\relax
+		\fi
+		\endgroup
+	\endpgfinterruptboundingbox
+	\begingroup
+		\pgfplots@glob@TMPa
+	\endgroup
+}%
+
+% TICK LABEL DIMENSION CONTROL
+%
+% The following framework is supposed to accumulate the value for
+% \pgfplotsvalueoflargesttickdimen.
+%
+% It works like this:
+%
+% \pgfplots@ticklabel@maxtickdimen@reset{#1}
+% \pgfplots@ticklabel@maxtickdimen@prepare@for@normalvec{#1}{<normal vector>}
+% ...
+% \pgfplots@ticklabel@maxtickdimen@updateforcurrentpath{#1}
+% \pgfplots@ticklabel@maxtickdimen@updateforcurrentpath{#1}
+% \pgfplots@ticklabel@maxtickdimen@updateforcurrentpath{#1}
+% ...
+% \pgfplots@ticklabel@maxtickdimen@finish{#1}
+%
+% -> then, \pgfplotsvalueoflargesttickdimen expands to the largest
+%  distance from a tick's coordinate to its tick label bounding box in
+%  direction of the outer normal vector.
+\def\pgfplots@ticklabel@maxtickdimen@reset#1{%
+	\expandafter\gdef\csname pgfplots@maxtickdimen@#1\endcsname{0pt}%
+	\expandafter\gdef\csname pgfplots@maxtickdimen@extrashift@#1\endcsname{0pt}%
+}%
+
+% Adds the extra shift '#2' along the normal vector for axis '#1'.
+\def\pgfplots@ticklabel@maxtickdimen@extrashift#1#2{%
+	\begingroup
+	\afterassignment\pgfplots@gobble@until@relax
+	\pgf@xa=#2pt\relax
+	\advance\pgf@xa by\csname pgfplots@maxtickdimen@extrashift@#1\endcsname\relax
+	\expandafter\xdef\csname pgfplots@maxtickdimen@extrashift@#1\endcsname{\the\pgf@xa}%
+	\endgroup
+}
+
+% Finalizes the maxtickdimen computation.
+%
+% This applies any extra shifts (including '#1ticklabel shift').
+\def\pgfplots@ticklabel@maxtickdimen@finish#1{%
+	\pgfkeysgetvalue{/pgfplots/#1ticklabel shift}\pgfmathresult
+	\ifx\pgfmathresult\pgfutil@empty
+	\else
+		\pgfplots@ticklabel@maxtickdimen@extrashift{#1}{\pgfkeysvalueof{/pgfplots/#1ticklabel shift}}%
+	\fi
+	\begingroup
+	\pgf@xa=\csname pgfplots@maxtickdimen@extrashift@#1\endcsname\relax
+	\advance\pgf@xa by \csname pgfplots@maxtickdimen@#1\endcsname\relax
+	\expandafter\xdef\csname pgfplots@maxtickdimen@#1\endcsname{\the\pgf@xa}%
+	\endgroup
+}%
+%
+% #1: the axis (x,y or z)
+% #2: the normal vector
+\def\pgfplots@ticklabel@maxtickdimen@prepare@for@normalvec#1#2{%
+	\pgf@process{#2}%
+	\edef\pgfplots@loc@vector@to@outside{\pgf@x=\the\pgf@x\space\pgf@y=\the\pgf@y\space}%
+	% Identify the corner point of the ticklabel bounding box which
+	% shall be used:
+	\ifdim\pgf@x>0sp
+		\ifdim\pgf@y>0sp
+			% NORTH EAST
+			\def\pgfplots@loc@ticklabel@bb@corner{\pgf@x=\pgf@picmaxx\pgf@y=\pgf@picmaxy}%
+		\else
+			% SOUTH EAST
+			\def\pgfplots@loc@ticklabel@bb@corner{\pgf@x=\pgf@picmaxx\pgf@y=\pgf@picminy}%
+		\fi
+	\else
+		\ifdim\pgf@y>0sp
+			% NORTH WEST
+			\def\pgfplots@loc@ticklabel@bb@corner{\pgf@x=\pgf@picminx\pgf@y=\pgf@picmaxy}%
+		\else
+			% SOUTH WEST
+			\def\pgfplots@loc@ticklabel@bb@corner{\pgf@x=\pgf@picminx\pgf@y=\pgf@picminy}%
+		\fi
+	\fi
+}%
+
+% #1: the axis (x,y or z)
+% #2: the point from where dimension shall be computed (the 'at'
+% argument of the tick label)
+\def\pgfplots@ticklabel@maxtickdimen@updateforcurrentpath#1#2{%
+	\pgfplotsscalarproductofvectors
+		{\pgfplots@loc@vector@to@outside}%
+		{\pgfpointdiff
+			{#2\pgf@pos@transform\pgf@x\pgf@y}%
+			{\pgfplots@loc@ticklabel@bb@corner}}%
+	\ifdim\pgf@x>\csname pgfplots@maxtickdimen@#1\endcsname\relax
+		\expandafter\xdef\csname pgfplots@maxtickdimen@#1\endcsname{\the\pgf@x}%
+	\fi
+}%
+
+% Expands to the largest distance of a tick position to its tick label
+% bounding box in direction of the outer unit normal vector.
+%
+% It does also include the value of the 'ticklabel shift' key.
+%
+% This function assumes that
+% \pgfplots@ticklabel@maxtickdimen@reset{#1}
+% \pgfplots@ticklabel@maxtickdimen@prepare@for@normalvec{#1}{<normal vector>}
+% ...
+% \pgfplots@ticklabel@maxtickdimen@updateforcurrentpath{#1}
+% \pgfplots@ticklabel@maxtickdimen@updateforcurrentpath{#1}
+% \pgfplots@ticklabel@maxtickdimen@updateforcurrentpath{#1}
+% ...
+% \pgfplots@ticklabel@maxtickdimen@finish{#1}
+% has been invoked completely.
+%
+% #1: either x,y or z. It denotes the axis for which the ticks are
+% requested.
+\def\pgfplotsvalueoflargesttickdimen#1{%
+	\csname pgfplots@maxtickdimen@#1\endcsname
+}%
+
+% Just like \pgfplotsqpointoutsideofaxis, but this one here uses the
+% axis on which tick labels will be drawn.
+%
+% #1: one of x, y or z.
+% #2: the coordinate on the tick axis designated by '#1'.
+% #3: a scale (a dimen) in which the point is moved in direction of
+% the outward normal vector of the axis.
+%
+% @see \pgfplotsqpointoutsideofaxis
+\def\pgfplotsqpointoutsideofticklabelaxis#1#2#3{%
+	\pgfplotsqpointoutsideofaxis{\csname pgfplots@#1ticklabelaxisspec\endcsname}{#2}{#3}%
+}%
+\def\pgfplotsqpointoutsideofticklabelaxisrel#1#2#3{%
+%\message{using \string\pgfplotsqpointoutsideofaxisrel{\csname pgfplots@#1ticklabelaxisspec\endcsname}{#2}{#3}^^J}%
+	\pgfplotsqpointoutsideofaxisrel{\csname pgfplots@#1ticklabelaxisspec\endcsname}{#2}{#3}%
+}%
+\def\pgfplotsqpointoutsideofticklabelaxistransformed#1#2#3{%
+	\pgfplotsqpointoutsideofaxistransformed{\csname pgfplots@#1ticklabelaxisspec\endcsname}{#2}{#3}%
+}%
+
+% Expands to the three-character-identification for the axis
+% containing tick labels for axis #1.
+%
+% #1: either x, y or z.
+\def\pgfplotsticklabelaxisspec#1{\csname pgfplots@#1ticklabelaxisspec\endcsname}%
+
+% The unit outer normal vector for axis #1.
+% #1: one of x, y or z.
+\def\pgfplotspointouternormalvectorofticklabelaxis#1{%
+	\pgfplotspointouternormalvectorofaxis{\csname pgfplots@#1ticklabelaxisspec\endcsname}%
+}
+
+% Defines \tick by applying any necessary math to the (possibly
+% transformed) tick value #2.
+%
+% #1: axis (x or y)
+% #2: tick value.
+\def\pgfplots@show@ticklabel@#1#2{%
+	\csname ifpgfplots@apply@datatrafo@#1\endcsname
+		\pgfplotscoordmath{#1}{datascaletrafo inverse}{#2}%
+		\ifcase\csname pgfplots@scaled@ticks@#1@choice\endcsname
+		\or
+			\expandafter\pgfmathfloatshift@\expandafter{\pgfmathresult}{\csname pgfplots@tick@scale@#1\endcsname}%
+		\or
+			\expandafter\pgfmathfloatshift@\expandafter{\pgfmathresult}{\csname pgfplots@tick@scale@#1\endcsname}%
+		\or
+			\expandafter\pgfmathfloatdivide@\expandafter{\pgfmathresult}{\csname pgfplots@tick@scale@#1\endcsname}%
+		\or
+			% scaled #1 ticks=manual. Invoke manual tick scaling code:
+			\expandafter\let\expandafter\pgfplots@loc@TMPa\csname pgfplots@tick@scale@#1\endcsname
+			\begingroup
+			\pgfkeys{/pgf/fpu=true}%
+			\expandafter\pgfplots@loc@TMPa\expandafter{\pgfmathresult}%
+			\pgfmath@smuggleone\pgfmathresult
+			\endgroup
+			\pgfmathfloatparsenumber\pgfmathresult%
+		\fi
+		% .. and this here provides \tick as fixed point repr:
+		\expandafter\pgfmathfloattofixed\expandafter{\pgfmathresult}%
+		\let\tick=\pgfmathresult
+	\else
+		\edef\tick{#2}%
+		\pgfplots@if{pgfplots@#1islinear}{%
+			\ifnum\csname pgfplots@scaled@ticks@#1@choice\endcsname=0
+			\else
+				\pgfmathfloatparsenumber{#2}%
+				\ifnum\csname pgfplots@scaled@ticks@#1@choice\endcsname=3
+					\expandafter\pgfmathfloatdivide@\expandafter{\pgfmathresult}{\csname pgfplots@tick@scale@#1\endcsname}%
+				\else
+					\expandafter\pgfmathfloatshift@\expandafter{\pgfmathresult}{\csname pgfplots@tick@scale@#1\endcsname}%
+				\fi
+				\expandafter\pgfmathfloattofixed\expandafter{\pgfmathresult}%
+				\let\tick=\pgfmathresult
+			\fi
+		}{}%
+	\fi
+	\pgfplots@coord@inv@trafo@apply{#1}{\tick}%
+	\let\tick=\pgfmathresult
+}%
+
+\def\pgfplots@user@ticklabel@list@x{%
+	\pgfplotslistselectorempty\ticknum\of\pgfplots@xticklabels\to\tick
+	\pgfplots@ticklabel@typeset@arg\tick
+}
+\def\pgfplots@user@ticklabel@list@y{%
+	\pgfplotslistselectorempty\ticknum\of\pgfplots@yticklabels\to\tick
+	\pgfplots@ticklabel@typeset@arg\tick
+}
+\def\pgfplots@user@ticklabel@list@z{%
+	\pgfplotslistselectorempty\ticknum\of\pgfplots@zticklabels\to\tick
+	\pgfplots@ticklabel@typeset@arg\tick
+}
+\def\pgfplots@user@extra@ticklabel@list@x{%
+	\pgfplotslistselectorempty\ticknum\of\pgfplots@extra@xticklabels\to\tick
+	\pgfplots@ticklabel@typeset@arg\tick
+}
+\def\pgfplots@user@extra@ticklabel@list@y{%
+	\pgfplotslistselectorempty\ticknum\of\pgfplots@extra@yticklabels\to\tick
+	\pgfplots@ticklabel@typeset@arg\tick
+}
+\def\pgfplots@user@extra@ticklabel@list@z{%
+	\pgfplotslistselectorempty\ticknum\of\pgfplots@extra@zticklabels\to\tick
+	\pgfplots@ticklabel@typeset@arg\tick
+}
+
+\def\pgfplots@limit@max@reg#1{%
+	\if\pgfkeysvalueof{/pgfplots/#1 dir/value}n%
+		\csname pgfplots@#1max@reg\endcsname
+	\else
+		\csname pgfplots@#1min@reg\endcsname
+	\fi
+}%
+\def\pgfplots@limit@min@reg#1{%
+	\if\pgfkeysvalueof{/pgfplots/#1 dir/value}n%
+		\csname pgfplots@#1min@reg\endcsname
+	\else
+		\csname pgfplots@#1max@reg\endcsname
+	\fi
+}%
+
+% Check if a label does not cross the x-axis
+\def\pgfplots@ytick@check@tickshow{%
+	\pgfplots@tickshowtrue
+	%
+	\pgfplots@if@has@axis@shift x{%
+		\pgfplots@hide@obscured@ytickfalse
+	}{}%
+	%
+	\ifpgfplots@hide@obscured@ytick
+		\if\pgfplots@yaxislinesnum2% center
+			\ifcase\pgfplots@xaxislinesnum\relax
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@min@reg{y}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@max@reg{y}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\or
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@min@reg{y}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{%
+					}%
+			\or
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@logical@ZERO@y pt}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\or
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@max@reg{y}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\fi
+		\fi
+	\fi
+}
+\def\pgfplots@ztick@check@tickshow{%
+	\pgfplots@tickshowtrue
+	%
+	\pgfplots@if@has@axis@shift x{%
+		\pgfplots@if@has@axis@shift y{%
+			\pgfplots@hide@obscured@ztickfalse
+		}{}%
+	}{}%
+	%
+	\ifpgfplots@hide@obscured@ztick
+		\if\pgfplots@zaxislinesnum2% center
+			\pgfplotsmath@ifapproxequal@dim
+				{\pgfplots@tmpa}{\pgfplots@logical@ZERO@z pt}
+				{\pgfplots@loc@tick@placement@tolerance}
+				{%
+					\pgfplots@tickshowfalse
+				}{}%
+		\fi
+	\fi
+}%
+
+% Fills the macros
+%   \pgfplots@tick@LOWER@b \pgfplots@tick@end@a
+%   \pgfplots@tick@UPPER@b \pgfplots@tick@end@b
+% with coordinates such that
+%   (\pgfplots@tick@LOWER@b,\pgfplots@tmpa) -- (\pgfplots@tick@end@a,\pgfplots@tmpa)
+% produces a correct tick line.
+%
+% The '@b' variant is only used in case of \pgfplots@ytickposnum = 0
+%
+% #1 : the current axis (x or y).
+% #2 : the current tick width
+%
+\def\pgfplots@prepare@tick@offsets@for@#1#2{%
+	%
+	% FIXME : this stuff was ok for 2D.
+	% For 3D, it works only for the cases of boxed axes or centered
+	% axis lines.
+	\ifcase\csname pgfplots@#1tickposnum\endcsname\relax
+		% both
+		%(\pgfplots@xcoordminTEX-\pgfplots@tick@offset, \pgfplots@tmpa) -- ++( #2, 0pt)
+		\edef\pgfplots@tick@LOWER@b{\csname pgfplots@\pgfplotspointonorientedsurfaceB min\endcsname}%
+		%
+		%(\pgfplots@xcoordmaxTEX+\pgfplots@tick@offset, \pgfplots@tmpa)	-- ++(-#2, 0pt)
+		\edef\pgfplots@tick@UPPER@b{\csname pgfplots@\pgfplotspointonorientedsurfaceB max\endcsname}%
+	\or
+		% left
+		% (\pgfplots@xcoordminTEX-\pgfplots@tick@offset, \pgfplots@tmpa) -- ++( #2, 0pt);
+		\edef\pgfplots@tick@LOWER@b{\csname pgfplots@\pgfplotspointonorientedsurfaceB min\endcsname}%
+	\or
+		% center
+		% (\pgfplots@ZERO@x      -\pgfplots@tick@offset, \pgfplots@tmpa) -- ++( #2, 0pt);
+		\edef\pgfplots@tick@LOWER@b{\csname pgfplots@logical@ZERO@\pgfplotspointonorientedsurfaceB \endcsname}%
+	\or
+		% right
+		% (\pgfplots@xcoordmaxTEX+\pgfplots@tick@offset, \pgfplots@tmpa)	-- ++(-#2, 0pt);
+		\edef\pgfplots@tick@UPPER@b{\csname pgfplots@\pgfplotspointonorientedsurfaceB max\endcsname}%
+	\else
+		% FALL BACK. never used, I guess?
+		% (\pgfplots@xcoordminTEX-\pgfplots@tick@offset, \pgfplots@tmpa) -- ++( #2, 0pt);
+		\edef\pgfplots@tick@LOWER@b{\csname pgfplots@\pgfplotspointonorientedsurfaceB min\endcsname}%
+	\fi
+	%
+	\ifcase\csname pgfplots@#1tickalignnum\endcsname\relax
+		\def\pgfplots@tick@offset{0}%
+	\or
+		\edef\pgfplots@tick@offset{#2}%
+	\or
+		\pgfmathmultiply@{0.5}{#2}%
+		\let\pgfplots@tick@offset=\pgfmathresult%
+	\fi
+	%
+	\edef\pgfplots@tick@LOWER@shiftbeg{-\pgfplots@tick@offset}%
+	\pgfmathadd@{\pgfplots@tick@LOWER@shiftbeg}{#2}%
+	\let\pgfplots@tick@LOWER@shiftend=\pgfmathresult
+	%
+	\edef\pgfplots@tick@UPPER@shiftbeg{\pgfplots@tick@offset}%
+	\pgfmathsubtract@{\pgfplots@tick@UPPER@shiftbeg}{#2}%
+	\let\pgfplots@tick@UPPER@shiftend=\pgfmathresult
+	%
+	% Assemble the \pgfplots@drawticklines@for@placecomputedtick
+	% command.
+	\def\pgfplots@drawticklines@for@placecomputedtick@LOWER{%
+		\pgfpathmoveto{\pgfplotspointonorientedsurfaceabwithbshift{\pgfplots@curtickpos}{\pgfplots@tick@LOWER@b}{\pgfplots@tick@LOWER@shiftbeg pt}}%
+		\pgfpathlineto{\pgfplotspointonorientedsurfaceabwithbshift{\pgfplots@curtickpos}{\pgfplots@tick@LOWER@b}{\pgfplots@tick@LOWER@shiftend pt}}%
+	}%
+	\def\pgfplots@drawticklines@for@placecomputedtick@UPPER{%
+		\pgfpathmoveto{\pgfplotspointonorientedsurfaceabwithbshift{\pgfplots@curtickpos}{\pgfplots@tick@UPPER@b}{\pgfplots@tick@UPPER@shiftbeg pt}}%
+		\pgfpathlineto{\pgfplotspointonorientedsurfaceabwithbshift{\pgfplots@curtickpos}{\pgfplots@tick@UPPER@b}{\pgfplots@tick@UPPER@shiftend pt}}%
+	}%
+%\message{place computed tick: LOWEROK=\pgfplots@drawticklines@for@placecomputedtick@LOWEROK; UPPEROK=\pgfplots@drawticklines@for@placecomputedtick@UPPEROK.}%
+}%
+
+
+\newif\ifpgfplots@needsminorloop
+
+\def\pgfplots@draw@tick@scale@label@for#1{%
+	\csname ifpgfplots@#1islinear\endcsname
+		\begingroup
+			\def\pgfplots@temp@isbaseten{0}%
+			\ifcase\csname pgfplots@scaled@ticks@#1@choice\endcsname
+				\global\let\pgfplots@glob@TMPa=\pgfutil@empty
+			\or
+				\xdef\pgfplots@glob@TMPa{\csname pgfplots@tick@scale@#1\endcsname}%
+				\def\pgfplots@temp@isbaseten{1}%
+			\or
+				\xdef\pgfplots@glob@TMPa{\csname pgfplots@tick@scale@#1\endcsname}%
+				\def\pgfplots@temp@isbaseten{1}%
+			\or
+				\xdef\pgfplots@glob@TMPa{\csname pgfplots@tick@scale@#1\endcsname}%
+				% real:
+			\or
+				% manual:
+				\global\def\pgfplots@glob@TMPa{dummyargument}%
+			\fi
+			\if1\pgfplots@temp@isbaseten
+				\expandafter\c@pgf@counta\pgfplots@glob@TMPa\relax
+				\multiply\c@pgf@counta by-1
+				\ifnum\c@pgf@counta=0\relax
+					\global\let\pgfplots@glob@TMPa=\pgfutil@empty
+				\else
+					\xdef\pgfplots@glob@TMPa{\the\c@pgf@counta}%
+				\fi
+			\fi
+		\endgroup
+		\ifx\pgfplots@glob@TMPa\pgfutil@empty
+		\else
+			\begingroup
+			\pgfkeysgetvalue{/pgfplots/#1tick scale label code/.@cmd}\pgfplots@loc@TMPa
+			\ifx\pgfplots@loc@TMPa\pgfplots@empty@command@key
+			\else
+				\edef\pgfplots@tick@scale@labels{\noexpand\pgfplots@invoke@pgfkeyscode{/pgfplots/#1tick scale label code/.@cmd}{\pgfplots@glob@TMPa}}%
+				%
+				\pgfplots@change@pgfpoints@to@descriptioncs
+				%
+				\node[%
+					/pgfplots/every tick label,%
+					/pgfplots/every #1 tick label,%
+					/pgfplots/every #1 tick scale label]
+				{\pgfplots@tick@scale@labels};
+			\fi
+			\endgroup
+		\fi
+	\fi
+}
+
+% Check if the current tick position, stored in \pgfplots@tmpa,
+% does not cross the y-axis.
+%
+% This is just a special case for centered axis lines.
+\def\pgfplots@xtick@check@tickshow{%
+	\pgfplots@tickshowtrue
+	%
+	\pgfplots@if@has@axis@shift y{%
+		\pgfplots@hide@obscured@xtickfalse
+	}{}%
+	%
+	\ifpgfplots@hide@obscured@xtick
+		\if\pgfplots@xaxislinesnum2% center
+			\ifcase\pgfplots@yaxislinesnum\relax
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@min@reg{x}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@max@reg{x}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\or
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@min@reg{x}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\or
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@logical@ZERO@x pt}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\or
+				\pgfplotsmath@ifapproxequal@dim
+					{\pgfplots@tmpa}{\pgfplots@limit@max@reg{x}}
+					{\pgfplots@loc@tick@placement@tolerance}
+					{%
+						\pgfplots@tickshowfalse
+					}{}%
+			\fi
+		\fi
+	\fi
+}
+
+% Draws extra ticks including grid lines, tick lines and tick labels
+% along the current oriented surface.
+%
+% See \pgfplots@drawticklines@onorientedsurf@ for a description of the
+% oriented surface.
+%
+% #1 : tick position list
+\def\pgfplots@draw@extra@ticks@onorientedsurf{%
+	\expandafter\pgfplots@draw@extra@ticks@onorientedsurf@\pgfplotspointonorientedsurfaceA
+}%
+% #1: axis (x or y)
+% #2: tick position list
+\def\pgfplots@draw@extra@ticks@onorientedsurf@#1#2{%
+	\begingroup
+	\def\pgfplots@scaled@ticks@x@choice{0}%
+	\def\pgfplots@scaled@ticks@y@choice{0}%
+	\def\pgfplots@scaled@ticks@z@choice{0}%
+	\csname pgfplots@#1minorticksfalse\endcsname
+	\csname pgfplots@#1minorgridsfalse\endcsname
+	\expandafter\let\expandafter\pgfplots@ticklabel@pos@orig\csname pgfplots@#1ticklabel@pos\endcsname%
+	\expandafter\let\expandafter\axis@TMP\csname pgfplots@extra@#1ticklabel\endcsname
+	\expandafter\let\csname pgfplots@#1ticklabel\endcsname=\axis@TMP
+%	\pgfplotsset{/pgfplots/every extra #1 tick}%
+	% use a scope here such that line width and draw color can be set.
+	\scope[/pgfplots/.cd,/pgfplots/every extra #1 tick]
+	\expandafter\let\expandafter\pgfplots@ticklabel@pos@\csname pgfplots@#1ticklabel@pos\endcsname
+	\ifx\pgfplots@ticklabel@pos@\pgfplots@ticklabel@pos@orig
+	\else
+		\pgfplots@init@ticklabelaxisspec
+	\fi
+	\pgfplots@prepare@tick@coordlists@for{#1}{#2}%
+	\pgfplots@drawgridlines@onorientedsurf%
+	\pgfplots@drawticklines@onorientedsurf%
+	\pgfplots@drawticklabels@onorientedsurf%
+	\endscope
+	\endgroup
+}
+
+% Computes final major and minor tick positions into global lists
+% \pgfplots@prepared@tick@positions@major@x
+% and
+% \pgfplots@prepared@tick@positions@minor@x.
+%
+% Both lists contain entries of the form {<index>}{<logical position>}
+% @see \pgfplots@prepared@tick@pos@unpack
+%
+% #1: the axis
+% #2: the tick list.
+%
+% PRECONDITION:
+% 	- \pgfplots@determinedefaultvalues has been executed.
+% 		That means particularly that \pgfplots@[xy][min,max] are available in TeX point
+% 		range (after datascaling and logs).
+% POSTCONDITION:
+% 	- the lists
+% 		\pgfplots@prepared@tick@positions@major@x
+% 		\pgfplots@prepared@tick@positions@major@tickindices@x
+% 		\pgfplots@prepared@tick@positions@minor@x
+% 	  are ready.
+\def\pgfplots@prepare@tick@coordlists@for#1#2{%
+	\begingroup
+	\expandafter\let\expandafter\ifpgfplots@islinear\csname ifpgfplots@#1islinear\endcsname
+	\expandafter\let\expandafter\ifpgfplots@minorticks\csname ifpgfplots@#1minorticks\endcsname
+	\expandafter\let\expandafter\ifpgfplots@minorgrids\csname ifpgfplots@#1minorgrids\endcsname
+	% these lists need to be global such that I can fill them inside
+	% of \foreach statements. And, yes: I have also added a TeX group
+	% on my own (but that's not the problem).
+	\global\pgfplotslistnewempty\pgfplots@prepared@tick@positions@major
+	%\global\pgfplotslistnewempty\pgfplots@prepared@tick@positions@major@tickindices
+	\global\pgfplotslistnewempty\pgfplots@prepared@tick@positions@minor
+	%
+	\pgfplots@prepare@tick@coordlists@for@handletolerance#1%
+	%
+	\edef\pgfplots@loc@TMPa{#2}%
+	\ifx\pgfplots@loc@TMPa\pgfutil@empty
+	\else
+		\ifpgfplots@minorticks
+			\pgfplots@needsminorlooptrue
+		\else
+			\ifpgfplots@minorgrids
+				\pgfplots@needsminorlooptrue
+			\else
+				\pgfplots@needsminorloopfalse
+			\fi
+		\fi
+		\pgfkeysgetvalue{/pgfplots/minor #1tick}\pgfplots@minor@tick@list
+		\ifx\pgfplots@minor@tick@list\pgfutil@empty
+		\else
+			\pgfplots@needsminorloopfalse
+		\fi
+		%
+		\ifpgfplots@needsminorloop
+			\ifpgfplots@islinear
+				\pgfkeysgetvalue{/pgfplots/minor #1 tick num}\pgfplots@minor@tick@num
+				\begingroup
+				\c@pgf@counta=\pgfplots@minor@tick@num\relax
+				\advance\c@pgf@counta by1\relax
+				\pgfplots@tmpa=\csname pgfplots@tick@distance@#1\endcsname pt %
+				\divide\pgfplots@tmpa by\c@pgf@counta
+				\edef\pgfmathresult{\pgf@sys@tonumber{\pgfplots@tmpa}}%
+				\pgfmath@smuggleone\pgfmathresult
+				\endgroup
+				\let\pgfplots@minor@tick@dist=\pgfmathresult
+			\else
+				\def\pgfplots@minor@tick@num{9}%
+			\fi
+		\fi
+		%
+		% Prepare the [xy]tick[min|max] key processing:
+		\let\pgfplots@checktickminmax=\pgfutil@empty
+		\expandafter\ifx\csname pgfplots@#1tickmin\endcsname\pgfutil@empty
+		\else
+			\expandafter\def\expandafter\pgfplots@checktickminmax\expandafter{%
+				\pgfplots@checktickminmax
+				\pgfplots@prepare@tick@coordlists@for@checktickmin#1%
+			}%
+		\fi
+		\expandafter\ifx\csname pgfplots@#1tickmax\endcsname\pgfutil@empty
+		\else
+			\expandafter\def\expandafter\pgfplots@checktickminmax\expandafter{%
+				\pgfplots@checktickminmax
+				\pgfplots@prepare@tick@coordlists@for@checktickmax#1%
+			}%
+		\fi
+		%
+		% compute the last index because it won't receive minor ticks:
+		\gdef\pgfplots@glob@TMPd{-1}%
+		\foreach \x in {#2} {%
+			\pgfplotsutil@advancestringcounter@global\pgfplots@glob@TMPd
+		}%
+		\let\c@pgfplots@ticknum@last=\pgfplots@glob@TMPd
+		%
+		\gdef\pgfplots@glob@TMPa{0}%
+		\foreach \x in {#2} {%
+			\let\pgfplots@ticknum=\pgfplots@glob@TMPa
+			%
+			\pgfplots@prepare@tick@coordlists@for@assign\pgfplots@tmpa=\x
+			\csname pgfplots@#1tick@check@tickshow\endcsname
+			\pgfplots@prepare@tick@coordlists@for@checkdatalimits#1%
+			\pgfplots@checktickminmax
+			%
+			\ifpgfplots@tickshow
+				\edef\x{{\pgfplots@ticknum}{\x}}%
+				\expandafter\pgfplotslistpushbackglobal\x\to\pgfplots@prepared@tick@positions@major
+				%\expandafter\pgfplotslistpushbackglobal\pgfplots@ticknum\to\pgfplots@prepared@tick@positions@major@tickindices
+			\fi
+			% X-Axis ticks bottom and top
+			\ifpgfplots@needsminorloop
+				\ifnum\pgfplots@ticknum=\c@pgfplots@ticknum@last\relax
+				\else
+					% SEE BELOW for the 'minor #1tick' feature -- it has a
+					% separate loop.
+					\foreach \pgfplots@i in {1,...,\pgfplots@minor@tick@num} {%
+						% XXX : bug 165 [minor ticks] minor tick drawn after the last xtick 
+						% seems to be here. Does that harm?
+						\ifpgfplots@islinear
+							\pgfmathmultiply@{\pgfplots@i}{\pgfplots@minor@tick@dist}%
+						\else
+							% in log:
+							%  log( i*10^k ) = log\pgfplots@i + k\log10 -> draw ticks for i=1..9
+							\pgfplotscoordmath{#1}{log unsigned int}{\pgfplots@i}%
+							\pgfplotscoordmath{#1}{tofixed}{\pgfmathresult}%
+						\fi
+						\pgfplots@prepare@tick@coordlists@for@advance\pgfplots@tmpa by\pgfmathresult
+						\pgfplots@tickshowtrue
+						\pgfplots@prepare@tick@coordlists@for@checkdatalimits#1%
+						\pgfplots@checktickminmax
+						\ifpgfplots@tickshow
+							\pgfplots@prepare@tick@coordlists@for@tofixed\pgfplots@tmpa
+							\c@pgf@counta=\pgfplots@ticknum\relax
+							\advance\c@pgf@counta by\pgfplots@i\relax
+							\edef\x{{\the\c@pgf@counta}{\pgfmathresult}}%
+							\expandafter\pgfplotslistpushbackglobal\x\to\pgfplots@prepared@tick@positions@minor
+						\fi
+					}%
+				\fi
+			\fi
+			\pgfplotsutil@advancestringcounter\pgfplots@ticknum
+			% carry \ticknum outside of this scope:
+			\global\let\pgfplots@glob@TMPa=\pgfplots@ticknum
+		}%
+		%
+		\ifx\pgfplots@minor@tick@list\pgfutil@empty
+		\else
+			% handle the 'minor #1tick' feature:
+			\def\pgfplots@loc@TMPa{\foreach \x in }%
+			\expandafter\pgfplots@loc@TMPa\expandafter{\pgfplots@minor@tick@list} {%
+				\let\pgfplots@ticknum=\pgfplots@glob@TMPa
+				%
+				\pgfplots@prepare@tick@coordlists@for@assign\pgfplots@tmpa=\x
+				\pgfplots@tickshowtrue
+				\pgfplots@prepare@tick@coordlists@for@checkdatalimits#1%
+				\pgfplots@checktickminmax
+				\ifpgfplots@tickshow
+					\edef\x{{\pgfplots@ticknum}{\x}}%
+					\expandafter\pgfplotslistpushbackglobal\x\to\pgfplots@prepared@tick@positions@minor
+				\fi
+				\pgfplotsutil@advancestringcounter\pgfplots@ticknum
+				% carry \ticknum outside of this scope:
+				\global\let\pgfplots@glob@TMPa=\pgfplots@ticknum
+			}%
+		\fi
+	\fi
+	\endgroup
+	\expandafter\let\csname pgfplots@prepared@tick@positions@minor@#1\endcsname=\pgfplots@prepared@tick@positions@minor
+	\expandafter\let\csname pgfplots@prepared@tick@positions@major@#1\endcsname=\pgfplots@prepared@tick@positions@major
+	%\expandafter\let\csname pgfplots@prepared@tick@positions@major@tickindices@#1\endcsname=\pgfplots@prepared@tick@positions@major@tickindices
+	\global\let\pgfplots@prepared@tick@positions@major=\relax
+	%\global\let\pgfplots@prepared@tick@positions@major@tickindices=\relax
+	\global\let\pgfplots@prepared@tick@positions@minor=\relax
+}%
+
+% The following set of macros can be used to replace the TeX register
+% arithmetics used to speed up the computations inside of
+% \pgfplots@prepare@tick@coordlists@for by something different.
+%
+% FIXME : that is no clean programming! Perhaps a new math class
+% should be used here, and implemented during the complete tick
+% preparation!?
+%
+% Keep in mind that these tick positions are in "transformed
+% range", i.e. they are expected to be in the range -16000...16000. At
+% the time of this modification, only the smithchart lib needs special
+% handling here... does this justify a re-design?
+\def\pgfplots@prepare@tick@coordlists@for@assign#1=#2{%
+	#1=#2pt
+}%
+\def\pgfplots@prepare@tick@coordlists@for@advance#1by#2{%
+	\advance#1 by#2 pt %
+}%
+\def\pgfplots@prepare@tick@coordlists@for@tofixed#1{%
+	\edef\pgfmathresult{\pgf@sys@tonumber{#1}}%
+}%
+\def\pgfplots@prepare@tick@coordlists@for@handletolerance#1{%
+	\afterassignment\pgfplots@gobble@until@relax
+	\pgfplots@tmpa=\pgfkeysvalueof{/pgfplots/#1tick placement tolerance}pt\relax
+	\pgfplots@tmpa=\csname pgfplots@#1@inverseveclength\endcsname\pgfplots@tmpa
+	\edef\pgfplots@loc@tick@placement@tolerance{\the\pgfplots@tmpa}%
+	%
+	\advance\csname pgfplots@#1min@reg\endcsname by-\pgfplots@tmpa
+	\advance\csname pgfplots@#1max@reg\endcsname by\pgfplots@tmpa
+}%
+\def\pgfplots@prepare@tick@coordlists@for@checktickmin#1{%
+	\ifdim\pgfplots@tmpa<\csname pgfplots@#1tickmin\endcsname pt
+		\pgfplots@tickshowfalse
+	\fi
+}%
+\def\pgfplots@prepare@tick@coordlists@for@checktickmax#1{%
+	\ifdim\pgfplots@tmpa>\csname pgfplots@#1tickmax\endcsname pt
+		\pgfplots@tickshowfalse
+	\fi
+}%
+\def\pgfplots@prepare@tick@coordlists@for@checkdatalimits#1{%
+	\ifdim\pgfplots@tmpa<\csname pgfplots@#1min@reg\endcsname
+		\pgfplots@tickshowfalse
+	\else
+		\ifdim\pgfplots@tmpa>\csname pgfplots@#1max@reg\endcsname
+			\pgfplots@tickshowfalse
+		\fi
+	\fi
+}%
+
+% Unpacks entries of the \pgfplots@prepared@tick@positions@* lists.
+%
+% Defines \pgfplots@tick to be the actual tick position and
+% \pgfplots@ticknum to be its index.
+%
+% Usage:
+% \expandafter\pgfplots@prepared@tick@pos@unpack\entry
+\def\pgfplots@prepared@tick@pos@unpack#1#2{%
+	\def\pgfplots@tick{#2}%
+	\def\pgfplots@ticknum{#1}%
+}
+
+% Draws grid lines at the a-positions of the currently set oriented
+% surface.
+%
+% Tick positions are taken out of the already precomputed list
+% \pgfplots@prepared@tick@positions@major@...
+%
+% See \pgfplots@drawticklines@onorientedsurf@ for a description of the
+% oriented surface.
+%
+% #1 : the verbatim axis name (either 'x' or 'y')
+% #2 : the index of the axis  (either 0 or 1)
+\def\pgfplots@drawgridlines@onorientedsurf{%
+	\expandafter\pgfplots@drawgridlines@onorientedsurf@\pgfplotspointonorientedsurfaceA
+}%
+\def\pgfplots@drawgridlines@onorientedsurf@#1{%
+	\pgfplots@if{pgfplots@shownothingof@\pgfplotspointonorientedsurfaceB}{%
+		\relax
+	}{%
+		\begingroup
+		\pgfplots@ifgridlines@onorientedsurf@should@be@drawn{%
+			\expandafter\let\expandafter\pgfplots@prepared@tick@positions@major@\csname pgfplots@prepared@tick@positions@major@#1\endcsname
+			\expandafter\let\expandafter\pgfplots@prepared@tick@positions@minor@\csname pgfplots@prepared@tick@positions@minor@#1\endcsname
+			\pgfplots@loop@CONTINUEfalse
+			\pgfplots@if{pgfplots@#1majorgrids}{\pgfplots@loop@CONTINUEtrue}{}%
+			\pgfplots@if{pgfplots@#1minorgrids}{\pgfplots@loop@CONTINUEtrue}{}%
+			\ifpgfplots@loop@CONTINUE
+				% I support only ONE layer for both, minor and major
+				% grid lines -- no distinction! I am lazy... FIXME
+				\pgfplotsgetlayerforstyle{%
+					every axis grid,%
+					every minor grid,%
+					every axis #1 grid,%
+					every major grid,%
+					every minor #1 grid,%
+					every major #1 grid%
+				}%
+				\pgfplotsonlayer{\pgfplotsretval}{#1 grid style}%
+				\scope
+				\pgfplots@drawgridlines@INSTALLCLIP@onorientedsurf#1%
+				%
+				\pgfplots@if{pgfplots@#1minorgrids}{%
+					\draw[%
+						/pgfplots/every axis grid,
+						/pgfplots/every minor grid,
+						/pgfplots/every axis #1 grid,
+						/pgfplots/every minor #1 grid]%
+					\pgfextra
+					\pgfplotslistforeach\pgfplots@prepared@tick@positions@minor@\as\pgfplots@curgridpos{%
+						\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curgridpos
+						\pgfplots@drawgridlines@onorientedsurf@fromto\pgfplots@tick
+					}%
+					\endpgfextra;
+				}{}%
+				%
+				\pgfplots@if{pgfplots@#1majorgrids}{%
+					\draw[%
+						/pgfplots/every axis grid,
+						/pgfplots/every major grid,
+						/pgfplots/every axis #1 grid,
+						/pgfplots/every major #1 grid]%
+					\pgfextra
+					\pgfplotslistforeach\pgfplots@prepared@tick@positions@major@\as\pgfplots@curgridpos{%
+						\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curgridpos
+						\pgfplots@drawgridlines@onorientedsurf@fromto\pgfplots@tick
+					}%
+					\endpgfextra;
+				}{}%
+				%
+				\endscope
+				\endpgfplotsonlayer
+			\fi
+		}{}%
+		\endgroup
+	}%
+}
+
+% Should draw a single grid line on the actual oriented surface.
+% #1 the value of the grid line.
+%
+% PRECONDITION
+% 	\pgfplots@ticknum contains the index of the current tick.
+\def\pgfplots@drawgridlines@onorientedsurf@fromto#1{%
+	\pgfpathmoveto{\pgfplotspointonorientedsurfaceab{#1}{\csname pgfplots@\pgfplotspointonorientedsurfaceB min\endcsname}}%
+	\pgfpathlineto{\pgfplotspointonorientedsurfaceab{#1}{\csname pgfplots@\pgfplotspointonorientedsurfaceB max\endcsname}}%
+}%
+
+% Implements 'axis x line shift' and its friends.
+%
+% It is called by grid line drawing instructions, tick lines, and tick
+% labels and installs a common shift. The purpose is to shift the
+% _entire_ axis along the outer normal.
+%
+% The operation is supposed to be used when an oriented surf is
+% installed.
+%
+%  #1 used to determine the axis for which the outer normal is to be
+%  determined. #1 is the 'a' value of the current oriented surf, i.e.
+%  one of 0, 1, 2, or v.
+%
+%  #2 same as '#1', but this here determines the 'b' value of the
+%  current oriented surf.
+%
+% exactly one of '#1' or '#2' must be 'v' such that a unique line can
+% be identified.
+\def\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf#1#2{%
+	\if v#1%
+		\let\pgfplots@loc@TMPa=\pgfplotspointonorientedsurfaceA%
+		\if v#2%
+			\pgfplots@error{Exactly one of '#1' or '#2' must be 'v', not both}%
+		\fi
+	\else
+		\if v#2%
+			\let\pgfplots@loc@TMPa=\pgfplotspointonorientedsurfaceB%
+		\else
+			\pgfplots@error{One of '#1' or '#2' must be 'v'}%
+		\fi
+	\fi
+	%
+	\pgfkeysgetvalue{/pgfplots/axis \pgfplots@loc@TMPa\space line shift}\pgfplots@loc@TMPb
+	\ifx\pgfplots@loc@TMPb\pgfutil@empty
+	\else
+		% in this case, we KNOW that is is
+		%   (a) parsed and
+		%   (b) nonzero and
+		%   (c) a dimension WITHOUT unit
+		% See \pgfplots@init@axis@shift
+		\ifpgfplots@separate@axis@lines
+		\else
+			\pgfplots@error{Internal error encountered: separate axis lines=false but axis shift found}%
+		\fi
+		\pgftransformshift{%
+			\expandafter\pgfqpointscale\expandafter{\pgfplots@loc@TMPb}{%
+				\begingroup
+				\pgf@process{%
+					\pgfplotspointonorientedsurfaceabtolinespec{#1}{#2}%
+					\expandafter\pgfplotspointouternormalvectorofaxis\expandafter{\pgfplotsretval}%
+				}%
+				\endgroup
+			}%
+		}%
+	\fi
+}%
+\let\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf@orig=\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf
+
+% Draws ticks on the currently active "oriented surface".
+%
+% The oriented surface is two dimensional and has been initialised
+% with \pgfplotspointonorientedsurfaceabsetupfor*** somehow.
+%
+% The idea is now the following:
+% - the tick positions change along the FIRST coordinate of this
+%   surface:
+%
+%   x ---- x ---- x ---- x
+%   --> FIRST -->
+%
+% - the tick lines are drawn along the SECOND coordinate of this
+%   surface:
+%
+%   | ---- | ---- | ---- |   | SECOND
+%   |      |      |      |   v
+%
+% for example,
+% \pgfplotspointonorientedsurfaceab@setupfor@xyZ{1}
+% \pgfplots@drawticklines@onorientedsurf
+% will draw ticks at x-positions designated by \pgfplots@xtick. The
+% small tick lines will be drawn along the y axis.  For each processed
+% point, the z coordinate will be fixed to '1'.
+%
+% Another example:
+% \pgfplotspointonorientedsurfaceab@setupfor@yxZ{-1}
+% \pgfplots@drawticklines@onorientedsurf
+% will draw ticks at y-positions designated by \pgfplots@ytick. The
+% small tick lines will be drawn along the x axis.  For each processed
+% point, the z coordinate will be fixed to '-1'.
+%
+% Tick positions are taken out of the already precomputed list
+% \pgfplots@prepared@tick@positions@major@...
+\def\pgfplots@drawticklines@onorientedsurf{%
+	\expandafter\pgfplots@drawticklines@onorientedsurf@\pgfplotspointonorientedsurfaceA
+}%
+
+\def\pgfplots@drawticklines@INSTALLCLIP@onorientedsurf#1{%
+	\pgfplots@drawtickgridlines@INSTALLCLIP@onorientedsurf{#1}%
+}%
+\def\pgfplots@drawgridlines@INSTALLCLIP@onorientedsurf#1{%
+	\pgfplots@drawtickgridlines@INSTALLCLIP@onorientedsurf{#1}%
+}%
+% Avoids tick lines which are too thick by introducing a clipping
+% region. Tick lines (and grid lines) won't extend to the left or
+% right of axis #1.
+\def\pgfplots@drawtickgridlines@INSTALLCLIP@onorientedsurf#1{%
+	\ifpgfplots@enable@tick@line@clipping
+		% this is LEGACY SUPPORT only. I did not want to change
+		% existing bounding boxes.
+		% Starting with 'compat=1.11', this is OFF.
+		\pgfplots@drawtickgridlines@INSTALLCLIP@onorientedsurf@{#1}%
+	\fi
+}
+\def\pgfplots@drawtickgridlines@INSTALLCLIP@onorientedsurf@#1{%
+	\pgfinterruptboundingbox%
+	\begingroup
+		% the case ||e_b|| == 0 should never happen here! Should be
+		% caught before entering this routine.
+		\let\pgfplots@loc@LENGTH=\pgfmathresult
+		\expandafter\let\expandafter\pgfplots@loc@MIN\csname pgfplots@\pgfplotspointonorientedsurfaceB min\endcsname
+		\expandafter\let\expandafter\pgfplots@loc@MAX\csname pgfplots@\pgfplotspointonorientedsurfaceB max\endcsname
+		\pgfpathmoveto{\pgfplotspointonorientedsurfaceabwithbshift{\csname pgfplots@#1min\endcsname}{\pgfplots@loc@MIN}{-5cm}}%
+		\pgfpathlineto{\pgfplotspointonorientedsurfaceabwithbshift{\csname pgfplots@#1max\endcsname}{\pgfplots@loc@MIN}{-5cm}}%
+		\pgfpathlineto{\pgfplotspointonorientedsurfaceabwithbshift{\csname pgfplots@#1max\endcsname}{\pgfplots@loc@MAX}{5cm}}%
+		\pgfpathlineto{\pgfplotspointonorientedsurfaceabwithbshift{\csname pgfplots@#1min\endcsname}{\pgfplots@loc@MAX}{5cm}}%
+		\pgfusepath{clip}%
+	\endgroup
+	\endpgfinterruptboundingbox%
+}%
+
+\def\pgfplots@drawticklines@onorientedsurf@#1{%
+	\pgfplots@if{pgfplots@shownothingof@\pgfplotspointonorientedsurfaceB}{%
+		\relax
+	}{%
+		\begingroup
+		% FIXME : the case of centered axis lines is missing here.
+		% 0 = lower limit,
+		% 1 = upper limit
+		% 2 = centered...
+		% Currently, "0" is also used for centered axis lines.
+		\if 2\csname pgfplots@\pgfplotspointonorientedsurfaceA axislinesnum\endcsname
+			% centered axis lines need no complicated visibility
+			% checks. They need tick lines in *any* case.
+			% So, enabled them here!
+			\def\pgfplots@drawticklines@for@placecomputedtick@LOWEROK{1}%
+			%
+			% UPPEROK is unnecessary in this context - for centered
+			% axis lines, it will be disabled anyway because the
+			% styles like 'axis x line*' set \pgfplots@xticknum to
+			% center - and that means that no tick lines will be drawn
+			% (see below).
+			\def\pgfplots@drawticklines@for@placecomputedtick@UPPEROK{0}%
+		\else
+			\pgfplots@ifaxisline@B@onorientedsurf@should@be@drawn{0}{%
+				\def\pgfplots@drawticklines@for@placecomputedtick@LOWEROK{1}%
+			}{%
+				\def\pgfplots@drawticklines@for@placecomputedtick@LOWEROK{0}%
+			}%
+			\pgfplots@ifaxisline@B@onorientedsurf@should@be@drawn{1}{%
+				\def\pgfplots@drawticklines@for@placecomputedtick@UPPEROK{1}%
+			}{%
+				\def\pgfplots@drawticklines@for@placecomputedtick@UPPEROK{0}%
+			}%
+		\fi
+		\if\pgfkeysvalueof{/pgfplots/\pgfplotspointonorientedsurfaceB\space dir/value}r%
+			% local special handling for reversed axes: exchange
+			% meaning of 'left' and 'right' here.
+			%
+			% the rest of the pgfplots code does that automatically because
+			% there, tickposnum is relevant to determine the axes
+			% which contains tick labels. And this algorithm checks
+			% for reversed axes implicitly.
+			%
+			\if1\csname pgfplots@\pgfplotspointonorientedsurfaceA tickposnum\endcsname
+				\expandafter\def\csname pgfplots@\pgfplotspointonorientedsurfaceA tickposnum\endcsname{3}%
+			\else
+				\if3\csname pgfplots@\pgfplotspointonorientedsurfaceA tickposnum\endcsname
+					\expandafter\def\csname pgfplots@\pgfplotspointonorientedsurfaceA tickposnum\endcsname{1}%
+				\fi
+			\fi
+		\fi
+		\ifcase\csname pgfplots@\pgfplotspointonorientedsurfaceA tickposnum\endcsname\relax
+			% both
+		\or
+			% lower
+			\def\pgfplots@drawticklines@for@placecomputedtick@UPPEROK{0}%
+		\or
+			% center
+			\def\pgfplots@drawticklines@for@placecomputedtick@UPPEROK{0}%
+		\or
+			% upper
+			\def\pgfplots@drawticklines@for@placecomputedtick@LOWEROK{0}%
+		\else
+			% never used?
+			\def\pgfplots@drawticklines@for@placecomputedtick@UPPEROK{0}%
+		\fi
+		\expandafter\let\expandafter\pgfplots@prepared@tick@positions@major@\csname pgfplots@prepared@tick@positions@major@#1\endcsname
+		\expandafter\let\expandafter\pgfplots@prepared@tick@positions@minor@\csname pgfplots@prepared@tick@positions@minor@#1\endcsname
+		%
+		% There is only ONE layer for both, minor and major
+		% tick lines -- no distinction!
+		\pgfplotsgetlayerforstyle{%
+			every tick,%
+			every minor tick,%
+			every #1 tick,%
+			every major tick,%
+			every minor #1 tick,%
+			every major #1 tick%
+		}%
+		\pgfplotsonlayer{\pgfplotsretval}{#1tick style}%
+		\scope
+		\pgfplots@drawticklines@INSTALLCLIP@onorientedsurf#1
+		%
+		\pgfplots@if{pgfplots@#1minorticks}{%
+			\draw[%
+				/pgfplots/every tick,
+				/pgfplots/every minor tick,
+				/pgfplots/every #1 tick,
+				/pgfplots/every minor #1 tick]%
+			\pgfextra
+			\pgfmathparse{\pgfplots@subtickwidth}%
+			\let\pgfplots@subtickwidth@=\pgfmathresult
+			\let\pgfplots@subtickwidth@=\pgfmathresult
+			\let\pgfplots@subtickwidth=\pgfmathresult
+			\pgfplots@prepare@tick@offsets@for@{#1}{\pgfplots@subtickwidth@}%
+			\if\pgfplots@drawticklines@for@placecomputedtick@LOWEROK1%
+				\begingroup
+				\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf v0%
+				\pgfplotslistforeach\pgfplots@prepared@tick@positions@minor@\as\pgfplots@curtickpos{%
+					\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curtickpos
+					\let\pgfplots@curtickpos=\pgfplots@tick
+					\pgfplotspointouternormalvectorofaxissetv{}{\pgfplots@curtickpos}%
+					\pgfplots@drawticklines@for@placecomputedtick@LOWER
+				}%
+				\endgroup
+			\fi
+			%
+			\if\pgfplots@drawticklines@for@placecomputedtick@UPPEROK1%
+				\begingroup
+				\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf v1%
+				\pgfplotslistforeach\pgfplots@prepared@tick@positions@minor@\as\pgfplots@curtickpos{%
+					\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curtickpos
+					\let\pgfplots@curtickpos=\pgfplots@tick
+					\pgfplotspointouternormalvectorofaxissetv{}{\pgfplots@curtickpos}%
+					\pgfplots@drawticklines@for@placecomputedtick@UPPER
+				}%
+				\endgroup
+			\fi
+			\endpgfextra;
+		}{}%
+		%
+		\pgfplots@if{pgfplots@#1majorticks}{%
+			\draw[%
+				/pgfplots/every tick,
+				/pgfplots/every major tick,
+				/pgfplots/every #1 tick,
+				/pgfplots/every major #1 tick]%
+			\pgfextra
+			\pgfmathparse{\pgfplots@tickwidth}%
+			\let\pgfplots@tickwidth@=\pgfmathresult
+			\let\pgfplots@tickwidth@=\pgfmathresult
+			\let\pgfplots@tickwidth=\pgfmathresult
+			\pgfplots@prepare@tick@offsets@for@{#1}{\pgfplots@tickwidth@}%
+			\if\pgfplots@drawticklines@for@placecomputedtick@LOWEROK1%
+				\begingroup
+				\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf v0%
+				\pgfplotslistforeach\pgfplots@prepared@tick@positions@major@\as\pgfplots@curtickpos{%
+					\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curtickpos
+					\let\pgfplots@curtickpos=\pgfplots@tick
+					\pgfplotspointouternormalvectorofaxissetv{}{\pgfplots@curtickpos}%
+					\pgfplots@drawticklines@for@placecomputedtick@LOWER
+				}%
+				\endgroup
+			\fi
+			%
+			\if\pgfplots@drawticklines@for@placecomputedtick@UPPEROK1%
+				\begingroup
+				\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf v1%
+				\pgfplotslistforeach\pgfplots@prepared@tick@positions@major@\as\pgfplots@curtickpos{%
+					\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curtickpos
+					\let\pgfplots@curtickpos=\pgfplots@tick
+					\pgfplotspointouternormalvectorofaxissetv{}{\pgfplots@curtickpos}%
+					\pgfplots@drawticklines@for@placecomputedtick@UPPER
+				}%
+				\endgroup
+			\fi
+			\endpgfextra;
+		}{}%
+		%
+		\endscope
+		\endpgfplotsonlayer
+		\endgroup
+	}%
+}
+
+% Draws tick labels at the positions of the currently set oriented
+% surface.
+%
+% Tick positions are taken out of the already precomputed list
+% \pgfplots@prepared@tick@positions@major@...
+%
+% For 2D axes, this task is relatively simple:
+% we iterate through every prepared major tick position and place a
+% tick label. Open points, however, are the question whether to use
+% the RIGHT or the LEFT axis line on the current oriented surface:
+%
+%        direction 'b' (second oriented)
+% |-------------------------|
+% |                         |
+% |                         |
+% |                         |direction 'a' (first oriented)
+% |                         |
+% |                         |
+% |                         |
+% |-------------------------|
+% Left                  Right
+%
+% Given the axis line which shall contain the labels, we have to
+% decide how to align the tick label nodes: on the left or on the
+% right? Of course, we want to align them such that they are "outside"
+% of the figure! That's simple as well: "Left axis line => outside
+% means left", "Right axis line => outside means right".
+%
+% That's all, basically.
+%
+% For 3D axes, all these points are basically ... the same!
+% Now it can happen that the current oriented surface shall
+% not contain ANY tick label. In that case, we do nothing.
+% Furthermore, the "outside" direction (i.e. the anchoring of the
+% label nodes) is a little bit more difficult.
+%
+%
+% See \pgfplots@drawticklines@onorientedsurf@ for a description of the
+% oriented surface.
+\def\pgfplots@drawticklabels@onorientedsurf{%
+	\expandafter\pgfplots@drawticklabels@onorientedsurf@\pgfplotspointonorientedsurfaceA
+}
+\def\pgfplots@drawticklabels@onorientedsurf@#1{%
+	\begingroup
+	\expandafter\let\expandafter\pgfplots@prepared@tick@positions@major@\csname pgfplots@prepared@tick@positions@major@#1\endcsname
+	% check whether
+	% - we need to place tick labels on the LEFT side,
+	% - we need to place tick labels on the RIGHT side,
+	% - we don't need tick labels for the current surface at all:
+	\pgfplotspointonorientedsurfaceabmatchaxisline{\csname pgfplots@#1ticklabelaxisspec\endcsname}{\pgfplots@ticklabelside}%
+	\ifx\pgfplots@ticklabelside\pgfutil@empty
+		% SKIP. The current oriented surface shall not get tick labels
+		% for #1.
+	\else
+		\pgfplots@if{pgfplots@#1majorticks}{%
+			\pgfplots@if{pgfplots@#1islinear}{%
+				\pgfplots@init@scaled@tick@for{#1}%
+			}{\relax}%
+			\begingroup
+			\pgfplotsgetlayerforstyle{every tick label,every #1 tick label}%
+			\pgfplotsonlayer\pgfplotsretval{#1tick label style}%
+			\pgfkeysalso{/tikz/every node/.append style={/pgfplots/every tick label,/pgfplots/every #1 tick label}}%
+			\pgfplots@drawticklabels@onorientedsurf@prepareanchor#1%
+			%
+			\pgfplotsmath@ifzero{\csname pgfplots@\pgfplotspointonorientedsurfaceB @veclength\endcsname}{%
+				\def\pgfplots@tick@offset{0}%
+			}{%
+				\ifcase\csname pgfplots@#1tickalignnum\endcsname\relax
+					\def\pgfmathresult{0}%
+				\or
+					\pgfmathparse{\pgfplots@tickwidth}%
+				\or
+					\pgfmathmultiply{0.5}{\pgfplots@tickwidth}%
+				\fi
+				\let\pgfplots@tick@offset=\pgfmathresult
+				\pgfplots@ticklabel@maxtickdimen@extrashift{#1}{\pgfplots@tick@offset}%
+				%\pgfmathmultiply@{\pgfplots@tick@offset}{\csname pgfplots@\pgfplotspointonorientedsurfaceB @inverseveclength\endcsname}%
+				\let\pgfplots@tick@offset=\pgfmathresult
+			}%
+			%
+			\if2\csname pgfplots@#1axislinesnum\endcsname % Centered axis lines?
+				\expandafter\let\expandafter\pgfplots@tick@origin\csname pgfplots@logical@ZERO@\pgfplotspointonorientedsurfaceB\endcsname%
+				% FIXME : that stuff here does not respect
+				% '[xyz]tickpos num' keys!
+				\pgfplots@tickposchoiceb%<-- backw. compat, is usually empty.
+				\if r\pgfkeysvalueof{/pgfplots/\pgfplotspointonorientedsurfaceB\space dir/value}%
+					% special handling for reversed axes.
+					\pgfmathmultiply{-1}{\pgfplots@tick@offset}%
+					\let\pgfplots@tick@offset=\pgfmathresult
+				\fi
+				\edef\pgfplots@tick@offset{-\pgfplots@tick@offset}%
+				%\pgfmathsubtract@{\pgfplots@tick@origin}{\pgfplots@tick@offset}%
+			\else
+				\if0\pgfplots@ticklabelside
+				 	\expandafter\let\expandafter\pgfplots@tick@origin\csname pgfplots@\pgfplotspointonorientedsurfaceB min\endcsname%
+					\pgfplots@tickposchoiceb%<-- backw. compat, is usually empty.
+					%\pgfmathsubtract@{\pgfplots@tick@origin}{\pgfplots@tick@offset}%
+					\edef\pgfplots@tick@offset{-\pgfplots@tick@offset}%
+				\else
+					\if1\pgfplots@ticklabelside
+						\expandafter\let\expandafter\pgfplots@tick@origin\csname pgfplots@\pgfplotspointonorientedsurfaceB max\endcsname%
+						\pgfplots@tickposchoicea%<-- backw. compat, is usually empty.
+						%\pgfmathadd@{\pgfplots@tick@origin}{\pgfplots@tick@offset}%
+					\else
+						% FIXME : ticklabelside == 2 is, in principle,
+						% a valid choice. It is the case handled with
+						% "if 2 == pgfplots@#1axislinesnum" above,
+						% isn't it!?
+						\expandafter\let\expandafter\pgfplots@tick@origin\csname pgfplots@logical@ZERO@\pgfplotspointonorientedsurfaceB\endcsname%
+						% FIXME : is that correct!?
+						\pgfplots@tickposchoiceb%<-- backw. compat, is usually empty.
+						%\pgfmathsubtract@{\pgfplots@tick@origin}{\pgfplots@tick@offset}%
+						%\edef\pgfplots@tick@offset{-\pgfplots@tick@offset}%
+						%
+						\if2\pgfplots@ticklabelside
+						\else
+							% Should never happen.
+							\pgfplots@error{Internal logic error during tick label placement (got placement character '\pgfplots@ticklabelside').
+									Please report this as a bug or verify your input arguments to #1ticklabel pos.}%
+						\fi
+					\fi
+				\fi
+			\fi
+			%\let\pgfplots@tick@origin=\pgfmathresult%
+			\edef\pgfplots@tick@offset{\pgfplots@tick@offset pt}%
+			%
+			% make sure the \pgfmathlogtologten method works even for
+			% non-standard 'log basis #1':
+			\def\pgfmathlogtologten@{\pgfplotscoordmath{#1}{log to log 10}}%
+			%
+			\xdef\pgfplots@show@ticklabel@LASTTICK{}%
+			\pgfplots@transformshift@along@outer@normal@on@line@of@oriented@surf{v}{\pgfplots@ticklabelside}%
+			\pgfplotslistforeachungrouped\pgfplots@prepared@tick@positions@major@\as\pgfplots@curtickpos{%
+				\expandafter\pgfplots@prepared@tick@pos@unpack\pgfplots@curtickpos
+				\let\pgfplots@curtickpos=\pgfplots@tick
+				%\expandafter\pgfplotslistpopfront\csname pgfplots@prepared@tick@positions@major@tickindices@#1\endcsname\to\pgfplots@ticknum
+				\pgfplotspointouternormalvectorofaxissetv{}{\pgfplots@curtickpos}%
+				\pgfplots@show@ticklabel
+					{#1}{\pgfplots@curtickpos}(\pgfplots@curtickpos,\pgfplots@tick@origin+\pgfplots@tick@offset)%
+					{\pgfplots@ticknum}%
+			}%
+			\pgfplots@ticklabel@maxtickdimen@finish{#1}%
+			\endpgfplotsonlayer
+			\endgroup
+			\pgfplots@draw@tick@scale@label@for #1%
+		}%
+		{\relax}% if...@major==false
+	\fi
+	\endgroup
+}
+
+% This here does the main work for any tick label ANCHORING.
+%
+% FIXME : the new 'near ticklabel' anchors are now the default method
+% to place tick labels.
+%
+% This feature here can be used to disable this anchoring, i.e. set
+% 'ticklabel anchor=tikz' to use a stupid heuristics.
+%
+% There are actually two choices:
+% Choice 1: near ticklabel
+% 	This choice places tick labels fully automatic outside of the
+% 	figure, all on a line which is parallel to the axis which
+% 	contains the corresponding tick labels.
+%
+% 	Since we are currently working on a restricted surface, we have
+% 	three direction related to that surface:
+% 		'a': this is the direction in which tick positions are known.
+% 		'b': the 'orthogonal' axis to 'a' which is also in the surface.
+% 		'n': the surface normal.
+% 	Now, the idea for tick labels is to place them at
+% 		SCALE_b * b + SCALE_n * n,
+% 	where the SCALE_[bn] numbers are choosen such that the label is
+% 	outside of the axis.
+%
+% 	The offset is simply added to the transformation matrix (as a
+% 	shift).
+%
+% Choice 2: tikz.
+%   This is more or less a backwards compatibility feature. It does
+%   not change the transformation matrix. It simply sets the 'at' key
+%   of each node to the tick position and prepares the correct anchors
+%   for the TikZ '\node' commands.
+%   That's all here.
+%
+% PRECONDITION:
+%   - called inside of  \pgfplots@drawticklabels@onorientedsurf@
+% POSTCONDITION:
+% 	- defines
+% 		- \pgfplots@tickposchoicea
+% 			if called, sets keys such that tick labels are RIGHT (TOP) of
+% 			the axis,
+%
+% 		- \pgfplots@tickposchoiceb
+% 			if called, sets keys such that tick labels are LEFT
+% 			(BOTTOM) of the axis,
+\def\pgfplots@drawticklabels@onorientedsurf@prepareanchor#1{%
+	\if\csname pgfplots@ticklabel@anchor@#1\endcsname0%
+		% auto is the same as 'near ticklabel':
+		\expandafter\def\csname pgfplots@ticklabel@anchor@#1\endcsname{1}%
+	\fi
+	\ifcase\csname pgfplots@ticklabel@anchor@#1\endcsname%
+		% 0: doesn't happen, see above.
+	\or
+		% 1: near ticklabel.
+		% The following code contains automatically
+		% aligned tick labels (especially for 3D axes).
+		% see the manual for the |near ticklabel| anchors.
+		\pgfkeys{/tikz/anchor=near #1ticklabel}%
+		%
+		% process the (optional) ticklabel distance:
+		\begingroup
+			\pgfkeysgetvalue{/pgfplots/#1ticklabel shift}\pgfmathresult
+			\ifx\pgfmathresult\pgfutil@empty
+			\else
+				\afterassignment\pgfplots@gobble@until@relax
+				\pgf@xa=\pgfkeysvalueof{/pgfplots/#1ticklabel shift}pt\relax
+				\multiply\pgf@xa by-1 % the direction vector points to the INSIDE. the shift should have opposite sign.
+				\edef\pgfmathresult{\pgf@sys@tonumber\pgf@xa}%
+			\fi
+			\pgfmath@smuggleone\pgfmathresult
+		\endgroup
+		\let\pgfplots@loc@TMPc=\pgfmathresult
+		\ifx\pgfplots@loc@TMPc\pgfutil@empty
+		\else
+			\pgfplotspointouternormalvectorofaxis@ifdependson@v{\pgfplotsticklabelaxisspec{#1}}{%
+				\tikzset{every node/.append code={%
+						\tikz@addtransform{%
+							\pgftransformshift{%
+								\pgfplotspointouternormalvectorofticklabelaxis{#1}%
+								\pgfqpointscale{-\pgfplots@loc@TMPc}{}%
+							}%
+						}%
+					}%
+				}%
+			}{%
+				\pgftransformshift{%
+					\pgfplotspointouternormalvectorofticklabelaxis{#1}%
+					\pgfqpointscale{-\pgfplots@loc@TMPc}{}%
+				}%
+			}%
+		\fi
+		%
+		% these things are irrelevant here:
+		\let\pgfplots@tickposchoicea=\pgfutil@empty
+		\let\pgfplots@tickposchoiceb=\pgfutil@empty
+	\or
+		% 2: tikz.
+		% We simply prepare the default anchor.
+		% Actually, this code is just for backwards compatibility -
+		% there may be people who prefer to set anchors. The
+		% 'near ticklabel' implementation is much more general, however.
+		\if\pgfplotspointonorientedsurfaceB x
+			\def\pgfplots@tickposchoicea{\tikzset{right}}%
+			\def\pgfplots@tickposchoiceb{\tikzset{left}}%
+		\else
+			\if\pgfplotspointonorientedsurfaceB y
+				\def\pgfplots@tickposchoicea{\tikzset{above}}%
+				\def\pgfplots@tickposchoiceb{\tikzset{below}}%
+			\else
+				\def\pgfplots@tickposchoicea{\tikzset{anchor=north east}}%
+				\def\pgfplots@tickposchoiceb{\tikzset{anchor=south east}}%
+			\fi
+		\fi
+	\fi
+	%
+	%
+	%
+	% suppress this warning. We cannot avoid it currently.
+	\pgfkeysgetvalue{/pgfplots/warning/missing near ticklabel at/.@cmd}\pgfplots@loc@neartickllabel@at
+	\pgfkeysdef{/pgfplots/warning/missing near ticklabel at}{}%
+	%
+	\pgfplots@ticklabel@maxtickdimen@prepare@for@normalvec
+		{#1}%
+		{\pgfplotspointouternormalvectorofticklabelaxis{#1}}%
+	%
+	\pgfkeyslet{/pgfplots/warning/missing near ticklabel at/.@cmd}\pgfplots@loc@neartickllabel@at%
+}%
+
+\newif\ifpgfplots@checkuniform@isfirst
+% Checks whether the argument to xtick or ytick is a UNIFORM tick
+% sequence.
+%
+% A uniform tick sequence is 0,...,10 and 3,4,5 and -5,-4,-2 but
+% NOT 0,2,4 or 4,10.
+%
+% Furthermore, any NON-integer tick arguments are also assumed to be
+% NOT uniform.
+%
+% INPUT:
+% #1: a tick argument (i.e. something which can be put to
+%     \foreach \x in {#1})
+%
+% OUTPUT:
+%    \pgfplots@isuniformticktrue
+%    or
+%    \pgfplots@isuniformtickfalse
+% depending on the check.
+% This variable will be set globally.
+\def\pgfplots@checkisuniformLOGtick#1{%
+	\begingroup
+	\global\pgfplots@isuniformticktrue
+	\pgfplots@checkuniform@isfirsttrue
+	\foreach \x in {#1}{%
+		\pgfmathmultiply@\x\reciproclogten
+		\let\cur=\pgfmathresult
+		% check whether
+		%  \cur - last == 1 (last = \pgfplots@glob@TMPb)
+		\ifpgfplots@checkuniform@isfirst
+			\global\pgfplots@checkuniform@isfirstfalse
+		\else
+			\pgfmathsubtract@\cur\pgfplots@glob@TMPb%
+			\pgfmathapproxequalto@\pgfmathresult{1.0}%
+			\ifpgfmathcomparison
+			\else
+				\global\pgfplots@isuniformtickfalse
+				\breakforeach
+			\fi
+		\fi
+		\global\let\pgfplots@glob@TMPb=\cur
+	}%
+	\endgroup
+}
+
+% Checks whether the linear tick sequence #1 is a uniform tick.
+%
+% It also assigns pgfplots@tick@distance@#1  as the distance.
+%
+% see \pgfplots@checkisuniformLOGtick for details.
+%
+% #1: a tick sequence (expanded)
+% #2: a macro which will be filled with the tick distance. This is
+% only valid if \pgfplots@isuniformticktrue.
+\def\pgfplots@checkisuniformLINEARtick#1#2{%
+	\begingroup
+	\global\pgfplots@isuniformticktrue
+	\pgfplots@checkuniform@isfirsttrue
+	\global\let\pgfplots@glob@TMPb=\pgfutil@empty
+	\global\def\pgfplots@glob@TMPa{1}%
+	\foreach \x in {#1}{%
+		\ifx\pgfplots@glob@TMPb\pgfutil@empty
+		\else
+			\pgfmathsubtract@\x\pgfplots@glob@TMPb
+			\ifpgfplots@checkuniform@isfirst
+				% remember first distance h = x_1 - x_0
+				\global\let\pgfplots@glob@TMPa=\pgfmathresult
+				\global\pgfplots@checkuniform@isfirstfalse
+			\else
+				% check whether x_i - x_{i-1} = h
+				\pgfmathapproxequalto@\pgfmathresult\pgfplots@glob@TMPa%
+				\ifpgfmathcomparison
+				\else
+					\global\pgfplots@isuniformtickfalse
+					\breakforeach
+				\fi
+			\fi
+		\fi
+		\global\let\pgfplots@glob@TMPb=\x%
+	}%
+	\endgroup
+	\let#2=\pgfplots@glob@TMPa
+}
+
+% helper method which computes log10*\x foreach \x in {#1}.
+% The result will be \xdef'ed into #2.
+%
+% #1: the ticks
+% #2: the output macro
+% #3: the axis
+\def\pgfplots@compute@tick@times@logten#1\to#2#3{%
+	\global\let#2=\pgfutil@empty
+	\foreach \pgfplots@loc@TMPb in {#1} {%
+		\pgfplotscoordmath{#3}{log to display log}{\pgfplots@loc@TMPb}%
+		\pgfplotscoordmath{#3}{tofixed}{\pgfmathresult}%
+		\ifx#2\pgfutil@empty
+			\xdef#2{\pgfmathresult}%
+		\else
+			\xdef#2{#2,\pgfmathresult}%
+		\fi
+	}%
+}
+
+% This is part of \pgfplots@assign@default@tick@foraxis and relies on
+% its temporary variables.
+%
+% INPUT:
+%   \MIN : the lower axis limit of #1 (a TeX register, in transformed range)
+%   \MAX : same with upper axis limit of #1
+%   \desirednumticks: the number of ticks to use
+% OUTPUT:
+%   \H   :  will contain the (transformed) distance between adjacent ticks
+%   \aftergroup\pgfplots@isuniformticktrue set if it applies
+%
+%
+\def\pgfplots@ticks@compute@tick@distance#1{%
+	% compute step size 'H':
+	\H=\MAX
+	\advance\H by-\MIN
+	\ifdim\H<0pt \H=-1\H \fi
+%\message{Axis limit #1: [\the\MIN:\the\MAX], diff = \the\H.^^J}%
+	\c@pgf@counta=\desirednumticks
+	\advance\c@pgf@counta by-1 %
+	\divide\H by\c@pgf@counta
+%\message{determining ticks for #1-axis: Wr := (width/max space between ticks) = \the\Wr, desirednumticks=max(\axisdefaulttryminticks, trunc(Wr)) = \the\desirednumticks, H#1=(axis range/(desirednumticks-1)) = \the\H^^J}%
+	%
+	% SEARCH for the NEXT FEASABLE H.
+	\edef\Hmacro{\pgf@sys@tonumber\H}%
+	\ifpgfplots@cur@is@linear
+		% CASE LINEAR AXIS
+		\ifpgfplots@is@datascaled
+			% This here works if the scaling trafo is linear.
+			\pgfplotscoordmath{#1}{datascaletrafo noshift inverse}{\Hmacro}%
+			\let\Hmacro=\pgfmathresult
+		\else
+			\pgfmathfloatparsenumber{\Hmacro}%
+			\let\Hmacro=\pgfmathresult
+		\fi
+%\message{Got T^{-1}(H#1) = \Hmacro^^J}%
+		%
+		\pgfplots@assign@default@tick@foraxis@normalizetickdist#1\Hmacro
+		\let\Hmacro=\pgfmathresult
+		%
+		% Ok, we are ready.
+		% Now, convert everything into the fixed point data
+		% range:
+		\ifpgfplots@is@datascaled
+			\pgfplotscoordmath{#1}{datascaletrafo noshift}{\Hmacro}%
+			\H=\pgfmathresult pt
+		\else
+			\pgfmathfloattofixed\Hmacro
+			\H=\pgfmathresult pt
+		\fi
+		%
+%\message{snapped-to-nicest = \Hmacro^^J}%
+		\aftergroup\pgfplots@isuniformticktrue
+	\else
+		% CASE LOG AXIS
+		%
+		% search for the "best" H= j* log(10),  j an integer.
+		%
+		% And prefer j=1 if that is possible (otherwise minor
+		% ticks are not useful).
+		\pgfmath@basic@multiply@{\Hmacro}{\pgfplots@loc@log@from@display@log@scale}%
+		\let\Hmacrobaseten=\pgfmathresult
+		\expandafter\H\pgfmathresult pt
+%\message{ [ H / log(10) = \pgfmathresult ]}%
+		\ifdim\H<2pt
+			\H=1pt
+		\else
+			\ifnum\H<1pt
+				\H=1pt
+			\else
+				\expandafter\pgfmathfloor\expandafter{\pgfmathresult}%
+				\expandafter\H\pgfmathresult pt
+			\fi
+		\fi
+		\ifdim\H=1pt
+			\aftergroup\pgfplots@isuniformticktrue
+			\pgfplots@isuniformticktrue
+		\else
+			\aftergroup\pgfplots@isuniformtickfalse
+			\pgfplots@isuniformtickfalse
+		\fi
+%\message{final H=\pgf@sys@tonumber{\H} * log(10)}%
+		\H=\pgfplots@loc@log@to@display@log@scale\H
+	\fi
+	%
+}
+
+% This is part of \pgfplots@assign@default@tick@foraxis
+%
+% It overwrites the internal macros of
+% \pgfplots@assign@default@tick@foraxis
+\def\pgfplots@assign@default@tick@foraxis@autoadjust@result#1{%
+	\pgfplots@if{pgfplots@#1islinear}{%
+		\begingroup
+			\def\pgfplots@tick@returnval@ready{1}%
+			\pgfutil@ifundefined{pgfplots@assign@default@tick@foraxis@recurselevel}{%
+				\def\pgfplots@assign@default@tick@foraxis@recurselevel{1}%
+			}{%
+				\pgfplotsutil@advancestringcounter\pgfplots@assign@default@tick@foraxis@recurselevel
+			}%
+			\ifnum\pgfplots@assign@default@tick@foraxis@recurselevel<15
+				\c@pgf@counta=\axisdefaulttryminticks\relax
+				\advance\c@pgf@counta by1
+				\edef\axisdefaulttryminticks{\the\c@pgf@counta}%
+%\message{**TOO FEW TICK LABELS FOR #1. RECURSION with try min ticks=\axisdefaulttryminticks.**^^J}%
+				% recurse.
+				\pgfplots@assign@default@tick@foraxis{#1}%
+				\expandafter\global\expandafter\let\expandafter\pgfplots@glob@TMPa\csname pgfplots@#1tick\endcsname
+				\expandafter\global\expandafter\let\expandafter\pgfplots@glob@TMPb\csname pgfplots@tick@distance@#1\endcsname
+			\else
+				\pgfplotswarning{tick computation failed}{#1}{\axisdefaulttryminticks}\pgfeov%
+				\def\pgfplots@tick@returnval@ready{0}%
+			\fi
+			\pgfmath@smuggleone\pgfplots@tick@returnval@ready
+		\endgroup
+	}{%
+		% Case logarithmic axes and too few ticks.
+		\aftergroup\pgfplots@isuniformtickfalse
+		% ok, do something special.
+		%
+		% The idea is now to place ticks at
+		% 10^{i*h} with properly choosen 'h'.
+		%
+		% So: apply basically the SAME code as above for linear
+		% axis, just everything log 10!  And keep in mind that all
+		% coordinates are actually given as natural logarithms.
+		\MIN\csname pgfplots@#1min\endcsname pt
+		\H=\MAX
+		\advance\H by-\MIN
+		\ifdim\H<0pt \H=-1\H \fi
+		\H=\pgfplots@loc@log@from@display@log@scale \H
+%\message{Axis limit #1: [\the\MIN:\the\MAX], diff/log(10) = \the\H.}%
+		\c@pgf@counta=\desirednumticks\relax
+		\advance\c@pgf@counta by-1
+		\ifnum\c@pgf@counta>2
+			% subtract one more. This algorithm here produces more
+			% ticks than the normal one which is designed for 10^i
+			\advance\c@pgf@counta by-1
+		\fi
+		\divide\H by\c@pgf@counta\relax
+%\message{determining ticks for #1-axis: Wr := (width/max space between ticks) = \the\Wr, desirednumticks=max(\axisdefaulttryminticks, trunc(Wr)) = \the\desirednumticks, H#1=(axis range/(desirednumticks-1)) = \the\H}%
+		%
+		% SEARCH for the NEXT FEASABLE H.
+		\edef\Hmacro{\pgf@sys@tonumber\H}%
+		\pgfmathfloatparsenumber{\Hmacro}%
+		\pgfplots@assign@default@tick@foraxis@normalizetickdist#1\pgfmathresult
+		%
+		\expandafter\pgfmathfloattofixed\expandafter{\pgfmathresult}%
+		\let\Hmacro=\pgfmathresult
+		\H=\Hmacro pt %
+		% Ok, our step size h for 10^{i*h} is ready!
+%\message{determined step size 10^{\Hmacro}}%
+		% Now, we want to activate the Tick set {10^{i*H}, i in \Z}
+		% compute I such that
+		%   10^{min} = 10^{I * H + rest};  |rest| < H
+		% -> I = round(xmin/H)
+		% -> MIN = I * H
+		% BUT EVERYTHING to log(10) basis!
+		\MIN=\pgfplots@loc@log@from@display@log@scale \MIN
+		\pgfmathlog@invoke@expanded\pgfmathdivide@{%
+			{\pgf@sys@tonumber\MIN}%
+			{\Hmacro}%
+		}%
+		\pgfmathsetcount{\c@pgf@counta}{\pgfmathresult}%
+		\ifdim\MIN<0pt
+			% the truncation rounds TOWARDS 0 which is not what I want.
+			\advance\c@pgf@counta by-1
+		\fi
+		\MIN=\H\relax
+		\multiply\MIN by\c@pgf@counta\relax
+		%
+		% convert back to basis 'e':
+		\MIN=\pgfplots@loc@log@to@display@log@scale\MIN\relax
+		\H=\pgfplots@loc@log@to@display@log@scale\H\relax
+		\MINH=\MIN\relax
+		\advance\MINH by\H\relax
+	}%
+}%
+
+% The key 'xtick distance' can have any value. This routine avoids
+% 'dimension too large' programmatically
+\def\pgfplots@assign@default@tick@foraxis@cap#1#2{%
+	\ifpgfplots@is@datascaled
+		\pgfplotscoordmath{#1}{datascaletrafo noshift inverse}{16000}%
+		\let\pgfplots@loc@TMPb=\pgfmathresult
+		\pgfplotscoordmath{#1}{min}{\pgfplots@loc@TMPb}{#2}%
+		\edef\pgfplots@loc@TMPc{#2}%
+		\ifx\pgfplots@loc@TMPc\pgfmathresult
+			\def\pgfplotsretval{0}%
+		\else
+			\def\pgfplotsretval{1}%
+		\fi
+	\else
+		\pgfplotscoordmath{#1}{parsenumber}{16000}%
+		\let\pgfplots@loc@TMPb=\pgfmathresult
+		\pgfplotscoordmath{#1}{parsenumber}{#2}%
+		\let\pgfplots@loc@TMPc=\pgfmathresult
+		\pgfplotscoordmath{#1}{min}{\pgfplots@loc@TMPb}{\pgfplots@loc@TMPc}%
+		%
+		\ifx\pgfplots@loc@TMPc\pgfmathresult
+			\def\pgfplotsretval{0}%
+		\else
+			\def\pgfplotsretval{1}%
+		\fi
+	\fi
+	%
+	\ifx1\pgfplotsretval
+		\pgfplotswarning{dimension too large in ticks}{#1}{#2}{\pgfmathresult}\pgfeov%
+	\fi
+}
+
+% This macro is NECESSARILY part of
+% \pgfplots@assign@default@tick@foraxis:
+%
+% OUTPUT:
+% \csname pgfplots@#1tick\endcsname
+% \csname pgfplots@tick@distance@#1\endcsname
+% FIXME : some \aftergroup trickery...
+\def\pgfplots@assign@default@tick@foraxis@compute#1{%
+	% Ok, we have either log or linear axis and need default
+	% ticks MIN,MIN+H,...,MAX.
+	\let\MINH=\pgf@xa
+	\let\H=\pgf@xb
+	\let\MAX=\pgf@ya
+	\let\MIN=\pgf@yb
+	%
+	\MAX=\csname pgfplots@#1max\endcsname pt %
+	\advance\MAX by0.001pt % avoid round errors
+	%\expandafter\MIN\the\c@pgf@counta pt
+	\MIN=\csname pgfplots@#1min\endcsname pt %
+	%
+	\pgfkeysgetvalue{/pgfplots/#1tick distance}\Hmacro
+	\ifx\Hmacro\pgfutil@empty
+		\pgfplots@ticks@compute@tick@distance#1%
+		\def\b@pgfplots@ticks@computed@tick@distance@is@final{0}%
+	\else
+		% ah - a user argument. OK, prefer that one over the default:
+		\pgfplotscoordmath{default}{parse}{abs(\Hmacro)}%
+		\let\Hmacro=\pgfmathresult
+		\pgfplots@assign@default@tick@foraxis@cap{#1}{\Hmacro}%
+		\let\Hmacro=\pgfmathresult
+		\ifpgfplots@cur@is@linear
+			\ifpgfplots@is@datascaled
+				\pgfplotscoordmath{#1}{datascaletrafo noshift}{\pgfmathresult}%
+			\else
+				\pgfplotscoordmath{#1}{tofixed}{\pgfmathresult}%
+			\fi
+			\H=\pgfmathresult pt
+			\pgfplots@isuniformticktrue
+			\aftergroup\pgfplots@isuniformticktrue
+		\else
+			\pgfplotscoordmath{#1}{log}{\pgfmathresult}%
+			\let\Hmacro=\pgfmathresult
+			\H=\pgfmathresult pt
+			%
+			% we don't know it. check it.
+			\pgfplots@isuniformtickfalse
+			\aftergroup\pgfplots@isuniformtickfalse
+		\fi
+		\def\b@pgfplots@ticks@computed@tick@distance@is@final{1}%
+	\fi
+%\message{Got H=\the\H\space(\Hmacro)^^J}%
+	%
+	% OK, now compute MIN/MAX :
+	\ifpgfplots@cur@is@linear
+		% The following code is carried out in floating point
+		% arithmetics because it requires large data ranges.
+		%
+		% I want to compute MIN@new := I*H  where I is chosen
+		% such that MIN = I*H + rest with rest < H.
+		% The problem is the possibly large range of MIN. I
+		% can't work completely in the transformed datarange,
+		% so numbers get too large.
+		%
+		% So, compute I := int( MIN / H )  (integer truncation)
+		% in float arithmetics and then MIN@new := I*H
+		\pgfmathfloatdivide@{\csname pgfplots@#1min@unscaled@as@float\endcsname}{\Hmacro}%
+		\pgfmathfloatint@{\pgfmathresult}%
+		\pgfmathfloatmultiply@{\pgfmathresult}{\Hmacro}%
+		\let\MIN@new=\pgfmathresult
+		%
+		% Now, convert everything into the fixed point data
+		% range:
+		\ifpgfplots@is@datascaled
+			\pgfplotscoordmath{#1}{datascaletrafo}{\MIN@new}%
+			\MIN=\pgfmathresult pt
+		\else
+			\pgfmathfloattofixed\MIN@new
+			\MIN=\pgfmathresult pt
+		\fi
+		%
+		% And, since we have used finite precision, I is most
+		% likely to be large. So: subtract one H. In the worst
+		% case, this produces one tick position too much (but
+		% it won't be printed).
+		\advance\MIN by-\H\relax
+	\else
+		% Now, we want to activate the Tick set
+		%   {lowest, lowest+H, ..., highest}
+		%
+		% Where
+		% 	lowest =  I * log(10) + rest, |rest| < log(10).
+		% this is conceptionally different from the approach for
+		% linear axes, because H = j*log(10).
+		%
+		% remember the original xmin in MINH:
+		\MINH=\MIN
+		%
+		% and compute I and I*log(10) here:
+		\MIN=\pgfplots@loc@log@from@display@log@scale \MIN
+		\edef\pgfmathresult{\pgf@sys@tonumber{\MIN}}%
+		\pgfmathsetcount{\c@pgf@counta}{\pgfmathresult}%
+		\ifdim\MIN<0pt
+			% the truncation rounds TOWARDS 0 which is not what I want.
+			\advance\c@pgf@counta by-1
+		\fi
+		\MIN=\pgfplots@loc@log@to@display@log@scale pt
+		\multiply\MIN by\c@pgf@counta
+		\ifpgfplots@isuniformtick
+		\else
+			% This here is a special case to move the first tick
+			% near the lower axis limit.
+			%
+			% "Near" means either directly above or directly below ymin.
+			%
+			% My application example is as follows:
+			% Let H = 2*log(10).
+			% Furthermore, ymin = 3e-6, ymax= 8e-2. That means we can choose either
+			%    10^{-5}, 10^{-3}, 10^{-1}
+			% or
+			%    10^{-4}, 10^{-2}
+			% as ticks. Well, I prefer the first one.
+			%
+			% HEURISTICS: start as near to ymin as possible!
+			%
+			% We check here if we can come nearer to ymin if we
+			% shift the current tick by log(10):
+			%  if( ymin - I * log(10) < 0.5*H ->  use I+1, that means add log(10).
+			%
+			% that's equivalent to
+			%  2*(ymin - I * log(10)) - H < 0.
+			\advance\MINH by-\MIN
+			\multiply\MINH by2
+			\advance\MINH by-\H
+			%
+			\ifdim\MINH<0pt
+				\advance\MIN \pgfplots@loc@log@to@display@log@scale pt
+			\fi
+		\fi
+	\fi
+	\MINH=\MIN
+	\advance\MINH by\H
+	% Ok, now it can happen that only ONE tick label is placed in
+	% this range.
+	% That's useless, so check for it.
+	%
+	% That's the case if
+	%    MIN < ORIGMIN && MAX < MIN+2 H
+	% MIN < ORIGMIN by construction (ok, MIN <= ORIGMIN by
+	% construction, but I don't care about this case).
+	% So: check only the second condition.
+%\message{Got MIN=\pgf@sys@tonumber\MIN; H=\pgf@sys@tonumber\H; MAX=\pgf@sys@tonumber\MAX.^^J}%
+	\def\pgfplots@tick@returnval@ready{0}%
+	%
+	\pgfplots@tmpa=\MINH
+	\advance\pgfplots@tmpa by\H
+	\ifdim\MAX<\pgfplots@tmpa
+		\if1\b@pgfplots@ticks@computed@tick@distance@is@final
+			% OK, we cannot auto adjust tick labels -- if if there are
+			% too few of them.
+		\else
+			\pgfplots@assign@default@tick@foraxis@autoadjust@result#1%
+		\fi
+	\fi
+	%
+	% assert that we have at least a minimal tick distance.
+	% FIXME :this minimum might actually be too small...
+	\ifdim\MINH=\MIN
+		\pgfplotsthrow{dimension too small in ticks}{#1}\pgfeov%
+		%
+		% make sure that the compilation succeeds:
+		\H=1pt %
+		\advance\MINH by\H
+	\fi
+	%
+	\if0\pgfplots@tick@returnval@ready
+		\xdef\pgfplots@glob@TMPb{\pgf@sys@tonumber{\H}}%
+		\advance\MAX by0.5\H % avoid rounding inaccuracies:
+		\xdef\pgfplots@glob@TMPa{\pgf@sys@tonumber{\MIN},\pgf@sys@tonumber{\MINH},...,\pgf@sys@tonumber{\MAX}}%
+	\fi
+%\message{final H=\the\H; returning \pgfplots@glob@TMPa.}%
+	\expandafter\let\csname pgfplots@#1tick\endcsname=\pgfplots@glob@TMPa
+	\expandafter\let\csname pgfplots@tick@distance@#1\endcsname=\pgfplots@glob@TMPb
+}%
+
+% Computes tick positions using the current axis limits.
+%
+% Parameters:
+% /pgfplots/max space between ticks
+%    Determines the maximum space which is not filled by at least one
+%    tick label (approximate, there is some rounding internally)
+% /pgfplots/try min ticks
+%    see manual
+%
+% Idea:
+% We want ticks at each
+%    { i*H, i in \Z }.
+% Of course, there shouldn't be TOO MUCH ticks.
+%
+% Our heuristics is to set
+%    desirednumticks = round(ACTUAL WIDTH / (max space between ticks) )
+% and generate H = (axis range) / (desirednumticks).
+%
+% Since not all step sizes H look well, restrict H to a set of allowed
+% step sizes such as
+%   { 1, 1/2, 1/5, 1/10 },
+% or, to be more precise:
+%   { 1*10^e, 2*10^e, 5*10^e }
+% -> round to the nearest matching number!
+% This yields H (for example as 2*10^e). Then, compute i*H, i \in \Z
+%
+% The data scaling transformation T(x) makes things more complicated.
+% Now, T(x) = q * x - p  and we need to check for problems with large
+% numbers:
+%   - q* H = ( T(Max) - T(Min) ) / desirednumticks = q * (Max - Min) / desirednumticks.
+%   - Using floating point arithmetics, (Max-Min)/desirednumticks (unscaled!)
+%   is analysed to restrict H to {1*10^e, 2*10^e, 5*10^e}.
+%   - So, we get  q * H  (we can't use the 'p' shift of the affine trafo here).
+%   - The next problem is to compute { I*H, I in \Z } because
+%   	I = trunc( Min / H )  = trunc(  ( T(Min) + p ) / (q*H) ).
+%     This can be seen by Min = I*H + rest   and thus T(Min) = I*q*H + q*rest -p.
+%     The Problem: (T(min)+p ) / (q*H) can be TOO BIG for pgfmath.
+%   -> for the data scaling case, I will use floating point
+%   arithmetics to compute that last step.
+%   I will acquire \pgfplots@[xy]min@unscaled@as@float here.
+%
+%
+%
+%
+%
+% For log-plots,
+% 	H in { j*log(10), j=1,2,3,... }
+% where the usual case should be j = 1.
+%
+% Then, the resulting tick is
+% TICK={MIN,MIN+H,...,MAX}
+% where
+%    MIN = I*H
+% is chosen such that
+%    axis minimum limit = I*H + rest; |rest| < H.
+%
+% Again, log plots follow a slightly different approach: here,
+%   MIN = I * log(10)
+% is chosen such that
+%    axis minimum limit = I*log(10) + rest; |rest| < log(10)
+% while H = j*log(10), j>=1.
+%
+%
+% PRECONDITION:
+% - limits are correct
+% - axis width/height is set correctly
+%
+% POSTCONDITION:
+% - Tick for axis #1 is assigned
+% - \ifpgfplots@determinedefaultvalues@needs@check@uniformtick is set
+%
+% REMARKS:
+% - this algorithms works also if the data range has been transformed
+%   with a LINEAR transformation.
+%   ATTENTION: as of 2008-05-15, the scaling trafo is AFFINE LINEAR.
+%   That means we have to eliminate the 'affine' shifting before the
+%   algorithms works correctly.
+\def\pgfplots@assign@default@tick@foraxis#1{%
+	\begingroup
+	% Shortcut-names:
+	\expandafter\let\expandafter\ifpgfplots@is@datascaled\csname ifpgfplots@apply@datatrafo@#1\endcsname
+	% Attention here: use UNSHIFTET scalings, see remark above
+	\expandafter\let\expandafter\ifpgfplots@cur@is@linear\csname ifpgfplots@#1islinear\endcsname
+	%
+	\let\desirednumticks=\c@pgf@countd
+	\let\Wr=\pgf@xc
+	\Wr=\csname pgfplotspoint#1axislength\endcsname\relax
+	% r = max place without ticks in pt -> choose desirednumticks >= W/r
+	\divide\Wr by\axisdefaulttickwidth\relax
+	\afterassignment\pgfplots@gobble@until@relax
+	\desirednumticks=\the\Wr\relax
+	\advance\desirednumticks by1
+	\ifpgfplots@cur@is@linear
+		\ifnum\axisdefaulttryminticks>\desirednumticks\relax
+			\desirednumticks=\axisdefaulttryminticks\relax
+		\fi
+	\else
+		\ifnum\pgfplots@default@try@minticks@log>\desirednumticks\relax
+			\desirednumticks=\pgfplots@default@try@minticks@log\relax
+		\fi
+		\expandafter\ifx\csname pgfplots@#1tickten\endcsname\pgfutil@empty
+		\else
+			% log plot and tickten-option: provide special processing.
+			\edef\pgfplots@loc@TMPa{\csname pgfplots@#1tickten\endcsname}%
+			\expandafter\pgfplots@compute@tick@times@logten\pgfplots@loc@TMPa\to\pgfplots@glob@TMPa{#1}%
+			\expandafter\let\csname pgfplots@#1tick\endcsname=\pgfplots@glob@TMPa
+		\fi
+	\fi
+	\ifpgfplots@cur@is@linear
+	\else
+		\pgfplotscoordmath{#1}{parsenumber}{1}%
+		\let\pgfplots@loc@TMPa=\pgfmathresult
+		\pgfplotscoordmath{#1}{log from display log}{\pgfplots@loc@TMPa}%
+		\pgfplotscoordmath{#1}{tofixed}{\pgfmathresult}%
+		\let\pgfplots@loc@log@from@display@log@scale=\pgfmathresult
+		%
+		\pgfplotscoordmath{#1}{log to display log}{\pgfplots@loc@TMPa}%
+		\pgfplotscoordmath{#1}{tofixed}{\pgfmathresult}%
+		\let\pgfplots@loc@log@to@display@log@scale=\pgfmathresult
+	\fi
+	%
+	\expandafter\ifx\csname pgfplots@#1tick\endcsname\pgfutil@empty
+		%
+		% OK, automatically compute ticks:
+		\pgfplots@assign@default@tick@foraxis@compute#1%
+		%
+		\expandafter\global\expandafter\let\expandafter\pgfplots@glob@TMPa\csname pgfplots@#1tick\endcsname
+		\expandafter\global\expandafter\let\expandafter\pgfplots@glob@TMPb\csname pgfplots@tick@distance@#1\endcsname
+		\ifpgfplots@isuniformtick
+			\aftergroup\pgfplots@determinedefaultvalues@needs@check@uniformtickfalse
+		\else
+			% in this case, we better check...
+			\aftergroup\pgfplots@determinedefaultvalues@needs@check@uniformticktrue
+		\fi
+	\else
+		\expandafter\global\expandafter\let\expandafter\pgfplots@glob@TMPa\csname pgfplots@#1tick\endcsname
+		\gdef\pgfplots@glob@TMPb{}% will be computed later, in 'check uniform tick'
+		\aftergroup\pgfplots@determinedefaultvalues@needs@check@uniformticktrue
+	\fi
+	\endgroup
+	\expandafter\let\csname pgfplots@#1tick\endcsname=\pgfplots@glob@TMPa
+	\expandafter\let\csname pgfplots@tick@distance@#1\endcsname=\pgfplots@glob@TMPb
+%\message{pgfplots.sty: #1tick set to \csname pgfplots@#1tick\endcsname [#1min=\csname pgfplots@#1min\endcsname, #1max=\csname pgfplots@#1max\endcsname].}%
+}
+
+% Takes the distance between adjacent ticks as floating point number
+% and returns a normalized tick distance.
+%
+%
+% The idea is to get "nice" (human readable) distances instead of
+% strange fractions or real numbers.
+%
+% The result will be assigned to \pgfmathresult (in float).
+%
+% #1 the axis (x or y or z)
+% #2 the unnormalized tick distance computed so far
+%
+% Example:
+% \pgfmathfloatparsenumber{x}{1234}
+% \pgfplots@assign@default@tick@foraxis@normalizetickdist{x}{\pgfmathresult}
+% \pgfmathfloatotfixed\pgfmathresult
+% -->
+%
+%  \pgfmathresult={1200}
+% or something like that.
+\def\pgfplots@assign@default@tick@foraxis@normalizetickdist#1#2{%
+	\begingroup
+	\let\H=\pgf@xb
+	\expandafter\pgfmathfloat@decompose#2\relax\pgfmathfloat@a@S\H\pgfmathfloat@a@E
+	% modify the mantisse:
+	\ifdim\H<2pt
+		\ifdim\H<1.5pt
+			\H=1.0pt
+		\else
+			\H=2.0pt
+		\fi
+	\else
+		\ifdim\H<4.9999pt
+			\ifdim\H<3.5pt
+				\H=2.0pt\relax
+			\else
+				\H=5.0pt\relax
+			\fi
+		\else
+			\ifdim\H<7.5pt
+				\H=5.0pt\relax
+			\else
+				\H=1.0pt\relax
+				\advance\pgfmathfloat@a@E by1
+			\fi
+		\fi
+	\fi
+	\pgfmathfloatcreate{\the\pgfmathfloat@a@S}{\pgf@sys@tonumber{\H}}{\the\pgfmathfloat@a@E}%
+	\pgfmath@smuggleone\pgfmathresult
+	\endgroup
+}%
+
+% Helper method for
+%  \pgfplots@apply@data@scale@trafo@to@options@for
+% #1: the ticks
+% #2: the trafo routine (not necessarily a single macro, but should
+% take one arg)
+% #3: the output macro name
+\long\def\pgfplots@apply@data@scale@trafo@to@user@ticks#1#2\to#3{%
+	\pgfplots@transform@csv@list{#1}\to{#3}\with{%
+		\pgfplotscoordmath{float}{parse}{\pgfmathresult}%
+		%\pgfmathfloatparsenumber{\pgfmathresult}%
+		#2{\pgfmathresult}%
+	}%
+}%
+
+% A helper macro which expects a \foreach-list '#1'. applies some transformation '#3' and stores results into a new csv list '#2'.
+% '#3' is some macro. Its input is '\pgfmathresult', its expected output is '\pgfmathresult'
+\long\def\pgfplots@transform@csv@list#1\to#2\with#3{%
+	\let#2=\pgfutil@empty
+	\begingroup
+	\pgfkeys{/pgf/fpu,/pgf/fpu/output format=float}%
+	% \pgfplotsforeachungrouped is more powerful than \foreach
+	\pgfplotsforeachungrouped \pgfplots@loc@TMPb in {#1} {%
+		\let\pgfmathresult=\pgfplots@loc@TMPb
+		#3%
+		\ifx#2\pgfutil@empty
+			\xdef#2{\pgfmathresult}%
+		\else
+			\xdef#2{#2,\pgfmathresult}%
+		\fi
+	}%
+	\endgroup
+	%
+}%
+
+% Helper method for
+%  \pgfplots@apply@data@scale@trafo@to@options@for
+% #1: the ticks ALREADY IN FLOAT FORMAT
+% #2: the trafo macro name
+% #3: the output macro name
+\long\def\pgfplots@apply@data@scale@trafo@to@user@ticks@isfloat#1#2\to#3{%
+	\let#3=\pgfutil@empty
+	\foreach \pgfplots@loc@TMPb in {#1} {%
+		#2{\pgfplots@loc@TMPb}%
+		\ifx#3\pgfutil@empty
+			\xdef#3{\pgfmathresult}%
+		\else
+			\xdef#3{#3,\pgfmathresult}%
+		\fi
+	}%
+	%
+}%
+
+% Adds a further, temporary anchor to every node which will be
+% processed. The anchor will be named '#2'. It is placed such that
+% 1. the node's center is on a line in direction of the inwards normal
+% vector of the #1 ticklabel axis and the 'at' position of the node,
+% 2. the node does not intrude the axis.
+%
+% In effect, one of the node's standard anchors (north, east, ... )
+% will be placed on the line
+%    \draw[blue,thick,->] (xticklabel cs:0,0) -- (xticklabel cs:1,0);
+%
+% This command is identical to calling
+% \pgfplotsdeclareborderanchorforaxis{#1}{<three-char-string-of-#1-ticklabel-axis>}{#2}
+%
+% @REMARKS:
+% 	- it is -by no means- necessary that any ticks or tick labels are
+% 	drawn or defined for this method.
+% 	- in fact, tick labels use such an anchor (the 'near #1ticklabel'
+% 	anchor is defined in this way)
+\def\pgfplotsdeclareborderanchorforticklabelaxis#1#2{%
+	\pgfplotsdeclareborderanchorforaxis{#1}{\pgfplotsticklabelaxisspec{#1}}{#2}%
+}