$|L'(1,\chi)/L(1,\chi)|$

You may use the local method of Landau with some bounds for L(s,chi) (expressing L'/L in terms of the local zeros, the approximation being controlled by an upper bound for |L(s,chi)| in a slightly larger region). Then, the convexity bound, as given by Rademacher for instance, should give you about what you get here. I would say, this would remove the +15.1, but both results should be rather close.

Tim (Trudgian) is the one who has the most precise explicit Landau formula (well, this is a combination of Caratheodory's inequality [bounds for Log f] with Koebe's inequality for derivatives).

That's all I have in my purse! Pintz-II (series of papers "Elementary methods in L-functions theory"] 1976 may be a good read.

Best, Olivier

Suppose that $\chi(-1)=1$ and that all non-trivial zeros $\beta+i\gamma$of $L(s,\chi)$ with $|\gamma|\le 1/2$ are on the critical line $\beta=1/2$. Recall the Hadamard factorization formula (see Davenport Chapter 12) which gives $$ \frac{L^{\prime}}{L}(s,\chi) = -\frac 12 \log \frac q\pi - \frac 12 \frac{\Gamma^{\prime}}{\Gamma}(s/2) + B(\chi) + \sum_{\rho} \Big( \frac{1}{s-\rho} +\frac{1}{\rho}\Big). \tag{1} $$ The constant $B(\chi)$ is a little tricky to bound, but its real part is well known to equal $-\sum_{\rho} \text{Re} (1/\rho)$. Thus we also have $$ \text{Re} \frac{L^{\prime}}{L}(s,\chi) = -\frac 12\log \frac{q}{\pi} -\frac 12 \frac{\Gamma^{\prime}}{\Gamma}(s/2) + \sum_{\rho} \text{Re} \Big(\frac{1}{s-\rho}\Big). \tag{2} $$

Apply (1) with $s=1$ and $s=3/2$ and subtract. This gives $$ \frac{L'}{L}(1,\chi) - \frac{L'}{L}(3/2,\chi) = \frac 12\Big(\frac{\Gamma^{\prime}}{\Gamma}(3/4) -\frac{\Gamma'}{\Gamma}(1/2) \Big) + \sum_{\rho} \frac{1/2}{(1-\rho)(3/2-\rho)}. $$ Therefore, by the triangle inequality, and a trivial bound for $|L'/L(3/2,\chi)|$ we find $$ \Big|\frac{L'}{L}(1,\chi)\Big| \le -\frac{\zeta'}{\zeta}(3/2) +\frac 12\Big| \frac{\Gamma'}{\Gamma}(3/4) -\frac{\Gamma'}{\Gamma}(1/2)\Big| +\sum_{\rho} \frac{1/2}{|(1-\rho)(3/2-\rho)|}. \tag{3} $$ By assumption $|\gamma|\le 1/2$ implies that $\beta=1/2$. This means that $|1-\rho|\ge 1/2$ always and that $$ |3/2-\rho| \le 1/2 + |1-\rho| \le 2|1-\rho|. $$ Therefore, the sum over zeros in (3) is bounded above by $$ \le \sum_{\rho} \frac{1}{|(3/2-\rho)|^2} \le 2 \sum_{\rho} \text{Re} \frac{1}{3/2-\rho} = \log \frac q{\pi} +\frac{\Gamma'}{\Gamma}(3/4) +2 \text{Re}\frac{L^{\prime}}{L}(3/2,\chi), $$ upon using (2) in the last estimate. Inserting this in (3), and again bounding $L'/L(3/2,\chi)$ trivially, we conclude that $$ \Big|\frac{L'}{L}(1,\chi)\Big| \le \log \frac{q}{\pi} -3 \frac{\zeta'}{\zeta}(3/2) + \frac{\Gamma'}{\Gamma}(3/4) + \frac 12 \Big|\frac{\Gamma'}{\Gamma}(3/4)-\frac{\Gamma'}{\Gamma}(1/2)\Big|. $$ Calculating these constants gave a bound $\le \log q + 2.75$ in this case.

The case when $\chi(-1)=-1$ is similar -- you only need to modify the $\Gamma$-factors. Obviously one can play with the argument with a different $\sigma$ than $3/2$ (chosen more or less arbitrarily). If you don't want to make an assumption on the low lying zeros, you could isolate the contribution of zeros near $1$, and then bound the rest of the zeros as above. Obviously some condition on zeros very near $1$ is necessary to give bounds for $L'/L(1,\chi)$, but as can be seen from (2), a general one-sided bound is given by $$ -\text{Re} \frac{L'}{L}(1,\chi) \le \frac 12 \log \frac{q}{\pi} +\frac 12\frac{\Gamma'}{\Gamma}(1/2). $$

Few remarks I already wrote privately to the OP.

First, my joint paper with Alessandro Zaccagnini mentioned by Pace Nielsen was published in a slightly different form in 2009 on Experimental Mathematics 19 (3), 279-284 (with an interesting final section by Karl K. Norton). The goal there was to get an extremely good estimate (at least 100 digits) for the Meissel-Mertens constants in arithmetic progressions p\equiv a \mod q, for every (a,q)=1 and q up to a certain bound (q\le 100, in this case). The results up to 20 digits and for every q\le 300 are available at my web page dedicated to this paper (the link is the one written in the paper). Please remark that such computations were performed about ten years ago and it is clear that now we can beat such results using the same algorithm on the more efficient PCs we have now at our disposal. Possible improvements there can be obtained using the more efficient implementations we have now to compute the needed values of Dirichlet L- functions at positive integers (for example the one we have now in Pari/gp). All these remarks essentially apply also to the computation performed for another paper of Zaccagnini and myself on the Mertens’ product in arithmetic progressions published in 2007 on Mathematics of Computation 78 (265), 315-326.

Second: I think that the faster approach we have to compute the value at 1 of the logarithmic derivative of Dirichlet L-functions mod q is, at least for q prime, the one described in a paper by Ford, Luca and Moree and improved in my arxiv preprint https://arxiv.org/abs/1903.05487 (see the references there for the paper of Ford, Luca and Moree). In such a preprint of mine I also tried to analyze the cost of directly compute the mentioned quantity using the Pari/gp function lfun; it seems that this approach is slower if compared with the other ones presented there (I never tried to compare the speediness of ARB and Pari/gp, though). The algorithms used to compute the log-derivative at 1 are based on classical formulae that link L(1,\chi) and L’(1, chi) to some special functions evaluated at some rational numbers in (0,1). For more details please refer to my preprint mentioned above. Final remark on this point: I have already precomputed and stored the special functions values needed to get the Euler-Kronecker constants for every q prime up to 5*10^5 and I am now performing the needed computations for reaching 10^6. This means that with suitable modifications to my programs I can get the data about the logarithmic derivative evaluated at 1 for every Dirichlet character mod q, q prime up to 5*10^5, by now, and in a month or so (I hope...) for every prime q up to 10^6. This clearly means to divert some computational power (and time to write and test the new programs) from the project I am now working on; if there is no hurry in getting such data, I will insert them in the next arxiv version of my work on the Euler-Kronecker constants (hoping that such a version will be the final one...)

Third: about computing the values of L(1,\chi), chi mod q, q prime: a similar approach used in my preprint mentioned in the second item of this list was used for performing the needed computation in an ongoing research project I am now working on (it is a collaboration with Pieter Moree, Sumaia Saad Eddin ad Alisa Sedunova). It is about the Kummer ratio for the class number of cyclotomic fields (please see our preprint https://arxiv.org/abs/1908.01152). Directly using the Pari/gp function lfun seems to be slower than computing L(1,\chi) using some classical formulae that make use of some special functions values. In this case too I did not compare the speediness of Pari/gp and ARB. From a computational point of view, this problem is simpler than the one on the Euler-Kronecker constants because the special function involved in the needed formulae here is the Psi function and, using a decimation in frequency strategy, we can in fact use the values of the cotangent function. An alternative algorithm uses the first chi-Bernoulli number instead. Again, please refer to our preprint mentioned above for more details.