From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 69676 invoked by alias); 8 May 2018 19:08:19 -0000 Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org Received: (qmail 69664 invoked by uid 89); 8 May 2018 19:08:18 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-1.9 required=5.0 tests=BAYES_00,RCVD_IN_DNSWL_NONE,SPF_HELO_PASS,SPF_PASS autolearn=ham version=3.3.2 spammy=Rev, siam, journal, Greetings X-HELO: NAM03-BY2-obe.outbound.protection.outlook.com Received: from mail-by2nam03on0120.outbound.protection.outlook.com (HELO NAM03-BY2-obe.outbound.protection.outlook.com) (104.47.42.120) by sourceware.org (qpsmtpd/0.93/v0.84-503-g423c35a) with ESMTP; Tue, 08 May 2018 19:08:17 +0000 Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=patrick.alken@colorado.edu; Received: from palken-co-ll.ngdc.noaa.gov (140.172.179.43) by DM2PR0301MB1247.namprd03.prod.outlook.com (2a01:111:e400:5020::24) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_256_CBC_SHA384_P256) id 15.20.735.18; Tue, 8 May 2018 19:08:14 +0000 Subject: Re: functions related to the Hurwitz zeta function To: gsl-discuss@sourceware.org References: From: Patrick Alken Message-ID: Date: Tue, 08 May 2018 19:08:00 -0000 User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.6.0 MIME-Version: 1.0 In-Reply-To: Content-Type: text/plain; charset=utf-8; format=flowed Content-Transfer-Encoding: 8bit X-ClientProxiedBy: DM5PR06CA0028.namprd06.prod.outlook.com (2603:10b6:3:5d::14) To DM2PR0301MB1247.namprd03.prod.outlook.com (2a01:111:e400:5020::24) X-MS-PublicTrafficType: Email X-Microsoft-Antispam: UriScan:;BCL:0;PCL:0;RULEID:(7020095)(4652020)(5600026)(4534165)(4627221)(201703031133081)(201702281549075)(2017052603328)(7153060)(7193020);SRVR:DM2PR0301MB1247; X-Microsoft-Exchange-Diagnostics: 1;DM2PR0301MB1247;3:SeoLyC4gNqHtTarug+uWugokUTOrGItbXAQlWRw0tITzE4WWhjV6FH5sSik07M2CJxgm5L+sFEIi5OebHLZd8xmSkH/qJh5e0/XS9XKAfFpaB9v9+PqBEDxQJnJWGE7j6++zosp/LsIVSto3PEjEMy2n6pNPzhSLT4SdHZDLVNBRnxdvbqWZI+FSUBG/7dT+rc40RjlQkSOjSWOMYza9xzKOVTJXyJ6rqh2JqUDMcAbrmGOGvWvo/iUPvgnFGyFn;25:SvrwWV31gRl/2CIRBoukHD1OrQp3uV05Aq+08lW62o4LOnGB0NwMTK9bY7GizIgtmaMrIykZFZ1jihYeREeIGAWJrUiZRAWQQxNMH40qbEsF0EAYMintHrBTTUBFHT5Q5OAiod0bquiBrSuj2nhRfw+7jcKQZagUStRemBcK01Ic+THZYVOHWFcZk0hLOViPs+vfbBaNroPF1qdIYgABIXfeu6ukiKfLY1I7eMoIqHMGBnNpKXH2p6J0gUdadyG66qXANY6wWNRIhmEVOgOlAerK8IKY2knysp3s/juEvcf3UnFTnqD44agWvguHvtTKRNhvn/JdJ+tL+K/P5/csCg==;31:0r0AVxu+sM1iKrcUvNH/aW3kG5tnaZQTO1ei4LOOeLwHDsBJFFcH2kPCwYWIJhdOmWANMgFpUAmckqhSdqdllHVNZFsf91PFpaKRwvaYpQp5jHvvvIP7Zg60NQQIKMhawEnekJf3F1/psG3T1WW9j8F0Fs6tmxM5UFN71ML2vkJJiot5f4giAivfbUvAWKz8YiPGtPCCz3SSiOIp3qqRXcF2e2PJwIlXIkGFU7hrT74= X-MS-TrafficTypeDiagnostic: DM2PR0301MB1247: X-Microsoft-Exchange-Diagnostics: 1;DM2PR0301MB1247;20:OOaKIDqgyYro+p6hjwzCFYpxOkEcnyCAAzqtEp00cZWc8Svl9Ls1SFBOSICb5GW1guQSv6VDct5Km5Cm1cd+zMJQs0vIrrheHNcOBR+4kxPQz7eX81H4srEBavOGpWSZoZQRg1pP/jTGGT/TEkzzODXmMidLsr/M5yotKYOCz+i1Lh+DTOAAPHme3+e7y9u90G9s5rQJLUC8KhnHTC/wFdxBsW65yInuyMDcOkDH5oJQkYxY2CGYOXyg+4vUJi6iyj9ZZOtILpHCMOdB7Y3p5TAUpnsasc9Cu9Qq12T7yjJZdZXZKzlzbYhfWEjCGdTuj2+NU5W1ue8oFrFojab5jaGedxI6ujs7B0XyU4BCkv5d89QCwfCNo5XOTS8mbXDXChKLA7GRW5Xq32h3JxOfDXJAhZh++nvTsOb6/c4UJkn0zEm/zfMED++LC6JdUZYsgi0UdkZfj3pOTXGDssOJOJMQWQ5synEw8018L9T6iFdsWt5D+XvfqLFagxSNqaax;4:eYawn5APPPt16CvdF8ONlxhfUdbpvoRKzS5bCaPvMjp8wo4A4TUbxCcZka5ouGMq+BrT8j6CXRm1ziocCNcNCMxMQ0lS8iXqfpTrkA3W76aDuEm6p5l1HYRP35Qa/CqBJIi+ewB52uGAr3z8WG5iLz8XS6Aa4rvmTNZN19W9Sb6sd2iBe1qYnFrqFJTvavGkXscqHOJaU1k6vTCWSHjsEN30Gtt3S5ht/mN47toBUIlt26ZgDAL9Zq9k5zptejYzc2fW5qypJMFdlv1+ufAiqg== X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-Test: UriScan:; X-MS-Exchange-SenderADCheck: 1 X-Exchange-Antispam-Report-CFA-Test: BCL:0;PCL:0;RULEID:(6040522)(2401047)(5005006)(8121501046)(10201501046)(93006095)(93001095)(3002001)(3231254)(944501410)(52105095)(149027)(150027)(6041310)(201703131423095)(201702281529075)(201702281528075)(20161123555045)(201703061421075)(201703061406153)(20161123560045)(20161123564045)(20161123558120)(20161123562045)(6072148)(201708071742011);SRVR:DM2PR0301MB1247;BCL:0;PCL:0;RULEID:;SRVR:DM2PR0301MB1247; X-Forefront-PRVS: 0666E15D35 X-Forefront-Antispam-Report: SFV:NSPM;SFS:(10019020)(39860400002)(396003)(376002)(366004)(346002)(39380400002)(199004)(189003)(51444003)(67846002)(2486003)(6306002)(6506007)(106356001)(2351001)(36756003)(52146003)(16526019)(47776003)(53546011)(2906002)(6916009)(229853002)(65806001)(53936002)(65956001)(76176011)(31686004)(66066001)(786003)(31696002)(105586002)(316002)(11346002)(6512007)(6666003)(386003)(75432002)(52116002)(2361001)(64126003)(97736004)(8936002)(5660300001)(476003)(6486002)(2616005)(478600001)(446003)(486006)(6116002)(7736002)(88552002)(58126008)(81166006)(3846002)(186003)(23676004)(305945005)(2870700001)(956004)(50466002)(6246003)(81156014)(68736007)(25786009)(8676002)(65826007)(26005);DIR:OUT;SFP:1102;SCL:1;SRVR:DM2PR0301MB1247;H:palken-co-ll.ngdc.noaa.gov;FPR:;SPF:None;LANG:en;PTR:InfoNoRecords;A:1;MX:1; Received-SPF: None (protection.outlook.com: colorado.edu does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: =?utf-8?B?MTtETTJQUjAzMDFNQjEyNDc7MjM6YitGOCt6Q2NHUnZyOFMvcG5BUlZUY2tw?= =?utf-8?B?bm1kYmFLdm5YaXJBOXRxNDFWY2tDbUREdzVEWksza25uNCsyS3BTamIzMEhK?= =?utf-8?B?QkZvSU9rQ1E3VWdSYmxnZ1FoRUt0dXdVOHR6Y2pzRWkxV085U0FnK2w2Ukdq?= =?utf-8?B?ajkzNXo0M3Z1dFA3MUlqeE1qT3M0cnZidTFDRjV4cTVnOVowNitHVVh5UlVl?= =?utf-8?B?bWNLNmdYdktMSWI5aWQzZHBCeDNIS2ZZK2J5dUoyeUJ0MURRVnFLc0txUzlU?= =?utf-8?B?UUNJNUZqN1VxTlVKVTBtcHUzVmhtMjE4Nk9wU3IzcVFvUEl4a0w4aUtTZzdr?= =?utf-8?B?Uld4YVRiN0lMUzQ5ckducW92TktRaXhyVEl4aXBuYXl2S1k1dkw4eGlESkhB?= =?utf-8?B?eTJkNTdnMFJ5aGd5VmNwR2tiNVg2WEVSYUVHQ2wrWVNpaSsyRGZ1VlJLc2Ry?= =?utf-8?B?OTNJOHJ2N0M2anZHTmlpdTgwM3JiRUFCR2VvSWp0TEFNWXpYNVBGc1JWSW5a?= =?utf-8?B?eFdaUTRpSWlYeU1rSnJBNloxYVhDYllNQjVoZTdUeEN6SUVwWHRPRXRLeEJV?= =?utf-8?B?Vk1yRUpVUEJJbk9NNDYrSzh4T2ZhaWtranNVRmFNNEV5VjhSOHVhTENTeTl6?= =?utf-8?B?SWU3MkZ5ZC9md2t2REYvQlhVMFlPejdNeExYLzRkdkpEQTZOZmwvKzB0UE45?= =?utf-8?B?akJnYUhNb2pzclpyb05vd3BmSHEwZVVDUUlCRlVDc1lxOWVtM2hMZ2NqWW9o?= =?utf-8?B?N1J5b3dmektrNGtRYyt1aGZ5ZnoraDRqOXFkbHZpOExFc0VwMG1ITy90YVJz?= =?utf-8?B?Z0hYUDNxUWhoTDN4MmZCSDZjT01Cd0c2OHkzV0lBQWdTZmRLNVc3OTloQkVD?= =?utf-8?B?WUZDVEdDUWhhRmhqNC9YMU5BQitBclJhZlZjcHNzVXV0MXVSK2xNdVhRbTRO?= =?utf-8?B?SUs4L1dVS0M5UEZxSDk5ZUkrMjVaRWxwNzZmY3hBWjR5R3p4VWdMeVhmNWVj?= =?utf-8?B?RmU0VHVCcEUrZlFKbjhadnVEcXJadFo4M1htbENUMGhtMEJva1QrcnE0YlpC?= =?utf-8?B?dlFPRE03RDViQzlBNGxFTGVDcFNiY25FdU5zWFlGMGx3cWg0MkVRT2dtWlpj?= =?utf-8?B?My9JRmFna3hXcGs2KzJQUXVyVGp3M3RGbXQ1K2U1VzVhRVVNelRPRFl3OC83?= =?utf-8?B?bk4rM2NrWWRLM3N2UzVzeGFLUEhhUUwxNkFZZVVnTjh2eGJiN0pqWTZRVWtj?= =?utf-8?B?NVA2cGo3ZXBvVzRpRzZxb1lONUdNcEZXNURoVmJ3OVVnZnVweWFvV2dlMTZL?= =?utf-8?B?ZWNWMDR5VXJSY1FXemxYdzdqVFRvSVRHem9saW5pZHhlOHJZc21uNUNyTW00?= =?utf-8?B?TzhrK00vVUo1bnZrNllVZmptM3NDU01qYjgvVHNpMjFhb2hna25sbjdnTFJD?= =?utf-8?B?cm1STXBkMExaS0FmTitjb0hYa3ZYdXFsNmRxUERNck5BU0VaalI2MHBheHNm?= =?utf-8?B?WTJIa2Z2NFIyYXJ0ZVVnWUM0clo3RjAvV25OSFg2cmxmZ2hJUlBCbXZZdW1Z?= =?utf-8?B?dHhXZFdGRzBwMnRTTndQMnpKeDdiZjA0WFhLVVFobnVucldVNDBNSDY2T0dh?= =?utf-8?B?ODh4cG0yamJOMWNmUjA4LzJ6endwRHRONGQ4Q1dpVXFiTWtCY3U5T0pwNC9C?= =?utf-8?B?Y1VvaldkdWt1OEFhOU0vbXkwL1JEYlVMam1DSlZqYlBhWmVWQkU5OHlTWE4r?= =?utf-8?B?RWZpeDhHY2k0bDRyRDd2b2FpMTY2Y1g2MGJmemJFY1FvenRwTEdiMGFUbkVQ?= =?utf-8?B?NG1xL3lwQWdqN05ScW1aUEhKb2N6Unp2ZlVHMGtzdWljNVJkcm5VN1lLaFl1?= =?utf-8?B?cm1FKzdrUjZSZUF1Q2x5a2lLM1Z4NWhlN3Z1Vm9RUUQzVWd1cVZqTTJXNys0?= =?utf-8?B?TDR0MU54dHd3djdRZnY2TVp2UWdsUmpjWit6UlNMT09iMWk5M0lyWVJCZ3k3?= =?utf-8?B?QnI3aTlReFNXWFQwcVRFMWZJeVFnLzRMU3NzZ2lnN1Q0Myt1YVJhUjBqeXEy?= =?utf-8?Q?mN6btg=3D?= X-Microsoft-Antispam-Message-Info: gKcpO4ymLP320At9fz0NfCdL7QKpjDIRrOX8tFDaX/ZWhyobLFvPGekp/Z22HaBvy5WZ0aPkpPPgBIcMA1eksVKDrm1QL4A9ruB1Io9B7lYekWBBmKYVpTwVqUN/CctYGz0BDtApRYG3eQk9fuAFBbGQlLXqrPgDuhV9SQ9SRmFBaELL6muYYkjM0oFRTHUk X-Microsoft-Exchange-Diagnostics: 1;DM2PR0301MB1247;6:FZ9DRISL1DWQgaR1lZAoT/bNCISo3V+AGyjuAhsPMp/9VFu+k1UENbcVz+AwmWRb/BfX3uGG2O943Cxh1dmP53bgLR3LVEg1kkYhTAGX02Ck1XCUtErqHY9/ISSuM/q69RBH3YbbKmZmQXkk/TQhN00J3f+SIxgFCaQRgVVLXhcxf5zFU1GfWOJN+ghZYmcYeM2Wzj8a7r2hds3WLV/0q/vxvxSoHTIOTPkNR6C+ZjrO9vD9SVsKKLAUN/xIS5oQpfsbNmO5C6GO1RL3Yja7NQ0E0GDbsYgV227hek/K/LgA7G9l/MIVTqYZ0xg8GMSX2uTRwZ4Pu+ECAGjn/8mHqG+S7hWDVYu5/YRShgthnNO72pO/ZcqL1ogBH8LPK6Oy6/lPp0PxS2cbJxgAxAwDAhBkPKhuWabMlSGfUiH6ggZQ42aq97YNc0YU1BRP9tNaGpUERyYrzG5NsoKzvi2oVw==;5:eA89QVFUvBsAC2f4etXU9lPwJesCBFzGR8O00y9L54mOAYAH97ihuJVHmJc1CQnn/+UqS2KghSb0LlX+j0wYniwjqZH6V9lP9I1HFFNpSYxoW4rBU/Snc0mV6iI+ygMkxLtyKWd8iv9h1s1QDEMO9Vkv/am3d2V8YOjndEKwSgg=;24:Gwk+1UZsUMS1quFyFKnAnDN2EUU2HDHfQ1CCvyhYhgtf0PXh+PLY4aVLvud6/oYxgkP58KSa+s0Kcm/Tcidr0MVYtTMi04W380Sk7LPk4+Y= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-Microsoft-Exchange-Diagnostics: 1;DM2PR0301MB1247;7:LAt3voJeutT6Cem5m4kaAanYpIGko6Svrzo/lUoLfDaqgMwd9Fd7vXiUvZtzcvpOVmSZXFOpTRk2rgTaBPzwzo9h7GMuIt+LI8JAHA+HgAb40+cqozlAWY4mESrJiwY1gjiwhKYekQfYmwwTTDl0h7xV4+B8D7aX7Jj/eL7bMXmI3lAVRLija3ETcGYHp/6hS0KIgX2cZp0n2cXSiUyp/YgwaL8Ddhn76+LGr7QdIFA7URwMdCxfEsyjqQZPeIo/ X-MS-Office365-Filtering-Correlation-Id: f522cb46-c6f2-4732-366d-08d5b5170ce5 X-OriginatorOrg: colorado.edu X-MS-Exchange-CrossTenant-OriginalArrivalTime: 08 May 2018 19:08:14.6923 (UTC) X-MS-Exchange-CrossTenant-Network-Message-Id: f522cb46-c6f2-4732-366d-08d5b5170ce5 X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-CrossTenant-Id: 3ded8b1b-070d-4629-82e4-c0b019f46057 X-MS-Exchange-Transport-CrossTenantHeadersStamped: DM2PR0301MB1247 X-SW-Source: 2018-q2/txt/msg00005.txt.bz2 Hello Jerome,   Sorry for the late reply. Could you prepare a patch against the latest master branch of the GSL git with your additions? Then I can take a look to see how this might fit in. It would be best to follow the API of the existing special functions in GSL as closely as possible. Thanks, Patrick On 04/03/2018 09:35 PM, Jerome BENOIT wrote: > Greetings, > > I have write within the GSL framework, > some functions related to the Hurwitz zeta function: > first and second derivatives, its logarithm and its logarithmic derivative. > > These functions are needed to fit properly discrete power law distributions [1]. > > I think that there is a place for them (along a plfit material) in GSL: > what is the procedure to follow ton bring them in GSL ? > > > Best regards, > Jerome > > > [1] > @ARTICLE{AClausetCRShaliziMEJNewman2009, > author = {Aaron Clauset and Cosma Rohilla Shalizi and Mark E. J. Newman}, > title = {Power-Law Distributions in Empirical Data}, > journal = "SIAM Rev.", > year = {2009}, > volume = {51}, > number = {4}, > pages = {661--703}, > numpages = {43}, > archivePrefix = ARXIV, > eprint = {0706.1062}, > doi = {10.1137/070710111}, > website = {http://tuvalu.santafe.edu/~aaronc/powerlaws/}, > }