Andrew Cooke | Contents | RSS | Previous

C[omp]ute

Welcome to my blog, which was once a mailing list of the same name and is still generated by mail. Please reply via the "comment" links.

Always interested in offers/projects/new ideas. Eclectic experience in fields like: numerical computing; Python web; Java enterprise; functional languages; GPGPU; SQL databases; etc. Based in Santiago, Chile; telecommute worldwide. CV; email.

© 2006-2015 Andrew Cooke (site) / post authors (content).

Estimating Power On Climbs

From: andrew cooke <andrew@...>

Date: Thu, 20 Jan 2022 19:31:49 -0300

I was slowly pedaling up a hill today, staring at my cycle computer
to avoid looking at the endless climb in front of me, when I realized
that you can easily estimate your power from the grade and speed
(km/h).

Me + my bike weight about 72kg.  Actually a bit more (I need to lose a
kilo or two), but 72 turns out to be a nice number in what follows.

72kg ascending vertically at 1m/s requires a power of 720W (because
g=10, more or less).

A speed of 1km/h is 1000m/3600s or 1/3.6 m/s.

Small angle approx means that if you are riding at a speed S(km/h) on
a gradient of G then your vertical speed is SG/100 km/h or SG/360 m/s.

So the power required (at 72kg) is SGx720/360 W or 2SG W.

For example, at 10km/h on a 5% grade your power is 2x10x5 = 100W.

Now this is incorrect if you're moving fast enough, or riding up a
shallow enough slope, that air resistance is important.  But for
climbing it's fine.

Obviously if you+bike weighs more or less than 72kg you need to
adjust accordingly.

Andrew

Permalink | Comment on this post

Previous Entries

For comments, see relevant pages (permalinks).

Applying Azure B2C Authentication To Function Apps

From: andrew cooke <andrew@...>

Date: Mon, 17 Jan 2022 20:09:08 -0300

Solving this has taken a long time (months) going back + forth with
Microsoft support.  Since I was stuck because I couldn't find this
info on the web I thought I should put something out there.


This is not a general explanation of how to enable B2C authentication
between a SPA (Single Page Application - JavaScript in a web page) and
a Function App.  Instead, what I will focus on is converting the
project at
https://github.com/Azure-Samples/ms-identity-b2c-javascript-spa.git to
work with a Function App that you have deployed.  If you can get that
working, you can carry everything across to your application.


So this is what you need to do:

1 - Clone
https://github.com/Azure-Samples/ms-identity-b2c-javascript-spa.git
and get it working as is.  Understand it.

2 - Deploy your application that uses a Function App (currently
without authentication), make sure it works, and note what the result
is if you call your Function App directly from a browser (in my case,
for example, I had a Function App that provided an access token to
Azure Maps - if I called it by hand, I got a token).

3 - Configure B2C.  There are three parts that need to be connected:

  3a - Define a User Flow.  Go to B2C in the Portal and select "User
  Flows" (left menu) then "New user flow" (top bar), "Sign up and sign
  in" (option box), "Recommended" (popup) and "Create" (button).
  Provide a name (eg "demoflow"), select "Email Signup" (checkbox) and
  "Create" (button).  Leave everything else as defaults.

  3b - Register the Function App.  Go to B2C in the Portal and select
  "App registrations" (left menu), "New registration" (top bar).
  Provide a name (eg "myfnapp"), set "Select a platform" as "Web" then
  select "Register".  Leave everything else as defaults.

  3c - Register the SPA.  Go to B2C in the Portal and select "App
  registrations" (left menu), "New registration" (top bar).  Provide a
  name (eg "msdemo"), set "Select a platform" as "SPA" with a redirect
  URI of "http://localhost:6420" then select "Register".  Leave
  everything else as defaults.

  3d - Connect things.  Go to B2C in the Portal and select "App
  registrations" (left menu), select "myfnapp", select "Expose an API"
  (left menu), select "Add a scope", save the default URI (popup on
  right), provide a name for the scope (eg "access") and provide some
  text for display (anything will do), then select "Add Scope"
  (button).

  Go to B2C in the Portal and select "App registrations" (left menu),
  select "msdemo", select "Add a permission", then "My APIs" (top
  menu), then "myfnapp", select "access" and then "Add permissions".
  In the "Configured permissions" page (displayed after the above),
  select "Grant admin consent" and "Yes".

4 - Enable authentication for the Function App.  Find the Function App
in the Portal, select "Authentication" (left menu), select "Add
identity provider", select "Microsoft", select "Pick an existing app
registration in this directory" and then select "myfnapp".  Choose
HTTP 401 for "Unauthenticated requests" and select "Add".

CORS also needs to be enabled for the Function App.  Select "CORS"
(left menu) and add "http://localhost:6420" then select "Save".

If you call your Function App from a browser now it will give an
error.

5 - Modify the MS demo.  All the above is standard, really; it's this
section that is critical.  The files that need altering are all in the
App directory and are all pure data structures (no code - the app is
nicely structured).

  5a - authConfig.js

  Change "clientID" to the "Application (client) ID" for the "msdemo"
  registration you defined above (it's displayed in the Portal in the
  App registration for that name).

  5b - policies.js

  Remove the "editProfile" entries because we're not using them.

  Change the "authorityDomain" to be "sts.windows.net" (more exactly,
  use the domain in the "Issuer URL" if you go to "Authentication" for
  the Function App and then select "Edit" for the provider).

  Change "fabrikamb2c" in "authentication" to match whatever you call
  your B2C.

  Change "B2C_1_susi_reset_v2" in the "authority" to be the user flow
  defined above ("B2C_1_demoflow")

  5c - apiConfig.js

  Change "b2cScopes" to be the default URI created for the "access"
  scope above (go to B2C, "App registrations", "myfnapp", "Expose an
  API" and copy the URI)

  Change "webApi" to be the URL of your Function App.

6 - Run the demo.

  > npm install
  > npm start

Open http://localhost:6420 and enable the dev console so you can see
network requests.  Sign-in and call the API.  In a perfect world, it
should work (should make successful calls; the web app itself won't do
much because the function being called is not what is expected).  You
should be able to see the expected return value from the API in the
network panel.

If you get an HTTP 401 error when calling the API then you may be able
to view the actual error from the MS auth system by doing the
following (this is deep magic, but it works):

Select the Function App in the Portal.  Select "New Support Request"
(left menu) (don't worry, we are not going to actually submit a
request).  Enter "authentication error" as "Summary", select
"Authentication and Authorization" as "Problem Type", select "Azure
App Service built-in authentication" as "Problem subtype".  Select
"Next:Solutions".  When something is displayed, scroll down and see if
there are any errors.

If you don't see any errors, wait a little and try again - it seems to
take 30min or so for the errors to propagate.


The critical step in all the above was learning how to see the errors.
From those I managed to infer setting the authorityDomain in 5b above.

Have fun,
Andrew

Permalink

[Bike] Gearing On The Back Of An Envelope

From: andrew cooke <andrew@...>

Date: Mon, 4 Oct 2021 21:28:20 -0300

Uphill

If you + bike weigh 100kg and you can put out a steady 100W then you can
ascend vertically at 0.1m/s (work against gravity if g=10m/s^2).

On a 10% incline (steep, but possible where I ride in the Andes foothills)
that means you can move forwards at 1m/s.

The circumference of a 700C wheel is about 2m and the slowest comfortable
cadence while seated is around 60rpm.  So with 1:1 gearing that's 2m/s.

So you need 1:2 gearing for the climb.

That kind of gearing is pushing the low end of available MTB gears.  In other
words, pretty much no-one has gearing that low.  Pretty much all entry-level
riders would have real problems climbing a 10% grade for any distance.


Downhill

Air resistance increases as v^3.  That has a sharp "bend" - either it's not
important or it's a brick wall.  So basically everyone (whatever their power
output) is limited to a similar top speed.  From experience it's roughly
60km/h.

Faster / steeper than that you need to get seriously aero and pedalling
efficiently isn't really an option.

You can probably manage 100rpm if you have to, in your top gear, all out.
With a 2m circumference wheel that's 200m/min.

60km/h is 1km/min.  So you need 5:1 gearing for the descent (to go from
200m/min to 1,000 m/min).

This is consistent with the top end of traditional "race" gearing (maybe 53
teeth at the front and 11 at the rear).


Gear Range

The above implies a gear range from 1:2 to 5:1.  That's a factor of 10, or
1,000%.

In practice, pro riders can put out around 4x the power / weight used above.
Compared to inexperienced riders they are absolute monsters.  That shifts the
lower gears needed up by a factor of 4.  So instead of a 1,000% range they
need a 250% range, which is well within typical "race bike" gearing.

Non-pro riders (and people with loaded touring bikes) that have the technical
understanding compromise at both ends.  At the top, 4:1 is generally enough,
and at the bottom they can live with maybe 1:1.5.  That gives a 600% range
which is possile with a Rohloff IG hub, or a frankenstien mix of road and MTB
components (as on my own road bike).

Currently fashionable 1x road / gravel can reach maybe 500% range.  Exactly
where that is "centred" depends on the front ring - there will likely be
significant compromise at both the high and low ends.

Andrew

Permalink

[Computing] Okular and Postscript in OpenSuse

From: andrew cooke <andrew@...>

Date: Thu, 9 Sep 2021 10:13:20 -0300

Wow.  PS is so old-fashioned now that support in okular is in a separate
package (okular-spectre) that is not installed by default.

Andrew

Permalink

[Computing] Fail2Ban on OpenSuse Leap 15.3 (NFTables)

From: andrew cooke <andrew@...>

Date: Fri, 30 Jul 2021 09:58:50 -0400

Couldn't find this documented anywhere, but it works just fine.

  localhost:/etc/fail2ban # cat jail.local

  [DEFAULT]
  bantime = 1h
  maxretry = 3

  [sshd]
  enabled = true
  backend = systemd
  banaction = nftables

  localhost:/etc/fail2ban # cat fail2ban.local

  [DEFAULT]
  loglevel = INFO

Typical log output:

  localhost:/var/log # tail fail2ban.log
  2021-07-30 08:19:02,184 fail2ban.filter         [16653]: INFO    [sshd] Found 49.88.112.69 - 2021-07-30 08:19:01
  2021-07-30 08:19:02,249 fail2ban.actions        [16653]: NOTICE  [sshd] Ban 49.88.112.69
  2021-07-30 08:28:00,927 fail2ban.actions        [16653]: NOTICE  [sshd] Unban 205.185.125.109
  2021-07-30 09:05:49,815 fail2ban.actions        [16653]: NOTICE  [sshd] Unban 77.227.158.129
  2021-07-30 09:19:02,048 fail2ban.actions        [16653]: NOTICE  [sshd] Unban 49.88.112.69
  2021-07-30 09:23:52,935 fail2ban.filter         [16653]: INFO    [sshd] Found 205.185.127.25 - 2021-07-30 09:23:52
  2021-07-30 09:23:54,934 fail2ban.filter         [16653]: INFO    [sshd] Found 205.185.127.25 - 2021-07-30 09:23:54
  2021-07-30 09:24:52,933 fail2ban.filter         [16653]: INFO    [sshd] Found 205.185.127.25 - 2021-07-30 09:24:52
  2021-07-30 09:24:53,125 fail2ban.actions        [16653]: NOTICE  [sshd] Ban 205.185.127.25
  2021-07-30 09:24:54,434 fail2ban.filter         [16653]: INFO    [sshd] Found 205.185.127.25 - 2021-07-30 09:24:54

All under systemd control.

Andrew

Permalink

[Cycling, Computing] Power Calculation and Brakes

From: andrew cooke <andrew@...>

Date: Sat, 24 Apr 2021 16:20:34 -0400

It's common (eg Strava) to estimate a cyclist's power output from their GPS
and elevation (barometric) data.

The basic idea is fairly simple: we can calculate the power required for the
bicycle to behave as recorded at each point in time using mass / gravity /
height, CdA / speed, and rolling resistance / speed.

This ignores wind, which makes the effective speed for wind resistance (CdA)
different from the measured speed, and so introduces errors.

Also ignored, but not commonly discussed, is braking.

If the cyclist does not use the brakes then all we have discussed so far is
correct.

If the cyclist brakes *and* pedals then we will underestimate their power
output because they are expending additional energy we do not "see" in the
data (assuming that there is no power meter data).

If the cyclist brakes when not pedalling then we *should* (if we have
sufficient resolution, perfect parameters, etc) see a negative input power.

In practice, when you estimate power, you do see negative spikes.  I had
assumed that these were noise, but the above suggests that it is also
reasonable to interpret them as braking.

In conclusion, then, braking may be visible in the output from modelling and
is not inherently a source of uncertainty *unless* the cyclist brakes when
also pedalling.  That exception is more likely than it sounds since "when
also" does not, in practice, mean (only) concurrently, but also at any time
within the "time step" used in the calculation (eg 10s).

Andrew

Permalink

[Hardware, Computing] Amazing Pockit Computer

From: andrew cooke <andrew@...>

Date: Thu, 18 Mar 2021 21:02:15 -0300

I am not sure I really believe this is real.  https://pockit.ai/

Watch the top video in the timeline.

Andrew

Permalink

Bullying

From: andrew cooke <andrew@...>

Date: Thu, 18 Mar 2021 04:59:59 -0300

I was buliied recently.  I am a 53 year old, white, male, successful
professional.  It was a bit weird.

I don't want to give details, because they would identify the people involved,
which would be complicated socially.  But I was harassed by email over a
couple of weeks by two people, who also contacted my partner, Paulina.

At first it just seemed like odd behaviour, but when I asked them to stop, and
they did not, I recognised the similarities with episodes in childhood.  Like
then I was confused but resistant (so you might call it an attempt at
bullying, rather than the real thing, because I never felt intimidated).  Like
then, the people involved were - I think (I can find no other explanation) -
resentful that I have a better life than they do.

There's really no word other than weird.  It was damn weird.

They apologised later and we haven't communicated since.

Andrew

Permalink

How I Am - 3 Years Post Accident, 8+ Years With MS

From: andrew cooke <andrew@...>

Date: Fri, 26 Feb 2021 21:28:47 -0300

Previously

  https://acooke.org/cute/MyBikeAcci0.html
  https://acooke.org/cute/HowIAm80.html
  ...
  http://acooke.org/cute/HowIAm0.html

My back was fixed by yoga and a harder mattress.  The legs (no idea what was
wrong with the right) have now healed enough that I am generally able to sleep
on my side (although not always after exercise - sometimes it is more
comfortable to switch to the softer sofa).  There can be some discomfort (not
pain) on the sides of the legs when walking and the tendon at the back of my
left leg (which interfered with the rod, I think) can get a little sore.

Recovering fitness post-accident has taken longer than expected.  Covid didn't
help - there was a time when we were allowed out for just a few hours twice a
week (time spent shopping, mainly).  I am now, almost three years on, capable
of climbing to Yerba Loca, but not to Farellones.  My power and endurance are
still not as they were - I suspect it will be another 6 months work.  On the
other hand, Yoga on days off the bike seems to have improved my core.

MS has been pretty steady for years.  In fact, a few months ago I was feeling
like I had 6 months almost symptom-free.  Unfortunately that ended with
intermittent "white noise" in my right leg, the muscle above the knee.  This
has been coming and going a month or more now - it is disturbing, but seems to
have no effect on strength or balance.

In a few weeks I may be due for the Covid vaccine.  This will likely be
Sinovac and my doctor is a little worried (something to do with the underlying
carrier and my borked immune system, not the Covid part).  We will see what
happens.

We've spent, what, maybe a year pretty much isolated.  Paulina is working from
home.  We go out only to buy supplies or exercise.

Hello to future self.  Good luck.

Andrew

Permalink

[USA Politics] In America's Uncivil War Republicans Are The Aggressors

From: andrew cooke <andrew@...>

Date: Mon, 8 Feb 2021 16:48:13 -0300

An important article from Perry Bacon Jr.

https://fivethirtyeight.com/features/in-americas-uncivil-war-republicans-are-the-aggressors/

Andrew

Permalink

[Programming] Selenium and Python

From: andrew cooke <andrew@...>

Date: Sun, 24 Jan 2021 09:14:51 -0300

I am testing user registration (including email confirmation) and it's this
easy:


    def test_registration_flow(self):

        with self.record_error():

            self.driver.get(self.url)
            self.driver.find_element_by_link_text("register").click()

            with self.smtpd():
                self.driver.find_element_by_id('username').send_keys('andrew')
                self.driver.find_element_by_id('email').send_keys('andrew@...')
                self.driver.find_element_by_id('password').send_keys(self.PASSWORD)
                self.driver.find_element_by_id('retype_password').send_keys(self.PASSWORD)
                self.driver.find_element_by_css_selector('input.btn').click()

            match = re.compile(r'(http[:/a-zA-Z0-9\-_]+)', re.MULTILINE).search(self.email)
            self.driver.get(match.group(1))

            self.assertTrue('Your email has been confirmed.' in self.driver.page_source)
            self.assertTrue('You have signed in successfully.' in self.driver.page_source)

            with self.app_context():
                users = User.query.all()
                self.assertEqual(1, len(users))
                user = users[0]
                self.assertEqual('andrew', user.username)
                self.assertEqual('andrew@...', user.email)
                self.assertEqual(dt.date.today(), user.email_confirmed_at)
                self.assertTrue(len(user.password) > 20)
                self.assertTrue(user.password.startswith('$'))
                self.assertTrue(user.active)

Obvs there's some infrastructure in there (smtpd, etc), but I still find it
amazing that a browser starts, email is sent, and the user is registered.

Also, if an error occurs, I leave the database on disk and the browser open,
to make debugging easier.

Andrew

Permalink

[Bike] How Fast Before Walking More Efficient Than Cycling?

From: andrew cooke <andrew@...>

Date: Thu, 21 Jan 2021 22:52:02 -0300

Cycling is generally more efficient, per km travelled, than walking.  However,
someone just pointed out to me that if you go fast enough this is not the case
(since losses to air resistance increase rapidly with velocity).

So what velocities are we talking about?

A typical walking speed is 3.5 mph.  That's roughly 1.5 m/s.  Looking at
http://sprott.physics.wisc.edu/technote/walkrun.htm that correponds to roughly
300 W.

So walking burns 300 / 1.5 J/m = 200 J/m.

Now let's consider cycling.  A typical CdA on a road bike is 0.3.  Force is
CdA x p x v^2 where p is air density (1.2kg/m^3).  So F = 0.36 v^2 N.  The
energy required to travel 1m is then 0.36 v^2 J (force x distance).  Setting
this equal to 200 (above) we can find the break-even velocity: 200 = 0.36 v^2
so v = 23 m/s (50 mph).

I may have made a mistake, but that ties in with my intuition, so I'll go with
that for now.

Andrew

Permalink